Changeset 12509
- Timestamp:
- 2020-03-04T18:56:01+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl/src/ICE/icethd_pnd.F90
r12500 r12509 238 238 239 239 ! Lid growing 240 F ( ln_pnd_lids ) lh_ip_1d(ji) = lh_ip_1d(ji) + zdh_frz240 IF ( ln_pnd_lids ) lh_ip_1d(ji) = lh_ip_1d(ji) + zdh_frz 241 241 ELSE 242 242 zdh_frz = 0._wp … … 262 262 IF ( ln_pnd_overflow ) THEN 263 263 264 ! If pond area exceeds a_pnd_avail_1d(ji)* a_i_1d(ji) then reduce the pond volume265 IF ( a_ip_1d(ji) > a_pnd_avail_1d(ji)* a_i_1d(ji) ) THEN264 ! If pond area exceeds zfr_mlt * a_i_1d(ji) then reduce the pond volume 265 IF ( a_ip_1d(ji) > zfr_mlt * a_i_1d(ji) ) THEN 266 266 v_ip_old = v_ip_1d(ji) ! Save original volume before leak for future use 267 dh_ip_over = zpnd_aspect * a_pnd_avail_1d(ji)- h_ip_1d(ji) ! This will be a negative number267 dh_ip_over = zpnd_aspect * zfr_mlt - h_ip_1d(ji) ! This will be a negative number 268 268 dh_ip_over = MAX(dh_ip_over,max_h_diff_ts) ! Apply a limit 269 269 h_ip_1d(ji) = MAX(0._wp, h_ip_1d(ji) + dh_ip_over)
Note: See TracChangeset
for help on using the changeset viewer.