--- trunk/Sources/phylmd/fisrtilp.f 2016/09/01 10:30:53 207 +++ trunk/Sources/phylmd/fisrtilp.f 2017/02/27 15:44:55 213 @@ -8,9 +8,8 @@ d_ql, rneb, radliq, rain, snow, pfrac_impa, pfrac_nucl, pfrac_1nucl, & frac_impa, frac_nucl, prfl, psfl, rhcl) - ! From phylmd/fisrtilp.F, version 1.2 2004/11/09 16:55:40 + ! From phylmd/fisrtilp.F, version 1.2, 2004/11/09 16:55:40 ! First author: Z. X. Li (LMD/CNRS), 20 mars 1995 - ! Other authors: Olivier, AA, IM, YM, MAF ! Objet : condensation et pr\'ecipitation stratiforme, sch\'ema de ! nuage, sch\'ema de condensation \`a grande \'echelle (pluie). @@ -23,8 +22,6 @@ USE suphec_m, ONLY: rcpd, rd, retv, rg, rlstt, rlvtt, rtt USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2 - ! Arguments: - REAL, INTENT (IN):: dtime ! intervalle du temps (s) REAL, INTENT (IN):: paprs(klon, klev+1) ! pression a inter-couche REAL, INTENT (IN):: pplay(klon, klev) ! pression au milieu de couche @@ -362,10 +359,10 @@ zcl(i) = cld_lc_lsc zct(i) = 1./cld_tau_lsc END IF - ! quantit\'e d'eau \`a \'elminier. + ! quantit\'e d'eau \`a \'eliminer zchau(i) = zct(i)*dtime/real(ninter)*zoliq(i)* & (1.0-exp(-(zoliq(i)/zneb(i)/zcl(i))**2))*(1.-zfice(i)) - ! meme chose pour la glace. + ! m\^eme chose pour la glace IF (ptconv(i, k)) THEN zfroi(i) = dtime/real(ninter)/zdz(i)*zoliq(i)* & fallvc(zrhol(i))*zfice(i) @@ -395,7 +392,7 @@ END IF END DO - ! Calculer les tendances de q et de t: + ! Calculer les tendances de q et de t : DO i = 1, klon d_q(i, k) = zq(i) - q(i, k) d_t(i, k) = zt(i) - t(i, k) @@ -403,8 +400,8 @@ ! Calcul du lessivage stratiforme DO i = 1, klon - zprec_cond(i) = max(zcond(i)-zoliq(i), 0.0)* & - (paprs(i, k)-paprs(i, k+1))/rg + zprec_cond(i) = max(zcond(i) - zoliq(i), 0.0) & + * (paprs(i, k)-paprs(i, k+1))/rg IF (rneb(i, k)>0.0 .AND. zprec_cond(i)>0.) THEN ! lessivage nucleation LMD5 dans la couche elle-meme IF (t(i, k)>=ztglace) THEN