Changeset 8453 for branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRD

Timestamp:

2017-08-22T18:54:47+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM --22-08-2017-- repair trends for passive tracers

Location:

branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRD

Files:

• trd_oce.F90 (modified) (1 diff)
• trdtra.F90 (modified) (2 diffs)

branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRD/trd_oce.F90

@@ -57 +57 @@
    !
    !                                                  !!!* Passive tracers trends indices (use if "key_top" defined)
-   INTEGER, PUBLIC, PARAMETER ::   jptra_sms  = 19     !: sources m. sinks
-   INTEGER, PUBLIC, PARAMETER ::   jptra_radn = 20     !: corr. trn<0 in trcrad
-   INTEGER, PUBLIC, PARAMETER ::   jptra_radb = 21     !: corr. trb<0 in trcrad (like atf)
+   INTEGER, PUBLIC, PARAMETER ::   jptra_sms  = 21     !: sources m. sinks
+   INTEGER, PUBLIC, PARAMETER ::   jptra_radn = 22     !: corr. trn<0 in trcrad
+   INTEGER, PUBLIC, PARAMETER ::   jptra_radb = 23     !: corr. trb<0 in trcrad (like atf)
    !
    !                                                  !!!* Momentum trends indices

branches/NERC/dev_r5518_GO6_package_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(:,:,:)