- Timestamp:
- 2017-12-22T13:26:37+01:00 (6 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRD
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRD/trd_oce.F90
r7179 r9163 57 57 ! 58 58 ! !!!* Passive tracers trends indices (use if "key_top" defined) 59 INTEGER, PUBLIC, PARAMETER :: jptra_sms = 19!: sources m. sinks60 INTEGER, PUBLIC, PARAMETER :: jptra_radn = 2 0!: corr. trn<0 in trcrad61 INTEGER, PUBLIC, PARAMETER :: jptra_radb = 2 1!: corr. trb<0 in trcrad (like atf)59 INTEGER, PUBLIC, PARAMETER :: jptra_sms = 21 !: sources m. sinks 60 INTEGER, PUBLIC, PARAMETER :: jptra_radn = 22 !: corr. trn<0 in trcrad 61 INTEGER, PUBLIC, PARAMETER :: jptra_radb = 23 !: corr. trb<0 in trcrad (like atf) 62 62 ! 63 63 ! !!!* Momentum trends indices -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90
r8104 r9163 20 20 USE trd_oce ! trends: ocean variables 21 21 USE trdtrc ! ocean passive mixed layer tracers trends 22 # if defined key_top 23 USE trc, ONLY: tra ! tracer definitions (trn, trb, tra, etc.) 24 # endif 22 25 USE trdglo ! trends: global domain averaged 23 26 USE trdpen ! trends: Potential ENergy … … 163 166 ENDIF 164 167 168 # if defined key_top 165 169 IF( ctype == 'TRC' ) THEN !== passive tracer trend ==! 166 170 ! … … 170 174 CASE( jptra_yad ) ; CALL trd_tra_adv( ptrd , pun , ptra, 'Y', ztrds ) 171 175 CASE( jptra_zad ) ; CALL trd_tra_adv( ptrd , pun , ptra, 'Z', ztrds ) 176 CASE( jptra_zdfp ) ! diagnose the "PURE" Kz trend (here: just before the swap) 177 ! ! iso-neutral diffusion case otherwise 178 ! jptra_zdf is "PURE" 179 CALL wrk_alloc( jpi, jpj, jpk, zws ) 180 ! 181 zws(:,:, 1 ) = 0._wp ! vertical diffusive fluxes 182 zws(:,:,jpk) = 0._wp 183 DO jk = 2, jpk 184 zws(:,:,jk) = avt(:,:,jk) * (tra(:,:,jk-1,ktra) - tra(:,:,jk,ktra) ) / fse3w(:,:,jk) * tmask(:,:,jk) 185 END DO 186 ! 187 ztrds(:,:,jpk) = 0._wp 188 DO jk = 1, jpkm1 189 ztrds(:,:,jk) = ( zws(:,:,jk) - zws(:,:,jk+1) ) / fse3t(:,:,jk) 190 END DO 191 CALL wrk_dealloc( jpi, jpj, jpk, zws ) 192 ! 172 193 CASE DEFAULT ! other trends: just masked 173 194 ztrds(:,:,:) = ptrd(:,:,:) * tmask(:,:,:) … … 177 198 ! 178 199 ENDIF 200 # endif 179 201 ! 180 202 CALL wrk_dealloc( jpi, jpj, jpk, ztrds )
Note: See TracChangeset
for help on using the changeset viewer.