Changeset 6140 for trunk/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90
- Timestamp:
- 2015-12-21T12:35:23+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90
r5836 r6140 31 31 USE dtatsd ! data: temperature & salinity 32 32 USE zdfmxl ! vertical physics: mixed layer depth 33 ! 33 34 USE in_out_manager ! I/O manager 34 35 USE lib_mpp ! MPP library … … 41 42 PRIVATE 42 43 43 PUBLIC tra_dmp ! routinecalled by step.F9044 PUBLIC tra_dmp_init ! routinecalled by nemogcm.F9044 PUBLIC tra_dmp ! called by step.F90 45 PUBLIC tra_dmp_init ! called by nemogcm.F90 45 46 46 47 ! !!* Namelist namtra_dmp : T & S newtonian damping * … … 52 53 53 54 !! * Substitutions 54 # include "domzgr_substitute.h90"55 55 # include "vectopt_loop_substitute.h90" 56 56 !!---------------------------------------------------------------------- … … 89 89 !! below the well mixed layer (nlmdmp=2) 90 90 !! 91 !! ** Action : - (ta,sa)tracer trends updated with the damping trend91 !! ** Action : - tsa: tracer trends updated with the damping trend 92 92 !!---------------------------------------------------------------------- 93 93 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 100 100 ! 101 101 CALL wrk_alloc( jpi,jpj,jpk,jpts, zts_dta ) 102 !103 102 IF( l_trdtra ) THEN !* Save ta and sa trends 104 103 CALL wrk_alloc( jpi,jpj,jpk,jpts, ztrdts ) … … 139 138 DO jj = 2, jpjm1 140 139 DO ji = fs_2, fs_jpim1 ! vector opt. 141 IF( fsdept(ji,jj,jk) >= hmlp (ji,jj) ) THEN140 IF( gdept_n(ji,jj,jk) >= hmlp (ji,jj) ) THEN 142 141 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & 143 142 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_tem) - tsb(ji,jj,jk,jp_tem) ) … … 177 176 !!---------------------------------------------------------------------- 178 177 INTEGER :: ios, imask ! local integers 179 ! !178 ! 180 179 NAMELIST/namtra_dmp/ ln_tradmp, nn_zdmp, cn_resto 181 180 !!---------------------------------------------------------------------- … … 229 228 END SUBROUTINE tra_dmp_init 230 229 230 !!====================================================================== 231 231 END MODULE tradmp
Note: See TracChangeset
for help on using the changeset viewer.