- Timestamp:
- 2017-05-29T16:58:42+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90
r8063 r8085 101 101 WRITE(numout,*)' Activate melt ponds ln_pnd = ', ln_pnd 102 102 WRITE(numout,*)' Type of melt pond scheme =0 presc, =1 empirical = 2 topo nn_pnd_scheme = ', nn_pnd_scheme 103 WRITE(numout,*)' Type of melt pond coupling =0 pass., =1 full, =2 rad, 3=fw nn_pnd_cpl = ', nn_pnd_cp 103 WRITE(numout,*)' Type of melt pond coupling =0 pass., =1 full, =2 rad, 3=fw nn_pnd_cpl = ', nn_pnd_cpl 104 104 WRITE(numout,*)' Prescribed pond fraction rn_apnd = ', rn_apnd 105 ENDIF 106 107 IF ( ln_pnd == .FALSE. ) THEN 108 WRITE(numout) ' Melt ponds are not activated ' 109 WRITE(numout) ' nn_pnd_scheme, nn_pnd_cpl and rn_apnd are set to zero ' 110 nn_pnd_scheme = 0 111 nn_pnd_cpl = 0 112 rn_apnd = 0._wp 105 113 ENDIF 106 114 ! … … 134 142 SELECT CASE ( nn_pnd_scheme ) 135 143 136 137 138 139 140 141 142 143 144 144 CASE (1) 145 146 CALL lim_mp_cesm ! empirical melt ponds 147 148 CASE (2) 149 150 CALL lim_mp_topo (at_i, a_i, & 151 & vt_i, v_i, v_s, t_i, s_i, a_ip_frac, & 152 & h_ip, t_su) 145 153 146 154 END SELECT … … 193 201 REAL(wp) :: z1_rhofw ! inverse freshwater density 194 202 REAL(wp) :: z1_zpnd_aspect ! inverse pond aspect ratio 195 REAL(wp) :: zv old! dummy pond volume203 REAL(wp) :: zvpold ! dummy pond volume 196 204 197 205 INTEGER :: ji, jj, jl, ij ! loop indices … … 261 269 ELSE !--- Update pond characteristics 262 270 263 !--- Add retained melt water 271 !--- Add retained melt water to melt ponds 264 272 v_ip(ji,jj,jl) = v_ip(ji,jj,jl) + zrfrac(ji,jj,jl) * z1_rhofw * zwfx_mlw(ji,jj) * a_i(ji,jj,jl) * rdt_ice 265 273
Note: See TracChangeset
for help on using the changeset viewer.