Difference between revisions of "History collections for chemistry and photolysis"
(→The RxnRates Collection) |
(→The RxnRates Collection) |
||
(10 intermediate revisions by 3 users not shown) | |||
Line 38: | Line 38: | ||
'O1DconcAfterChem ', 'GIGCchem', | 'O1DconcAfterChem ', 'GIGCchem', | ||
'O3PconcAfterChem ', 'GIGCchem', | 'O3PconcAfterChem ', 'GIGCchem', | ||
+ | 'O3concAfterChem ', 'GIGCchem', | ||
+ | 'RO2concAfterChem ', 'GIGCchem', | ||
:: | :: | ||
Line 68: | Line 70: | ||
|HO2concAfterChem | |HO2concAfterChem | ||
|HO2 concentration immediately after exiting the FlexChem solver | |HO2 concentration immediately after exiting the FlexChem solver | ||
− | | | + | |v/v |
| | | | ||
| | | | ||
Line 98: | Line 100: | ||
*[[List of diagnostics archived to bpch format|ND43]], #5 | *[[List of diagnostics archived to bpch format|ND43]], #5 | ||
− | | | + | |-valign="top" |
+ | |O3concAfterChem | ||
+ | |O3 concentration immediately after exiting the FlexChem solver | ||
+ | |molec/cm3 | ||
+ | | | ||
+ | | | ||
+ | *fullchem | ||
+ | | | ||
+ | | | ||
+ | *None | ||
+ | |-valign="top" | ||
+ | |RO2concAfterChem | ||
+ | |RO2 concentration immediately after exiting the FlexChem solver | ||
+ | |molec/cm3 | ||
+ | | | ||
+ | | | ||
+ | *fullchem | ||
+ | | | ||
+ | | | ||
+ | *None | ||
+ | |||
+ | |} | ||
== The JValues collection == | == The JValues collection == | ||
Line 759: | Line 782: | ||
== The RxnRates Collection == | == The RxnRates Collection == | ||
− | The '''RxnRates''' collection contains reaction rates | + | The '''RxnRates''' collection contains reaction rates from the chemical mechanism (as computed by the KPP-generated solver code). For example, in the case of the NO + O<sub>3</sub> → NO<sub>2</sub> + O<sub>2</sub> reaction the returned quantity is k[NO][O<sub>3</sub>] in molec/cm<sup>3</sup>/s. |
Here is a sample definition section for the RxnRates collection. | Here is a sample definition section for the RxnRates collection. | ||
Line 769: | Line 792: | ||
# It is best to list individual reactions to avoid using too much memory. | # It is best to list individual reactions to avoid using too much memory. | ||
# Reactions should be listed as "RxnRate_EQnnn", where nnn is the reaction | # Reactions should be listed as "RxnRate_EQnnn", where nnn is the reaction | ||
− | # index as listed in KPP/fullchem/gckpp_Monitor.F90 (pad zeroes as needed | + | # index as listed in KPP/fullchem/gckpp_Monitor.F90 or KPP/Hg/gckpp_monitor.F90 |
− | + | # (pad zeroes as needed) | |
− | + | ||
# | # | ||
RxnRates.template: '%y4%m2%d2_%h2%n2z.nc4', | RxnRates.template: '%y4%m2%d2_%h2%n2z.nc4', | ||
Line 797: | Line 819: | ||
|-valign=top" | |-valign=top" | ||
!colspan="7" bgcolor="#CCFFFF"|Diagnostics applicable only to the fullchem simulation with the aciduptake option | !colspan="7" bgcolor="#CCFFFF"|Diagnostics applicable only to the fullchem simulation with the aciduptake option | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
|-valign="top" | |-valign="top" | ||
|RxnRate_EQnnn | |RxnRate_EQnnn | ||
− | |Rate for reaction nnn (see <tt> | + | |Rate for reaction nnn (see <tt>gckpp_Monitor.F90</tt> for the reaction number. |
− | |1/ | + | |s<sup>-1</sup> |
| | | | ||
+ | *?RXN? | ||
| | | | ||
*fullchem | *fullchem | ||
*Hg (in 13.4.0 and later) | *Hg (in 13.4.0 and later) | ||
+ | *carboncycle (in 14.1.0 and later) | ||
| | | | ||
| | | |
Latest revision as of 15:32, 30 August 2023
Previous | Next | Guide to GEOS-Chem History diagnostics
- Overview of History diagnostics
- Legend for History diagnostics
- Collections for History diagnostics
- History collections for advection and mixing
- History collections for aerosols
- History collections for budgets
- History collections for chemistry and photolysis
- History collections for convection and wet deposition
- History collections for dry deposition
- History collections for mercury and POPs
- History collections for met fields and related quantities
- History collections for methane
- History collections for species concentrations
- History collections for TransportTracers
- History collections for CO2
- Replicating the behavior of the prior timeseries diagnostics
- Adding new History diagnostics
Contents
The ConcAfterChem collection
The ConcAfterChem collection contains diagnostics for OH, HO2, etc. species immediately upon exiting the FlexChem solver.
Here is a sample definition section for the ConcAfterChem collection.
- If this collection is not already present in the HISTORY.rc file in the GEOS-Chem run directory for your selected simulation, you can copy and paste this into your HISTORY.rc file and edit accordingly.
- To prevent an individual field from being included in the diagnostic output, place a comment character # in front of the field name.
- Please see our Legend for History diagnostics page for more information about each of the collection tags.
ConcAfterChem.template: '%y4%m2%d2_%h2%n2z.nc4', ConcAfterChem.format: 'CFIO', ConcAfterChem.frequency: 00000100 000000 ConcAfterChem.duration: 00000100 000000 ConcAfterChem.mode: 'time-averaged' ConcAfterChem.fields: 'OHconcAfterChem ', 'GIGCchem', 'HO2concAfterChem ', 'GIGCchem', 'O1DconcAfterChem ', 'GIGCchem', 'O3PconcAfterChem ', 'GIGCchem', 'O3concAfterChem ', 'GIGCchem', 'RO2concAfterChem ', 'GIGCchem', ::
This table describes the diagnostic quantities belonging to the ConcAfterChem collection:
Diagnostic name | Description | Units | Wildcards | Simulations | Notes | Bpch equiv. |
---|---|---|---|---|---|---|
OHconcAfterChem | OH concentration immediately after exiting the FlexChem solver | molec/cm3 |
|
| ||
HO2concAfterChem | HO2 concentration immediately after exiting the FlexChem solver | v/v |
|
| ||
O1DconcAfterChem | O1D concentration immediately after exiting the FlexChem solver | molec/cm3 |
|
| ||
O3PconcAfterChem | O3P concentration immediately after exiting the FlexChem solver | molec/cm3 |
|
| ||
O3concAfterChem | O3 concentration immediately after exiting the FlexChem solver | molec/cm3 |
|
| ||
RO2concAfterChem | RO2 concentration immediately after exiting the FlexChem solver | molec/cm3 |
|
|
The JValues collection
The JValues collection contains diagnostics for photolysis rates for various chemical species, obtained from the FAST-JX photolysis mechanism.
Here is a sample definition section for the JValues collection.
- If this collection is not already present in the HISTORY.rc file in the GEOS-Chem run directory for your selected simulation, you can copy and paste this into your HISTORY.rc file and edit accordingly.
- To prevent an individual field from being included in the diagnostic output, place a comment character # in front of the field name.
- Please see our Legend for History diagnostics page for more information about each of the collection tags.
JValues.template: '%y4%m2%d2_%h2%n2z.nc4', JValues.format: 'CFIO', JValues.frequency: 00000000 010000 JValues.duration: 00000000 010000 JValues.mode: 'instantaneous' JValues.fields: 'Jval_?PHO? ', 'GIGCchem', ::
This table describes the diagnostic quantities belonging to the JValues collection:
Diagnostic name | Description | Units | Wildcards | Simulations | Notes | Bpch equiv. |
---|---|---|---|---|---|---|
Jval_<spcname> | Photolysis rates | 1/s |
|
|
||
JvalO3O1D | Photolysis rate of O3 -> O1D | 1/s |
|
|
||
JvalO3O3P | Photolysis rate of O3 -> O3P | 1/s |
|
|
The KppDiags collection
The KppDiags collection contains KPP solver diagnostics. You can leave this collection disabled unless you are interested in assessing the solver's performance.
Here is a sample definition section for the KppDiags collection.
- If this collection is not already present in the HISTORY.rc file in the GEOS-Chem run directory for your selected simulation, you can copy and paste this into your HISTORY.rc file and edit accordingly.
- To prevent an individual field from being included in the diagnostic output, place a comment character # in front of the field name.
- Please see our Legend for History diagnostics page for more information about each of the collection tags.
KppDiags.template: '%y4%m2%d2_%h2%n2z.nc4', KppDiags.format: 'CFIO', KppDiags.frequency: 010000 KppDiags.duration: 010000 KppDiags.mode: 'time-averaged' KppDiags.fields: 'KppIntCounts ', 'GIGCchem', 'KppJacCounts ', 'GIGCchem', 'KppTotSteps ', 'GIGCchem', 'KppAccSteps ', 'GIGCchem', 'KppRejSteps ', 'GIGCchem', 'KppLuDecomps ', 'GIGCchem', 'KppSubsts ', 'GIGCchem', 'KppSmDecomps ', 'GIGCchem', ::
This table describes the diagnostic quantities belonging to the JValues collection:
- NOTE: fullchem refers to all simulations that use a full-chemistry mechanism (i.e. benchmark, complexSOA*, standard, tropchem. aciduptake, marinePOA, RRTMG, TOMAS).
Diagnostic name | Description | Units | Wildcards | Simulations | Notes | Bpch equiv. |
---|---|---|---|---|---|---|
KppIntCounts | Number of times the KPP integrator was called | count |
|
|
||
KppJacCounts | Number of times the Jacobian matrix was constructed | count |
|
|
||
KppTotSteps | Total number of integration timesteps | count |
|
|
||
KppAccSteps | Number of accepted integration timesteps | count |
|
|
||
KppRejSteps | Number of rejected integration timesteps | count |
|
|
||
KppLuDecomps | Number of LU decompositions performed | count |
|
|
||
KppSubsts | Number of matrix substitutions performed. Includes forward and backward substitutions. | count |
|
|
||
KppSmDecomps | Number of singular matrix decopositions performed.
|
count |
|
|
The Metrics collection
The Metrics collection contains diagnostics for computing OH metrics from a GEOS-Chem full chemistry simulation. To compute the OH metrics, you must run the Python script metrics.py that ships with each fullchem or CH4 run directory.
Here is a sample definition section for the JValues collection.
- If this collection is not already present in the HISTORY.rc file in the GEOS-Chem run directory for your selected simulation, you can copy and paste this into your HISTORY.rc file and edit accordingly.
- To prevent an individual field from being included in the diagnostic output, place a comment character # in front of the field name.
- Please see our Legend for History diagnostics page for more information about each of the collection tags.
Metrics.template: '%y4%m2%d2_%h2%n2z.nc4', Metrics.format: 'CFIO', Metrics.frequency: 'End', Metrics.duration: 'End', Metrics.mode: 'time-averaged' Metrics.fields: 'AirMassColumnFull ', 'GIGCchem', 'CH4emission ', 'GIGCchem', 'CH4massColumnFull ', 'GIGCchem' 'CH4massColumnTrop ', 'GIGCchem', 'LossOHbyCH4columnTrop ', 'GIGCchem', 'LossOHbyMCFcolumnTrop ', 'GIGCchem', 'OHwgtByAirMassColumnFull ', 'GIGCchem', ::
The table below describes diagnostic quantities belonging to the ProdLoss collection.
Diagnostic name | Description | Units | Wildcards | Simulations | Notes | Bpch equiv. |
---|---|---|---|---|---|---|
AirMassColumnFull | Air mass, full atmosphere column sums | kg |
|
|
||
LossOHbyCH4columnTrop | Loss rate of CH4 by OH, tropospheric column sums | molec cm-3 |
|
|
||
LossOHbyMCFcolumnTrop | Loss rate of CH3CCl3 (aka MCF) by OH, tropospheric column sums | molec cm-3 |
|
|
||
OHwgtByAirMassColumnFull | Airmass-weighted mean OH concentration, full-atmosphere column sums | kg air kg OH m-3 |
|
|
||
Diagnostics applicable only to the CH4 or tagCH4 simulations | ||||||
CH4emission | CH4 emission, used for computing OH metrics | kg s-1 |
|
|
||
CH4massColumnFull | Airmass-weighted CH4 concentration, full-atmosphere column sums | kg air kg CH4 m-3 |
|
|
||
CH4massColumnTrop | Airmass-weighted CH4 concentration, troposphere-only column sums | kg air kg CH4 m-3 |
|
|
The ProdLoss Collection
The ProdLoss collection contains chemical production and loss rates.
Here is a sample definition section for the ProdLoss collection.
- If this collection is not already present in the HISTORY.rc file in the GEOS-Chem run directory for your selected simulation, you can copy and paste this into your HISTORY.rc file and edit accordingly.
- To prevent an individual field from being included in the diagnostic output, place a comment character # in front of the field name.
- Please see our Legend for History diagnostics page for more information about each of the collection tags.
- NOTE: This example is for the benchmark simulation. Some quantities in this collection are not applicable to certain simulations.
ProdLoss.template: '%y4%m2%d2_%h2%n2z.nc4', ProdLoss.format: 'CFIO', ProdLoss.frequency: 00000100 000000 ProdLoss.duration: 00000100 000000 ProdLoss.mode: 'time-averaged' ProdLoss.fields: 'Prod_?PRD? ', 'GIGCchem', 'ProdBCPIfromBCPO ', 'GIGCchem', 'ProdOCPIfromOCPO ', 'GIGCchem', 'ProdSO4fromH2O2inCloud ', 'GIGCchem', 'ProdSO4fromO2inCloudMetal ', 'GIGCchem', 'ProdSO4fromO3inCloud ', 'GIGCchem', 'ProdSO4fromO3inSeaSalt ', 'GIGCchem', 'ProdSO4fromHOBrInCloud ', 'GIGCchem', 'ProdSO4fromSRO3 ', 'GIGCchem', 'ProdSO4fromSRHObr ', 'GIGCchem', 'ProdSO4fromO3s ', 'GIGCchem', 'Loss_?LOS? ', 'GIGCchem', 'LossHNO3onSeaSalt ', 'GIGCchem', 'ProdCOfromCH4 ', 'GIGCchem', 'ProdCOfromNMVOC ', 'GIGCchem', ::
The table below describes diagnostic quantities belonging to the ProdLoss collection.
- NOTE: fullchem refers to all simulations that use a full-chemistry mechanism (i.e. benchmark, complexSOA*, standard, tropchem. aciduptake, marinePOA, RRTMG, TOMAS).
Diagnostic name | Description | Units | Wildcards | Simulations | Notes | Bpch equiv. |
---|---|---|---|---|---|---|
Diagnostics applicable only to the fullchem simulation with the aciduptake option | ||||||
ProdSO4fromOxidationOnDust | Production of SO4 from oxidation on dust aerosols | kg S/s |
|
|||
ProdNITfromHNO3uptakeOnDust | Production of NIT from HNO3 uptake on dust aerosols | kg N/s |
|
|||
ProdSO4fromUptakeOfH2SO4g | Production of SO4 from uptake of H2SO4(g) | kg S/s |
|
|||
Diagnostics applicable only to the aerosol-only simulation | ||||||
ProdSO2fromDMSandOH | Production of SO2 from DMS + OH (in sulfate_mod.F) | kg S/s |
|
|||
ProdSO2fromDMSandNO3 | Production of SO2 from DMS + NO3R | kg S/s |
|
|||
ProdSO2fromDMS | Total P(SO2) from DMS | kg S/s |
|
|||
ProdMSAfromDMS | Production of MSA from DMS | kg S/s |
|
|||
ProdSO4fromGasPhase | Production of SO4 in the gas phase | kg S/s |
|
|||
Diagnostics applicable only to the aerosol-only and fullchem simulations | ||||||
ProdBCPIfromBCPO | Production of hydrophilic BC from hydrophobic BCs | kg |
|
|||
ProdOCPIfromOCPO | Production of hydrophilic BC from hydrophobic BCs | kg |
|
|||
ProdSO4fromH2O2inCloud | Production of SO4 from aqueous oxidation of H2O2 in clouds | kg S/s |
|
|||
ProdSO4fromO2inCloudMetal | Production of SO4 from aqueous oxidation of O2 from metals in cloud | kg S/s |
|
|
||
ProdSO4fromO3inCloud | Production of SO4 from aqueous oxidation of O3 in clouds | kg S/s |
|
|||
ProdSO4fromO3inSeaSalt | Production of SO4 from O3 in sea salt | kg S/s |
|
|||
ProdSO4fromO3s | Production of SO4 from aqueous phase SO3-- loss by OH | kg S/s |
|
|
||
ProdSO4fromSRO3 | Production of SO4 from sulfur production rate of O3 | kg S/s |
|
|
||
LossHNO3onSeaSalt | Loss of HNO3 on sea salt aerosols | kg S/s |
|
|||
Diagnostics applicable only to full-chemistry simulations | ||||||
ProdSO4fromHOBrinCloud | Production of SO4 from aqueous oxidation of HOBr in clouds | kg S/s |
|
|
||
ProdSO4fromSRHOBr | Production of SO4 from sulfur production rate of HOBr+O3 | kg S/s |
|
|
||
ProdCOfromCH4 | Production of CO from CH4 | molec cm-3 s-1 |
|
|
||
ProdCOfromNMVOC | Production of CO from NMVOCs | molec cm-3 s-1 |
|
|
||
Diagnostics for production and loss of species or chemical families (e.g. Ox) | ||||||
Prod_<spcname> | Chemical production for a given species or family | molec/cm3/s |
|
|
||
Loss_<spcname> | Chemical loss for a given species or family | molec/cm3/s |
|
|
The RxnRates Collection
The RxnRates collection contains reaction rates from the chemical mechanism (as computed by the KPP-generated solver code). For example, in the case of the NO + O3 → NO2 + O2 reaction the returned quantity is k[NO][O3] in molec/cm3/s.
Here is a sample definition section for the RxnRates collection.
- If this collection is not already present in the HISTORY.rc file in the GEOS-Chem run directory for your selected simulation, you can copy and paste this into your HISTORY.rc file and edit accordingly.
- To prevent an individual field from being included in the diagnostic output, place a comment character # in front of the field name.
- Please see our Legend for History diagnostics page for more information about each of the collection tags.
# # It is best to list individual reactions to avoid using too much memory. # Reactions should be listed as "RxnRate_EQnnn", where nnn is the reaction # index as listed in KPP/fullchem/gckpp_Monitor.F90 or KPP/Hg/gckpp_monitor.F90 # (pad zeroes as needed) # RxnRates.template: '%y4%m2%d2_%h2%n2z.nc4', RxnRates.frequency: 00000000 010000 RxnRates.duration: 00000000 010000 RxnRates.mode: 'time-averaged' RxnRates.fields: 'RxnRate_EQ001 ', 'RxnRate_EQ002 ',
The table below describes diagnostic quantities belonging to the ProdLoss collection.
- NOTE: fullchem refers to all simulations that use a full-chemistry mechanism (i.e. benchmark, complexSOA*, standard, tropchem. aciduptake, marinePOA, RRTMG, TOMAS).
Diagnostic name | Description | Units | Wildcards | Simulations | Notes | Bpch equiv. |
---|---|---|---|---|---|---|
Diagnostics applicable only to the fullchem simulation with the aciduptake option | ||||||
RxnRate_EQnnn | Rate for reaction nnn (see gckpp_Monitor.F90 for the reaction number. | s-1 |
|
|
The StateChm Collection
The StateChm collection contains quantities from State_Chm, the Chemistry state object (other than the species concentrations, which are stored in the SpeciesConc collection).
Here is a sample definition section for the StateChm collection.
- If this collection is not already present in the HISTORY.rc file in the GEOS-Chem run directory for your selected simulation, you can copy and paste this into your HISTORY.rc file and edit accordingly.
- To prevent an individual field from being included in the diagnostic output, place a comment character # in front of the field name.
- Other fields of the State_Chm object can be added to this collection by prefixing the field name with Chem_.
- Please see our Legend for History diagnostics page for more information about each of the collection tags.
StateChm.template: '%y4%m2%d2_%h2%n2z.nc4', StateChm.format: 'CFIO', StateChm.frequency: 00000100 000000 StateChm.duration: 00000100 000000 StateChm.mode: 'time-averaged' StateChm.fields: 'Chem_phSav ', 'GIGCchem', 'Chem_HplusSav ', 'GIGCchem', 'Chem_WaterSav ', 'GIGCchem', 'Chem_SulRatSav ', 'GIGCchem', 'Chem_NaRatSav ', 'GIGCchem', 'Chem_AcidPurSav ', 'GIGCchem', 'Chem_BiSulSav ', 'GIGCchem', 'Chem_pHCloud ', 'GIGCchem', 'Chem_SSAlk', ', 'GIGCchem', 'Chem_HSO3AQ ', 'GIGCchem', 'Chem_SO3AQ ', 'GIGCchem', 'Chem_fupdateHOBr ', 'GIGCchem', ::
The table below describes diagnostic quantities belonging to the StateChm collection.
- NOTE: fullchem refers to all simulations that use a full-chemistry mechanism (i.e. benchmark, complexSOA*, standard, tropchem. aciduptake, marinePOA, RRTMG, TOMAS).
Diagnostic name | Description | Units | Wildcards | Simulations | Notes | Bpch equiv. |
---|---|---|---|---|---|---|
Chem_phSav | ISORROPIA aerosol pH | 1 |
|
|||
Chem_HplusSav | ISORROPIA H+ concentration | M |
|
|||
Chem_WaterSav | ISORROPIA aerosol water | μg/m3 |
|
|||
Chem_SulRatSav | ISORROPIA sulfate concentration | M |
|
|||
Chem_NaRatSav | ISORROPIA Na+ concentration | M |
|
|||
Chem_AcidPurSav | ISORROPIA acidpur ?? concentration | M |
|
|||
Chem_BiSulSav | ISORROPIA bisulfate (general acid) concentration | M |
|
|||
Chem_phCloud | Cloud PH | 1 |
|
|
||
Chem_SSAlk | Sea salt alkalinity | 1 |
|
|
||
Chem_HSO3AQ | Cloud bisulfite concentration | mol/L |
|
|
||
Chem_SO3AQ | Cloud sulfite concentration | mol/L |
|
|
||
Chem_fupdateHOBr | Correction factor for HOBr removal by SO2 | mol/L |
|
|