--- trunk/phylmd/clouds_gno.f90 2013/11/15 18:45:49 76 +++ trunk/phylmd/clouds_gno.f 2018/02/05 10:39:38 254 @@ -6,7 +6,7 @@ SUBROUTINE CLOUDS_GNO(klon, ND, R, RS, QSUB, PTCONV, RATQSC, CLDF) - ! From LMDZ4/libf/phylmd/clouds_gno.F, version 1.2 2004/11/09 16:55:40 + ! From LMDZ4/libf/phylmd/clouds_gno.F, version 1.2, 2004/11/09 16:55:40 use numer_rec_95, only: nr_erf @@ -155,8 +155,8 @@ ! erfv for v large: erfcu(i) = 1.0-NR_ERF(u(i)) - ! !!! ATTENTION : rajout d'un seuil pour l'exponentiel - aux(i) = sqrtpi*erfcu(i)*EXP(min(v(i)*v(i), 100.)) + ! Attention : ajout d'un seuil pour l'exponentielle + aux(i) = sqrtpi*erfcu(i)*EXP(min(v(i)*v(i), 80.)) coeff(i) = 1.0 - 1./2./(v(i)**2.) + 3./4./(v(i)**4.) block(i) = coeff(i) * EXP(-v(i)*v(i)) / v(i) / sqrtpi dist(i) = v(i) * aux(i) / coeff(i) - beta(i)