/[lmdze]/trunk/Sources/phylmd/clouds_gno.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/clouds_gno.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 227 by guez, Thu Nov 2 15:47:03 2017 UTC
# Line 155  contains Line 155  contains
155                        ! erfv for v large:                        ! erfv for v large:
156    
157                        erfcu(i) = 1.0-NR_ERF(u(i))                        erfcu(i) = 1.0-NR_ERF(u(i))
158                        ! !!! ATTENTION : rajout d'un seuil pour l'exponentiel                        ! Attention : ajout d'un seuil pour l'exponentielle
159                        aux(i) = sqrtpi*erfcu(i)*EXP(min(v(i)*v(i), 100.))                        aux(i) = sqrtpi*erfcu(i)*EXP(min(v(i)*v(i), 80.))
160                        coeff(i) = 1.0 - 1./2./(v(i)**2.) + 3./4./(v(i)**4.)                        coeff(i) = 1.0 - 1./2./(v(i)**2.) + 3./4./(v(i)**4.)
161                        block(i) = coeff(i) * EXP(-v(i)*v(i)) / v(i) / sqrtpi                        block(i) = coeff(i) * EXP(-v(i)*v(i)) / v(i) / sqrtpi
162                        dist(i) = v(i) * aux(i) / coeff(i) - beta(i)                        dist(i) = v(i) * aux(i) / coeff(i) - beta(i)
Legend:
Removed from v.134  
changed lines
  Added in v.227
  ViewVC Help
Powered by ViewVC 1.1.21