Changeset 14983 for NEMO/trunk/src/OCE/ZDF/zdftke.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/zdftke.F90 (modified) (2 diffs)

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