Changeset 14345

Timestamp:

2021-01-26T14:25:33+01:00 (3 years ago)

Author:

cbricaud

Message:

commit changes for tickets #2604 (improve scaling of mixing under sea-ice)

Location:

NEMO/branches/2021/dev_r14318_MOI-02_MOI-GLS-ICE

Files:

• cfgs/SHARED/namelist_ref (modified) (1 diff)
• src/OCE/ZDF/zdfgls.F90 (modified) (5 diffs)

NEMO/branches/2021/dev_r14318_MOI-02_MOI-GLS-ICE/cfgs/SHARED/namelist_ref

@@ -1225 +1225 @@
    rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
    rn_hsro       =  0.02   !  Minimum surface roughness
+   nn_mxlice     =     0   ! type of scaling under sea-ice
    rn_hsri       =  0.03   !  Ice-ocean roughness
    rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met>1)

NEMO/branches/2021/dev_r14318_MOI-02_MOI-GLS-ICE/src/OCE/ZDF/zdfgls.F90

@@ -26 +26 @@
    USE zdfmxl         ! mixed layer
    USE sbcwave , ONLY : hsw   ! significant wave height
+#if defined key_si3
+   USE ice, ONLY: hm_i, h_i
+#endif
+#if defined key_cice
+   USE sbc_ice, ONLY: h_i
+#endif
    !
    USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
@@ -51 +57 @@
    LOGICAL  ::   ln_length_lim     ! use limit on the dissipation rate under stable stratification (Galperin et al. 1988)
    LOGICAL  ::   ln_sigpsi         ! Activate Burchard (2003) modification for k-eps closure & wave breaking mixing
+   INTEGER  ::   nn_mxlice         ! type of scaling under sea-ice (=0/1/2/3)
    INTEGER  ::   nn_bc_surf        ! surface boundary condition (=0/1)
    INTEGER  ::   nn_bc_bot         ! bottom boundary condition (=0/1)
@@ -217 +224 @@
       END SELECT
       !
-      ! adapt roughness where there is sea ice
-      zhsro(:,:) = ( (1._wp-zice_fra(:,:)) * zhsro(:,:) + zice_fra(:,:) * rn_hsri )*tmask(:,:,1)  + (1._wp - tmask(:,:,1))*rn_hsro
+      SELECT CASE( nn_mxlice )       ! Type of scaling under sea-ice
+      !
+      CASE( 0 )                      ! No scaling under sea-ice
+         zhsro(:,:) = zhsro(:,:)
+         !
+      CASE( 1 )                      ! scaling with constant sea-ice roughness
+         zhsro(:,:) = ( (1._wp-zice_fra(:,:)) * zhsro(:,:) + zice_fra(:,:) * rn_hsri )*tmask(:,:,1)  + (1._wp - tmask(:,:,1))*rn_hsro
+         !
+      CASE( 2 )                      ! scaling with mean sea-ice thickness
+#if defined key_si3
+         zhsro(:,:) = ( (1._wp-zice_fra(:,:)) * zhsro(:,:) + zice_fra(:,:) * hm_i(:,:) )*tmask(:,:,1)  + (1._wp - tmask(:,:,1))*rn_hsro
+#endif
+         !
+      CASE( 3 )                      ! scaling with max sea-ice thickness
+#if defined key_si3 || defined key_cice
+         zhsro(:,:) = ( (1._wp-zice_fra(:,:)) * zhsro(:,:) + zice_fra(:,:) * MAXVAL(h_i(:,:,:),3) )*tmask(:,:,1)  + (1._wp - tmask(:,:,1))*rn_hsro
+#endif
+         !
+      END SELECT
       !
       DO_3D( 0, 0, 0, 0, 2, jpkm1 )  !==  Compute dissipation rate  ==!
@@ -775 +799 @@
       NAMELIST/namzdf_gls/rn_emin, rn_epsmin, ln_length_lim,       &
          &            rn_clim_galp, ln_sigpsi, rn_hsro, rn_hsri,   &
-         &            rn_crban, rn_charn, rn_frac_hs,              &
+         &            nn_mxlice, rn_crban, rn_charn, rn_frac_hs,   &
          &            nn_bc_surf, nn_bc_bot, nn_z0_met, nn_z0_ice, &
          &            nn_stab_func, nn_clos
@@ -815 +839 @@
          WRITE(numout,*) '      Type of closure                               nn_clos        = ', nn_clos
          WRITE(numout,*) '      Surface roughness (m)                         rn_hsro        = ', rn_hsro
-         WRITE(numout,*) '      Ice-ocean roughness (used if nn_z0_ice/=0)    rn_hsri        = ', rn_hsri
+
+         WRITE(numout,*) '      type of scaling under sea-ice nn_mxlice = ', nn_mxlice
+         IF( nn_mxlice == 1 ) &
+            WRITE(numout,*) '      Ice-ocean roughness (used if nn_z0_ice/=0)    rn_hsri        = ', rn_hsri
+         SELECT CASE( nn_mxlice )             ! Type of scaling under sea-ice
+            CASE( 0 )   ;   WRITE(numout,*) '   ==>>>   No scaling under sea-ice'
+            CASE( 1 )   ;   WRITE(numout,*) '   ==>>>   scaling with constant sea-ice thickness'
+            CASE( 2 )   ;   WRITE(numout,*) '   ==>>>   scaling with mean sea-ice thickness'
+            CASE( 3 )   ;   WRITE(numout,*) '   ==>>>   scaling with max sea-ice thickness'
+            CASE DEFAULT
+               CALL ctl_stop( 'zdf_tke_init: wrong value for nn_mxlice, should be 0,1,2,3 ')
+         END SELECT
          WRITE(numout,*)
       ENDIF