New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10965 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DIA/diatmb.F90 – NEMO

Ignore:
Timestamp:
2019-05-10T18:02:51+02:00 (5 years ago)
Author:
davestorkey
Message:

2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps : DIA and stpctl.F90. Just testing in ORCA1 so far.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DIA/diatmb.F90

    r10499 r10965  
    9494 
    9595 
    96    SUBROUTINE dia_tmb 
     96   SUBROUTINE dia_tmb( Kmm ) 
    9797      !!---------------------------------------------------------------------- 
    9898      !!                 ***  ROUTINE dia_tmb  *** 
    … …  
    103103      !! 
    104104      !!-------------------------------------------------------------------- 
     105      INTEGER, INTENT(in) :: Kmm     ! time level index 
     106      ! 
    105107      REAL(wp) ::   zmdi =1.e+20     ! land value 
    106108      REAL(wp), DIMENSION(jpi,jpj,3) :: zwtmb    ! workspace  
    107109      !!-------------------------------------------------------------------- 
    108110      ! 
    109       CALL dia_calctmb( tsn(:,:,:,jp_tem), zwtmb ) 
     111      CALL dia_calctmb( ts(:,:,:,jp_tem,Kmm), zwtmb ) 
    110112      !ssh already output but here we output it masked 
    111113      IF( ll_wd ) THEN 
    112          CALL iom_put( "sshnmasked", (sshn(:,:)+ssh_ref)*tmask(:,:,1) + zmdi*(1.0 - tmask(:,:,1)) ) 
     114         CALL iom_put( "sshnmasked", (ssh(:,:,Kmm)+ssh_ref)*tmask(:,:,1) + zmdi*(1.0 - tmask(:,:,1)) ) 
    113115      ELSE 
    114          CALL iom_put( "sshnmasked", sshn(:,:)*tmask(:,:,1) + zmdi*(1.0 - tmask(:,:,1)) ) 
     116         CALL iom_put( "sshnmasked", ssh(:,:,Kmm)*tmask(:,:,1) + zmdi*(1.0 - tmask(:,:,1)) ) 
    115117      ENDIF 
    116118 
    … …  
    119121      CALL iom_put( "bot_temp"  , zwtmb(:,:,3) )    ! tmb Temperature 
    120122      ! 
    121       CALL dia_calctmb( tsn(:,:,:,jp_sal), zwtmb ) 
     123      CALL dia_calctmb( ts(:,:,:,jp_sal,Kmm), zwtmb ) 
    122124      CALL iom_put( "top_sal"   , zwtmb(:,:,1) )    ! tmb Salinity  
    123125      CALL iom_put( "mid_sal"   , zwtmb(:,:,2) )    ! tmb Salinity 
    124126      CALL iom_put( "bot_sal"   , zwtmb(:,:,3) )    ! tmb Salinity 
    125127      ! 
    126       CALL dia_calctmb( un(:,:,:), zwtmb ) 
     128      CALL dia_calctmb( uu(:,:,:,Kmm), zwtmb ) 
    127129      CALL iom_put( "top_u"     , zwtmb(:,:,1) )    ! tmb  U Velocity 
    128130      CALL iom_put( "mid_u"     , zwtmb(:,:,2) )    ! tmb  U Velocity 
    129131      CALL iom_put( "bot_u"     , zwtmb(:,:,3) )    ! tmb  U Velocity 
    130132      ! 
    131       CALL dia_calctmb( vn(:,:,:), zwtmb ) 
     133      CALL dia_calctmb( vv(:,:,:,Kmm), zwtmb ) 
    132134      CALL iom_put( "top_v"     , zwtmb(:,:,1) )    ! tmb  V Velocity 
    133135      CALL iom_put( "mid_v"     , zwtmb(:,:,2) )    ! tmb  V Velocity 
Note: See TracChangeset for help on using the changeset viewer.