Changeset 9173 for branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90

Timestamp:

2018-01-03T14:19:47+01:00 (6 years ago)

Author:

jpalmier

Message:

JPALM -- merge with GO6_package 9114-9163 -- help test up-to-date suite

File:

• branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90 (modified) (4 diffs)

branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90

@@ -20 +20 @@
    USE trd_oce        ! trends: ocean variables
    USE trdtrc         ! ocean passive mixed layer tracers trends 
+# if defined key_top
+   USE trc,           ONLY: tra    ! tracer definitions (trn, trb, tra, etc.) 
+# endif
    USE trdglo         ! trends: global domain averaged
    USE trdpen         ! trends: Potential ENergy
@@ -163 +166 @@
       ENDIF
 
+# if defined key_top
       IF( ctype == 'TRC' ) THEN                           !==  passive tracer trend  ==!
          !
@@ -170 +174 @@
          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(:,:,:)
@@ -177 +198 @@
          !
       ENDIF
+# endif
       !
       CALL wrk_dealloc( jpi, jpj, jpk, ztrds )