- Timestamp:
- 2018-07-12T17:55:41+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/ICE/icethd_pnd.F90
r9923 r9937 133 133 REAL(wp) :: zdv_mlt ! available meltwater for melt ponding 134 134 REAL(wp) :: z1_Tp ! inverse reference temperature 135 REAL(wp) :: z1_rhofw ! inverse freshwater density136 135 REAL(wp) :: z1_zpnd_aspect ! inverse pond aspect ratio 137 136 REAL(wp) :: zfac, zdum … … 139 138 INTEGER :: ji ! loop indices 140 139 !!------------------------------------------------------------------- 141 z1_rhofw = 1._wp / rhofw142 140 z1_zpnd_aspect = 1._wp / zpnd_aspect 143 141 z1_Tp = 1._wp / zTp … … 157 155 ! 158 156 ! available meltwater for melt ponding [m, >0] and fraction 159 zdv_mlt = -( dh_i_sum(ji)*rhoi c + dh_s_mlt(ji)*rhosn ) * z1_rhofw * a_i_1d(ji)157 zdv_mlt = -( dh_i_sum(ji)*rhoi + dh_s_mlt(ji)*rhos ) * r1_rhow * a_i_1d(ji) 160 158 zfr_mlt = zrmin + ( zrmax - zrmin ) * a_i_1d(ji) ! from CICE doc 161 159 !zfr_mlt = zrmin + zrmax * a_i_1d(ji) ! from Holland paper … … 168 166 ! melt pond mass flux (<0) 169 167 IF( ln_pnd_fwb .AND. zdv_mlt > 0._wp ) THEN 170 zfac = zfr_mlt * zdv_mlt * rho fw * r1_Dt_ice168 zfac = zfr_mlt * zdv_mlt * rhow * r1_Dt_ice 171 169 wfx_pnd_1d(ji) = wfx_pnd_1d(ji) - zfac 172 170 !
Note: See TracChangeset
for help on using the changeset viewer.