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

Timestamp:

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

Author:

edblockley

Message:

Fixing reproducibility issue in melt-ponds scheme

File:

• NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/src/ICE/icethd_pnd.F90 (modified) (2 diffs)

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