362 |
zcl(i) = cld_lc_lsc |
zcl(i) = cld_lc_lsc |
363 |
zct(i) = 1./cld_tau_lsc |
zct(i) = 1./cld_tau_lsc |
364 |
END IF |
END IF |
365 |
! quantit\'e d'eau \`a \'elminier. |
! quantit\'e d'eau \`a \'eliminer |
366 |
zchau(i) = zct(i)*dtime/real(ninter)*zoliq(i)* & |
zchau(i) = zct(i)*dtime/real(ninter)*zoliq(i)* & |
367 |
(1.0-exp(-(zoliq(i)/zneb(i)/zcl(i))**2))*(1.-zfice(i)) |
(1.0-exp(-(zoliq(i)/zneb(i)/zcl(i))**2))*(1.-zfice(i)) |
368 |
! meme chose pour la glace. |
! m\^eme chose pour la glace |
369 |
IF (ptconv(i, k)) THEN |
IF (ptconv(i, k)) THEN |
370 |
zfroi(i) = dtime/real(ninter)/zdz(i)*zoliq(i)* & |
zfroi(i) = dtime/real(ninter)/zdz(i)*zoliq(i)* & |
371 |
fallvc(zrhol(i))*zfice(i) |
fallvc(zrhol(i))*zfice(i) |
395 |
END IF |
END IF |
396 |
END DO |
END DO |
397 |
|
|
398 |
! Calculer les tendances de q et de t: |
! Calculer les tendances de q et de t : |
399 |
DO i = 1, klon |
DO i = 1, klon |
400 |
d_q(i, k) = zq(i) - q(i, k) |
d_q(i, k) = zq(i) - q(i, k) |
401 |
d_t(i, k) = zt(i) - t(i, k) |
d_t(i, k) = zt(i) - t(i, k) |