Changeset 12686 for NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl/src/ICE

Timestamp:

2020-04-06T11:56:20+02:00 (4 years ago)

Author:

dancopsey

Message:

• We don't need maximum outflow when reducing pond width or depth.
  • Move code changes in sbccpl to avoid conflicts.

File:

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

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

@@ -133 +133 @@
       REAL(wp), PARAMETER ::   zpnd_aspect = 0.8_wp   ! pond aspect ratio
       REAL(wp), PARAMETER ::   zTp         = -2._wp   ! reference temperature
-      REAL(wp), PARAMETER ::   max_h_diff_s = -1.0E-6 ! Maximum meltpond depth change due to leaking or overflow (m s-1)
       REAL(wp), PARAMETER ::   pnd_lid_max = 0.015_wp ! pond lid thickness above which the ponds disappear from the albedo calculation
       REAL(wp), PARAMETER ::   pnd_lid_min = 0.005_wp ! pond lid thickness below which the full pond area is used in the albedo calculation
@@ -149 +148 @@
       REAL(wp) ::   lh_ip_end        ! Lid thickness at end of timestep (temporary variable)
       REAL(wp) ::   zdh_frz          ! Amount of melt pond that freezes (m)
-      REAL(wp) ::   v_ip_old         ! Pond volume before leaking back to the ocean
-      REAL(wp) ::   dh_ip_over       ! Pond thickness change due to overflow or leaking
+      REAL(wp) ::   dh_ip_over       ! Pond thickness change due to leaking
       REAL(wp) ::   dh_i_pndleak     ! Grid box mean change in water depth due to leaking back to the ocean
       REAL(wp) ::   h_gravity_head   ! Height above sea level of the top of the melt pond
@@ -168 +166 @@
       z1_Tp          = 1._wp / zTp 
       z1_rhoi        = 1._wp / rhoi
-      max_h_diff_ts  = max_h_diff_s * rdt_ice
 
       ! Define time-independent field for use in refreezing
@@ -266 +263 @@
                ! If pond area exceeds zfr_mlt * a_i_1d(ji) then reduce the pond volume
                IF ( a_ip_1d(ji) > zfr_mlt * a_i_1d(ji) ) THEN
-                   v_ip_old = v_ip_1d(ji)   ! Save original volume before leak for future use
-                   dh_ip_over = zpnd_aspect * zfr_mlt - 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) = MAX(0._wp, h_ip_1d(ji) + dh_ip_over)
+                   h_ip_1d(ji) = zpnd_aspect * zfr_mlt
                    a_ip_frac_1d(ji) = h_ip_1d(ji) / zpnd_aspect
                    a_ip_1d(ji) = a_ip_frac_1d(ji) * a_i_1d(ji)
@@ -277 +271 @@
                ! If pond depth exceeds half the ice thickness then reduce the pond volume
                IF ( h_ip_1d(ji) > 0.5_wp * h_i_1d(ji) ) THEN
-                   v_ip_old = v_ip_1d(ji)   ! Save original volume before leak for future use
-                   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) = MAX(0._wp, h_ip_1d(ji) + dh_ip_over)
+                   h_ip_1d(ji) = 0.5_wp * h_i_1d(ji)
                    a_ip_frac_1d(ji) = h_ip_1d(ji) / zpnd_aspect
                    a_ip_1d(ji) = a_ip_frac_1d(ji) * a_i_1d(ji)
@@ -309 +300 @@
                    dh_ip_over = MIN(dh_ip_over, 0._wp)      ! Make sure it is negative
 
-                   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) = MAX(0._wp, h_ip_1d(ji) + dh_ip_over)
                    a_ip_frac_1d(ji) = h_ip_1d(ji) / zpnd_aspect