Changeset 14983 for NEMO/trunk/src/OCE/ZDF/zdfphy.F90

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)

File:

• NEMO/trunk/src/OCE/ZDF/zdfphy.F90 (modified) (1 diff)

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')
       !