chem/KPP/inc/radm2/kpp_mechd_ibu_radm2.inc

   1 !
   2 ! calculate rate constant for n2o5 + water in RADM2
   3
   4   es  = 1000.*0.6112*exp(17.67*(t_phy(i,k,j)-273.15)/(t_phy(i,k,j)- 29.65))
   5   qvs = es / ( p_phy(i,k,j) - es )
   6
   7
   8   rh =  moist(i,k,j,P_QV) / qvs
   9   rh = MIN ( MAX ( rh, 0.), 1.)
  10
  11
  12   rc_n2o5 =  REAL(1.0 / ( 3.6E4 * EXP( -( rh / 0.28 ) ** 2.8 ) + 300.0 ), KIND=dp)
  13
  14