|
|
| (229 intermediate revisions by 6 users not shown) |
| Line 1: |
Line 1: |
| − | On this page we provide information about the FlexChem chemical solver, which was introduced in [[GEOS-Chem v11-01#v11-01g|GEOS-Chem v11-01g]].
| + | This content has been migrated to: |
| | | | |
| − | == Overview ==
| + | * [https://kpp.readthedocs.io <tt>kpp.readthedocs.io</tt>]: User manual for The Kinetic Preprocessor (KPP) |
| − | | + | * [https://geos-chem.readthedocs.io/en/latest/geos-chem-shared-docs/supplemental-guides/using-kpp-with-gc.html '''Update chemical mechanisms with KPP''' at <tt>geos-chem.readthedocs.io</tt>] |
| − | The clean and flexible reimplementation of the Kinetic PreProcessor package (aka KPP)—known as '''FlexChem'''—is nearing full integration into GEOS-Chem. Most of the remaining FlexChem development work will focus on replacing legacy infrastructure—which prevents GEOS-Chem from operating efficiently in high-performance computing (HPC) environments—with newer, more efficient algorithms.
| + | |
| − | | + | |
| − | Check back soon for the latest information!
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 19:58, 19 May 2016 (UTC)
| + | |
| − | | + | |
| − | == Milestones ==
| + | |
| − | | + | |
| − | The following table shows several milestones that were achieved in the FlexChem implementation, development, as well as ongoing development tasks.
| + | |
| − | | + | |
| − | {| border=1 cellspacing=0 cellpadding=5
| + | |
| − | |-valign="top" bgcolor="#CCCCCC"
| + | |
| − | !width="600px"|Task
| + | |
| − | !width="150px"|Developer
| + | |
| − | !width="250px"|Status
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Add FlexChem into v11-01c
| + | |
| − | |Mike Long
| + | |
| − | |Completed 14 Dec 2015
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Enable the “tropchem” mechanism
| + | |
| − | |Mike Long
| + | |
| − | |Completed 14 Dec 2015
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Restore the OH and HO2 diagnostics (ND43)
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 18 Dec 2015
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove CSPEC array and replaced with State_Chm%Species
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 22 Dec 2015
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Enable a temporary workaround for family tracers (ISOPN, MMN)
| + | |
| − | |Mike Long
| + | |
| − | |Completed 25 Jan 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Enable FAST-JX photochemistry
| + | |
| − | |Mike Long
| + | |
| − | |Completed 25 Jan 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Enable the "benchmark" chemistry mechanism
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 29 Jan 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Fix HBr and HOBr hetchem rates for the "tropchem" mechanism
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 04 Feb 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Fix unit conversions for heterogeneous chemistry
| + | |
| − | |Mike Long
| + | |
| − | |Completed 18 Feb 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Restore the broken J-value diagnostic (ND22)
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 15 Mar 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Parallelize the main KPP driver loop; fixed other minor issues
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 30 Mar 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Enable the SOA and SOA-SVPOA mechanisms
| + | |
| − | |Lizzie Lundgren
| + | |
| − | |Completed 01 Apr 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Add more parallelization fixes
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 18 Apr 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Merge FlexChem with v11-01f
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 20 Apr 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Add KPP repository to Bitbucket (https://bitbucket.org/gcst/kpp)
| + | |
| − | |Mike Long
| + | |
| − | |Completed 19 Apr 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Create new gckpp* files from the updated KPP solver package
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 22 Apr 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Enable the UCX mechanism
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 26 Apr 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Introduce a new prod/loss diagnostic into the KPP solver package
| + | |
| − | |Mike Long
| + | |
| − | |Completed 29 Apr 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Fix various other bugs and cleaned up some leftover things
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 04 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Complete unit tests and 1-month benchmarks for the tropchem, benchmark, UCX, SOA, and SOA-SVPOA simulations
| + | |
| − | |Melissa Sulprizio<br>Lizzie Lundgren
| + | |
| − | |Completed 05 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove tracer indices of Rn, Pb, Be, Hg, and POPs simulations from <code>tracerid_mod.F</code>. (This paves the way for us to retire <code>tracerid_mod.F</code>.)
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 02 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Add a fast species name lookup algorithm
| + | |
| − | |Bob Yantosca
| + | |
| − | |Code is in place as of 04 May 2016, but has not been implemented throughout GEOS-Chem yet.
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Store the unique list of GEOS-Chem species in the GEOS-Chem species database object. (In other words, we combine the list of advected tracers with the list of KPP chemical species and remove duplicate entries).
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 09 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Store the indices of each species in the KPP chemical reaction matrix in the GEOS-Chem species database object.
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 09 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Merge FlexChem with other v11-01g updates
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 17 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove family tracer fields of the <code>Input_Opt</code> object, namely : <code>TRACER_N_CONST, TRACER_CONST, TRACER_COEFF, ID_EMITTED</code>. These were only used by SMVGEAR and are now obsolete.
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 18 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Create mapping vectors in the <code>State_Chm</code> object to store the ID numbers of species that are advected, dry-deposited, wet-deposited, and/or included in the KPP chemical mechanism
| + | |
| − | |Bob Yantosca
| + | |
| − | |Completed 19 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove all 1-D loop indexing variables (<code>JLOOP, KLOOP, KTLOOP, JLOP, IXSAVE, IYSAVE, IZSAVE</code>) that were used by SMVGEAR
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 19 May 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove SMVGEAR input files (<code>globchem.dat, mglob.dat</code>) and associated subroutines (<code>readchem.F, reader.F</code>)
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |Completed 30 Jun 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove <code>tracerid_mod.F</code> and related tracer ID flags (<code>IDTxxxx</code>, <code>IDxxxx</code>) from GEOS-Chem.
| + | |
| − | |GCST
| + | |
| − | |
| + | |
| − | *Completed 24 Jun 2016
| + | |
| − | *Validated 27 Jun 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Use the GEOS-Chem species database to retrieve species index from species name
| + | |
| − | |GCST
| + | |
| − | |
| + | |
| − | *Completed 24 Jun 2016
| + | |
| − | *Validated 27 Jun 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove tracers from restart files and read/write restart file species only
| + | |
| − | |Lizzie Lundgren
| + | |
| − | |
| + | |
| − | *Completed 12 Jul 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Create unit conversion routines for species that mimic existing routines for tracers or convert to/from molecules/cm3
| + | |
| − | |Lizzie Lundgren
| + | |
| − | |
| + | |
| − | *Completed 25 Jul 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Add State_Chm variable to track species units throughout GEOS-Chem (State_Chm%Spc_Units)
| + | |
| − | |Lizzie Lundgren
| + | |
| − | |
| + | |
| − | *Completed 25 Jul 2016 | + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove the remaining family tracers (ISOPN, MMN, CFCX, HCFCX). Henceforth, FlexChem will only work with individual species.
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |
| + | |
| − | *Completed 28 Jul 2016
| + | |
| − | *Validated 01 Aug 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Move the non-advected species initial unit conversion of background values from init_flexchem (called during chemistry) to read_gc_restart (called during restart file read)
| + | |
| − | |Lizzie Lundgren
| + | |
| − | |
| + | |
| − | *Completed 4 Aug 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Store species background concentrations values previously stored in globchem.dat in the species database
| + | |
| − | |Lizzie Lundgren
| + | |
| − | |
| + | |
| − | *Completed 11 Aug 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Attach the new KPP prod/loss output to the GEOS-Chem prod/loss diagnostic
| + | |
| − | |Melissa Sulprizio
| + | |
| − | |
| + | |
| − | *Completed 11 Aug 2016
| + | |
| − | *Validated 11 Aug 2016
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Remove <code>State_Chm%Tracers</code> and use species concentrations in <code>State_Chm%Species</code> instead
| + | |
| − | |GCST
| + | |
| − | |In progress
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Update the GEOS-Chem Makefiles so that a fresh version of KPP will be built (using a custom chemical mechanism that you specify) each time GEOS-Chem is compiled.
| + | |
| − | |GCST
| + | |
| − | |TBD
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Validate FlexChem with additional debugging, unit testing, and benchmarking
| + | |
| − | |GCST
| + | |
| − | |TBD
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Update GEOS-Chem documentation and user manual for v11-01
| + | |
| − | |GCST
| + | |
| − | |TBD
| + | |
| − | | + | |
| − | |}
| + | |
| − | | + | |
| − | --[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 21:41, 11 August 2016 (UTC)
| + | |
| − | | + | |
| − | == Adding chemical mechanisms in FlexChem ==
| + | |
| − | | + | |
| − | === KPP source code ===
| + | |
| − | | + | |
| − | The KPP source code for FlexChem is currently available as a Git repository on Bitbucket. It can be obtained using:
| + | |
| − | | + | |
| − | git clone https://bitbucket.org/gcst/kpp
| + | |
| − | | + | |
| − | Once you have KPP on your system, you will need to build it:
| + | |
| − | | + | |
| − | cd kpp/kpp-2.2.3_01
| + | |
| − | make clean
| + | |
| − | make
| + | |
| − | | + | |
| − | If the build completed successfully, you should see the executable <tt>kpp</tt> in the <tt>kpp/kpp-2.2.3_01/bin/</tt> directory. Next, add KPP to your PATH. For bash users, add the following lines the <tt>~/.bashrc</tt> file:
| + | |
| − | | + | |
| − | export PATH=$PATH:/PATH_TO_KPP/KPP/kpp-2.2.3_01/bin/
| + | |
| − | export KPP_HOME=/PATH_TO_KPP/KPP/kpp-2.2.3_01
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:18, 10 August 2016 (UTC) | + | |
| − | | + | |
| − | === Building a custom chemical mechanism ===
| + | |
| − | | + | |
| − | [[GEOS-Chem v11-01#v11-01g|GEOS-Chem v11-01g]] has five pre-built chemical mechanisms: standard, tropchem, SOA, SOA-SVPOA, and UCX. These mechanisms are described on our [[GEOS-Chem_chemistry_mechanisms#Mechanisms_in_GEOS-Chem_v10-01_and_later_versions|''GEOS-Chem chemistry mechanisms'' wiki page]].
| + | |
| − | | + | |
| − | GEOS-Chem will eventually include the option to create a custom chemistry mechanism and automatically build that mechanism with KPP when GEOS-Chem is compiled. For now, users still need to build custom chemistry mechanisms manually. The [[#KPP source code|KPP Git repository]] includes the following files:
| + | |
| − | | + | |
| − | geos2kpp_parser.pl
| + | |
| − | gckpp.kpp
| + | |
| − | rename_f90.sh
| + | |
| − | copy_gckpp.sh
| + | |
| − | | + | |
| − | These files can be used to build a custom chemistry mechanisms with KPP. Follow these steps:
| + | |
| − | | + | |
| − | #The <tt>geos2kpp_parser.pl</tt> script can be used to create <tt>globchem.eqn</tt>, <tt>globchem.spc</tt>, and <tt>globchem.def</tt> files from a <tt>globchem.dat</tt> file. You may also choose to start from the preexisting <tt>globchem.eqn</tt>, <tt>globchem.spc</tt>, and <tt>globchem.def</tt> files found in the <tt>KPP</tt> subdirectories of the GEOS-Chem source code (e.g. <tt>Code.v11-01/KPP/Standard/</tt>).
| + | |
| − | #If needed, edit <tt>globchem.eqn</tt> to add/modify reactions for your mechanism and edit <tt>globchem.spc</tt> to add/modify species to your mechanism. Put the final <tt>globchem.*</tt> files in the same directory as the <tt>gckpp.kpp</tt> file.
| + | |
| − | #Execute KPP using <code>kpp gckpp.kpp</code>. This will create the necessary <tt>gckpp_*.f90</tt> files.
| + | |
| − | #Run <code>rename_f90.sh</code> to rename <tt>gckpp_*.f90</tt> to <tt>gckpp_*.F90</tt>.
| + | |
| − | #Open <tt>copy_gckpp.sh</tt> and modify the lines for <tt>codedir</tt>. Save and run <code>copy_gckpp.sh</code>. This will copy the <tt>gckpp_*.F90</tt> files to <tt>codedir/Custom</tt>.
| + | |
| − | #Compile GEOS-Chem for the new mechanism by including <tt>CHEM=Custom</tt> in your make command.
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:18, 10 August 2016 (UTC)
| + | |
| − | | + | |
| − | === Production & loss diagnostic ===
| + | |
| − | | + | |
| − | Mike Long implemented prod/loss structure parsing in the KPP source code. This option can be turned on by adding the line <tt>#FLUX on</tt> to the <tt>gckpp.kpp</tt> file. For example:
| + | |
| − | | + | |
| − | #INTEGRATOR rosenbrock
| + | |
| − | #LANGUAGE Fortran90
| + | |
| − | #DRIVER none
| + | |
| − | #HESSIAN off
| + | |
| − | #MEX off
| + | |
| − | #STOICMAT off
| + | |
| − | <span style="color:red">#FLUX on</span>
| + | |
| − |
| + | |
| − | #INCLUDE globchem.spc
| + | |
| − | #INCLUDE globchem.eqn
| + | |
| − | | + | |
| − | For the FLUX option to work properly, MEX and HESSIAN must be turned off as shown above.
| + | |
| − | | + | |
| − | When changing any of the settings in <tt>gckpp.kpp</tt>, the chemistry mechanism will need to be rebuilt with KPP as [[#Building the chemical mechanism|described above]].
| + | |
| − | | + | |
| − | The pre-built chemistry mechanisms (Standard, Tropchem, SOA, SOA-SVPOA, and UCX) were built with the FLUX option turned on. The prod/loss rates from KPP can be saved out in GEOS-Chem via the [http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#PLMenu ND65 diagnostic].
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:18, 10 August 2016 (UTC) | + | |
| − | | + | |
| − | ==== Inefficiencies in prod/loss code ====
| + | |
| − | | + | |
| − | The current implementation of the prod/loss diagnostic output in KPP is somewhat inefficient. The code is set up so that KPP calculates the prod/loss of ALL species, even though a user may only want prod/loss rates for a few species. This likely adds considerable memory overhead because KPP has <tt>NREAC</tt> additional species to track. Making this code more efficient is a task for future development.
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:18, 10 August 2016 (UTC)
| + | |
| − | | + | |
| − | ==== Comparing KPP and SMVGEAR prod/loss output ====
| + | |
| − | | + | |
| − | During the implementation of FlexChem in [[GEOS-Chem v11-01#v11-01g|GEOS-Chem v11-01g]], KPP's new prod/loss diagnostics were validated against SMVGEAR's prod/loss diagnostics. We list the prod/loss rates [molec/cm3/s] for several families after a 20-minute run below:
| + | |
| − | | + | |
| − | {| border=1 cellspacing=0 cellpadding=5
| + | |
| − | |-bgcolor="#CCCCCC"
| + | |
| − | !width="100px"|Family Name
| + | |
| − | !width="100px"|v11-01f<br>(SMVGEAR)
| + | |
| − | !width="100px"|v11-01g<br>(KPP)
| + | |
| − | |-
| + | |
| − | |POX
| + | |
| − | |666.15
| + | |
| − | |194722.60
| + | |
| − | |-
| + | |
| − | |LOX
| + | |
| − | |15410.74
| + | |
| − | |208939.00
| + | |
| − | |-
| + | |
| − | |PCO
| + | |
| − | |5623.75
| + | |
| − | |5610.81
| + | |
| − | |-
| + | |
| − | |LCO
| + | |
| − | |9129.74
| + | |
| − | |8771.62
| + | |
| − | |-
| + | |
| − | |PBrOx
| + | |
| − | |5703.22
| + | |
| − | |58242.10
| + | |
| − | |-
| + | |
| − | |PBry
| + | |
| − | |2.04
| + | |
| − | |74904.60
| + | |
| − | |-
| + | |
| − | |LBry
| + | |
| − | |1.96E-04
| + | |
| − | |74902.59
| + | |
| − | |}
| + | |
| − | | + | |
| − | The differences in the prod/loss rates returned by SMVGEAR versus KPP are large for most of the families. Examination of the prod/loss code revealed that SMVGEAR is missing many of the reactions that KPP is flagging for inclusion in the prod/loss rate calculations. For example, the reactions that KPP uses for calculating PBrOx in the [[NOx-Ox-HC-Aer-Br_chemistry_mechanism|tropchem mechanism]] are listed below:
| + | |
| − | | + | |
| − | * HBr + OH --> RR240 + Br + H2O
| + | |
| − | 2 BrO --> RR241 + 2 Br + O2
| + | |
| − | BrO + NO --> RR243 + Br + NO2
| + | |
| − | * Br2 + OH --> RR245 + HOBr + Br
| + | |
| − | BrO + OH --> RR246 + Br + HO2
| + | |
| − | * CHBr3 + OH --> RR255 + 3 Br
| + | |
| − | * CH2Br2 + OH --> RR256 + 2 Br
| + | |
| − | * CH3Br + OH --> RR257 + Br
| + | |
| − | * Br2 --> RR382 + 2 Br
| + | |
| − | BrO --> RR383 + Br + O3
| + | |
| − | * HOBr --> RR384 + Br + OH
| + | |
| − | * BrNO3 --> RR385 + Br + NO3
| + | |
| − | * BrNO2 --> RR387 + Br + NO2
| + | |
| − | * CHBr3 --> RR388 + 3 Br
| + | |
| − | Br + O3 --> RR237 + BrO + O2
| + | |
| − | Br + NO3 --> RR247 + BrO + NO2
| + | |
| − | * BrNO3 --> RR386 + BrO + NO2
| + | |
| − | | + | |
| − | The starred reactions indicate reactions that SMVGEAR uses to compute PBrOx. KPP uses 17 total reactions for PBrOx, while SMVGEAR only uses 11 reactions. The number of reactions used in the prod/loss rate calculations for SMVGEAR and KPP (for the [[NOx-Ox-HC-Aer-Br_chemistry_mechanism|tropchem mechanism]]) are summarized below:
| + | |
| − | | + | |
| − | {| border=1 cellspacing=0 cellpadding=5
| + | |
| − | |-bgcolor="#CCCCCC"
| + | |
| − | !width="100px"|Family Name
| + | |
| − | !width="100px"|v11-01f<br>(SMVGEAR)
| + | |
| − | !width="100px"|v11-01g<br>(KPP)
| + | |
| − | |-
| + | |
| − | |POX
| + | |
| − | |92
| + | |
| − | |163
| + | |
| − | |-
| + | |
| − | |LOX
| + | |
| − | |45
| + | |
| − | |100
| + | |
| − | |-
| + | |
| − | |PCO
| + | |
| − | |57
| + | |
| − | |57
| + | |
| − | |-
| + | |
| − | |LCO
| + | |
| − | |7
| + | |
| − | |7
| + | |
| − | |-
| + | |
| − | |PBrOx
| + | |
| − | |11
| + | |
| − | |17
| + | |
| − | |-
| + | |
| − | |PBry
| + | |
| − | |4
| + | |
| − | |34
| + | |
| − | |-
| + | |
| − | |LBry
| + | |
| − | |0
| + | |
| − | |30
| + | |
| − | |}
| + | |
| − | | + | |
| − | The KPP generally use many more reactions in the prod/loss calculations than SMVGEAR. This may explain the large differences in the prod/loss rates returned by KPP as compared to SMVGEAR.
| + | |
| − | | + | |
| − | Note that SMVGEAR and KPP use the same number of reactions for calculating PCO and LCO (in fact, the reactions used are identical). The PCO and LCO rate values listed in the first table are also fairly similar between the SMVGEAR and KPP runs. This indicates that we can reasonably reproduce the old prod/loss rates from SMVGEAR when KPP uses the same reactions.
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:00, 11 August 2016 (UTC)
| + | |
| − | | + | |
| − | == Validation ==
| + | |
| − | | + | |
| − | TBD
| + | |
| − | | + | |
| − | == Known issues ==
| + | |
| − | | + | |
| − | TBD
| + | |
