Changeset 14866

Timestamp:

2021-05-14T17:21:36+02:00 (3 years ago)

Author:

edblockley

Message:

Fixing reproducibility issue in melt-ponds scheme

Location:

NEMO/branches/UKMO/NEMO_4.0.4_GO8_package

Files:

• cfgs/SHARED/namelist_ice_ref (modified) (1 diff)
• src/ICE/icestp.F90 (modified) (1 diff)
• src/ICE/icethd_pnd.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/cfgs/SHARED/namelist_ice_ref

@@ -199 +199 @@
          rn_apnd_min =   0.15         !     minimum meltwater fraction contributing to pond growth (TOPO and LEV)
          rn_apnd_max =   0.85         !     maximum meltwater fraction contributing to pond growth (TOPO and LEV)
-         rn_pnd_flush=   0.01         !     pond flushing efficiency (tuning parameter) (LEV)
+         rn_pnd_flush=   0.1          !     pond flushing efficiency (tuning parameter) (LEV)
       ln_pnd_CST     = .false.        !  constant  melt ponds
          rn_apnd     =   0.2          !     prescribed pond fraction, at Tsu=0 degC

NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/src/ICE/icestp.F90

@@ -269 +269 @@
       ENDIF
       CALL ice_var_glo2eqv
-      CALL ice_var_agg(1)
+      CALL ice_var_agg(2)
       !
       CALL ice_dyn_init                ! set ice dynamics parameters

NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/src/ICE/icethd_pnd.F90

@@ -582 +582 @@
       
       zvolp(:,:) = 0._wp
+      zvolp_res(:,:) = 0._wp
 
       DO jl = 1, jpl
@@ -647 +648 @@
                   zTavg = zTavg + t_su(ji,jj,jl)*a_i(ji,jj,jl)
                END DO
-!! EWB: bug here - should be divided by at_i(ji,jj) but zTavg is not used anywhere
-               zTavg = zTavg / a_i(ji,jj,jl) !!! could get a division by zero here
+               zTavg = zTavg / at_i(ji,jj) !!! could get a division by zero here
          
                DO jl = 1, jpl-1