Changeset 12394 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icethd_pnd.F90

Timestamp:

2020-02-18T10:44:02+01:00 (4 years ago)

Author:

dancopsey

Message:

Fix compile and XIOS errors.

File:

• NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icethd_pnd.F90 (modified) (5 diffs)

NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icethd_pnd.F90

@@ -151 +151 @@
 
       ! Define time-independent field for use in refreezing
-      omega_dt = 2.0_wp * rcnd_i * rdtice / (rLfus * rhow)
+      omega_dt = 2.0_wp * rcnd_i * rdt_ice / (rLfus * rhow)
 
       DO ji = 1, npti
@@ -167 +167 @@
             h_ip_1d(ji)      = 0._wp
             lh_ip_1d(ji)     = 0._wp
+
+            actual_mlt = 0._wp
+            actual_frz = 0._wp
             !                                                     !--------------------------------!
          ELSE                                                     ! Case ice thickness >= rn_himin !
@@ -190 +193 @@
                
                ! The following equation is a rearranged form of:
-               ! lid_thickness_end - lid_thickness_start = rcnd_i * t_grad * rdtice / (0.5*(lid_thickness_end + lid_thickness_start) * rLfus * rhow)
+               ! lid_thickness_end - lid_thickness_start = rcnd_i * t_grad * rdt_ice / (0.5*(lid_thickness_end + lid_thickness_start) * rLfus * rhow)
                ! where: lid_thickness_start = lh_ip_1d(ji)
                !        lid_thickness_end = lh_ip_end
@@ -209 +212 @@
             END IF
 
-            ! melt pond mass flux diagnostic (melt only)
-            zfac = actual_mlt * a_i_1d(ji * rhow * r1_rdtice
-            wfx_pnd_1d(ji) = wfx_pnd_1d(ji) - zfac
-            !
-            ! adjust ice/snow melting flux to balance melt pond flux (>0)
-            zdum = zfac / ( wfx_snw_sum_1d(ji) + wfx_sum_1d(ji) )
-            wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) * (1._wp + zdum)
-            wfx_sum_1d(ji)     = wfx_sum_1d(ji)     * (1._wp + zdum)
+            ! melt pond mass flux (<0)
+            IF( zdv_mlt > 0._wp ) THEN
+               zfac = actual_mlt * a_i_1d(ji) * rhow * r1_rdtice
+               wfx_pnd_1d(ji) = wfx_pnd_1d(ji) - zfac
+               !
+               ! adjust ice/snow melting flux to balance melt pond flux (>0)
+               zdum = zfac / ( wfx_snw_sum_1d(ji) + wfx_sum_1d(ji) )
+               wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) * (1._wp + zdum)
+               wfx_sum_1d(ji)     = wfx_sum_1d(ji)     * (1._wp + zdum)
+            ENDIF
             !
             ! Make sure pond volume or lid thickness has not gone negative
@@ -243 +248 @@
             write(numout,*)'icethd_pnd: a_i_1d(ji), v_ip_1d(ji), t_su_1d(ji), zfr_mlt, zdv_mlt = ',a_i_1d(ji), ' ', v_ip_1d(ji), ' ', t_su_1d(ji), ' ', zfr_mlt, ' ', zdv_mlt
             write(numout,*)'icethd_pnd: meltt = ', -( dh_i_sum(ji)*rhoi + dh_s_mlt(ji)*rhos ) / rhoi
+            write(numout,*)'icethd_pnd: lh_ip_1d(ji), actual_mlt, actual_frz, t_su_1d(ji) = ',lh_ip_1d(ji), ' ', actual_mlt, ' ', actual_frz, ' ', t_su_1d(ji)
          END IF