wetr change

[zywshoon]

Orginally, for mixed particle, we set the wet radius = dry radius when temperature< 190.
Now, we set the wet radius = wet radius(T=190) when temperature< 190.

wetr.F90

• subroutine getwetr(carma, igroup, rh, rdry, rwet, rhopdry, rhopwet, rc, h2o_mass, h2o_vp, temp, kappa)

• subroutine getwetr(carma, igroup, rh, rdry, rwet, rhopdry, rhopwet, rc, h2o_mass, h2o_vp, temp, kappa,cstate, iz)
• type(carmastate_type), intent(in),optional :: cstate !! the carma state object
• integer, intent(in),optional :: iz !! level index
• real(kind=f) :: rvap,gc_cgs,rwet190,rhopwet190,rh190,r_ratio190,pvap190h2o

• if (temp .le. 190._f) then

•   rwet = rdry
  

•   rhopwet = rhopdry
  

•   if (temp .le. 190._f) then
  

•     rvap = RGAS / gwtmol(igash2o)
  

•     gc_cgs = gc(iz,igash2o) / (zmet(iz)*xmet(iz)*ymet(iz))
  

• pvap190h2o = 10.0_f * exp(54.842763_f - (6763.22_f / 190._f) - (4.210_f * log(190._f)) + (0.000367_f * 190._f) + &

•     (tanh(0.0415_f * (190._f - 218.8_f)) * &
  

•     (53.878_f - (1331.22_f / 190._f) - (9.44523_f * log(190._f)) + 0.014025_f * 190._f)))
  

•     rh190 = gc_cgs * rvap * 190._f/pvap190h2o
  

•     rwet190 = rdry * (1._f + rh190*kappa/(1._f-rh190))**(1._f/3._f)
  

•     r_ratio190 = (rdry / rwet190)**3._f
  

•     rhopwet190 = r_ratio190 * rhopdry + (1._f - r_ratio190) * RHO_W
  

•     rwet = rwet190
  

•     rhopwet = rhopwet190
  

rhopart.F90

•           rhop(iz,ibin,igroup), rhop_wet(iz,ibin,igroup), rc, temp=t(iz), kappa=kappahygro(iz,ibin,igroup),cstate=cstate, iz=iz)
  

---

•           rhop(iz,ibin,igroup), rhop_wet(iz,ibin,igroup), rc, temp=t(iz), kappa=kappahygro(iz,ibin,igroup))
  

162c162

•           rhop(iz,ibin,igroup), rhop_wet(iz,ibin,igroup), rc, temp=t(iz), kappa=kappahygro(iz,ibin,igroup),cstate=cstate, iz=iz)
  

---

•           rhop(iz,ibin,igroup), rhop_wet(iz,ibin,igroup), rc, temp=t(iz), kappa=kappahygro(iz,ibin,igroup))
  

167c167

•           rhop(iz,ibin,igroup), rhop_wet(iz,ibin,igroup), rc, temp=t(iz), kappa=kappahygro(iz,ibin,igroup),cstate=cstate, iz=iz)
  

---

•           rhop(iz,ibin,igroup), rhop_wet(iz,ibin,igroup), rc, temp=t(iz), kappa=kappahygro(iz,ibin,igroup))