When I read the code file "module_mp_milbrandt2mom.F" , I find that "espQ" shows up so many times. However , I can't understand what the "epsQ" represent ? What is this value meaning?
code such as :
! Graupel:
if (QG(i,k)>epsQ) then
iQG = 1./QG(i,k)
iNG = 1./NG(i,k)
iLAMg = max( iLAMmin1, iLAMDA_x(DE(i,k),QG(i,k),iNG,iGG99,thrd) )
iLAMg2 = iLAMg**2
iLAMgB0= iLAMg**(bfg)
iLAMgB1= iLAMg**(bfg+1.)
iLAMgB2= iLAMg**(bfg+2.)
!No_g = (NG(i,k))*iGG31/iLAMg**(1.+alpha_g)
No_g= NG(i,k)*iGG31/iLAMg !optimized for alpha_g=0
vg0 = gamfact(i,k)*ckQg1*iLAMgB0
Dg = Dm_x(DE(i,k),QG(i,k),iNG,icmg,thrd)
else
iLAMg = 0.; vg0 = 0.; Dg = 0.; No_g = 0.
iLAMg2 = 0.; iLAMgB0= 0.; iLAMgB1= 0.; iLAMgB1= 0.
endif
code such as :
! Graupel:
if (QG(i,k)>epsQ) then
iQG = 1./QG(i,k)
iNG = 1./NG(i,k)
iLAMg = max( iLAMmin1, iLAMDA_x(DE(i,k),QG(i,k),iNG,iGG99,thrd) )
iLAMg2 = iLAMg**2
iLAMgB0= iLAMg**(bfg)
iLAMgB1= iLAMg**(bfg+1.)
iLAMgB2= iLAMg**(bfg+2.)
!No_g = (NG(i,k))*iGG31/iLAMg**(1.+alpha_g)
No_g= NG(i,k)*iGG31/iLAMg !optimized for alpha_g=0
vg0 = gamfact(i,k)*ckQg1*iLAMgB0
Dg = Dm_x(DE(i,k),QG(i,k),iNG,icmg,thrd)
else
iLAMg = 0.; vg0 = 0.; Dg = 0.; No_g = 0.
iLAMg2 = 0.; iLAMgB0= 0.; iLAMgB1= 0.; iLAMgB1= 0.
endif