Astrochem issue #64

h2o.chm

#H2O          + uv-photon                   -> OH           + H            3.28e-10  0.00e+00  1.63e+00 13 4314
H2O                                        -> H2O(ice)                    1.00e+00  1.80e+01  0.00e+00 20 6238
H2O(ice)                                   -> H2O                         0.00e+00  1.80e+01  5.70e+03 21 6239
#H2O(ice)     + cosmic-ray                  -> H2O                         4.40e-17  0.00e+00  0.00e+00 22 6240

input.ini

[files]
source = source.mdl
chem = h2o.chm
# Physical parameters
[phys]
chi = 1.0
cosmic = 1.3e-17
# Solver parameters
[solver]
ti = 1e-6
tf = 1e7
# Initial abundances
[abundances]
H2 = 0.5
H2O(ice) = 2.5e-4
H2O = 2.5e-4
He = 1.4e-1
CN = 6.6e-8
HCN = 1.00e-8
N = 5.10e-6
NH3(ice) =  9.9e-6
N2 = 1.1e-6
H3(+) = 1e-8
CS = 4e-9
SO = 5e-9
Si(+) = 1e-11
S(+) = 1e-11
Mg(+) = 1e-11
Fe(+) = 1e-11
C(+) = 1e-9
CH4 = 1e-7
CO = 1e-6
C = 7e-7
HCO(+) = 9e-9
C2H = 8e-9
e(-) = 2.004e-8
grain = 4.25106352511e-11
# Output
[output]
abundances = H2O,H2O(ice)
trace_routes = 0

source.mdl

# Source model file example
# shell number, Av [mag], nH [cm^-3], Tgas [K], Tdust [K]
#
0	0.0	1e+12	1289	1289