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 6060 for branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90 – NEMO

Ignore:
Timestamp:
2015-12-16T10:25:22+01:00 (8 years ago)
Author:
timgraham
Message:

Merged dev_r5836_noc2_VVL_BY_DEFAULT into branch

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90

    r5836 r6060  
    3131   USE dtatsd         ! data: temperature & salinity 
    3232   USE zdfmxl         ! vertical physics: mixed layer depth 
     33   ! 
    3334   USE in_out_manager ! I/O manager 
    3435   USE lib_mpp        ! MPP library 
     
    4142   PRIVATE 
    4243 
    43    PUBLIC   tra_dmp      ! routine called by step.F90 
    44    PUBLIC   tra_dmp_init ! routine called by nemogcm.F90 
     44   PUBLIC   tra_dmp        ! called by step.F90 
     45   PUBLIC   tra_dmp_init   ! called by nemogcm.F90 
    4546 
    4647   !                                           !!* Namelist namtra_dmp : T & S newtonian damping * 
     
    5253 
    5354   !! * Substitutions 
    54 #  include "domzgr_substitute.h90" 
    5555#  include "vectopt_loop_substitute.h90" 
    5656   !!---------------------------------------------------------------------- 
     
    8989      !!      below the well mixed layer (nlmdmp=2) 
    9090      !! 
    91       !! ** Action  : - (ta,sa)  tracer trends updated with the damping trend 
     91      !! ** Action  : - tsa: tracer trends updated with the damping trend 
    9292      !!---------------------------------------------------------------------- 
    9393      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
     
    139139            DO jj = 2, jpjm1 
    140140               DO ji = fs_2, fs_jpim1   ! vector opt. 
    141                   IF( fsdept(ji,jj,jk) >= hmlp (ji,jj) ) THEN 
     141                  IF( gdept_n(ji,jj,jk) >= hmlp (ji,jj) ) THEN 
    142142                     tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem)   & 
    143143                        &                 + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_tem) - tsb(ji,jj,jk,jp_tem) ) 
     
    177177      !!---------------------------------------------------------------------- 
    178178      INTEGER ::   ios, imask   ! local integers  
    179       !! 
     179      ! 
    180180      NAMELIST/namtra_dmp/ ln_tradmp, nn_zdmp, cn_resto 
    181181      !!---------------------------------------------------------------------- 
     
    229229   END SUBROUTINE tra_dmp_init 
    230230 
     231   !!====================================================================== 
    231232END MODULE tradmp 
Note: See TracChangeset for help on using the changeset viewer.