Changeset 14983

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

Files:

• cfgs/SHARED/field_def_nemo-oce.xml (modified) (1 diff)
• src/OCE/ZDF/zdfphy.F90 (modified) (1 diff)
• src/OCE/ZDF/zdftke.F90 (modified) (2 diffs)

NEMO/trunk/cfgs/SHARED/field_def_nemo-oce.xml

@@ -755 +755 @@
     <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                unit="s-2" />
 
+    <!-- dissipation diagnostics (note: ediss_k is only available with tke scheme) -->   
+    <field id="avt_k"        long_name="vertical eddy diffusivity from closure schemes" standard_name="ocean_vertical_eddy_diffusivity"       unit="m2/s" />
+    <field id="avm_k"        long_name="vertical eddy viscosity from closure schemes"   standard_name="ocean_vertical_eddy_viscosity"         unit="m2/s" />
+    <field id="ediss_k"      long_name="Kolmogorov energy dissipation (tke scheme)"     standard_name="Kolmogorov_energy_dissipation"         unit="W/kg" />
+    <field id="eshear_k"     long_name="energy source from vertical shear"              standard_name="energy_source_from_shear"              unit="W/kg" />
+    <field id="estrat_k"     long_name="energy sink from stratification"                standard_name="energy_sink_from_stratification"       unit="W/kg" />
+    
   </field_group>
 

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