DSMACC chemical box model
This page describes the DSMACC chemical box model (U. Leeds).
The Dynamically Simple Model of Atmospheric Chemical Complexity (DSMACC) is designed to be a simple box model which can easily change its chemical scheme and can be used both for free running and constrained simulations. It is written in FORTRAN and uses the KPP chemical pre-processor, the TUV photolysis scheme and can use the Master Chemical Mechanism or other chemical mechanisms.
DSMACC is a useful and flexible tool. We are using it to test new GEOS-Chem chemical mechanisms.
Instructions for Use
To use the DSMACC model with GEOS-Chem with the GEOS-Chem chemical mechanisms follow the instructions below after downloading the DSMACC code and KPP utilities from the DSMACC website:
- In the boxmodel directory type "make clean" to clear the directory.
- In the same directory type "Geos2kppbox_parser.pl globchem.dat" to generate the three kpp-compatible globchem.eqn, globchem.spc and globchem.def files. Ensure that the number of reaction lines in globchem.dat matches the number of lines read in with the perl script (under the section entitled "Evaluate Product Lines" in Geos2kppbox_parser.pl).
- Generate the kpp *.f90 files by typing "kpp model.kpp".
- In the tuv_old/ directory type "make clean" followed by "make". A warning message will appear ("no action performed for specified file(s)). This is anticipated.
- Go back to the boxmodel directory and compile the model by entering "make". You may get error messages about the incompatibility of the variable type with reaction rates of the order 1.0E-39. This should not affect the code from compiling.
- Run the code by entering "model". Or you can send the output straight to an output file by entering "model > model_output"
Authors and collaborators
Source code and data
The code is freely available although support is limited!
- Emmerson, K.M., M. J. Evans, Comparison of tropospheric gas-phase chemistry schemes for use within global models, Atmos. Chem. Phys., 9, 1831-1845, 2009. PDF
None at this time.