Changeset 14983 for NEMO/trunk/src/OCE/ZDF

Timestamp:

2021-06-11T17:11:20+02:00 (3 years ago)

Author:

clem

Message:

trunk: small sneaky enhancement related to ticket #2661 (adding dissipation outputs)

Location:

NEMO/trunk/src/OCE/ZDF

Files:

• zdfphy.F90 (modified) (1 diff)
• zdftke.F90 (modified) (2 diffs)

NEMO/trunk/src/OCE/ZDF/zdfphy.F90

@@ -375 +375 @@
       ENDIF
       !
+      ! diagnostics of energy dissipation
+      IF( iom_use('avt_k') .OR. iom_use('avm_k') .OR. iom_use('eshear_k') .OR. iom_use('estrat_k') ) THEN
+         IF( l_zdfsh2 ) THEN
+            zsh2(:,:,1  ) = 0._wp
+            zsh2(:,:,jpk) = 0._wp
+            CALL lbc_lnk( 'zdfphy', zsh2, 'W', 1._wp )
+            CALL iom_put( 'avt_k'   ,   avt_k       * wmask )
+            CALL iom_put( 'avm_k'   ,   avm_k       * wmask )
+            CALL iom_put( 'eshear_k',   zsh2        * wmask )
+            CALL iom_put( 'estrat_k', - avt_k * rn2 * wmask )
+         ENDIF
+      ENDIF
+      !
       IF( ln_timing )   CALL timing_stop('zdf_phy')
       !

NEMO/trunk/src/OCE/ZDF/zdftke.F90

@@ -219 +219 @@
       REAL(wp), DIMENSION(A2D(nn_hls))     ::   zice_fra, zhlc, zus3, zWlc2
       REAL(wp), DIMENSION(A2D(nn_hls),jpk) ::   zpelc, zdiag, zd_up, zd_lw
+      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztmp ! for diags
       !!--------------------------------------------------------------------
       !
@@ -446 +447 @@
          en(ji,jj,jk) = MAX( en(ji,jj,jk), rn_emin ) * wmask(ji,jj,jk)
       END_3D
+      !
+      ! Kolmogorov energy of dissipation (W/kg)
+      !    ediss = Ce*sqrt(en)/L*en
+      !    dissl = sqrt(en)/L
+      IF( iom_use('ediss_k') ) THEN
+         ALLOCATE( ztmp(A2D(nn_hls),jpk) )
+         ztmp(:,:,:) = zfact3 * dissl * en * wmask
+         CALL lbc_lnk( 'zdftke', ztmp, 'W', 1._wp )
+         CALL iom_put( 'ediss_k', ztmp )
+         DEALLOCATE( ztmp )
+      ENDIF
       !
       !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<