Changeset 12497

Timestamp:

2020-03-02T12:45:12+01:00 (3 years ago)

Author:

dancopsey

Message:

Fix volume and area calculations after melt ponds overflow to the ocean.

File:

• NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl/src/ICE/icethd_pnd.F90 (modified) (4 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl/src/ICE/icethd_pnd.F90

@@ -245 +245 @@
                 dh_ip_over = zpnd_aspect * a_pnd_avail_1d(ji) - h_ip_1d(ji)   ! This will be a negative number
                 dh_ip_over = MAX(dh_ip_over,max_h_diff_ts)                       ! Apply a limit
-                h_ip_1d(ji) = h_ip_1d(ji) + dh_ip_over
-                a_ip_1d(ji) = a_pnd_avail_1d(ji) * a_i_1d(ji)
+                h_ip_1d(ji) = MAX(0._wp, h_ip_1d(ji) + dh_ip_over)
+                a_ip_frac_1d(ji) = h_ip_1d(ji) / zpnd_aspect
+                a_ip_1d(ji) = a_ip_frac_1d(ji) * a_i_1d(ji)
                 v_ip_1d(ji) = h_ip_1d(ji) * a_ip_1d(ji)
-                a_ip_frac_1d(ji) = a_ip_1d(ji) / a_i_1d(ji)
 
             ENDIF
@@ -257 +257 @@
                 dh_ip_over = 0.5_wp * h_i_1d(ji) - h_ip_1d(ji)                ! This will be a negative number
                 dh_ip_over = MAX(dh_ip_over,max_h_diff_ts)                       ! Apply a limit
-                h_ip_1d(ji) = h_ip_1d(ji) + dh_ip_over
+                h_ip_1d(ji) = MAX(0._wp, h_ip_1d(ji) + dh_ip_over)
                 a_ip_frac_1d(ji) = h_ip_1d(ji) / zpnd_aspect
                 a_ip_1d(ji) = a_ip_frac_1d(ji) * a_i_1d(ji)
@@ -289 +289 @@
                 v_ip_old = v_ip_1d(ji)   ! Save original volume before leak for future use
                 dh_ip_over = MAX(dh_ip_over,max_h_diff_ts)                       ! Apply a limit
-                h_ip_1d(ji) = h_ip_1d(ji) + dh_ip_over
+                h_ip_1d(ji) = MAX(0._wp, h_ip_1d(ji) + dh_ip_over)
                 a_ip_frac_1d(ji) = h_ip_1d(ji) / zpnd_aspect
                 a_ip_1d(ji) = a_ip_frac_1d(ji) * a_i_1d(ji)
@@ -314 +314 @@
          ENDIF
 
-      END DO
       END DO
       !