Changeset 3281 for trunk/LMDZ.GENERIC
- Timestamp:
- Mar 20, 2024, 4:29:27 PM (9 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/rain_generic.F90
r3102 r3281 281 281 d_t(i,k) = - d_q(i,k) * RLVTT_generic/RCPD 282 282 else 283 ! zqev/zqevt are the maximum amount of vapour that we are allowed to add by evaporation of rain 283 284 zqev = MAX (0.0, (zqs(i,k)-q(i,k)))*dmass(i,k)/(ptimestep*(1.d0+dqsat(i,k))) 284 285 !max evaporation to reach saturation implictly accounting for temperature reduction … … 286 287 *sqrt(precip_rate(i))*dmass(i,k)/pplay(i,k)*zt(i,k)*R) ! BC modif here, is it R or r/(mu/1000) ? 287 288 zqev = MIN (zqev, zqevt) 288 zqev = MAX (zqev, 0.0) ! a priori inutile d'après les précédentes lignes 289 zqev = MAX (zqev, 0.0) ! not necessary according to previous lines 290 291 ! we withdraw the evaporated rain 289 292 precip_rate_tmp(i)= precip_rate(i) - zqev 290 293 precip_rate_tmp(i)= max(precip_rate_tmp(i),0.0)
Note: See TracChangeset
for help on using the changeset viewer.