Changeset 8471 for branches/2017/nemo_v3_6_STABLE_trdtrc/NEMOGCM/NEMO/OPA_SRC

Timestamp:

2017-08-29T18:12:42+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- #1933 -- add and correct passive tracer trends - needs improvements and testing

File:

• branches/2017/nemo_v3_6_STABLE_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90 (modified) (2 diffs)

branches/2017/nemo_v3_6_STABLE_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90

@@ -20 +20 @@
    USE trd_oce        ! trends: ocean variables
    USE trdtrc         ! ocean passive mixed layer tracers trends 
+   USE trc,           ONLY: tra    ! tracer definitions (trn, trb, tra, etc.) 
    USE trdglo         ! trends: global domain averaged
    USE trdpen         ! trends: Potential ENergy
@@ -170 +171 @@
          CASE( jptra_yad )   ;   CALL trd_tra_adv( ptrd , pun , ptra, 'Y', ztrds ) 
          CASE( jptra_zad )   ;   CALL trd_tra_adv( ptrd , pun , ptra, 'Z', ztrds ) 
+         CASE( jptra_zdfp )           ! diagnose the "PURE" Kz trend (here: just before the swap)
+            !                         ! iso-neutral diffusion case otherwise
+            !                         jptra_zdf is "PURE"
+            CALL wrk_alloc( jpi, jpj, jpk, zws )
+            !
+            zws(:,:, 1 ) = 0._wp                        ! vertical diffusive fluxes
+            zws(:,:,jpk) = 0._wp  
+            DO jk = 2, jpk
+               zws(:,:,jk) = avt(:,:,jk) * (tra(:,:,jk-1,ktra) - tra(:,:,jk,ktra) ) / fse3w(:,:,jk) * tmask(:,:,jk)
+            END DO
+            !
+            ztrds(:,:,jpk) = 0._wp   
+            DO jk = 1, jpkm1
+               ztrds(:,:,jk) = ( zws(:,:,jk) - zws(:,:,jk+1) ) / fse3t(:,:,jk)
+            END DO
+            CALL wrk_dealloc( jpi, jpj, jpk, zws )
+            !
          CASE DEFAULT                 ! other trends: just masked 
                                  ztrds(:,:,:) = ptrd(:,:,:) * tmask(:,:,:)