#include atoms_red
#include SAPRC99.spc
#include SAPRC99.eqn

#include ../boxmox/wrapper
#include ../boxmox/setup

#LOOKATALL

#INLINE F90_GLOBAL

  ! global variables used by special rate constants of this mechanism

  ! air/water concentration (molecules/cm3)
  REAL(KIND=dp)      :: c_m, c_h2o
  ! relative humidity
  REAL(KIND=dp)      :: rh

  REAL(KIND=dp)      :: rc_n2o5

#ENDINLINE

#INLINE F90_INIT

  ! 1E-03 * (1.0/28.97) * 6.022e23 * air density
  c_m    = 2.078702e+19 * 1.25

  ! 1E-03 * (1.0/18.02) * 6.022e23 * air density, then times kg(water)/kg(air)
  c_h2o  = 3.341842e+19 * 1.25 * 0.010

  ! relative humidity (not really used...)
  rh = 0.5

  ! n2o5 rate const
  rc_n2o5 = 1e-30

#ENDINLINE