Git Product home page Git Product logo

uclchem's Introduction

UCLCHEM

UCLCHEM is a gas-grain chemical code that propagates the abundances of chemical species through a network of user-defined reactions according to the physical conditions of the gas. We provide several physical models to enable the modelling of different astrophysical environments and a utility script MakeRates to help the user produce a chemical network from simple lists of reactions and species.


Installation Instructions


Full documentation is available from the website: uclchem.github.io

UCLCHEM is intended to be used as a python module but must be installed from source rather than an online index such as Pypi. This is because users are expected to modify the source code, at least by creating their own networks. To obtain and install the code simply run:

[!IMPORTANT]
From version 3.3 onwards, we need to install in editable mode. This means that the directory that you store UCLCHEM in will directly be used by Python instead of being copied to the site-packages. Keep in mind that code changes you make in the uclchem directory will now affect you directly!

git clone https://github.com/uclchem/UCLCHEM.git
cd UCLCHEM
pip install -e .

You can then import uclchem in any python script. You need to pip install -e . whenever you change your network.

To see the contents of this python module, check our Python API docs. To see some example notebooks, check the tutorial section of the docs or the notebooks in Tutorials/.

If you want to build an executable from the Fortran source, head to src/fortran_src and run make. You can then run the executable with ./uclchem CLOUD input_file.inp where there examples of input files in the examples/ directory. We do not suggest users use the code this way unt

Prerequisites

To build UCLCHEM, you'll need gfortran, make and python 3.9+.

To run the python module, you'll need the python modules listed in requirements.txt


Change Log


See change.log! We've made a large number of improvements for v3.0. The code has been restructured to be Python first in its intended use, different physical models can be accessed without recompilation, and MakeRates is more helpful than ever before.


Contributing


This is an open source science code for the community and are open to pull requests. We are also happy to work with you to produce a physics module if none of the models available in the python module uclchem.model suit the modelling work you wish to do. If you are contributing, please try to work with our current code style. Feel free to checkout the latest developments with git fetch; git checkout develop We have the following general guidelines:

Github

  • Work in a personal branch or fork to your own Github to develop features.
  • Make sure you base your new work on the develop branch.
  • Pull requests should be opened with the develop branch as target.
  • In principle, squash and merge is preferred over keeping the entire git commit history when merging into develop.
  • Releases will be coordinated by the core developers, this will consitute a push to main and creating a release.

Python

  • Use Black to format your code.
  • snake_case variables and functions with self-explanatory names
  • Docstrings for all functions, they're used to produce the online docs!

Fortran

  • camelCase variable and subroutines names that are self-explanatory where possible
  • CAPITALIZED fortran built in functions to make code structure apparent.
  • Modularization, related subroutines should be added as modules. Small tweaks should be inserted into relevant module

Citing UCLCHEM


If you use UCLCHEM for your research, please cite the following paper (https://ui.adsabs.harvard.edu/abs/2017AJ....154...38H/abstract):

@ARTICLE{2017AJ....154...38H,
       author = {{Holdship}, J. and {Viti}, S. and {Jim{\'e}nez-Serra}, I. and {Makrymallis}, A. and {Priestley}, F.},
        title = "{UCLCHEM: A Gas-grain Chemical Code for Clouds, Cores, and C-Shocks}",
      journal = {\aj},
     keywords = {astrochemistry, ISM: molecules, shock waves, Astrophysics - Astrophysics of Galaxies},
         year = 2017,
        month = jul,
       volume = {154},
       number = {1},
          eid = {38},
        pages = {38},
          doi = {10.3847/1538-3881/aa773f},
archivePrefix = {arXiv},
       eprint = {1705.10677},
 primaryClass = {astro-ph.GA},
       adsurl = {https://ui.adsabs.harvard.edu/abs/2017AJ....154...38H},
      adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}

uclchem's People

Contributors

bamash avatar gijsvermarien avatar jonholdship avatar marcus-keil avatar ross-odonoghue avatar tomasjames avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar  avatar

uclchem's Issues

Incorrect production and destruction in analysis.py

Analysis.py provides clearly incorrect destruction rates, they never seem to overtake the production, so the "derivative" never becomes smaller than 0. See the files below. Above plot shows the abundances for a static model like in the test runs, the bottom plots show the following three lines:

  • The sign of the first order forward difference of the abundance (terrible approximation, but the sign should be correct)
  • The derivative retrieved with get_odes (directly from UCLCHEM wrap), this has the correct sign but a much different magnitude than analysis.
  • The deriatives obtained by subtracting the destruction from production from analysis.py

This is done for three species, the deriatives for non abundant species are more terrible (fluctuate around zero), so I didnt even bother with those. Plots below are on UCLCHEM v3.2.0 release candidate.

Image

Below a snippet of the text output of analysis.py in UCLCHEM v3.1.0

All Reactions
************************
CH3OH + CRPHOT -> H2CO + H2
CH3OH + CRPHOT -> OH + CH3
#CH3OH + DESCR -> CH3OH
#CH3OH + DESOH2 -> CH3OH
#CH3OH + DEUVCR -> CH3OH
CH3OH + FREEZE -> #CH3OH
#H + #CH2OH + LHDES -> CH3OH
#H + #H3CO + LHDES -> CH3OH
#H3CO + #HNCO + LHDES -> CH3OH + OCN
@H + @CH2OH + LHDES -> CH3OH
@H + @H3CO + LHDES -> CH3OH
@H3CO + @HNCO + LHDES -> CH3OH + OCN
CH3OH + PHOTON -> H2CO + H2
CH3OH + PHOTON -> H3CO+ + H + E-
CH3OH + PHOTON -> OH + CH3
#CH3OH + THERM -> CH3OH
@CH3OH + THERM -> CH3OH
C+ + CH3OH -> H3CO+ + CH
C+ + CH3OH -> HCO + CH3+
CH + CH3OH -> CH3CHO + H
CH+ + CH3OH -> CH3OH2+ + C
CH+ + CH3OH -> H2CO + CH3+
CH+ + CH3OH -> H3CO+ + CH2
CH3+ + CH3OH -> H3CO+ + CH4
CH3OH + S2+ -> H2S2+ + H2CO
CH3OH2+ + E- -> CH3OH + H
CH3OH2+ + NH3 -> CH3OH + NH4+
CH4+ + CH3OH -> CH3OH2+ + CH3
H+ + CH3OH -> CH3+ + H2O
H+ + CH3OH -> H3CO+ + H2
H+ + CH3OH -> HCO+ + H2 + H2
H2CO+ + CH3OH -> CH3OH2+ + HCO
H3+ + CH3CHO -> CH3+ + CH3OH
H3+ + CH3OH -> CH3+ + H2O + H2
H3+ + CH3OH -> CH3OH2+ + H2
H3CO+ + CH3OH -> CH3OH2+ + H2CO
H3O+ + CH3OH -> CH3OH2+ + H2O
HCO+ + CH3OH -> CH3OH2+ + CO
HE+ + CH3OH -> OH+ + CH3 + HE
HE+ + CH3OH -> OH + CH3+ + HE
N+ + CH3OH -> H2CO+ + NH + H
N+ + CH3OH -> H3CO+ + NH
N+ + CH3OH -> NO+ + CH3 + H
N+ + CH3OH -> NO + CH3+ + H
O+ + CH3OH -> H2CO+ + H2O
O+ + CH3OH -> H3CO+ + OH
O2+ + CH3OH -> O2 + H3CO+ + H
SI+ + CH3OH -> SIOH+ + CH3


***************************
New Important Reactions At: 0.00e+00 years
Formation = 3.62e-37 from:
CH3OH2+ + E- -> CH3OH + H : 100.00%

Destruction = -2.54e-38 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e-07 years
Formation = 2.80e-32 from:
#H + #H3CO + LHDES -> CH3OH : 51.74%
#H + #CH2OH + LHDES -> CH3OH : 48.26%

Destruction = -2.54e-38 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e-05 years
Formation = 2.78e-30 from:
#H + #H3CO + LHDES -> CH3OH : 51.74%
#H + #CH2OH + LHDES -> CH3OH : 48.26%

Destruction = -2.54e-38 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e-03 years
Formation = 1.68e-24 from:
#H + #H3CO + LHDES -> CH3OH : 99.98%

Destruction = -2.24e-36 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e-02 years
Formation = 1.29e-20 from:
#H + #H3CO + LHDES -> CH3OH : 99.98%

Destruction = -2.60e-31 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e-01 years
Formation = 4.03e-19 from:
#H + #H3CO + LHDES -> CH3OH : 99.98%

Destruction = -1.49e-28 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e+00 years
Formation = 1.44e-17 from:
#H + #H3CO + LHDES -> CH3OH : 99.98%

Destruction = -3.45e-26 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e+01 years
Formation = 8.43e-16 from:
#H + #H3CO + LHDES -> CH3OH : 99.98%

Destruction = -6.82e-24 from:
C+ + CH3OH -> HCO + CH3+ : 79.49%
C+ + CH3OH -> H3CO+ + CH : 19.87%

***************************
New Important Reactions At: 1.00e+02 years
Formation = 4.99e-14 from:
#H + #H3CO + LHDES -> CH3OH : 99.98%

Destruction = -5.42e-22 from:
C+ + CH3OH -> HCO + CH3+ : 79.34%
C+ + CH3OH -> H3CO+ + CH : 19.83%

***************************
New Important Reactions At: 1.00e+03 years
Formation = 8.90e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.97%

Destruction = -1.60e-20 from:
C+ + CH3OH -> HCO + CH3+ : 77.43%
C+ + CH3OH -> H3CO+ + CH : 19.36%
SI+ + CH3OH -> SIOH+ + CH3 : 2.76%

***************************
New Important Reactions At: 1.00e+04 years
Formation = 2.93e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.93%

Destruction = -8.02e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.54%
H+ + CH3OH -> H3CO+ + H2 : 7.50%
C+ + CH3OH -> HCO + CH3+ : 6.58%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.73%
C+ + CH3OH -> H3CO+ + CH : 1.64%
H+ + CH3OH -> CH3+ + H2O : 1.15%

***************************
New Important Reactions At: 1.11e+04 years
Formation = 3.04e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.92%

Destruction = -8.30e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 78.92%
H+ + CH3OH -> H3CO+ + H2 : 8.27%
C+ + CH3OH -> HCO + CH3+ : 6.87%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.91%
C+ + CH3OH -> H3CO+ + CH : 1.72%
H+ + CH3OH -> CH3+ + H2O : 1.27%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.33%

***************************
New Important Reactions At: 1.41e+04 years
Formation = 3.33e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.92%

Destruction = -9.19e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 74.63%
H+ + CH3OH -> H3CO+ + H2 : 9.79%
C+ + CH3OH -> HCO + CH3+ : 8.22%
H+ + CH3OH -> HCO+ + H2 + H2 : 2.26%
C+ + CH3OH -> H3CO+ + CH : 2.06%
H+ + CH3OH -> CH3+ + H2O : 1.50%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.47%
HE+ + CH3OH -> OH + CH3+ + HE : 0.47%

***************************
New Important Reactions At: 1.91e+04 years
Formation = 3.82e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.92%

Destruction = -1.10e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 72.69%
H+ + CH3OH -> H3CO+ + H2 : 9.91%
C+ + CH3OH -> HCO + CH3+ : 9.19%
C+ + CH3OH -> H3CO+ + CH : 2.30%
H+ + CH3OH -> HCO+ + H2 + H2 : 2.28%
H+ + CH3OH -> CH3+ + H2O : 1.52%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.59%
HE+ + CH3OH -> OH + CH3+ + HE : 0.59%

***************************
New Important Reactions At: 2.11e+04 years
Formation = 4.01e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.91%

Destruction = -1.18e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 72.90%
H+ + CH3OH -> H3CO+ + H2 : 9.55%
C+ + CH3OH -> HCO + CH3+ : 9.30%
C+ + CH3OH -> H3CO+ + CH : 2.32%
H+ + CH3OH -> HCO+ + H2 + H2 : 2.20%
H+ + CH3OH -> CH3+ + H2O : 1.47%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.61%
HE+ + CH3OH -> OH + CH3+ + HE : 0.61%
CH3OH + FREEZE -> #CH3OH : 0.44%

***************************
New Important Reactions At: 2.31e+04 years
Formation = 4.19e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.91%

Destruction = -1.25e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 73.32%
C+ + CH3OH -> HCO + CH3+ : 9.33%
H+ + CH3OH -> H3CO+ + H2 : 9.12%
C+ + CH3OH -> H3CO+ + CH : 2.33%
H+ + CH3OH -> HCO+ + H2 + H2 : 2.10%
H+ + CH3OH -> CH3+ + H2O : 1.40%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.63%
HE+ + CH3OH -> OH + CH3+ + HE : 0.63%
CH3OH + FREEZE -> #CH3OH : 0.47%

***************************
New Important Reactions At: 3.81e+04 years
Formation = 5.50e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.91%

Destruction = -1.75e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 77.30%
C+ + CH3OH -> HCO + CH3+ : 9.04%
H+ + CH3OH -> H3CO+ + H2 : 6.05%
C+ + CH3OH -> H3CO+ + CH : 2.26%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.39%
H+ + CH3OH -> CH3+ + H2O : 0.93%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.67%
HE+ + CH3OH -> OH + CH3+ + HE : 0.67%
CH3OH + FREEZE -> #CH3OH : 0.67%
CH3OH + CRPHOT -> H2CO + H2 : 0.41%

***************************
New Important Reactions At: 3.91e+04 years
Formation = 5.58e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.91%

Destruction = -1.79e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 77.51%
C+ + CH3OH -> HCO + CH3+ : 9.02%
H+ + CH3OH -> H3CO+ + H2 : 5.88%
C+ + CH3OH -> H3CO+ + CH : 2.25%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.36%
H+ + CH3OH -> CH3+ + H2O : 0.90%
CH3OH + FREEZE -> #CH3OH : 0.68%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.67%
HE+ + CH3OH -> OH + CH3+ + HE : 0.67%
CH3OH + CRPHOT -> H2CO + H2 : 0.42%

***************************
New Important Reactions At: 4.61e+04 years
Formation = 6.14e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.00e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 78.75%
C+ + CH3OH -> HCO + CH3+ : 8.92%
H+ + CH3OH -> H3CO+ + H2 : 4.89%
C+ + CH3OH -> H3CO+ + CH : 2.23%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.13%
CH3OH + FREEZE -> #CH3OH : 0.76%
H+ + CH3OH -> CH3+ + H2O : 0.75%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.68%
HE+ + CH3OH -> OH + CH3+ + HE : 0.68%
CH3OH + CRPHOT -> H2CO + H2 : 0.47%

***************************
New Important Reactions At: 5.01e+04 years
Formation = 6.44e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.11e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 79.27%
C+ + CH3OH -> HCO + CH3+ : 8.90%
H+ + CH3OH -> H3CO+ + H2 : 4.43%
C+ + CH3OH -> H3CO+ + CH : 2.23%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.02%
CH3OH + FREEZE -> #CH3OH : 0.80%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.69%
HE+ + CH3OH -> OH + CH3+ + HE : 0.69%
H+ + CH3OH -> CH3+ + H2O : 0.68%
CH3OH + CRPHOT -> H2CO + H2 : 0.49%

***************************
New Important Reactions At: 5.71e+04 years
Formation = 6.96e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.31e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 79.91%
C+ + CH3OH -> HCO + CH3+ : 8.93%
H+ + CH3OH -> H3CO+ + H2 : 3.75%
C+ + CH3OH -> H3CO+ + CH : 2.23%
CH3OH + FREEZE -> #CH3OH : 0.87%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.86%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.70%
HE+ + CH3OH -> OH + CH3+ + HE : 0.70%
H+ + CH3OH -> CH3+ + H2O : 0.58%
CH3OH + CRPHOT -> H2CO + H2 : 0.54%

***************************
New Important Reactions At: 6.01e+04 years
Formation = 7.17e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.39e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.09%
C+ + CH3OH -> HCO + CH3+ : 8.96%
H+ + CH3OH -> H3CO+ + H2 : 3.51%
C+ + CH3OH -> H3CO+ + CH : 2.24%
CH3OH + FREEZE -> #CH3OH : 0.90%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.81%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.70%
HE+ + CH3OH -> OH + CH3+ + HE : 0.70%
CH3OH + CRPHOT -> H2CO + H2 : 0.56%
H3+ + CH3OH -> CH3+ + H2O + H2 : 0.54%

***************************
New Important Reactions At: 6.31e+04 years
Formation = 7.38e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.47e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.23%
C+ + CH3OH -> HCO + CH3+ : 8.99%
H+ + CH3OH -> H3CO+ + H2 : 3.29%
C+ + CH3OH -> H3CO+ + CH : 2.25%
CH3OH + FREEZE -> #CH3OH : 0.93%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.76%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.71%
HE+ + CH3OH -> OH + CH3+ + HE : 0.71%
H3+ + CH3OH -> CH3+ + H2O + H2 : 0.58%
CH3OH + CRPHOT -> H2CO + H2 : 0.58%

***************************
New Important Reactions At: 6.61e+04 years
Formation = 7.59e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.55e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.32%
C+ + CH3OH -> HCO + CH3+ : 9.04%
H+ + CH3OH -> H3CO+ + H2 : 3.09%
C+ + CH3OH -> H3CO+ + CH : 2.26%
CH3OH + FREEZE -> #CH3OH : 0.96%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.72%
HE+ + CH3OH -> OH + CH3+ + HE : 0.72%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.71%
H3+ + CH3OH -> CH3+ + H2O + H2 : 0.62%
CH3OH + CRPHOT -> H2CO + H2 : 0.59%

***************************
New Important Reactions At: 7.01e+04 years
Formation = 7.86e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.65e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.39%
C+ + CH3OH -> HCO + CH3+ : 9.12%
H+ + CH3OH -> H3CO+ + H2 : 2.85%
C+ + CH3OH -> H3CO+ + CH : 2.28%
CH3OH + FREEZE -> #CH3OH : 1.00%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.72%
HE+ + CH3OH -> OH + CH3+ + HE : 0.72%
H3+ + CH3OH -> CH3+ + H2O + H2 : 0.67%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.66%
CH3OH + CRPHOT -> H2CO + H2 : 0.62%

***************************
New Important Reactions At: 7.31e+04 years
Formation = 8.05e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.73e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.40%
C+ + CH3OH -> HCO + CH3+ : 9.19%
H+ + CH3OH -> H3CO+ + H2 : 2.69%
C+ + CH3OH -> H3CO+ + CH : 2.30%
CH3OH + FREEZE -> #CH3OH : 1.03%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.73%
HE+ + CH3OH -> OH + CH3+ + HE : 0.73%
H3+ + CH3OH -> CH3+ + H2O + H2 : 0.71%
CH3OH + CRPHOT -> H2CO + H2 : 0.64%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.62%

***************************
New Important Reactions At: 7.51e+04 years
Formation = 8.18e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.78e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.40%
C+ + CH3OH -> HCO + CH3+ : 9.24%
H+ + CH3OH -> H3CO+ + H2 : 2.59%
C+ + CH3OH -> H3CO+ + CH : 2.31%
CH3OH + FREEZE -> #CH3OH : 1.05%
H3+ + CH3OH -> CH3+ + H2O + H2 : 0.74%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.74%
HE+ + CH3OH -> OH + CH3+ + HE : 0.74%
CH3OH + CRPHOT -> H2CO + H2 : 0.65%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.60%

***************************
New Important Reactions At: 8.11e+04 years
Formation = 8.56e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.90%

Destruction = -2.92e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 80.30%
C+ + CH3OH -> HCO + CH3+ : 9.40%
C+ + CH3OH -> H3CO+ + CH : 2.35%
H+ + CH3OH -> H3CO+ + H2 : 2.32%
CH3OH + FREEZE -> #CH3OH : 1.11%
H3+ + CH3OH -> CH3+ + H2O + H2 : 0.83%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.75%
HE+ + CH3OH -> OH + CH3+ + HE : 0.75%
CH3OH + CRPHOT -> H2CO + H2 : 0.69%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.54%

***************************
New Important Reactions At: 9.41e+04 years
Formation = 9.32e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.89%

Destruction = -3.22e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 79.81%
C+ + CH3OH -> HCO + CH3+ : 9.82%
C+ + CH3OH -> H3CO+ + CH : 2.46%
H+ + CH3OH -> H3CO+ + H2 : 1.89%
CH3OH + FREEZE -> #CH3OH : 1.23%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.03%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.78%
HE+ + CH3OH -> OH + CH3+ + HE : 0.78%
CH3OH + CRPHOT -> H2CO + H2 : 0.76%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.43%
CH3OH + CRPHOT -> OH + CH3 : 0.36%

***************************
New Important Reactions At: 1.10e+05 years
Formation = 1.02e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.89%

Destruction = -3.55e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 78.87%
C+ + CH3OH -> HCO + CH3+ : 10.41%
C+ + CH3OH -> H3CO+ + CH : 2.60%
H+ + CH3OH -> H3CO+ + H2 : 1.54%
CH3OH + FREEZE -> #CH3OH : 1.36%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.30%
CH3OH + CRPHOT -> H2CO + H2 : 0.84%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.82%
HE+ + CH3OH -> OH + CH3+ + HE : 0.82%
CH3OH + CRPHOT -> OH + CH3 : 0.40%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.35%

***************************
New Important Reactions At: 1.20e+05 years
Formation = 1.06e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.89%

Destruction = -3.75e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 78.09%
C+ + CH3OH -> HCO + CH3+ : 10.84%
C+ + CH3OH -> H3CO+ + CH : 2.71%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.47%
CH3OH + FREEZE -> #CH3OH : 1.44%
H+ + CH3OH -> H3CO+ + H2 : 1.40%
CH3OH + CRPHOT -> H2CO + H2 : 0.89%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.84%
HE+ + CH3OH -> OH + CH3+ + HE : 0.84%
CH3OH + CRPHOT -> OH + CH3 : 0.42%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.33%

***************************
New Important Reactions At: 1.60e+05 years
Formation = 1.21e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.88%

Destruction = -4.47e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 74.13%
C+ + CH3OH -> HCO + CH3+ : 12.91%
C+ + CH3OH -> H3CO+ + CH : 3.23%
H3+ + CH3OH -> CH3+ + H2O + H2 : 2.22%
CH3OH + FREEZE -> #CH3OH : 1.67%
H+ + CH3OH -> H3CO+ + H2 : 1.17%
CH3OH + CRPHOT -> H2CO + H2 : 1.03%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.93%
HE+ + CH3OH -> OH + CH3+ + HE : 0.93%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.50%
CH3OH + CRPHOT -> OH + CH3 : 0.49%

***************************
New Important Reactions At: 1.90e+05 years
Formation = 1.27e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.87%

Destruction = -4.96e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 70.38%
C+ + CH3OH -> HCO + CH3+ : 14.72%
C+ + CH3OH -> H3CO+ + CH : 3.68%
H3+ + CH3OH -> CH3+ + H2O + H2 : 2.87%
CH3OH + FREEZE -> #CH3OH : 1.82%
H+ + CH3OH -> H3CO+ + H2 : 1.21%
CH3OH + CRPHOT -> H2CO + H2 : 1.12%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.00%
HE+ + CH3OH -> OH + CH3+ + HE : 1.00%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.65%
CH3OH + CRPHOT -> OH + CH3 : 0.53%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.28%

***************************
New Important Reactions At: 2.00e+05 years
Formation = 1.29e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.87%

Destruction = -5.12e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 69.00%
C+ + CH3OH -> HCO + CH3+ : 15.35%
C+ + CH3OH -> H3CO+ + CH : 3.84%
H3+ + CH3OH -> CH3+ + H2O + H2 : 3.12%
CH3OH + FREEZE -> #CH3OH : 1.87%
H+ + CH3OH -> H3CO+ + H2 : 1.25%
CH3OH + CRPHOT -> H2CO + H2 : 1.15%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.02%
HE+ + CH3OH -> OH + CH3+ + HE : 1.02%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.71%
CH3OH + CRPHOT -> OH + CH3 : 0.55%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.29%

***************************
New Important Reactions At: 2.30e+05 years
Formation = 1.32e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.86%

Destruction = -5.56e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 64.53%
C+ + CH3OH -> HCO + CH3+ : 17.24%
C+ + CH3OH -> H3CO+ + CH : 4.31%
H3+ + CH3OH -> CH3+ + H2O + H2 : 3.94%
CH3OH + FREEZE -> #CH3OH : 2.03%
H+ + CH3OH -> H3CO+ + H2 : 1.43%
CH3OH + CRPHOT -> H2CO + H2 : 1.26%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.12%
HE+ + CH3OH -> OH + CH3+ + HE : 1.12%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.89%
CH3OH + CRPHOT -> OH + CH3 : 0.60%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.49%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.33%

***************************
New Important Reactions At: 2.40e+05 years
Formation = 1.33e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.86%

Destruction = -5.70e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 62.95%
C+ + CH3OH -> HCO + CH3+ : 17.85%
C+ + CH3OH -> H3CO+ + CH : 4.46%
H3+ + CH3OH -> CH3+ + H2O + H2 : 4.25%
CH3OH + FREEZE -> #CH3OH : 2.10%
H+ + CH3OH -> H3CO+ + H2 : 1.50%
CH3OH + CRPHOT -> H2CO + H2 : 1.30%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.15%
HE+ + CH3OH -> OH + CH3+ + HE : 1.15%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.96%
CH3OH + CRPHOT -> OH + CH3 : 0.62%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.58%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.35%

***************************
New Important Reactions At: 2.50e+05 years
Formation = 1.02e-11 from:
#H + #H3CO + LHDES -> CH3OH : 99.83%

Destruction = -5.71e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 61.51%
C+ + CH3OH -> HCO + CH3+ : 18.42%
C+ + CH3OH -> H3CO+ + CH : 4.60%
H3+ + CH3OH -> CH3+ + H2O + H2 : 4.52%
CH3OH + FREEZE -> #CH3OH : 2.14%
H+ + CH3OH -> H3CO+ + H2 : 1.57%
CH3OH + CRPHOT -> H2CO + H2 : 1.32%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.18%
HE+ + CH3OH -> OH + CH3+ + HE : 1.18%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.02%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.68%
CH3OH + CRPHOT -> OH + CH3 : 0.63%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.39%

***************************
New Important Reactions At: 2.80e+05 years
Formation = 2.55e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.81%

Destruction = -2.83e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 66.37%
C+ + CH3OH -> HCO + CH3+ : 17.76%
C+ + CH3OH -> H3CO+ + CH : 4.44%
H3+ + CH3OH -> CH3+ + H2O + H2 : 3.33%
H+ + CH3OH -> H3CO+ + H2 : 1.44%
CH3OH + FREEZE -> #CH3OH : 1.39%
CH3OH + CRPHOT -> H2CO + H2 : 0.86%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.80%
HE+ + CH3OH -> OH + CH3+ + HE : 0.80%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.75%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.62%
CH3OH + CRPHOT -> OH + CH3 : 0.41%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.33%

***************************
New Important Reactions At: 2.90e+05 years
Formation = 1.63e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.82%

Destruction = -1.84e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 68.29%
C+ + CH3OH -> HCO + CH3+ : 17.72%
C+ + CH3OH -> H3CO+ + CH : 4.43%
H3+ + CH3OH -> CH3+ + H2O + H2 : 2.74%
H+ + CH3OH -> H3CO+ + H2 : 1.38%
CH3OH + FREEZE -> #CH3OH : 1.10%
CH3OH + CRPHOT -> H2CO + H2 : 0.68%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.65%
HE+ + CH3OH -> OH + CH3+ + HE : 0.65%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.62%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.52%
CH3OH + CRPHOT -> OH + CH3 : 0.32%

***************************
New Important Reactions At: 3.00e+05 years
Formation = 1.07e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.83%

Destruction = -1.15e-19 from:
SI+ + CH3OH -> SIOH+ + CH3 : 69.11%
C+ + CH3OH -> HCO + CH3+ : 18.22%
C+ + CH3OH -> H3CO+ + CH : 4.55%
H3+ + CH3OH -> CH3+ + H2O + H2 : 2.26%
H+ + CH3OH -> H3CO+ + H2 : 1.36%
CH3OH + FREEZE -> #CH3OH : 0.88%
CH3OH + CRPHOT -> H2CO + H2 : 0.55%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.53%
HE+ + CH3OH -> OH + CH3+ + HE : 0.53%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.51%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.43%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.31%

***************************
New Important Reactions At: 3.10e+05 years
Formation = 7.17e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.83%

Destruction = -7.33e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 68.76%
C+ + CH3OH -> HCO + CH3+ : 19.29%
C+ + CH3OH -> H3CO+ + CH : 4.82%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.91%
H+ + CH3OH -> H3CO+ + H2 : 1.38%
CH3OH + FREEZE -> #CH3OH : 0.73%
CH3OH + CRPHOT -> H2CO + H2 : 0.45%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.45%
HE+ + CH3OH -> OH + CH3+ + HE : 0.45%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.43%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.36%

***************************
New Important Reactions At: 3.20e+05 years
Formation = 4.98e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.84%

Destruction = -4.86e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 67.38%
C+ + CH3OH -> HCO + CH3+ : 20.85%
C+ + CH3OH -> H3CO+ + CH : 5.21%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.68%
H+ + CH3OH -> H3CO+ + H2 : 1.45%
CH3OH + FREEZE -> #CH3OH : 0.63%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.39%
HE+ + CH3OH -> OH + CH3+ + HE : 0.39%
CH3OH + CRPHOT -> H2CO + H2 : 0.39%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.38%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.33%

***************************
New Important Reactions At: 3.30e+05 years
Formation = 3.59e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.85%

Destruction = -3.37e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 65.13%
C+ + CH3OH -> HCO + CH3+ : 22.82%
C+ + CH3OH -> H3CO+ + CH : 5.71%
H+ + CH3OH -> H3CO+ + H2 : 1.56%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.55%
CH3OH + FREEZE -> #CH3OH : 0.57%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.36%
HE+ + CH3OH -> OH + CH3+ + HE : 0.36%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.36%
CH3OH + CRPHOT -> H2CO + H2 : 0.35%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.35%

***************************
New Important Reactions At: 3.40e+05 years
Formation = 2.70e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.86%

Destruction = -2.45e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 62.14%
C+ + CH3OH -> HCO + CH3+ : 25.17%
C+ + CH3OH -> H3CO+ + CH : 6.29%
H+ + CH3OH -> H3CO+ + H2 : 1.71%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.50%
CH3OH + FREEZE -> #CH3OH : 0.55%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.39%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.35%
HE+ + CH3OH -> OH + CH3+ + HE : 0.35%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.34%
CH3OH + CRPHOT -> H2CO + H2 : 0.34%

***************************
New Important Reactions At: 3.70e+05 years
Formation = 1.50e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.88%

Destruction = -1.25e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 50.06%
C+ + CH3OH -> HCO + CH3+ : 33.88%
C+ + CH3OH -> H3CO+ + CH : 8.47%
H+ + CH3OH -> H3CO+ + H2 : 2.44%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.60%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.56%
CH3OH + FREEZE -> #CH3OH : 0.55%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.38%
HE+ + CH3OH -> OH + CH3+ + HE : 0.38%
H+ + CH3OH -> CH3+ + H2O : 0.37%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.36%

***************************
New Important Reactions At: 3.80e+05 years
Formation = 1.33e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.88%

Destruction = -1.08e-20 from:
SI+ + CH3OH -> SIOH+ + CH3 : 45.49%
C+ + CH3OH -> HCO + CH3+ : 37.05%
C+ + CH3OH -> H3CO+ + CH : 9.26%
H+ + CH3OH -> H3CO+ + H2 : 2.75%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.69%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.63%
CH3OH + FREEZE -> #CH3OH : 0.57%
H+ + CH3OH -> CH3+ + H2O : 0.42%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.40%
HE+ + CH3OH -> OH + CH3+ + HE : 0.40%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.38%

***************************
New Important Reactions At: 3.90e+05 years
Formation = 1.22e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.88%

Destruction = -9.75e-21 from:
SI+ + CH3OH -> SIOH+ + CH3 : 40.90%
C+ + CH3OH -> HCO + CH3+ : 40.18%
C+ + CH3OH -> H3CO+ + CH : 10.04%
H+ + CH3OH -> H3CO+ + H2 : 3.08%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.78%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.71%
CH3OH + FREEZE -> #CH3OH : 0.59%
H+ + CH3OH -> CH3+ + H2O : 0.47%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.42%
HE+ + CH3OH -> OH + CH3+ + HE : 0.42%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.40%
CH3OH + CRPHOT -> H2CO + H2 : 0.37%

***************************
New Important Reactions At: 4.00e+05 years
Formation = 1.15e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.89%

Destruction = -9.07e-21 from:
C+ + CH3OH -> HCO + CH3+ : 43.19%
SI+ + CH3OH -> SIOH+ + CH3 : 36.41%
C+ + CH3OH -> H3CO+ + CH : 10.80%
H+ + CH3OH -> H3CO+ + H2 : 3.42%
H3+ + CH3OH -> CH3+ + H2O + H2 : 1.89%
H+ + CH3OH -> HCO+ + H2 + H2 : 0.79%
CH3OH + FREEZE -> #CH3OH : 0.62%
H+ + CH3OH -> CH3+ + H2O : 0.53%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.44%
HE+ + CH3OH -> OH + CH3+ + HE : 0.44%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.43%
CH3OH + CRPHOT -> H2CO + H2 : 0.38%

***************************
New Important Reactions At: 4.40e+05 years
Formation = 1.09e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.89%

Destruction = -8.50e-21 from:
C+ + CH3OH -> HCO + CH3+ : 52.96%
SI+ + CH3OH -> SIOH+ + CH3 : 21.22%
C+ + CH3OH -> H3CO+ + CH : 13.24%
H+ + CH3OH -> H3CO+ + H2 : 4.77%
H3+ + CH3OH -> CH3+ + H2O + H2 : 2.37%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.10%
H+ + CH3OH -> CH3+ + H2O : 0.73%
CH3OH + FREEZE -> #CH3OH : 0.72%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.54%
HE+ + CH3OH -> OH + CH3+ + HE : 0.54%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.54%
CH3OH + CRPHOT -> H2CO + H2 : 0.45%

***************************
New Important Reactions At: 4.70e+05 years
Formation = 1.16e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.88%

Destruction = -9.28e-21 from:
C+ + CH3OH -> HCO + CH3+ : 57.46%
C+ + CH3OH -> H3CO+ + CH : 14.36%
SI+ + CH3OH -> SIOH+ + CH3 : 13.52%
H+ + CH3OH -> H3CO+ + H2 : 5.61%
H3+ + CH3OH -> CH3+ + H2O + H2 : 2.78%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.29%
H+ + CH3OH -> CH3+ + H2O : 0.86%
CH3OH + FREEZE -> #CH3OH : 0.80%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.63%
HE+ + CH3OH -> OH + CH3+ + HE : 0.63%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.63%
CH3OH + CRPHOT -> H2CO + H2 : 0.49%

***************************
New Important Reactions At: 4.80e+05 years
Formation = 1.20e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.88%

Destruction = -9.67e-21 from:
C+ + CH3OH -> HCO + CH3+ : 58.48%
C+ + CH3OH -> H3CO+ + CH : 14.62%
SI+ + CH3OH -> SIOH+ + CH3 : 11.60%
H+ + CH3OH -> H3CO+ + H2 : 5.85%
H3+ + CH3OH -> CH3+ + H2O + H2 : 2.92%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.35%
H+ + CH3OH -> CH3+ + H2O : 0.90%
CH3OH + FREEZE -> #CH3OH : 0.82%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.66%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.66%
HE+ + CH3OH -> OH + CH3+ + HE : 0.66%
CH3OH + CRPHOT -> H2CO + H2 : 0.51%

***************************
New Important Reactions At: 4.90e+05 years
Formation = 1.24e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.87%

Destruction = -1.01e-20 from:
C+ + CH3OH -> HCO + CH3+ : 59.30%
C+ + CH3OH -> H3CO+ + CH : 14.83%
SI+ + CH3OH -> SIOH+ + CH3 : 9.95%
H+ + CH3OH -> H3CO+ + H2 : 6.06%
H3+ + CH3OH -> CH3+ + H2O + H2 : 3.07%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.40%
H+ + CH3OH -> CH3+ + H2O : 0.93%
CH3OH + FREEZE -> #CH3OH : 0.85%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.70%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.69%
HE+ + CH3OH -> OH + CH3+ + HE : 0.69%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.52%
CH3OH + CRPHOT -> H2CO + H2 : 0.52%

***************************
New Important Reactions At: 5.20e+05 years
Formation = 1.38e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.87%

Destruction = -1.18e-20 from:
C+ + CH3OH -> HCO + CH3+ : 60.79%
C+ + CH3OH -> H3CO+ + CH : 15.20%
H+ + CH3OH -> H3CO+ + H2 : 6.57%
SI+ + CH3OH -> SIOH+ + CH3 : 6.30%
H3+ + CH3OH -> CH3+ + H2O + H2 : 3.56%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.51%
H+ + CH3OH -> CH3+ + H2O : 1.01%
CH3OH + FREEZE -> #CH3OH : 0.93%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.81%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.79%
HE+ + CH3OH -> OH + CH3+ + HE : 0.79%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.62%
CH3OH + CRPHOT -> H2CO + H2 : 0.57%

***************************
New Important Reactions At: 5.50e+05 years
Formation = 1.54e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.85%

Destruction = -1.38e-20 from:
C+ + CH3OH -> HCO + CH3+ : 61.21%
C+ + CH3OH -> H3CO+ + CH : 15.30%
H+ + CH3OH -> H3CO+ + H2 : 6.93%
H3+ + CH3OH -> CH3+ + H2O + H2 : 4.13%
SI+ + CH3OH -> SIOH+ + CH3 : 4.06%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.60%
H+ + CH3OH -> CH3+ + H2O : 1.06%
CH3OH + FREEZE -> #CH3OH : 1.01%
H3+ + CH3OH -> CH3OH2+ + H2 : 0.93%
HE+ + CH3OH -> OH+ + CH3 + HE : 0.89%
HE+ + CH3OH -> OH + CH3+ + HE : 0.89%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.73%
CH3OH + CRPHOT -> H2CO + H2 : 0.62%

***************************
New Important Reactions At: 5.90e+05 years
Formation = 1.79e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.84%

Destruction = -1.69e-20 from:
C+ + CH3OH -> HCO + CH3+ : 60.75%
C+ + CH3OH -> H3CO+ + CH : 15.19%
H+ + CH3OH -> H3CO+ + H2 : 7.24%
H3+ + CH3OH -> CH3+ + H2O + H2 : 4.99%
SI+ + CH3OH -> SIOH+ + CH3 : 2.34%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.67%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.13%
CH3OH + FREEZE -> #CH3OH : 1.13%
H+ + CH3OH -> CH3+ + H2O : 1.11%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.06%
HE+ + CH3OH -> OH + CH3+ + HE : 1.06%
HCO+ + CH3OH -> CH3OH2+ + CO : 0.92%
CH3OH + CRPHOT -> H2CO + H2 : 0.70%

***************************
New Important Reactions At: 6.10e+05 years
Formation = 1.92e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.83%

Destruction = -1.87e-20 from:
C+ + CH3OH -> HCO + CH3+ : 60.27%
C+ + CH3OH -> H3CO+ + CH : 15.07%
H+ + CH3OH -> H3CO+ + H2 : 7.34%
H3+ + CH3OH -> CH3+ + H2O + H2 : 5.47%
SI+ + CH3OH -> SIOH+ + CH3 : 1.80%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.69%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.24%
CH3OH + FREEZE -> #CH3OH : 1.18%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.14%
HE+ + CH3OH -> OH + CH3+ + HE : 1.14%
H+ + CH3OH -> CH3+ + H2O : 1.13%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.03%
CH3OH + CRPHOT -> H2CO + H2 : 0.73%

***************************
New Important Reactions At: 6.20e+05 years
Formation = 1.99e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.82%

Destruction = -1.97e-20 from:
C+ + CH3OH -> HCO + CH3+ : 59.99%
C+ + CH3OH -> H3CO+ + CH : 15.00%
H+ + CH3OH -> H3CO+ + H2 : 7.37%
H3+ + CH3OH -> CH3+ + H2O + H2 : 5.71%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.70%
SI+ + CH3OH -> SIOH+ + CH3 : 1.59%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.29%
CH3OH + FREEZE -> #CH3OH : 1.21%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.19%
HE+ + CH3OH -> OH + CH3+ + HE : 1.19%
H+ + CH3OH -> CH3+ + H2O : 1.13%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.08%
CH3OH + CRPHOT -> H2CO + H2 : 0.75%

***************************
New Important Reactions At: 6.30e+05 years
Formation = 2.06e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.82%

Destruction = -2.07e-20 from:
C+ + CH3OH -> HCO + CH3+ : 59.69%
C+ + CH3OH -> H3CO+ + CH : 14.92%
H+ + CH3OH -> H3CO+ + H2 : 7.40%
H3+ + CH3OH -> CH3+ + H2O + H2 : 5.96%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.71%
SI+ + CH3OH -> SIOH+ + CH3 : 1.40%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.35%
CH3OH + FREEZE -> #CH3OH : 1.24%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.23%
HE+ + CH3OH -> OH + CH3+ + HE : 1.23%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.14%
H+ + CH3OH -> CH3+ + H2O : 1.14%
CH3OH + CRPHOT -> H2CO + H2 : 0.77%

***************************
New Important Reactions At: 6.40e+05 years
Formation = 2.13e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.81%

Destruction = -2.17e-20 from:
C+ + CH3OH -> HCO + CH3+ : 59.39%
C+ + CH3OH -> H3CO+ + CH : 14.85%
H+ + CH3OH -> H3CO+ + H2 : 7.42%
H3+ + CH3OH -> CH3+ + H2O + H2 : 6.20%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.71%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.40%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.28%
HE+ + CH3OH -> OH + CH3+ + HE : 1.28%
CH3OH + FREEZE -> #CH3OH : 1.26%
SI+ + CH3OH -> SIOH+ + CH3 : 1.24%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.20%
H+ + CH3OH -> CH3+ + H2O : 1.14%
CH3OH + CRPHOT -> H2CO + H2 : 0.78%

***************************
New Important Reactions At: 6.50e+05 years
Formation = 2.20e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.81%

Destruction = -2.28e-20 from:
C+ + CH3OH -> HCO + CH3+ : 59.08%
C+ + CH3OH -> H3CO+ + CH : 14.77%
H+ + CH3OH -> H3CO+ + H2 : 7.44%
H3+ + CH3OH -> CH3+ + H2O + H2 : 6.44%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.71%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.46%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.32%
HE+ + CH3OH -> OH + CH3+ + HE : 1.32%
CH3OH + FREEZE -> #CH3OH : 1.29%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.26%
H+ + CH3OH -> CH3+ + H2O : 1.14%
SI+ + CH3OH -> SIOH+ + CH3 : 1.10%
CH3OH + CRPHOT -> H2CO + H2 : 0.80%

***************************
New Important Reactions At: 6.60e+05 years
Formation = 2.27e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.80%

Destruction = -2.39e-20 from:
C+ + CH3OH -> HCO + CH3+ : 58.76%
C+ + CH3OH -> H3CO+ + CH : 14.69%
H+ + CH3OH -> H3CO+ + H2 : 7.45%
H3+ + CH3OH -> CH3+ + H2O + H2 : 6.68%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.72%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.51%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.36%
HE+ + CH3OH -> OH + CH3+ + HE : 1.36%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.31%
CH3OH + FREEZE -> #CH3OH : 1.31%
H+ + CH3OH -> CH3+ + H2O : 1.14%
SI+ + CH3OH -> SIOH+ + CH3 : 0.98%
CH3OH + CRPHOT -> H2CO + H2 : 0.81%

***************************
New Important Reactions At: 6.80e+05 years
Formation = 2.42e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.79%

Destruction = -2.63e-20 from:
C+ + CH3OH -> HCO + CH3+ : 58.13%
C+ + CH3OH -> H3CO+ + CH : 14.53%
H+ + CH3OH -> H3CO+ + H2 : 7.45%
H3+ + CH3OH -> CH3+ + H2O + H2 : 7.16%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.72%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.62%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.45%
HE+ + CH3OH -> OH + CH3+ + HE : 1.45%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.43%
CH3OH + FREEZE -> #CH3OH : 1.35%
H+ + CH3OH -> CH3+ + H2O : 1.15%
CH3OH + CRPHOT -> H2CO + H2 : 0.84%
SI+ + CH3OH -> SIOH+ + CH3 : 0.78%

***************************
New Important Reactions At: 7.00e+05 years
Formation = 2.58e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.78%

Destruction = -2.88e-20 from:
C+ + CH3OH -> HCO + CH3+ : 57.50%
C+ + CH3OH -> H3CO+ + CH : 14.37%
H3+ + CH3OH -> CH3+ + H2O + H2 : 7.62%
H+ + CH3OH -> H3CO+ + H2 : 7.44%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.73%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.72%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.54%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.53%
HE+ + CH3OH -> OH + CH3+ + HE : 1.53%
CH3OH + FREEZE -> #CH3OH : 1.39%
H+ + CH3OH -> CH3+ + H2O : 1.14%
CH3OH + CRPHOT -> H2CO + H2 : 0.86%
SI+ + CH3OH -> SIOH+ + CH3 : 0.63%

***************************
New Important Reactions At: 7.10e+05 years
Formation = 2.66e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.78%

Destruction = -3.02e-20 from:
C+ + CH3OH -> HCO + CH3+ : 57.19%
C+ + CH3OH -> H3CO+ + CH : 14.30%
H3+ + CH3OH -> CH3+ + H2O + H2 : 7.85%
H+ + CH3OH -> H3CO+ + H2 : 7.43%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.78%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.71%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.59%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.57%
HE+ + CH3OH -> OH + CH3+ + HE : 1.57%
CH3OH + FREEZE -> #CH3OH : 1.41%
H+ + CH3OH -> CH3+ + H2O : 1.14%
CH3OH + CRPHOT -> H2CO + H2 : 0.87%
SI+ + CH3OH -> SIOH+ + CH3 : 0.57%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.43%

***************************
New Important Reactions At: 7.40e+05 years
Formation = 2.91e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.76%

Destruction = -3.45e-20 from:
C+ + CH3OH -> HCO + CH3+ : 56.26%
C+ + CH3OH -> H3CO+ + CH : 14.07%
H3+ + CH3OH -> CH3+ + H2O + H2 : 8.53%
H+ + CH3OH -> H3CO+ + H2 : 7.40%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.93%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.74%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.70%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.69%
HE+ + CH3OH -> OH + CH3+ + HE : 1.69%
CH3OH + FREEZE -> #CH3OH : 1.46%
H+ + CH3OH -> CH3+ + H2O : 1.14%
CH3OH + CRPHOT -> H2CO + H2 : 0.90%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.47%
CH3OH + CRPHOT -> OH + CH3 : 0.43%

***************************
New Important Reactions At: 7.50e+05 years
Formation = 2.99e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.75%

Destruction = -3.60e-20 from:
C+ + CH3OH -> HCO + CH3+ : 55.96%
C+ + CH3OH -> H3CO+ + CH : 13.99%
H3+ + CH3OH -> CH3+ + H2O + H2 : 8.75%
H+ + CH3OH -> H3CO+ + H2 : 7.38%
H3+ + CH3OH -> CH3OH2+ + H2 : 1.98%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.79%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.73%
HE+ + CH3OH -> OH + CH3+ + HE : 1.73%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.70%
CH3OH + FREEZE -> #CH3OH : 1.47%
H+ + CH3OH -> CH3+ + H2O : 1.13%
CH3OH + CRPHOT -> H2CO + H2 : 0.91%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.49%
CH3OH + CRPHOT -> OH + CH3 : 0.43%

***************************
New Important Reactions At: 7.60e+05 years
Formation = 3.08e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.75%

Destruction = -3.76e-20 from:
C+ + CH3OH -> HCO + CH3+ : 55.65%
C+ + CH3OH -> H3CO+ + CH : 13.91%
H3+ + CH3OH -> CH3+ + H2O + H2 : 8.97%
H+ + CH3OH -> H3CO+ + H2 : 7.36%
H3+ + CH3OH -> CH3OH2+ + H2 : 2.03%
HCO+ + CH3OH -> CH3OH2+ + CO : 1.84%
HE+ + CH3OH -> OH+ + CH3 + HE : 1.76%
HE+ + CH3OH -> OH + CH3+ + HE : 1.76%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.70%
CH3OH + FREEZE -> #CH3OH : 1.48%
H+ + CH3OH -> CH3+ + H2O : 1.13%
CH3OH + CRPHOT -> H2CO + H2 : 0.92%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.50%

***************************
New Important Reactions At: 8.90e+05 years
Formation = 4.31e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.66%

Destruction = -6.29e-20 from:
C+ + CH3OH -> HCO + CH3+ : 51.61%
C+ + CH3OH -> H3CO+ + CH : 12.90%
H3+ + CH3OH -> CH3+ + H2O + H2 : 11.87%
H+ + CH3OH -> H3CO+ + H2 : 7.08%
H3+ + CH3OH -> CH3OH2+ + H2 : 2.69%
HCO+ + CH3OH -> CH3OH2+ + CO : 2.39%
HE+ + CH3OH -> OH+ + CH3 + HE : 2.29%
HE+ + CH3OH -> OH + CH3+ + HE : 2.29%
CH3OH + FREEZE -> #CH3OH : 1.64%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.63%
H+ + CH3OH -> CH3+ + H2O : 1.09%
CH3OH + CRPHOT -> H2CO + H2 : 1.01%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.66%

***************************
New Important Reactions At: 9.30e+05 years
Formation = 4.72e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.63%

Destruction = -7.24e-20 from:
C+ + CH3OH -> HCO + CH3+ : 50.34%
H3+ + CH3OH -> CH3+ + H2O + H2 : 12.79%
C+ + CH3OH -> H3CO+ + CH : 12.58%
H+ + CH3OH -> H3CO+ + H2 : 6.98%
H3+ + CH3OH -> CH3OH2+ + H2 : 2.90%
HCO+ + CH3OH -> CH3OH2+ + CO : 2.54%
HE+ + CH3OH -> OH+ + CH3 + HE : 2.45%
HE+ + CH3OH -> OH + CH3+ + HE : 2.45%
CH3OH + FREEZE -> #CH3OH : 1.68%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.61%
H+ + CH3OH -> CH3+ + H2O : 1.07%
CH3OH + CRPHOT -> H2CO + H2 : 1.04%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.71%

***************************
New Important Reactions At: 9.70e+05 years
Formation = 5.13e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.60%

Destruction = -8.28e-20 from:
C+ + CH3OH -> HCO + CH3+ : 49.06%
H3+ + CH3OH -> CH3+ + H2O + H2 : 13.72%
C+ + CH3OH -> H3CO+ + CH : 12.26%
H+ + CH3OH -> H3CO+ + H2 : 6.88%
H3+ + CH3OH -> CH3OH2+ + H2 : 3.11%
HCO+ + CH3OH -> CH3OH2+ + CO : 2.68%
HE+ + CH3OH -> OH+ + CH3 + HE : 2.62%
HE+ + CH3OH -> OH + CH3+ + HE : 2.62%
CH3OH + FREEZE -> #CH3OH : 1.73%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.59%
CH3OH + CRPHOT -> H2CO + H2 : 1.07%
H+ + CH3OH -> CH3+ + H2O : 1.06%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.75%

***************************
New Important Reactions At: 1.10e+06 years
Formation = 6.49e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.50%

Destruction = -1.21e-19 from:
C+ + CH3OH -> HCO + CH3+ : 44.90%
H3+ + CH3OH -> CH3+ + H2O + H2 : 16.75%
C+ + CH3OH -> H3CO+ + CH : 11.22%
H+ + CH3OH -> H3CO+ + H2 : 6.55%
H3+ + CH3OH -> CH3OH2+ + H2 : 3.79%
HE+ + CH3OH -> OH+ + CH3 + HE : 3.17%
HE+ + CH3OH -> OH + CH3+ + HE : 3.17%
HCO+ + CH3OH -> CH3OH2+ + CO : 3.09%
CH3OH + FREEZE -> #CH3OH : 1.86%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.51%
CH3OH + CRPHOT -> H2CO + H2 : 1.15%
H+ + CH3OH -> CH3+ + H2O : 1.01%
H3O+ + CH3OH -> CH3OH2+ + H2O : 0.90%

***************************
New Important Reactions At: 1.20e+06 years
Formation = 7.50e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.42%

Destruction = -1.55e-19 from:
C+ + CH3OH -> HCO + CH3+ : 41.81%
H3+ + CH3OH -> CH3+ + H2O + H2 : 19.03%
C+ + CH3OH -> H3CO+ + CH : 10.45%
H+ + CH3OH -> H3CO+ + H2 : 6.29%
H3+ + CH3OH -> CH3OH2+ + H2 : 4.31%
HE+ + CH3OH -> OH+ + CH3 + HE : 3.58%
HE+ + CH3OH -> OH + CH3+ + HE : 3.58%
HCO+ + CH3OH -> CH3OH2+ + CO : 3.35%
CH3OH + FREEZE -> #CH3OH : 1.95%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.45%
CH3OH + CRPHOT -> H2CO + H2 : 1.21%
H3O+ + CH3OH -> CH3OH2+ + H2O : 1.02%
H+ + CH3OH -> CH3+ + H2O : 0.97%

***************************
New Important Reactions At: 1.30e+06 years
Formation = 8.43e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.33%

Destruction = -1.90e-19 from:
C+ + CH3OH -> HCO + CH3+ : 38.92%
H3+ + CH3OH -> CH3+ + H2O + H2 : 21.18%
C+ + CH3OH -> H3CO+ + CH : 9.73%
H+ + CH3OH -> H3CO+ + H2 : 6.04%
H3+ + CH3OH -> CH3OH2+ + H2 : 4.80%
HE+ + CH3OH -> OH+ + CH3 + HE : 3.96%
HE+ + CH3OH -> OH + CH3+ + HE : 3.96%
HCO+ + CH3OH -> CH3OH2+ + CO : 3.57%
CH3OH + FREEZE -> #CH3OH : 2.04%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.39%
CH3OH + CRPHOT -> H2CO + H2 : 1.26%
H3O+ + CH3OH -> CH3OH2+ + H2O : 1.14%
H+ + CH3OH -> CH3+ + H2O : 0.93%
CH3OH + CRPHOT -> OH + CH3 : 0.60%

***************************
New Important Reactions At: 1.50e+06 years
Formation = 9.98e-13 from:
#H + #H3CO + LHDES -> CH3OH : 99.18%

Destruction = -2.59e-19 from:
C+ + CH3OH -> HCO + CH3+ : 33.99%
H3+ + CH3OH -> CH3+ + H2O + H2 : 24.89%
C+ + CH3OH -> H3CO+ + CH : 8.50%
H3+ + CH3OH -> CH3OH2+ + H2 : 5.64%
H+ + CH3OH -> H3CO+ + H2 : 5.60%
HE+ + CH3OH -> OH+ + CH3 + HE : 4.61%
HE+ + CH3OH -> OH + CH3+ + HE : 4.61%
HCO+ + CH3OH -> CH3OH2+ + CO : 3.91%
CH3OH + FREEZE -> #CH3OH : 2.18%
H3O+ + CH3OH -> CH3OH2+ + H2O : 1.37%
CH3OH + CRPHOT -> H2CO + H2 : 1.35%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.29%
H+ + CH3OH -> CH3+ + H2O : 0.86%
CH3OH + CRPHOT -> OH + CH3 : 0.64%

***************************
New Important Reactions At: 1.80e+06 years
Formation = 1.15e-12 from:
#H + #H3CO + LHDES -> CH3OH : 99.01%

Destruction = -3.43e-19 from:
C+ + CH3OH -> HCO + CH3+ : 28.97%
H3+ + CH3OH -> CH3+ + H2O + H2 : 28.72%
C+ + CH3OH -> H3CO+ + CH : 7.24%
H3+ + CH3OH -> CH3OH2+ + H2 : 6.50%
HE+ + CH3OH -> OH+ + CH3 + HE : 5.24%
HE+ + CH3OH -> OH + CH3+ + HE : 5.24%
H+ + CH3OH -> H3CO+ + H2 : 5.13%
HCO+ + CH3OH -> CH3OH2+ + CO : 4.20%
CH3OH + FREEZE -> #CH3OH : 2.32%
H3O+ + CH3OH -> CH3OH2+ + H2O : 1.66%
CH3OH + CRPHOT -> H2CO + H2 : 1.43%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.18%
H+ + CH3OH -> CH3+ + H2O : 0.79%
CH3OH + CRPHOT -> OH + CH3 : 0.68%

***************************
New Important Reactions At: 1.90e+06 years
Formation = 1.19e-12 from:
#H + #H3CO + LHDES -> CH3OH : 98.96%
#H + #CH2OH + LHDES -> CH3OH : 1.04%

Destruction = -3.65e-19 from:
H3+ + CH3OH -> CH3+ + H2O + H2 : 29.59%
C+ + CH3OH -> HCO + CH3+ : 27.84%
C+ + CH3OH -> H3CO+ + CH : 6.96%
H3+ + CH3OH -> CH3OH2+ + H2 : 6.70%
HE+ + CH3OH -> OH+ + CH3 + HE : 5.38%
HE+ + CH3OH -> OH + CH3+ + HE : 5.38%
H+ + CH3OH -> H3CO+ + H2 : 5.01%
HCO+ + CH3OH -> CH3OH2+ + CO : 4.25%
CH3OH + FREEZE -> #CH3OH : 2.35%
H3O+ + CH3OH -> CH3OH2+ + H2O : 1.73%
CH3OH + CRPHOT -> H2CO + H2 : 1.45%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.16%
H+ + CH3OH -> CH3+ + H2O : 0.77%
CH3OH + CRPHOT -> OH + CH3 : 0.69%

***************************
New Important Reactions At: 2.00e+06 years
Formation = 1.22e-12 from:
#H + #H3CO + LHDES -> CH3OH : 98.92%
#H + #CH2OH + LHDES -> CH3OH : 1.08%

Destruction = -3.83e-19 from:
H3+ + CH3OH -> CH3+ + H2O + H2 : 30.29%
C+ + CH3OH -> HCO + CH3+ : 26.93%
H3+ + CH3OH -> CH3OH2+ + H2 : 6.86%
C+ + CH3OH -> H3CO+ + CH : 6.73%
HE+ + CH3OH -> OH+ + CH3 + HE : 5.49%
HE+ + CH3OH -> OH + CH3+ + HE : 5.49%
H+ + CH3OH -> H3CO+ + H2 : 4.92%
HCO+ + CH3OH -> CH3OH2+ + CO : 4.30%
CH3OH + FREEZE -> #CH3OH : 2.38%
H3O+ + CH3OH -> CH3OH2+ + H2O : 1.80%
CH3OH + CRPHOT -> H2CO + H2 : 1.47%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.13%
H+ + CH3OH -> CH3+ + H2O : 0.76%
CH3OH + CRPHOT -> OH + CH3 : 0.70%

***************************
New Important Reactions At: 2.60e+06 years
Formation = 1.29e-12 from:
#H + #H3CO + LHDES -> CH3OH : 98.81%
#H + #CH2OH + LHDES -> CH3OH : 1.19%

Destruction = -4.39e-19 from:
H3+ + CH3OH -> CH3+ + H2O + H2 : 32.30%
C+ + CH3OH -> HCO + CH3+ : 24.31%
H3+ + CH3OH -> CH3OH2+ + H2 : 7.31%
C+ + CH3OH -> H3CO+ + CH : 6.08%
HE+ + CH3OH -> OH+ + CH3 + HE : 5.79%
HE+ + CH3OH -> OH + CH3+ + HE : 5.79%
H+ + CH3OH -> H3CO+ + H2 : 4.65%
HCO+ + CH3OH -> CH3OH2+ + CO : 4.42%
CH3OH + FREEZE -> #CH3OH : 2.45%
H3O+ + CH3OH -> CH3OH2+ + H2O : 2.01%
CH3OH + CRPHOT -> H2CO + H2 : 1.52%
H+ + CH3OH -> HCO+ + H2 + H2 : 1.07%
CH3OH + CRPHOT -> OH + CH3 : 0.72%
H+ + CH3OH -> CH3+ + H2O : 0.71%

Setup.py calls

Currently the code is compiled 3 times instead of just one, causing the installer to run much slower. Confirmed with the following snippet added to setup.py

from datetime import datetime
with open("log.txt", 'a') as logfile:
    logfile.write(str(datetime.now()) + '\n')

Which creates the file (only running once):

  14:36:54 with uclchem_osx as gijsv in ~/uclchem/UCLCHEM at main(!?)
πŸ‘¨πŸ»β€πŸ’» >  cat log.txt
2023-09-01 14:34:40.932011
2023-09-01 14:35:16.858450
2023-09-01 14:35:50.135776

Models have occassional final time step truncation

Occassional models, when run with grid.py, have their final time step truncated such that the abundance of some molecules is not written out.

This is leading to errors when plotting these abundances as a function of time, as the array dimensions do not match (i.e. python is able to find the time step for that entry, but not the molecular abundance as it has been truncated).

Move to main and dev

We still use master which is outdated, this should be updated to main. With this change we should do the following:

  • Move to main as the default branch, ensure it has all the release tags in it
  • Create a dev branch, which will be the primary developer branch
  • Update the documentation to reflect the improved workflow for both users and developers.
    • documentation
    • README.md
    • ...?

Precision

We're using high precision real variables for most arrays.
However, is there an issue with ODE calculations where a high rate of change term results in later terms being rounded off as the are added?

Makerates species masses

Added a function find_constituents() to functions.py that checks mass and warns if expected mass (sum of elements) is different to input mass in species file.

Several improvements to be made:

  1. Add warning for unknown species - currently get a python error
  2. Move function to the species class
  3. expand elements list to include a few more elements (Fe) and isotopes (13C, 18O)

Increasing memory footprint

Current version of UCLCHEM footprint increases when the python wrapper is used in a script to perform multiple models in a row, which is the default scenario for our grids: https://uclchem.github.io/docs/running_a_grid.

It was determined that this was probably due to the interplay of Fortran and f2py. Since we cannot really address it directly, the best solution is to make sure workers are restarted every so often. This can be done using the joblib workers, i.e. with Parallel and delayed:

from joblib import Parallel, delayed
results = Parallel(n_jobs=6)(
        delayed(pool_func)(row) for row in model_table.iterrows()
    )

Fix comparenetworks.py

At some points we added user_defined_bulk to the return of read_species_file, so we cannot directly pipe its input into read_reaction_file. Currently this breaks comparenetworks.py.

Add mac M1 installation instructions

There currently are new installation instructions for the M1 macs, we should provide these in the documentation so people can use their new shiny macbooks without needing Rosetta 2.

Large exponents overwrite exponent sign

When trying to analyze some of the files I ran into a value that was marked as "non-numeric" by pandas:

0          0.00000E+00
1          3.05698-105
2          4.50159E-67
3          7.77345E-60

It appears that the large exponent overwrites "E". Do we want to do anything about this, maybe just format the exponent with one more padded 0 just in case?

Vdiff indexing

Most recent update to network.f90 broke the indexing of the vdiff array in the thermal desorption branch.

Vdiff has length nspecies but we only need it for surface species.

Currently, if a grain species is at position i in the grain array and j in the species array, it's diffusion rate is placed at position i in the vdiff array.

Then, in rates.f90, it is accessed at position j so we get the wrong value.

Move from setup.py

With pip 22.3.1 we get the following deprecation warning, announcing the deprecation of setup.py. We need to find some solution that still allows us to compile the fortran as we do in setup.py

DEPRECATION: uclchem is being installed using the legacy 'setup.py install' method, because it does not have a 'pyproject.toml' and the 'wheel' package is not installed. pip 23.1 will enforce this behaviour change. A possible replacement is to enable the '--use-pep517' option. Discussion can be found at https://github.com/pypa/pip/issues/8559`

Makerates Auxiliary Scripts Improvement

We have several scripts (comparenetworks.py, comparespecies.py and speciesnetwork.py) that help analyze reaction networks and check for differences between networks. They're fairly basic and could do with updating.

Comparenetworks.py should read two reaction networks (either reactions.csv output of makerates or any valid input file) and print out all reactions that are only in one network or have different rate coefficients between networks.

Comparespecies.py might not be useful, it just lists differences between two lists of species.

Speciesnetwork.py should list all formation and destruction routes for a single species.

Contributor should make these more robust and user-friendly, add missing functionality (comparenetworks does not currently check coefficients) and consider whether other scripts would be useful.

Grain variables

The grain surface area per hydrogen atom (GRAIN_AREA=2.4d-22) is precalculated

but we have everything we need to do this consistently by setting average grain radius (GRAIN_RADIUS) and calculating surface area of a dust grain from that. Then dividing byGAS_DUST_DENSITY_RATIO which is 1/n_d (the fractional abunance of dust grains).

Without it, our grain treatment is not entirely consistent.

fail to recompile with toggle "three_phase" off

Hi,

After installing, UCLCHEM with everything default works fine.

I want to exclude bulk ice chemistry for starters. Setting the three_phase toggle to False in usersettings.yaml and executing MakeRates.py spits some warnings, but works. Then trying to recompile I get an error. Compiling from src/fortran_src as recommended in the docs, the error looks like this:


(uclcenv) CN09:/vol/thchem/jsalaris/UCLCHEM/src/fortran_src $ make python
gfortran -O3 -fPIC -ffree-line-length-0 -c chemistry.f90
chemistry.f90:176:49:

176 | abund(nBulk,dstep)=sum(abund(bulkList,dstep))
| 1
Error: Symbol β€˜bulklist’ at (1) has no IMPLICIT type
make: *** [Makefile:40: chemistry.o] Error 1

Any idea why I get this error and how I can fix it?

I am using:

uclchem-3.2.0
python-3.9.16
pip 22.2.2
numpy 1.22.3

Move to result code 0 on success

Currently we use 1 with negative error codes, I think it would be best to move to 0 as the success code to be in line with common practices.

Deprecation of pip temporary directories

We need to address the following deprecation warning:

DEPRECATION: A future pip version will change local packages to be built in-place without first copying to a temporary directory. We recommend you use --use-feature=in-tree-build to test your packages with this new behavior before it becomes the default.
   pip 21.3 will remove support for this functionality. You can find discussion regarding this at https://github.com/pypa/pip/issues/7555.

failed to install in Ubuntu

I was trying to install UCLCHEM in my Ubuntu 23.04. After I created the environment, installed all the requirements, and "pip install -e .", there was a long error which was ended by:

      AttributeError: module 'distutils' has no attribute 'core'
      make: *** [Makefile:37: python] Error 1
      [end of output]
  
  note: This error originates from a subprocess, and is likely not a problem with pip.
  ERROR: Failed building editable for uclchem
Failed to build uclchem
ERROR: Could not build wheels for uclchem, which is required to install pyproject.toml-based projects

Besides, in the Docs of UCLCHEM, it is required to "conda install gfortran". But this package seems not available, and I have had installed gfortran with apt on my computer. Does it matter?

Deprecate or migrate main.f90

There is still main.f90 as a legacy interface to run models, we should either move the command line interface to Python, or include main.f90 in the scope of code we maintain.

analysis.py

Analysis.py estimates the rate of change of a species wrt to each reaction in the network and then for each time step, lists the most important destruction and formation routes.

However, since it's not used often it hasn't been maintained. Several things need to be improved:

  1. User interface. There are almost no comments, not even to explain inputs
  2. Rate calculations. rates.f90 is effectively rewritten in python in plot_functions.py. It needs to be updated to account for new reaction rates OR better yet, can we make it use a compile python version of UCLCHEM?

Ionpol reactions

The KIDA network has two types of reactions that are not in UMIST and so are not coded for in UCLCHEM: ionpol1 and ionpol2. There's details in this document.

Adding a section to rates.f90 to deal with ionpol reactions would be trivial, it's just a case of letting UCLCHEM know which reactions are ionpols in a way that doesn't break the way MakeRates currently reads in UMIST files.

Rates Efficiency Improvements

This should be a fairly small update that will have major improvements. Probably of interest to @Marcus-Keil and Johannes if they're running UCLCHEM a lot.

We should use Makerates to massively optimize rates.f90. We can sort the reaction list by type and then store the indices of each reaction type. Then the rates can be calculated together by index rather than using the massive loop and case statement in rates.f90.

Eg make a list of two body gas reactions called gasTwoBodyReacs then calculate all gas rates together:
rates(gasTwoBodyReacs)=alpha(gasTwoBodyReacs)*((temp/300)**beta(gasTwoBodyReacs))*exp(-gamma(gasTwoBodyReacs)/temp)

That would be a massive speed up because array operations are much better optimized than loops. We'd have one list (or just start/end indices for sorted reactions) per reaction type and can calculate all of the reactions for a given type as one array operation.

It's worth branching the code, doing this and then testing the run time over 10 runs or so. I think we'd see a huge improvement.

Refractory Species

Currently every species is either in the gas phase or in the ices and all frozen species are treated the same way. However, certain species are refractory and will not sublimate at high temperatures.

It would be useful to introduce a flag that prevents a species from being removed from the dust through thermal sublimation.

One could then edit the shock modules so that we include both mantle and core sputtering using the expressions from Jimenez-Serra 2008. Currently, we assume mantle sputtering for all species and simply deplete abundances of core species.

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    πŸ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. πŸ“ŠπŸ“ˆπŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❀️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.