Archive output with the History diagnostics
Introduction
GEOS-Chem Classic and GCHP allow you to save various diagnostic fields to netCDF files via the History component. In the following sections you will learn how to schedule diagnostics for output using History.
Note
HEMCO has its own diagnostic archiving capability. Please see the
HEMCO diagnostics
chapter of The HEMCO User’s Guide
more information. Note that HEMCO diagnostics are configured
from configuration file HEMCO_Diagn.rc
for both GCHP and GEOS-Chem Classic
but diagnostics are also specified in HISTORY.rc
for GCHP only.
The HISTORY.rc configuration file
You can request GEOS-Chem diagnostics by editing the
HISTORY.rc
configuration file. This file lists the groups of
diagnostic outputs (called collections), and the parameters for
each collection (frequency of archival, mode of archiving, fields per
collection, etc.). GEOS-Chem Classic will write netCDF files for each
collection at the output intervals that you specify in
HISTORY.rc
. If using GCHP this is done by the MAPL History component
also using configuration file HISTORY.rc
.
The following is a stripped-down HISTORY.rc
file example for
illustration purposes. The HISTORY.rc
file that you will find
in your GEOS-Chem run directory will contain more collections than
what is shown below. This example is for GEOS-Chem Classic.
#============================================================================
# EXPID allows you to specify the beginning of the file path corresponding
# to each diagnostic collection. For example:
#
# EXPID: ./GEOSChem
# Will create netCDF files whose names begin "GEOSChem",
# in this run directory.
#
# EXPID: ./OutputDir/GEOSChem
# Will create netCDF files whose names begin with "GEOSChem"
# in the OutputDir sub-folder of this run directory.
#
#============================================================================
EXPID: ./OutputDir/GEOSChem
#==============================================================================
# %%%%% COLLECTION NAME DECLARATIONS %%%%%
#
# To disable a collection, place a "#" character in front of its name.
#==============================================================================
COLLECTIONS: 'SpeciesConc',
'SpeciesConcSubset',
'ConcAfterChem',
::
#==============================================================================
# %%%%% THE SpeciesConc COLLECTION %%%%%
#
# GEOS-Chem species concentrations (default = advected species)
#
# Available for all simulations
#==============================================================================
SpeciesConc.template: '%y4%m2%d2_%h2%n2z.nc4',
SpeciesConc.frequency: 00000000 060000
SpeciesConc.duration : 00000001 000000
SpeciesConc.format: 'CFIO',
SpeciesConc.mode: 'instantaneous',
SpeciesConc.fields: 'SpeciesConcVV_?ADV?',
'SpeciesConcMND_?ALL?',
::
#==============================================================================
# %%%%% THE SpeciesConcSubset COLLECTION %%%%%
#
# Same as the SpeciesConc collection, but will subset data in the horizontal
# and vertical dimensions so that the netCDF diagnostic files will cover
# a smaller region of the globe. This can save disk space and memory.
#
# NOTE: This capability will be available in GEOS-Chem "Classic" 12.5.0
# and later versions.
#
# Available for all simulations
#==============================================================================
SpeciesConcSubset.template: '%y4%m2%d2_%h2%n2z.nc4',
SpeciesConcSubset.frequency: 00000000 060000
SpeciesConcSubset.duration: 00000001 000000
SpeciesConcSubset.format: 'CFIO',
SpeciesConcSubset.mode: 'instantaneous',
SpeciesConcSubset.LON_RANGE: -40.0 60.0,
SpeciesConcSubset.LAT_RANGE: -10.0 50.0,
SpeciesConcSubset.levels: 1 2 3 4 5,
SpeciesConcSubset.fields: 'SpeciesConcVV_?ADV?',
::
#==============================================================================
# %%%%% THE ConcAfterChem COLLECTION %%%%%
#
# Concentrations of OH, HO2, O1D, O3P immediately after exiting the KPP solver
# or OH after the CH4 specialty-simulation chemistry routine.
#
# OH: Available for all full-chemistry simulations and CH4 specialty sim
# HO2: Available for all full-chemistry simulations
# O1D, O3P: Availalbe for full-chemistry simulations using UCX mechanism
#==============================================================================
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',
'HO2concAfterChem',
'O1DconcAfterChem',
'O3PconcAfterChem',
::
In this HISTORY.rc
file, we are requesting three collections, or
types of netCDF file output. The table below explains in more detail
parameters shown in the HISTORY.rc
file above.
This parameter controls the filename prefix. In this example,,
EXPID
is set to./OutputDir/GEOSChem
by default. This means that all diagnostic files will be written to the./OutputDir
subdirectory of the GEOS-Chem run directory, and will start with the prefixGEOSChem
.Note
Restart files are placed in the
./Restarts
subdirectory of the run directory instead of./OutputDir
, which only contains diagnostic files.
COLLECTIONS
The
COLLECTIONS:
tag specifies all of the diagnostic collections that you wish to activate during a GEOS-Chem simulation. Each collection represents a group of diagnostic quantities that will be written to disk in netCDF file format. The collection name will be automatically added to the netCDF file name along with the date/or time.Each GEOS-Chem run directory will ship with its own customized
HISTORY.rc
configuration file. Only the diagnostic collections pertaining to a particular GEOS-Chem simulation will be included in the correspondingHISTORY.rc
file.Each collection name must be bracketed by single quotes, and be followed by a comma.
To disable an entire diagnostic collection, simply put a
#
comment character in front of the collection name in theCOLLECTIONS:
section.GEOS-Chem will expect to find a collection definition section for each of the activated collections listed under the
COLLECTIONS:
section. In other words, if you haveSpeciesConc
listed underCOLLECTIONS:
, but there is no further information provided about theSpeciesConc
collection, then GEOS-Chem will halt with an error message.Note
You are not limited to the collections that are pre-defined in the
HISTORY.rc
configuration file. You may create additional diagnostic collections to suit your research purposes.
SpeciesConc
Name of the first collection in this
HISTORY.rc
file.
SpeciesConc.template
Determines the date and time format for the
SpeciesConc
collection filename suffix, as described below:
%y4%m2%d2_%h2%n2z.nc4
printsYYYYMMDD_hhmmz.nc4
to the end of each netCDF filename.
YYYYMMDD
is the date in year/month/day format;
hhmm
is the time inhour:minutes
format.
z
denotes “Zulu”, which is an abbreviation for UTC time.
.nc4
denotes that the data file is in the netCDF-4 format.Note
For example, the complete file path for the
SpeciesConc
collection at 00:00 UTC on 2020/01/01 will be./OutputDir/GEOSChem.SpeciesConc.20200101_0000z.nc4
, where:
EXPID
specifies the filename prefix (./OutputDir/GEOSChem
)
SpeciesConc.template
specifies the format of the filename suffix (.20200101_0000z.nc4
).
SpeciesConc.frequency
Determines how often the diagnostic quantities belonging to the
SpeciesConc
collection will be saved to a netCDF file. This can be specified as eitherhhmmss
orYYYYMMDD hhmmss
.In the above example, data belonging to the collection will be written to the file every 6 hours. Because
SpeciesConc
is an instantaneous collection, no time-averaging will be performed.
SpeciesConc.format
For GCHP simulations only. This tag indicates the I/O library that will be used.
SpeciesConc.duration
Determines how often a new
SpeciesConc
netCDF file will be created. Uses the same format asSpeciesConc.frequency
.
SpeciesConc.mode
Determines the averaging method for the
SpeciesConc
collection. Allowable values are:
instantaneous
: Archives instantaneous values at the interval specified bySpeciesConc.frequency
.
time-averaged
: Archives values averaged over the interval specified bySpeciesConc.frequency
.
SpeciesConc.fields
Lists the diagnostic quantities to be archived in the
SpeciesConc
collection. Some diagnostic quantities (e.g. concentrations, fluxes, masses) may also have an extra dimension, which represents species, size bins, reaction numbers, etc.For example, to request the ozone species concentration (in mixing ratio units) you may use the field name
SpeciesConcVV_O3
. The species name is separated from the quantity name by a single underscore.Note
For GCHP, each diagnostic field must be followed by the name of the ESMF gridded component that it is associated with. For arrays in GEOS-Chem objects State_Met, State_Chm, and State_Diag this is
'GCHPchem',
. A few diagnostics may also be output from the advection component of GCHP which has ESMF gridded component nameDYNAMICS
.If you are using GEOS-Chem Classic,you may also use a wildcard to specify a given category of species. In the above example
SpeciesConcVV_?ADV?
refers to all advected species andSpeciesConcVV_?ALL?
refers to all species (both advected and non-advected).Note
GCHP does not allow the use of wildcards. Each diagnostic quantity must be listed individually.
::
Signifies the end of the SpeciesConc definition section.
::
may be placed at any column.
SpeciesConc.subset
Name of the second diagnostic collection specified in this sample
HISTORY.rc
configuration file. In this collection we will request output to be restricted to a subset of the horizontal grid.The
.template
,.frequency
,.duration
,:mode
, and.fields
are described for theSpeciesConc
collection above, so we will not repeat them here.
SpeciesConcSubset.LON_RANGE
Defines the longitude range (
min max
) where diagnostic data will be archived. Data outside of this range will be ignored. If this option is omitted, values at all longitudes (-180 180
) will be included.
SpeciesConcSubset.LAT_RANGE
Defines the latitude range (
min max
) where diagnostic data will be archived. Data outside of this range will be ignored. If this option is omitted, values at all latitudes (-90 90
) will be included.
SpeciesConcSubset.levels
Specifies the levels that you wish to be included in the diagnostic archiving. If omitted, data at all levels will be included.
Note
In GEOS-Chem Classic, all levels between the minimum and maximum level specified will be included in the diagnostic archival. This differs from the behavior in GCHP, which archives only the specified levels.
ConcAfterChem
Name of the third collection specified in this sample
HISTORY.rc
configuration file.The
.template
,.frequency
,.duration
,:mode
, and.fields
are described for theSpeciesConc
collection above, so we will not repeat them here.
ConcAftercChem.mode
In this example, the
ConcAfterChem.mode
setting indicates that theConcAfterChem
collection will contain time-averaged data. The averaging interval is set in the
Wildcards (GEOS-Chem Classic only)
For GEOS-Chem Classic diagnostic output, you can use the following wildcards with diagnostic quantities that have a species/bin/reaction dimension:
Wildcard |
Description |
Example |
---|---|---|
?ADV? |
Advected species |
SpeciesConcVV_?ADV? |
?AER? |
Aerosol species |
SpeciesConcVV_?AER? |
?ALL? |
All species |
SpeciesConcVV_?ALL? |
?DRY? |
Dry-deposited species |
SpeciesConcVV_?DRY? |
?DRYALT? |
Species for the histguide-concabovechem collection |
SpeciesConcVV_?DRYALT |
?DUSTBIN? |
Dust bin number |
AODdust550nm_?DUSTBIN? |
?FIX? |
Fixed species in the KPP chemistry mechanism |
SpeciesConcVV_?FIX? |
?GAS? |
Gas-phase species |
SpeciesConcVV_?GAS? |
?HYG? |
Aerosol species that undergo hygroscopic growth (e.g. black carbon) |
AODhyg550nm_?HYG? |
?KPP? |
All species (fixed + variable) in the KPP chemistry mechanism |
SpeciesConcVV_?KPP? |
?LOS? |
Chemical loss species or families |
SpeciesConcVV_?LOS? |
?PHO? |
Photolyzed species |
SpeciesConcVV_?PHO? |
?PRD? |
Chemical production species or families |
SpeciesConcVV_?PRD? |
?RRTMG? |
RRTMG-computed fluxes |
RadAllSkywSurf_?RRTMG? |
?RXN? |
KPP reaction rates |
RxnRate_?RXN? |
?TOMASBIN? |
TOMAS size bins |
TomasH2SO4Mass_?TOMASBIN? |
?UVFLX? |
UV flux bins |
UVFluxDiffuse_?UVFLX? |
?VAR? |
Variable species in the KPP mechanism |
SpeciesConcVV_?VAR? |
?WET? |
Wet-deposited species |
SpeciesConcVV_?WET |
Prefixes
You may add any field from the State_Met
and State_Chm
objects to any diagnostic collection as well. These fields must be
prefixed as described below:
Wildcard |
Description |
Example |
---|---|---|
Chem_ |
Request diagnostic output
from the |
Chem_pHCloud |
Met_ |
Request diagnostic output
from the |
Met_SPHU |
File naming convention
As mentioned above, SpeciesConc.template
, GEOS-Chem History
files adhere to the following naming convention:
EXPID.collection-name.collection-template
e.g.
../OutputDir/GEOSChem.SpeciesConc.20200101_0000z.nc4
The duration tag of each collection in HISTORY.rc
controls how
often a new file will be written to disk, as we saw above:
SpeciesConc.duration: 00000001 000000 # Write a new file each day
SpeciesConcSubset.duration: 00000001 000000 # Write a new file each day
ConcAfterChem.duration: 00000100 000000 # Write a new file each month
Therefore, based on all of these settings in our example
HISTORY.rc
file, GEOS-Chem will write the following netCDF
files to disk in the current run directory:
GEOSChem.SpeciesConc.20160101_0000z.nc4 GEOSChem.SpeciesConcSubset.20160101_0000z.nc4
GEOSChem.SpeciesConc.20160102_0000z.nc4 GEOSChem.SpeciesConcSubset.20160102_0000z.nc4
GEOSChem.SpeciesConc.20160103_0000z.nc4 GEOSChem.SpeciesConcSubset.20160102_0000z.nc4
GEOSChem.SpeciesConc.20160104_0000z.nc4 GEOSChem.SpeciesConcSubset.20160104_0000z.nc4
... etc ...
GEOSChem.ConcAfterChem.20160101_0000z.nc4
GEOSChem.ConcAfterChem.20160201_0000z.nc4
GEOSChem.ConcAfterChem.20160301_0000z.nc4
GEOSChem.ConcAfterChem.20160401_0000z.nc4
... etc ...
Vertical coordinates in netCDF files
All netCDF files produced by GEOS-Chem (i.e. diagnostic files and restart files) adhere to the the COARDS netCDF convention. for the lon, lat, and time dimensions.
For the vertical dimension, we have chosen to use the following coordinate variables to include in GEOS-Chem Classic output files, emulating the file format of the NCAR Community Earth System Model (CESM):
variables:
double lev(lev) ;
lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;
lev:units = "level" ;
lev:positive = "down" ;\
lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;
double hyam(lev) ;
hyam:long_name = "hybrid A coefficient at layer midpoints" ;
double hybm(lev) ;
hybm:long_name = "hybrid B coefficient at layer midpoints" ;
double ilev(ilev) ;
ilev:long_name = "hybrid level at interfaces (1000*(A+B))" ;
ilev:units = "level" ;
ilev:positive = "down" ;
ilev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
ilev:formula_terms = "a: hyai b: hybi p0: P0 ps: PS" ;
double hyai(ilev) ;
hyai:long_name = "hybrid A coefficient at layer interfaces" ;
double hybi(ilev) ;
hybi:long_name = "hybrid B coefficient at layer interfaces" ;
double P0 ;
P0:long_name = "reference pressure" ;
The lev variable is used for data that is placed on the midpoints between vertical levels. This is an “approximate” eta coordinate, which is close to 1 at the surface and close to zero at the atmosphere top.
lev = 0.99250002413, 0.97749990013, 0.962499776, 0.947499955, 0.93250006,
0.91749991, 0.90249991, 0.88749996, 0.87249996, 0.85750006, 0.842500125,
0.82750016, 0.8100002, 0.78750002, 0.762499965, 0.737500105, 0.7125001,
0.6875001, 0.65625015, 0.6187502, 0.58125015, 0.5437501, 0.5062501,
0.4687501, 0.4312501, 0.3937501, 0.3562501, 0.31279158, 0.26647905,
0.2265135325, 0.192541016587707, 0.163661504087706, 0.139115, 0.11825,
0.10051436, 0.085439015, 0.07255786, 0.06149566, 0.05201591, 0.04390966,
0.03699271, 0.03108891, 0.02604911, 0.021761005, 0.01812435, 0.01505025,
0.01246015, 0.010284921, 0.008456392, 0.0069183215, 0.005631801,
0.004561686, 0.003676501, 0.002948321, 0.0023525905, 0.00186788,
0.00147565, 0.001159975, 0.00090728705, 0.0007059566, 0.0005462926,
0.0004204236, 0.0003217836, 0.00024493755, 0.000185422, 0.000139599,
0.00010452401, 7.7672515e-05, 5.679251e-05, 4.0142505e-05, 2.635e-05,
1.5e-05 ;
The lev variable may be used for quick plotting. To compute the actual pressure at the midpoint of the grid box (I,J,L), you will need to supply your own 2-D surface pressure field (e.g. saved from another diagnostic file):
Pmid = ( hyam(L) * PS(I,J) ) + hybm(L)
The ilev variable is used for data that is placed on vertical level edges or “interfaces” (hence the “i” in ilev). This is also an “approximate” eta coordinate.
ilev = 1, 0.98500004826, 0.969999752, 0.9549998, 0.94000011, 0.92500001,
0.90999981, 0.89500001, 0.87999991, 0.86500001, 0.85000011, 0.83500014,
0.82000018, 0.80000022, 0.77499982, 0.75000011, 0.7250001, 0.7000001,
0.6750001, 0.6375002, 0.6000002, 0.5625001, 0.5250001, 0.4875001,
0.4500001, 0.4125001, 0.3750001, 0.3375001, 0.28808306, 0.24487504,
0.208152025, 0.176930008175413, 0.150393, 0.127837, 0.108663, 0.09236572,
0.07851231, 0.06660341, 0.05638791, 0.04764391, 0.04017541, 0.03381001,
0.02836781, 0.02373041, 0.0197916, 0.0164571, 0.0136434, 0.0112769,
0.009292942, 0.007619842, 0.006216801, 0.005046801, 0.004076571,
0.003276431, 0.002620211, 0.00208497, 0.00165079, 0.00130051, 0.00101944,
0.0007951341, 0.0006167791, 0.0004758061, 0.0003650411, 0.0002785261,
0.000211349, 0.000159495, 0.000119703, 8.934502e-05, 6.600001e-05,
4.758501e-05, 3.27e-05, 2e-05, 1e-05 ;
To compute the actual pressure at the bottom and top edges of the grid box (I,J,L), you will need to supply your own 2-D surface pressure field (e.g. saved from another diagnostic file):
Pbot = ( hyai(L ) * PS(I,J) ) + hybi(L )
Ptop = ( hyai(L+1) * PS(I,J) ) + hybi(L+1)
GCHP history files contain simply the lev coordinate which is the level index of the model.
Level 1 corresponds to surface for all diagnostics retrieved from GEOS-Chem
arrays State_Met, State_Diag, and State_Chm. These include all diagnostics
that have prefix Met_
or Chm_
as well as all other
GCHPchem
diagnostics that do not begin with Emis
.
Diagnostics that begin with Emis
or that have gridded component
name other than GCHPchem
have Level 1 correspond to top-of-atmosphere.
Diagnostic collections
The diagnostic collections described in the sections below are used by
default in GEOS-Chem simulations. You can create your own customized
collections by modifying the HISTORY.rc
file.
The only restriction is that you cannot mix data that is placed on grid box layer edges in the same collection as data placed on grid box layer centers. This violates the netCDF convention that all data variables have to be defined with the same vertical dimension.
Note
For diagnostic quantities that have a species/bin/reaction
dimension, we will use the abbreviation <name|wc>
to
indicate a species/bin/reaction name or wildcard.
For example, SpeciesConcVV_<name|wc>
means that the
diagnostic quantity can be a single species
(SpeciesConcVV_O3
) or a wildcarded subset
of species (SpeciesConcVV_?ADV?
).
AdvFluxVert
The AdvFluxVert collection contains diagnostics for vertical transport in GEOS-Chem Classic advection. In practice, this collection is only used to obtain the vertical flux of O3 from GEOS-Chem Classic fullchem benchmark simulations. Most GEOS-Chem users will not need to activate this collection.
Sample definition section for HISTORY.rc
AdvFluxVert.template: '%y4%m2%d2_%h2%n2z.nc4',
AdvFluxVert.frequency: 00000100 000000
AdvFluxVert.duration: 00000100 000000
AdvFluxVert.mode: 'time-averaged'
AdvFluxVert.fields: 'AdvFluxVert_O3',
::
List of diagnostic fields in the AdvFluxVert collection
Diagnostic field |
Description |
Units |
Wildcards |
---|---|---|---|
AdvFluxVert_<name|wc> |
Vertical flux of species in advection |
kg/s |
?ADV? |
AerosolMass
The AerosolMass collection contains diagnostics for aerosol mass and particulate matter from full-chemistry simulations.
Sample definition section for HISTORY.rc
AerosolMass.template: '%y4%m2%d2_%h2%n2z.nc4',
AerosolMass.frequency: 00000100 000000
AerosolMass.duration: 00000100 000000
AerosolMass.mode: 'time-averaged'
AerosolMass.fields: 'AerMassASOA ',
'AerMassBC ',
'AerMassINDIOL ',
'AerMassISN1OA ',
'AerMassISOA ',
'AerMassLVOCOA ',
'AerMassNH4 ',
'AerMassNIT ',
'AerMassOPOA ',
'AerMassPOA ',
'AerMassSAL ',
'AerMassSO4 ',
'AerMassSOAGX ',
'AerMassSOAIE ',
'AerMassSOAME ',
'AerMassSOAMG ',
'AerMassTSOA ',
'BetaNO ',
'PM25 ',
'TotalBiogenicOA',
'TotalOA ',
'TotalOC ',
::
List of diagnostic fields in the AerosolMass collection
Diagnostic field |
Description |
Units |
---|---|---|
AerMassASOA[1] |
Aerosol products of light aromatics + IVOC oxidation |
\({\mu}g/m^3\) |
AerMassBC |
Aerosol products of light aromatics + IVOC oxidation |
\({\mu}g/m^3\) |
AerMassINDIOL[1] |
Generic aerosol-phase organonitrate hydrolysis product |
\({\mu}g/m^3\) |
AerMassISN10A[1] |
Aerosol phase 2nd generation hydroxynitrates formed from ISOP + NO3 rxn pathway |
\({\mu}g/m^3\) |
AerMassISOA[1] |
Aerosol products of isoprene oxidation |
\({\mu}g/m^3\) |
AerMassLVOCOA[1] |
Aerosol-phase low-volatility non-IEPOX product of ISOPOOH (RIP) oxidation |
\({\mu}g/m^3\) |
AerMassNH4 |
Ammonium |
\({\mu}g/m^3\) |
AerMassNIT |
Inorganic nitrate aerosol |
\({\mu}g/m^3\) |
AerMassPOA[1] |
Aerosols from SVOCs |
\({\mu}g/m^3\) |
AerMassOPOA[1] |
Aerosols products of POG oxidation |
\({\mu}g/m^3\) |
AerMassSAL |
Sea salt aerosol (SALA+SALC) |
\({\mu}g/m^3\) |
AerMassSO4 |
Sulfate |
\({\mu}g/m^3\) |
AerMassSOAGX[1] |
Aerosol phase glyoxal |
\({\mu}g/m^3\) |
AerMassSOAIE[1] |
Aerosol phase IEPOX |
\({\mu}g/m^3\) |
AerMassSOAME[1] |
Aerosol phase IMAE |
\({\mu}g/m^3\) |
AerMassSOAMG[1] |
Aerosol phase methylglyoxal |
\({\mu}g/m^3\) |
AerMassTSOA[1] |
Aerosol products of terpene oxidation |
\({\mu}g/m^3\) |
BetaNO[1] |
NO branching ratio |
\({\mu}g/m^3\) |
PM25 |
Particulate matter (d < 2.5 :math:{mu}m`) |
\({\mu}g/m^3\) |
TotalBiogenicOA[2] |
Sum of all organic aerosol |
\({\mu}g/m^3\) |
TotalOA[1] |
Sum of all organic aerosol |
\({\mu}g/m^3\) |
TotalOC |
Sum of all organic carbon |
\({\mu}g/m^3\) |
Notes for the AerosolMass collection
Aerosols
The Aerosols collection contains diagnostics for aerosol optical depth and related quantities from full-chemistry simulations.
Note
Some diagnostic fields in the Aerosols collection may be computed
at up to 3 wavelengths (WL1
, WL2
,
WL3
) as specified in this menu of the
geoschem_config.yml
file:
rrtmg_rad_transfer_model:
... etc ...
aod_wavelengths_in_nm:
- 550
For GEOS-Chem simulations that do not use the RRTMG radiative
transfer model, you may specify only one wavelength WL1
,
which is set to a default value of 550 nm. For GEOS-Chem
simulations using RRTMG, you may specify up to 2 more additional
wavelengths (WL2
and WL3
). GEOS-Chem will
replace the tokens WL1
, WL2
, WL3
in diagnostic field names with the corresponding wavelength.
For example, these diagnostic fields:
AODHygWL1_BCPI
AODDustWL1_DST1
AODStratLiquidAerWL1
AODPolarStratCloudWL1
AODSOAfromAqIsopreneWL1
AODStratLiquidAerWL1
AODPolarStratCloudWL1
will be saved to the GEOSChem.Aerosols.YYYYMMDD_hhmmz.nc4
file(s) under these names:
AODHyg550nm_BCPI
AODDust550nm_DST1
AODStratLiquidAer550nm
AODPolarStratCloud550nm
AODSOAfromAqIsoprene550nm
AODStratLiquidAer550nm
AODPolarStratCloud550nm
Sample definition section for HISTORY.rc
Aerosols.template: '%y4%m2%d2_%h2%n2z.nc4',
Aerosols.frequency: 00000100 000000
Aerosols.duration: 00000100 000000
Aerosols.mode: 'time-averaged'
Aerosols.fields: 'AODDust ',
'AODDustWL1_?DUSTBIN? ',
'AODHygWL1_?HYG? ',
'AODSOAfromAqIsopreneWL1 ',
'AODStratLiquidAerWL1 ',
'AODPolarStratCloudWL1 ',
'AerHygroscopicGrowth_?HYG? ',
'AerNumDensityStratLiquid ',
'AerNumDensityStratParticulate',
'AerAqueousVolume ',
'AerSurfAreaDust ',
'AerSurfAreaHyg_?HYG? ',
'AerSurfAreaStratLiquid ',
'AerSurfAreaPolarStratCloud ',
'Chem_AeroAreaMDUST1 ',
'Chem_AeroAreaMDUST2 ',
'Chem_AeroAreaMDUST3 ',
'Chem_AeroAreaMDUST4 ',
'Chem_AeroAreaMDUST5 ',
'Chem_AeroAreaMDUST6 ',
'Chem_AeroAreaMDUST7 ',
'Chem_AeroAreaSULF ',
'Chem_AeroAreaBC ',
'Chem_AeroAreaOC ',
'Chem_AeroAreaSSA ',
'Chem_AeroAreaSSC ',
'Chem_AeroAreaBGSULF ',
'Chem_AeroAreaICEI ',
'Chem_AeroRadiMDUST1 ',
'Chem_AeroRadiMDUST2 ',
'Chem_AeroRadiMDUST3 ',
'Chem_AeroRadiMDUST4 ',
'Chem_AeroRadiMDUST5 ',
'Chem_AeroRadiMDUST6 ',
'Chem_AeroRadiMDUST7 ',
'Chem_AeroRadiSULF ',
'Chem_AeroRadiBC ',
'Chem_AeroRadiOC ',
'Chem_AeroRadiSSA ',
'Chem_AeroRadiSSC ',
'Chem_AeroRadiBGSULF ',
'Chem_AeroRadiICEI ',
'Chem_WetAeroAreaMDUST1 ',
'Chem_WetAeroAreaMDUST2 ',
'Chem_WetAeroAreaMDUST3 ',
'Chem_WetAeroAreaMDUST4 ',
'Chem_WetAeroAreaMDUST5 ',
'Chem_WetAeroAreaMDUST6 ',
'Chem_WetAeroAreaMDUST7 ',
'Chem_WetAeroAreaSULF ',
'Chem_WetAeroAreaBC ',
'Chem_WetAeroAreaOC ',
'Chem_WetAeroAreaSSA ',
'Chem_WetAeroAreaSSC ',
'Chem_WetAeroAreaBGSULF ',
'Chem_WetAeroAreaICEI ',
'Chem_WetAeroRadiMDUST1 ',
'Chem_WetAeroRadiMDUST2 ',
'Chem_WetAeroRadiMDUST3 ',
'Chem_WetAeroRadiMDUST4 ',
'Chem_WetAeroRadiMDUST5 ',
'Chem_WetAeroRadiMDUST6 ',
'Chem_WetAeroRadiMDUST7 ',
'Chem_WetAeroRadiSULF ',
'Chem_WetAeroRadiBC ',
'Chem_WetAeroRadiOC ',
'Chem_WetAeroRadiSSA ',
'Chem_WetAeroRadiSSC ',
'Chem_WetAeroRadiBGSULF ',
'Chem_WetAeroRadiICEI ',
'Chem_StatePSC ',
'Chem_KhetiSLAN2O5H2O ',
'Chem_KhetiSLAN2O5HCl ',
'Chem_KhetiSLAClNO3H2O ',
'Chem_KhetiSLAClNO3HCl ',
'Chem_KhetiSLAClNO3HBr ',
'Chem_KhetiSLABrNO3H2O ',
'Chem_KhetiSLABrNO3HCl ',
'Chem_KhetiSLAHOClHCl ',
'Chem_KhetiSLAHOClHBr ',
'Chem_KhetiSLAHOBrHCl ',
'Chem_KhetiSLAHOBrHBr ',
::
List of diagnostic fields in the Aerosols collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
AODDust[3] |
Mineral dust optical depth @ WL1 |
1 |
|
AODDustWL1_<name|wc>[3] |
AOD for each dust bin @ WL1 |
1 |
?DUSTBIN? |
AODHygWL1_<name|wc>[3] |
AOD @ WL1 for aerosol species |
1 |
?HYG? |
AODSOAfromAqIsopreneWL1[4] |
Optical depth of SOA from aqueous isoprene @ WL1 |
1 |
|
AODStratLiquidAerWL1 |
Stratospheric liquid optical depth @ WL1 |
1 |
|
AODPolarStratCloudWL1 |
Polar stratospheric cloud type 1a/2 optical depth @ WL1 |
1 |
|
AerHygroscopicGrowth_<name|wc>[3] |
Hygroscopic growth of aerosol species |
1 |
?HYG? |
AerNumDensityStratLiquid |
Stratospheric liquid aerosol number density |
1/cm3 |
|
AerNumDensityStratParticulate |
Stratospheric particulate aerosol number density |
1/cm3 |
|
AerAqueousVolume |
Aqueous aerosol volume |
cm2/cm3 |
|
AerSurfAreaDust |
Surface area of mineral dust |
cm2/cm3 |
|
AerSurfAreaHyg_<name|wc> |
Surface area of aerosol species |
cm2/cm3 |
?HYG? |
AerSurfAreaStratLiquid |
Stratospheric liquid surface area |
cm2/cm3 |
|
AerSurfaceAreaPolarStratCloud |
Polar stratospheric cloud type 1a/2 surface area |
cm2/cm3 |
|
Chem_AeroAreaMDUST1[3] |
Dry aerosol area for mineral dust (0.15 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaMDUST2[3] |
Dry aerosol area for mineral dust (0.25 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaMDUST3[3] |
Dry aerosol area for mineral dust (0.4 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaMDUST4[3] |
Dry aerosol area for mineral dust (0.8 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaMDUST5[3] |
Dry aerosol area for mineral dust (1.5 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaMDUST6[3] |
Dry aerosol area for mineral dust (2.5 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaMDUST7[3] |
Dry aerosol area for mineral dust (4.0 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaSULF[3] |
Dry aerosol area for sulfate aerosol |
cm2/cm3 |
|
Chem_AeroAreaBC[3] |
Dry aerosol area for black carbon |
cm2/cm3 |
|
Chem_AeroAreaOC[3] |
Dry aerosol area for organic carbon |
cm2/cm3 |
|
Chem_AeroAreaSSA[3] |
Dry aerosol area for sea salt aerosol, accumulation mode |
cm2/cm3 |
|
Chem_AeroAreaSSC[3] |
Dry aerosol area for sea salt aerosol, coarse mode |
cm2/cm3 |
|
Chem_AeroAreaBGSULF |
Dry aerosol area for background stratospheric sulfate |
cm2/cm3 |
|
Chem_AeroAreaICEI |
Dry aerosol area for irregular ice cloud |
cm2/cm3 |
|
Chem_AeroRadiMDUST1[3] |
Dry aerosol radius for mineral dust (0.15 \({\mu}m\)) |
cm |
|
Chem_AeroRadiMDUST2[3] |
Dry aerosol radius for mineral dust (0.25 \({\mu}m\)) |
cm |
|
Chem_AeroRadiMDUST3[3] |
Dry aerosol radius for mineral dust (0.4 \({\mu}m\)) |
cm |
|
Chem_AeroRadiMDUST4[3] |
Dry aerosol radius for mineral dust (0.8 \({\mu}m\)) |
cm |
|
Chem_AeroRadiMDUST5[3] |
Dry aerosol radius for mineral dust (1.5 \({\mu}m\)) |
cm |
|
Chem_AeroRAdiMDUST6[3] |
Dry aerosol radius for mineral dust (2.5 \({\mu}m\)) |
cm |
|
Chem_AeroRadiMDUST7[3] |
Dry aerosol radius for mineral dust (4.0 \({\mu}m\)) |
cm |
|
Chem_AeroRadiSULF[3] |
Dry aerosol radius for sulfate aerosol |
cm |
|
Chem_AeroRadiBC[3] |
Dry aerosol radius for black carbon |
cm |
|
Chem_AeroRadiOC[3] |
Dry aerosol radius for organic carbon |
cm |
|
Chem_AeroRadiSSA[3] |
Dry aerosol radius for sea salt aerosol, accumulation mode |
cm |
|
Chem_AeroRadiSSC[3] |
Dry aerosol radius for sea salt aerosol, coarse mode |
cm |
|
Chem_AeroRadiBGSULF |
Dry aerosol radius for background stratospheric sulfate |
cm |
|
Chem_AeroRadiICEI |
Dry aerosol Radius for irregular ice cloud |
cm |
|
Chem_WetAeroAreaMDUST1[3] |
Wet aerosol area for mineral dust (0.15 \({\mu}m\)) |
cm2/cm3 |
|
Chem_WetAeroAreaMDUST2[3] |
Wet aerosol area for mineral dust (0.25 \({\mu}m\)) |
cm2/cm3 |
|
Chem_WetAeroAreaMDUST3[3] |
Wet aerosol area for mineral dust (0.4 \({\mu}m\)) |
cm2/cm3 |
|
Chem_WetAeroAreaMDUST4[3] |
Wet aerosol area for mineral dust (0.8 \({\mu}m\)) |
cm2/cm3 |
|
Chem_WetAeroAreaMDUST5[3] |
Wet aerosol area for mineral dust (1.5 \({\mu}m\)) |
cm2/cm3 |
|
Chem_WetAeroAreaMDUST6[3] |
Wet aerosol area for mineral dust (2.5 \({\mu}m\)) |
cm2/cm3 |
|
Chem_AeroAreaMDUST7[3] |
Dry aerosol area for mineral dust (4.0 \({\mu}m\)) |
cm2/cm3 |
|
Chem_WetAeroAreaSULF[3] |
Wet aerosol area for sulfate aerosol |
cm2/cm3 |
|
Chem_WetAeroAreaBC[3] |
Wet aerosol area for black carbon |
cm2/cm3 |
|
Chem_WetAeroAreaOC[3] |
Wet aerosol area for organic carbon |
cm2/cm3 |
|
Chem_WetAeroAreaSSA[3] |
Wet aerosol area for sea salt aerosol, accumulation mode |
cm2/cm3 |
|
Chem_WetAeroAreaSSC[3] |
Wet aerosol area for sea salt aerosol, coarse mode |
cm2/cm3 |
|
Chem_WetAeroAreaBGSULF |
Wet aerosol area for background stratospheric sulfate |
cm2/cm3 |
|
Chem_WetAeroAreaICEI |
Wet aerosol area for irregular ice cloud |
cm2/cm3 |
|
Chem_WetAeroRadiMDUST1[3] |
Wet aerosol radius for mineral dust (0.15 \({\mu}m\)) |
cm |
|
Chem_WetAeroRadiMDUST2[3] |
Wet aerosol radius for mineral dust (0.25 \({\mu}m\)) |
cm |
|
Chem_WetAeroRadiMDUST3[3] |
Wet aerosol radius for mineral dust (0.4 \({\mu}m\)) |
cm |
|
Chem_WetAeroRadiMDUST4[3] |
Wet aerosol radius for mineral dust (0.8 \({\mu}m\)) |
cm |
|
Chem_WetAeroRadiMDUST5 [3] |
Wet aerosol radius for mineral dust (1.5 \({\mu}m\)) |
cm |
|
Chem_WetAeroRAdiMDUST6[3] |
Wet aerosol radius for mineral dust (2.5 \({\mu}m\)) |
cm |
|
Chem_WetAeroRadiMDUST7[3] |
Wet aerosol radius for mineral dust (4.0 \({\mu}m\)) |
cm |
|
Chem_WetAeroRadiSULF[3] |
Wet aerosol radius for sulfate aerosol |
cm |
|
Chem_WetAeroRadiBC[3] |
Wet aerosol radius for black carbon |
cm |
|
Chem_WetAeroRadiOC[3] |
Wet aerosol radius for organic carbon |
cm |
|
Chem_WetAeroRadiSSA[3] |
Wet aerosol radius for sea salt aerosol, accumulation mode |
cm |
|
Chem_WetAeroRadiSSC[3] |
Wet aerosol radius for sea salt aerosol, coarse mode |
cm |
|
Chem_WetAeroRadiBGSULF |
Wet aerosol radius for background stratospheric sulfate |
cm |
|
Chem_WetAeroRadiICEI |
Wet aerosol Radius for irregular ice cloud |
cm |
|
Chem_KhetiSLAN2O5H2O |
Sticking coefficient for N2O5 + H2O rxn |
1 |
|
Chem_KhetiSLAN2O5HCl |
Sticking coefficient for N2O5 + HCl rxn |
1 |
|
Chem_KhetiSLACLNO3H2O |
Sticking coefficient for ClNO3 + H2O rxn |
1 |
|
Chem_KhetiSLACLNO3HCL |
Sticking coefficient for ClNO3 + HCl rxn |
1 |
|
Chem_KhetiSLACLNO3HBR |
Sticking coefficient for ClNO3 + HBr rxn |
1 |
|
Chem_KhetiSLABRNO3H2O |
Sticking coefficient for BrNO3 + H2O rxn |
1 |
|
Chem_KhetiSLABRNO3HCL |
Sticking coefficient for BrNO3 + HCl rxn |
1 |
|
Chem_KhetiSLAHOCLHCL |
Sticking coefficient for HOCl + HCl rxn |
1 |
|
Chem_KhetiSLAHOCLHBR |
Sticking coefficient for HOCl + HBr rxn |
1 |
|
Chem_KhetiSLAHOBRHCL |
Sticking coefficient for HOBr + HCl rxn |
1 |
|
Chem_KhetiSLAHOBRHBR |
Sticking coefficient for HOBr + HBr rxn |
1 |
Notes for the Aerosols colletion
Defined for aerosol-only and fullchem simulations.
Only defined for fullchem simulation with complex SOA species.
BoundaryConditions
The BoundaryConditions diagnostic collection contains advected species concentrations (archived from a global simulation) that will be used by GEOS-Chen Classic nested-grid simulations as transport boundary conditions.
Sample definition section for HISTORY.rc
BoundaryConditions.template: '%y4%m2%d2_%h2%n2z.nc4',
BoundaryConditions.frequency: 00000000 030000
BoundaryConditions.duration: 00000100 000000
BoundaryConditions.mode: 'instantaneous'
BoundaryConditions.fields: 'SpeciesBC_?ADV?',
::
List of diagnostic fields in the BoundaryConditions collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
SpeciesBC_<name|wc>[5] |
Advected species concentrations used as boundary conditions GEOS-Chem Classic nested-grid simulations |
mol/mol dry air |
?ADV? |
Footnotes
This diagnostic is only for use with GEOS-Chem Classic.
Budget
The Budget diagnostic collection is a 2D diagnostic containing the mass tendencies per grid cell, in kg/s, for each species within a region of the column and across each GEOS-Chem component. The diagnostic is calculated by taking the difference in vertically summed column mass before and after major GEOS-Chem components.
There are three pre-defined column regions defined for this diagnostic: troposphere-only, PBL-only, and full column, as well as a user-defined column region. By post-processing this diagnostic you can calculate global mass change or mass change across regions by summing the diagnostic values for the relevant grid cells. You can also retrieve the mass change across a longer chunk of time by multiplying the time-averaged output by the number of seconds in the averaging period.
While there are seven major components in GEOS-Chem, there are only six implemented for the budget diagnostics: chemistry, mixing, convection, transport (GEOS-Chem Classic only), wet deposition, and combined emissions and dry deposition.
Attention
Emissions and dry deposition components are combined together for this diagnostic because of the way they are applied at the same time. Furthermore, if using non-local PBL mixing then the emissions and dry deposition budget diagnostic will not capture all fluxes from these sources and sinks. This is because emissions and dry deposition tendencies below the PBL are applied within mixing instead. When using full mixing, however, mixing and emissions/dry deposition budget diagnostics are fully separated.
Sample definition section for HISTORY.rc
Budget.template: '%y4%m2%d2_%h2%n2z.nc4',
Budget.frequency: 00000100 000000
Budget.duration: 00000100 000000
Budget.mode: 'time-averaged'
Budget.fields: 'BudgetChemistryFull_?ADV? ',
'BudgetChemistryPBL_?ADV? ',
'BudgetChemistryTrop_?ADV? ',
'BudgetEmisDepFull_?ADV? ',
'BudgetEmisDepTrop_?ADV? ',
'BudgetEmisDepPBL_?ADV? ',
'BudgetTransportFull_?ADV? ',
'BudgetTransportTrop_?ADV? ',
'BudgetTransportPBL_?ADV? ',
'BudgetMixingFull_?ADV? ',
'BudgetMixingTrop_?ADV? ',
'BudgetMixingPBL_?ADV? ',
'BudgetConvectionFull_?ADV?',
'BudgetConvectionTrop_?ADV?',
'BudgetConvectionPBL_?ADV? ',
'BudgetWetDepFull_?WET? ',
'BudgetWetDepTrop_?WET? ',
'BudgetWetDepPBL_?WET? ',
::
List of diagnostic fields in the Budget collection
Diagnostic field |
Mass tendency (kg/s) across … |
Wildcard |
---|---|---|
BudgetChemistryFull_<name|wc> |
Chemistry (full atmosphere) |
?ADV? |
BudgetChemistryLevs1to35_<name|wc>[6] |
Chemistry (fixed level range) |
?ADV? |
BudgetChemistryPBL_<name|wc> |
Chemistry (PBL only) |
?ADV? |
BudgetChemistryTrop_<name|wc> |
Chemistry (troposphere only) |
?ADV? |
BudgetConvectionFull_<name|wc> |
Convection (full atmosphere) |
?ADV? |
BudgetConvectionLevs1to35_<name|wc> [6] |
Convection (fixed level range) |
?ADV? |
BudgetConvectionPBL_<name|wc> |
Convection (PBL only) |
?ADV? |
BudgetConvectionTrop_<name|wc> |
Convection (troposphere only) |
?ADV? |
BudgetEmisDepFull_<name|wc>[7] |
Emissions & dry deposition (full atmosphere) |
?ADV? |
Emissions & dry deposition (fixed level range) |
?ADV? |
|
BudgetEmisDepPBL_<name|wc>[7] |
Emissions & dry deposition (PBL only) |
?ADV? |
BudgetEmisDepTrop_<name|wc>[7] |
Emissions & dry deposition (troposphere only) |
?ADV? |
BudgetMixingFull_<name|wc> [8] |
PBL mixing (full atmosphere) |
?ADV? |
PBL mixing (full atmosphere) (fixed level range) |
?ADV? |
|
BudgetMixingPBL_<name|wc>[8] |
PBL mixing (PBL only) |
?ADV? |
BudgetMixingTrop_<name|wc>[8] |
PBL mixing (troposphere only) |
?ADV? |
BudgetTransportFull_<name|wc>[9] |
Transport (full attmosphere) |
?ADV? |
Transport (fixed level range) |
?ADV? |
|
BudgetTransportPBL_<name|wc>[9] |
Transport (PBL only) |
?ADV? |
BudgetTransportTrop_<name|wc>[9] |
Transport (troposphere only) |
?ADV? |
BudgetWetDepFull_<name|wc> |
Wet deposition (full atmosphere) |
?WET? |
BudgetWetDepLevs1to35_<name|wc>[6] |
Wet deposition (fixed level range) |
?WET? |
BudgetWetDepPBL_<name|wc> |
Wet deposition (PBL only) |
?WET |
Notes for the Budget collection
These diagnostic quantities allow you to compute mass
tendencies in a fixed level range. The lower level and upper level of
the range is specified in the diagnostic name
(LevsXtoY
). Levels 1 to 35 (surface to
approximately the tropopause) are the default settings.
The emissions and dry deposition budget diagnostics will not
capture all fluxes if using the non-local PBL mixing scheme since
these tendencies are applied within mixing in vdiff_mod.F90
below the PBL. When using full mixing, however, mixing and emissions/dry
deposition are fully separated.
The mixing budget diagnostics includes the application of emissions and dry deposition below the PBL if using the non-local PBL mixing scheme (vdiff).
GEOS-Chem Classic only
Carbon
The Carbon collection contains diagnostic fields specific to the GEOS-Chem carbon gases simulation.
Sample definition section for HISTORY.rc
Carbon.template: '%y4%m2%d2_%h2%n2z.nc4',
Carbon.frequency: 00000100 000000
Carbon.duration: 00000100 000000
Carbon.mode: 'time-averaged'
Carbon.fields: 'OHconcAfterChem',
'ProdCOfromCH4 ',
'ProdCOfromNMVOC',
'ProdCO2fromCO ',
::
List of diagnostic fields in the Carbon collection
Diagnostic field |
Description |
Units |
---|---|---|
OHconcAfterChem |
OH concentration immediately after chemistry |
molec/cm3 |
ProdCOfromCH4 |
Production of CO from CH4 |
molec/cm3 |
ProdCOfromNMVOC |
Production of CO from non-methane VOCs |
molec/cm3 |
ProdCO2fromCO |
Production of CO2 from CO oxidation |
molec/cm3 |
CH4
The CH4 collection contains diagnostics for loss of CH4 and OH concentration for the CH4 simulation.
Attention
This simulation is slated to be replaced by the GEOS-Chem carbon gases simulation. when this happens, the CH4 collection will be replaced by the Carbon collection.
Sample definition section for HISTORY.rc
CH4.template: '%y4%m2%d2_%h2%n2z.nc4',
CH4.frequency: 00000100 000000
CH4.duration: 00000100 000000
CH4.mode: 'time-averaged'
CH4.fields: 'OHconcAfterChem ',
'LossCH4byClinTrop ',
'LossCH4byOHinTrop ',
'LossCH4inStrat ',
::
List of diagnostic fields in the CH4 collection
Diagnostic field |
Description |
Units |
---|---|---|
LossCH4byClinTrop |
Loss of CH4 by raction with Cl in the troposphere |
kg/s |
LossCH4byOHinTrop |
Loss of CH4 by raction with OH in the troposphere |
kg/s |
LossCH4inStrat |
Loss of CH4 in the stratosphere |
kg/s |
OHconcAfterChem |
OH concentration after chemistry |
kg/s |
CloudConvFlux
The CloudConvFlux collection contains diagnostics for mass fluxes in cloud convection.
Sample definition section for HISTORY.rc
CloudConvFlux.template: '%y4%m2%d2_%h2%n2z.nc4',
CloudConvFlux.frequency: 00000100 000000
CloudConvFlux.duration: 00000100 000000
CloudConvFlux.mode: 'time-averaged'
CloudConvFlux.fields: 'CloudConvFlux_?ADV?',
::
List of diagnostic fields in the CloudConvFlux collection
Diagnostic field |
Description |
Units |
Wildcards |
---|---|---|---|
CloudConvFlux_<name|wc> |
Mass change due to cloud convection |
kg/s |
?ADV? |
CO
The CO collection contains diagnostic fields for the GEOS-Chem tagged CO simulation.
Attention
The tagged CO simulation is slated to be replaced by the GEOS-Chem carbon gases simulation. When this happens, the CO collection will be replaced with the Carbon collection.
CO.template: '%y4%m2%d2_%h2%n2z.nc4',
CO.frequency: 00000100 000000
CO.duration: 00000100 000000
CO.mode: 'time-averaged'
CO.fields: 'ProdCOfromCH4 ',
'ProdCOfromNMVOC',
::
List of diagnostic fields in the CO collection
Diagnostic field |
Description |
Units |
---|---|---|
ProdCOfromCH4 |
Production of CO from CH4 |
kg/s |
ProdCOfromNMVOC |
Production of CO from non-methane VOCs |
kg/s |
CO2
The CO2 collection contains diagnostic outputs from the GEOS-Chem CO2 simulation.
Attention
The CO2 simulation is slated to be replaced by the new GEOS-Chem carbon gases simulation. When this happens, the CO2 collection will be replaced with the Carbon collection.
Note
Several other diagnostics for the CO2 simulation are archived via HEMCO diagnostics.
Sample definition section for HISTORY.rc
CO2.template: '%y4%m2%d2_%h2%n2z.nc4',
CO2.frequency: 00000100 000000
CO2.duration: 00000100 000000
CO2.mode: 'time-averaged'
CO2.fields: 'ProdCO2fromCO',
::
List of diagnostic fields in the CO2 collection
Diagnostic field |
Description |
Units |
---|---|---|
ProdCO2fromCO |
Chemical production of CO2 from CO oxidation |
kg/m2/s |
ConcAboveSfc
The ConcAboveSfc diagnostic collection uses dry deposition quantities (surface resistance, dry deposition velocity) to compute the species concentration of O3 and HNO3 at a given altitude (such as 10 m) above the surface. This will facilitate comparison between GEOS-Chem and observational networks (e.g. CASTNET), which often place instruments above the canopy at approx. 10m height.
Attention
If dry deposition is turned off in your simulation, then you must disable this collection, or else your run will stop with an error.
Sample definition section for HISTORY.rc
ConcAboveSfc.template: '%y4%m2%d2_%h2%n2z.nc4',
ConcAboveSfc.mode: 'instantaneous'
ConcAboveSfc.fields: 'DryDepRaALT1 ',
'DryDepVelForALT1_?DRYALT?',
'SpeciesConcALT1_?DRYALT? ',
::
List of diagnostic fields in the ConcAboveSfc collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
DryDepRaALT1[10] |
Dry deposition aerodynamic resistance at ALT1 meters above the surface |
s/cm |
|
Dry deposition velocity of species tagged with the ?DRYALT?wildcard |
cm/s |
?DRYALT? |
|
Concentration of species tagged with the ?DRYALT? wildcard |
mol/mol dry air |
?DRYALT? |
Notes about the ConcAboveSfc collection
Replace ALT1
with the altitude in meters above the
surface at which you would like these quantities computed.
For example: DryDepVelFor10m_?DRYALT?
, etc.
Currently the ?DRYALT?
species are O3 and HNO3.
ConcAfterChem
The ConcAfterChem collection contains diagnostics for OH, HO2, etc. species immediately upon exiting the chemical solver.
Sample definition section for HISTORY.rc
ConcAfterChem.template: '%y4%m2%d2_%h2%n2z.nc4',
ConcAfterChem.frequency: 00000100 000000
ConcAfterChem.duration: 00000100 000000
ConcAfterChem.mode: 'time-averaged'
ConcAfterChem.fields: 'OHconcAfterChem ',
'HO2concAfterChem',
'O1DconcAfterChem',
'O3PconcAfterChem',
'O3concAfterChem ',
'RO2concAfterChem',
::
List of diagnostic fields in the ConcAfterChem collection
Diagnostic field |
Description |
Units |
---|---|---|
HO2concAfterChem |
HO2 immediately after exiting the chemical solver |
mol/mol |
O1DconcAfterChem |
O1D immediately after exiting the chemical solver |
molec/cm3 |
O3concAfterChem |
O3 immediately after exiting the chemical solver |
molec/cm3 |
O3PconcAfterChem |
O3P immediately after exiting the chemical solver |
molec/cm3 |
OHconcAfterChem |
OH immediately after exiting the chemical solver |
molec/cm3 |
RO2concAfterChem |
RO2 immediately after exiting the chemical solver |
molec/cm3 |
DryDep
The DryDep collection contains diagnostics for the flux and velocity of each species lost to dry-deposition.
Sample definition section for HISTORY.rc
DryDep.template: '%y4%m2%d2_%h2%n2z.nc4',
DryDep.frequency: 00000100 000000
DryDep.duration: 00000100 000000
DryDep.mode: 'time-averaged'
DryDep.fields: 'DryDepVel_?DRY?',
'DryDepMix_?DRY?',
'DryDepChm_?DRY?',
'DryDep_?DRY? ',
::
List of diagnostic fields in the DryDep collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
DryDep_<name|wc> |
Total dry deposition flux |
molec/cm2/s |
?DRY? |
DryDepChm_<name|wc> |
Dry deposition flux (computed in chemistry) |
molec/cm2/s |
?DRY? |
DryDepMix_<name|wc> |
Dry deposition flux (computed in the PBL) |
molec/cm2/s |
?DRY? |
DryDepVel_<name|wc> |
Dry deposition velocity |
cm/s |
?DRY? |
JValues
The JValues collection contains diagnostics for photolysis rates for various chemical species, obtained from the photolysis mechanism.
Sample definition section for HISTORY.rc
JValues.template: '%y4%m2%d2_%h2%n2z.nc4',
JValues.frequency: 00000100 000000
JValues.duration: 00000100 000000
JValues.mode: 'time-averaged'
JValues.fields: 'Jval_?PHO?',
'JvalO3O1D ',
'JvalO3O3P ',
::
List of diagnostic fields in the JValues collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
Jval_<name|wc> |
Photolysis rates |
1/s |
?PHO? |
JvalO3O1D |
Photolysis rate of O3 \(\rightarrow\) O1D |
1/s |
|
JvalO3O3P |
Photolysis rate of O3 \(\rightarrow\) O3P |
1/s |
KppARDiags
The KppARDiags collection contains diagnostics for the KPP Rosenbrock solver with mechanism auto-reduction. You may leave this collection disabled unless you are interested in assessing the solver’s performance.
Sample definition section for HISTORY.rc
KppARDiags.template: '%y4%m2%d2_%h2%n2z.nc4',
KppARDiags.frequency: 00000100 000000
KppARDiags.duration: 00000100 000000
KppARDiags.mode: 'time-averaged'
KppARDiags.fields: 'KppAutoReducerNVAR',
'KppAutoReduceThres',
'KppcNONZERO ',
::
List of diagnostic fields in the KppARDiags collection
Diagnostic field |
Description |
Units |
---|---|---|
KppAutoReducerNVAR |
Number of species ( |
count |
KppAutoReduceThres |
Auto-reduction threshold |
molec/cm3/s |
KppcNONZERO |
Number of nonzero elements
( |
count |
KppDiags
The KppDiags collection contains KPP solver diagnostics. You may leave this collection disabled unless you are interested in assessing the solver’s performance.
Sample definition section for HISTORY.rc
KppDiags.template: '%y4%m2%d2_%h2%n2z.nc4',
KppDiags.frequency: 00000100 000000
KppDiags.duration: 00000100 000000
KppDiags.mode: 'time-averaged'
KppDiags.fields: 'KppIntCounts',
'KppJacCounts',
'KppTotSteps ',
'KppAccSteps ',
'KppRejSteps ',
'KppLuDecomps',
'KppSubsts ',
'KppSmDecomps',
::
List of diagnostic fields in the KppDiags collection
Diagnostic field |
Description |
Units |
---|---|---|
KppAccSteps |
Number of accepted integration timesteps |
count |
KppIntCounts |
Number of times the KPP integrator was called |
count |
KppJacCounts |
Number of times the KPP Jacobian matrix was constructed |
count |
KppLuDecomps |
Number of LU decompositions performed |
count |
KppSmDecomps[12] |
Number of singular matrix decompositions performed |
count |
KppSubsts |
Number of matrix substitutions performed (both forward & backward substitutions) |
count |
KppRejSteps |
Number of rejected integration timesteps |
count |
KppTotSteps |
Total number of integration timesteps |
count |
Footnotes
For Rosenbrock solvers, KppSmDecomps will be zero everywhere, because the Rosenbrock method utilizes LU decomposition instead of singular matrix decomposition.
LevelEdgeDiags
The LevelEdgeDiags collection contains diagnostics for quantities (mostly met fields) that are defined on the vertical edges of each grid box. According to the COARDS convention, all of the data variables in a netCDF file must be defined with the same vertical dimension.
Sample definition section for HISTORY.rc
LevelEdgeDiags.template: '%y4%m2%d2_%h2%n2z.nc4',
LevelEdgeDiags.frequency: 00000100 000000
LevelEdgeDiags.duration: 00000100 000000
LevelEdgeDiags.mode: 'time-averaged'
LevelEdgeDiags.fields: 'Met_CMFMC ',
'Met_PEDGE ',
'Met_PEDGEDRY',
'Met_PFICU ',
'Met_PFILSAN ',
'Met_PFLCU ',
'Met_PFLLSAN ',
::
List of diagnostic fields in the LevelEdgeDiags collection
Diagnostic field |
Description |
Units |
---|---|---|
Met_CMFMC |
Upward moist convective mass flux |
kg/m2/s |
Met_PEDGE |
Surface pressure at level edges (based on moist air) |
hPa |
Met_PEDGEDRY |
Surface pressure at level edges (based on dry air) |
hPa |
Met_PFICU |
3d flux of ice convective precipitation |
kg/m2/s |
Met_PFILSAN |
3d flux of ice non-convective precipitation |
kg/m2/s |
Met_PFLCU |
3d flux of liquid convective precipitation |
kg/m2/s |
Met_PFLLSAN |
3d flux of liquid non-convective precipitation |
kg/m2/s |
MercuryChem
The MercuryChem collection contains concentrations and prod/loss diagnostic outputs for the GEOS-Chem mercury simulation.
Sample definition section for HISTORY.rc
MercuryChem.template: '%y4%m2%d2_%h2%n2z.nc4',
MercuryChem.frequency: ${RUNDIR_HIST_TIME_AVG_FREQ}
MercuryChem.duration: ${RUNDIR_HIST_TIME_AVG_DUR}
MercuryChem.mode: 'time-averaged'
MercuryChem.fields: 'HgBrAfterChem ',
'HgClAfterChem ',
'HgOHAfterChem ',
'HgBrOAfterChem ',
'HgClOAfterChem ',
'HgOHOAfterChem ',
'Hg2GToHg2P ',
'Hg2PToHg2G ',
'Hg2GasToHg2StrP',
'Hg2GasToSSA ',
::
List of diagnostic fields in the MercuryChem collection
Diagnostic field |
Description |
Units |
---|---|---|
Hg2GToHg2P |
Hg2 gas transferred to Hg2P |
molec/cm3/s |
Hg2GasToHg2StrP |
Hg2 gas transferred to Hg2StrP |
molec/cm3/s |
Hg2GasToSSA |
Hg2 gas transferred to sea salt aerosol |
molec/cm3/s |
Hg2PToHg2G |
Hg2P transferred to Hg2 gas |
molec/cm3/s |
HgBrAfterChem |
HgBr concentration immediately after chemistry |
mol/mol |
HgBrOAfterChem |
HgBrO concentration immediately after chemistry |
mol/mol |
HgClAfterChem |
HgCl concentration immediately after chemistry |
mol/mol |
HgClOAfterChem |
HgClO concentration immediately after chemistry |
mol/mol |
HgOHAfterChem |
HgOH concentration immediately after chemistry |
mol/mol |
HgOHOAfterChem |
HgOHO concentration immediately after chemistry |
mol/mol |
MercuryEmis
The MercuryEmis collection contains emission diagnostics for the GEOS-Chem mercury simulation.
Note
Several other mercury emission diagnostics are archived via HEMCO diagnostics.
Sample definition section for HISTORY.rc
MercuryEmis.template: '%y4%m2%d2_%h2%n2z.nc4',
MercuryEmis.frequency: 00000100 000000
MercuryEmis.duration: 00000100 000000
MercuryEmis.mode: 'time-averaged'
MercuryEmis.fields: 'EmisHg0land ',
'EmisHg0ocean',
'EmisHg0soil ',
'EmisHg0snow ',
::
List of diagnostic fields in the MercuryEmis collection
Diagnostic field |
Description |
Units |
---|---|---|
EmisHg0land |
Re-emission of Hg0 from land |
kg/s |
EmisHg0ocean |
Emissions of Hg0 from oceans |
kg/s |
EmisHg0snow |
Emission of Hg0 from snowpack |
kg/s |
EmisHg0soil |
Emissions of Hg0 from soils |
kg/s |
MercuryOcean
The MercuryOcean collection contains diagnostics from the mercury ocean model, used in the GEOS-Chem mercury simulation.
Sample definition section for HISTORY.rc
MercuryOcean.template: '%y4%m2%d2_%h2%n2z.nc4',
MercuryOcean.frequency: 00000000 040000
MercuryOcean.duration: 00000000 040000
MercuryOcean.mode: 'time-averaged'
MercuryOcean.fields: 'FluxHg0fromAirToOcean ',
'FluxHg0fromOceanToAir ',
'FluxHg2HgPfromAirToOcean',
'FluxHg2toDeepOcean ',
'FluxOCtoDeepOcean ',
'MassHg0inOcean ',
'MassHg2inOcean ',
'MassHgPinOcean ',
'MassHgTotalInOcean ',
::
List of diagnostic fields in the MercuryOcean collection
Diagnostic field |
Description |
Units |
---|---|---|
FluxHg0fromAirToOcean |
Deposition flux of Hg0 from the atmosphere to the ocean |
kg/s |
FluxHg0fromOceanToAir |
Volatization flux of Hg0 from the ocean to the atmosphere |
kg/s |
FluxHg2HgPfromAirToOcean |
Deposition flux of Hg2 + HgP from atmosphere to ocean |
kg/s |
FluxHg2toDeepOcean |
Flux of Hg2 sunk to the deep ocean |
kg/s |
MassHg0inOcean |
Total mass of oceanic Hg0 |
kg |
MassHg2inOcean |
Total mass of oceanic Hg2 |
kg |
MassHgPinOcean |
Total mass of oceanic HgP |
kg |
MassHgTotalInOcean |
Total mass of all organic mercury |
kg |
Metrics
The Metrics collection contains diagnostics for computing OH metrics from a GEOS-Chem full chemistry simulation (needed for benchmarking).
To compute the OH metrics, you must run the Python script
metrics.py
that ships with each fullchem or CH4 run directory.
Sample definition section for HISTORY.rc
Metrics.template: '%y4%m2%d2_%h2%n2z.nc4',
Metrics.frequency: 'End',
Metrics.duration: 'End',
Metrics.mode: 'time-averaged'
Metrics.fields: 'AirMassColumnFull ',
'LossOHbyCH4columnTrop ',
'LossOHbyMCFcolumnTrop ',
'OHwgtByAirMassColumnFull',
::
List of diagnostic fields in the Metrics collection
Diagnostic field |
Description |
Units |
---|---|---|
AirMassColumnFull |
Air mass column (full atmosphere) |
kg |
LossOHbyCH4columnTrop |
Loss rate of CH4 by OH (tropospheric column sums) |
molec/cm3 |
LossOHbyMCFcolumnTrop |
Loss rate of CH4 by CH3CCl3 aka MCF (tropospheric column sums) |
molec/cm3 |
OHwgtByAirMassColumnFull |
Airmass-weighted OH concentration (full atmosphere column sums) |
kg air/kg OH/m3 |
ProdLoss
The ProdLoss collection contains chemical production and loss rates.
Sample definition section for HISTORY.rc
ProdLoss.template: '%y4%m2%d2_%h2%n2z.nc4',
ProdLoss.frequency: 00000100 000000
ProdLoss.duration: 00000100 000000
ProdLoss.mode: 'time-averaged'
ProdLoss.fields: 'Prod_?PRD? ',
'ProdBCPIfromBCPO ',
'ProdOCPIfromOCPO ',
'ProdHMSfromSO2andHCHOinCloud',
'ProdSO2andHCHOfromHMSinCloud',
'ProdSO4fromHMSinCloud ',
'ProdSO4fromH2O2inCloud ',
'ProdSO4fromO2inCloudMetal ',
'ProdSO4fromO3inCloud ',
'ProdSO4fromO3inSeaSalt ',
'ProdSO4fromHOBrInCloud ',
'ProdSO4fromSRO3 ',
'ProdSO4fromSRHObr ',
'ProdSO4fromO3s ',
'Loss_?LOS? ',
'LossHNO3onSeaSalt ',
'ProdCOfromCH4 ',
'ProdCOfromNMVOC ',
::
List of diagnostic fields in the ProdLoss collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
Loss_<name|wc> |
Chemical loss for a given species or family |
molec/cm3 |
?LOS? |
LossHNO3onSeaSalt[15] |
L(HNO3) on sea salt aerosols |
kg S/s |
|
Prod_<name|wc> |
Chemical production for a given species or family |
molec/cm3 |
?PRD? |
ProdBCPIfromBCPO[15] |
Production of hydrophilic BC from hydrophobic BC |
kg |
|
ProdCOfromCH4[16] |
P(CO) from CH4 |
molec/cm3 |
|
ProdCOfromNMVOC[16] |
P(CO) from NMVOCs SO3– loss by OH |
molec/cm3 |
|
ProdOCPIfromOCPO[15] |
Production of hydrophilic OC from hydrophobic OC |
kg |
|
ProdMSAfromDMS[14] |
P(MSA) from DMS |
kg S/s |
|
ProdNITfromHNO3uptakeOnDust[13] |
P(NIT) from HNO3 uptake on dust aerosols |
kg N/s |
|
ProdSO2fromDMS[14] |
Total P(SO2) from DMS |
kg S/s |
|
ProdSO2fromDMSandNO3[14] |
P(SO2) from DMS + NO3 |
kg S/s |
|
ProdSO2fromDMSandOH[14] |
P(SO2) from DMS + OH |
kg S/s |
|
ProdSO2fromOxidationOnDust[13] |
P(SO2) from DMS+OH on dust aerosols |
kg S/s |
|
ProdSO4fromGasPhase[14] |
P(SO4) in gas phase |
kg S/s |
|
ProdSO4fromH2O2inCloud[15] |
P(SO4) from aqueous oxidation of H2O2 in clouds |
kg S/s |
|
ProdSO4fromHOBrinCloud[16] |
P(SO4) from aqueous oxidation of HOBr in clouds |
kg S/s |
|
ProdSO4fromO2inCloudMetal[15] |
P(SO4) from aqueous oxidation of O2 from metals in cloud |
kg S/s |
|
ProdSO4fromO3inCloud[15] |
P(SO4) from aqueous oxidation of O3 in clouds |
kg S/s |
|
ProdSO4fromO3inSeaSalt[15] |
P(SO4) from O3 in sea salt |
kg S/s |
|
ProdSO4fromO3s[15] |
P(SO4) from aqueous phase SO3– loss by OH |
kg S/s |
|
ProdSO4fromSRHOBR[16] |
P(SO4) from sulfur production rate of HOBr + O3 |
kg S/s |
|
ProdSO4fromSRO3[15] |
P(SO4) from sulfur production rate of O3 |
kg S/s |
|
ProdSO4fromUptakeOfH2SO4g[13] |
P(SO4) from H2SO4 uptake on dust aerosols |
kg S/s |
Notes for the ProdLoss collection
Only defined for fullchem simulation with aciduptake on dust.
Only defined for the aerosol-only simulation.
Defined for aerosol-only and fullchem simulations.
Only defined for fullchem simulations.
RadioNuclide
The RadioNuclide collection contains diagnostic outputs for radionuclide species in the GEOS-Chem TransportTracers simulation.
Note
Emissions of Rn222, Be7, and Be10 species are archived to diagnostic output by HEMCO diagnostics, and are thus not contained in this collection.
Sample definition section for HISTORY.rc
RadioNuclide.template: '%y4%m2%d2_%h2%n2z.nc4',
RadioNuclide.format: 'CFIO',
RadioNuclide.frequency: 00000100 000000
RadioNuclide.duration: 00000100 000000
RadioNuclide.mode: 'time-averaged'
RadioNuclide.fields: 'PbFromRnDecay ',
'RadDecay_Rn222 ',
'RadDecay_Pb210 ',
'RadDecay_Pb210s',
'RadDecay_Be7 ',
'RadDecay_Be7s ',
'RadDecay_Be10 ',
'RadDecay_Be10s ',
::
List of diagnostic fields in the RadioNuclide collection
Diagnostic field |
Description |
Units |
---|---|---|
PbFromRnDecay |
Pb210 created from radioactive decay |
kg/s |
RadDecay_Be7 |
Loss of Be7 due to radioactive decay |
kg/s |
RadDecay_Be7s |
Loss of Be7 (produced in the stratosphere) due to radioactive decay |
kg/s |
RadDecay_Be10 |
Loss of Be10 due to radioactive decay |
kg/s |
RadDecay_Be10s |
Loss of Be10 (produced in the stratosphere) due to radioactive decay |
kg/s |
RadDecay_Pb210 |
Loss of Pb210 due to radioactive decay |
kg/s |
RadDecay_Pb210s |
Loss of Pb210 (produced in the stratosphere) due to radioactive decay |
kg/s |
RadDecay_Rn222 |
Loss of Rn222 due to readioactive decay |
kg/s |
Restart
The Restart diagnostic collection contains fields for saving out to
the GEOS-Chem Classic restart file. This type of diagnostic output is used in
all GEOS-Chem Classic simulations; therefore, we have listed Restart first in
the HISTORY.rc
files that ship with each GEOS-Chem run directory. Note
that GCHP restart files are not handled by the MAPL History component and therefore
do not appear in HISTORY.rc
.
Note
The restart file will be created in the Restarts/
subdirectory of the run directory, not in OutputDir/
.
Sample definition section for HISTORY.rc
Restart.filename: './GEOSChem.Restart.%y4%m2%d2_%h2%n2z.nc4',
Restart.frequency: 'End',
Restart.duration: 'End',
Restart.mode: 'instantaneous'
Restart.fields: 'SpeciesRst_?ALL? ',
'Chem_H2O2AfterChem ',
'Chem_SO2AfterChem ',
'Chem_DryDepNitrogen',
'Chem_WetDepNitrogen',
'Chem_KPPHvalue ',
'Met_DELPDRY ',
'Met_PS1WET ',
'Met_PS1DRY ',
'Met_SPHU1 ',
'Met_TMPU1 ',
::
List of diagnostic fields in the Restart collection
Diagnostic field |
Description |
Units |
---|---|---|
Chm_DryDepNitrogen |
Dry deposited nitrogen |
molec/cm2/s |
Chm_H2O2AfterChem |
Concentration of H2O2 after sulfate chemistry |
v/v |
Chm_KPPHvalue |
H-value for Rosenbrock solver |
unitless |
Chm_SO2AfterChem |
Concentration of SO2 after sulfate chemistry |
v/v |
Chm_StatePSC |
Polar stratospheric cloud type |
count |
Chm_WetDepNitrogen |
Wet deposited nitrogen |
molec/cm2/s |
Met_DELPDRY |
Delta-pressure across grid box (dry air) |
hPa |
Met_PS1WET |
Wet surface pressure at dt start |
hPa |
Met_PS1DRY |
Dry surface pressure at dt start |
hPa |
Met_SPHU1 |
Instantaneous specific humidity at time=T |
g kg-1 |
Met_TMPU1 |
Instantaneous temperature at time=T |
K |
SpeciesRst_?ALL? |
Instantaneous chemical species concentrations for use in starting subsequent GEOS-Chem simulations |
mol/mol dry air |
RRTMG
The RRTMG collection contains radiative flux diagnostics computed by the RRTMG radiative transfer model. You can leave this collection disabled unless your simulation uses RRTMG.
Note
You may compute RRTMG diagnostic quantities at up to 3 wavelengths
(WL1
, WL2
, WL3
). Specify the
wavelengths in this menu of the geoschem_config.yml
file:
rrtmg_rad_transfer_model:
activate: true
aod_wavelengths_in_nm:
- 550 700 1000
... etc ...
GEOS-Chem will replace the tokens WL1
, WL2
,
WL3
in diagnostic field names with the corresponding
wavelength. For example:
RadAODWL1_SU
RadAODWL2_SU
RadAODWL3_SU
will be saved to the GEOSChem.RRTMG.YYYYMMDD_hhmmz.nc4
file(s) under these names:
RadAOD550nm_SU
RadAOD700nm_SU
RadAOD1000nm_SU
Sample definition section for HISTORY.rc
#==============================================================================
# %%%%% THE RRTMG COLLECTION %%%%%
#
# Outputs for different species from the RRTMG radiative transfer model:
# (See http://wiki.geos-chem.org/Coupling_GEOS-Chem_with_RRTMG)
#
# 0=BA (Baseline ) 1=O3 (Ozone ) 2=ME (Methane )
# 3=SU (Sulfate ) 4=NI (Nitrate ) 5=AM (Ammonium )
# 6=BC (Black carbon) 7=OA (Organic aerosol) 8=SS (Sea Salt )
# 9=DU (Mineral dust) 10=PM (All part. matter) 12=ST (Strat aer., UCX only)
#
# NOTES:
# (1) Only request diagnostics you need to reduce the overall run time.
# (2) The ?RRTMG? wildcard includes all output except ST (strat aerosols).
# However, if ST is included explicitly for one diagnostic then it
# will be included for all others that use the wildcard.
# (3) Only enable ST if running with UCX.
# (4) Optics diagnostics have a reduced set of output species (no BASE, O3, ME)
#==============================================================================
RRTMG.template: %y4%m2%d2_%h2%n2z.nc4',
RRTMG.frequency: 00000100 000000
RRTMG.duration: 00000100 000000
RRTMG.mode: time-averaged'
RRTMG.fields: 'RadClrSkyLWSurf_?RRTMG?',
'RadAllSkyLWSurf_?RRTMG?',
'RadClrSkySWSurf_?RRTMG?',
'RadAllSkySWSurf_?RRTMG?',
'RadClrSkyLWTOA_?RRTMG? ',
'RadAllSkyLWTOA_?RRTMG? ',
'RadClrSkySWTOA_?RRTMG? ',
'RadAllSkySWTOA_?RRTMG? ',
'RadAODWL1_?RRTMG? ',
'RadAsymWL1_?RRTMG? ',
::
List of diagnostic fields in the RRTMG collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
DynHeating |
Dynamic heating rate in baseline simulation |
K/day |
?RRTMG? |
DtRad |
Temperature change due to radiative heating |
K |
?RRTMG? |
RadAODWL1_<name|wc> |
Aerosol optical depth computed @ WL1 |
1 |
?RRTMG? |
RadAODWL2_<name|wc> |
Aerosol optical depth computed @ WL2 |
1 |
?RRTMG? |
RadAODWL3_<name|wc> |
Aerosol optical depth computed @ WL3 |
1 |
?RRTMG? |
RadAsymWL1_<name|wc> |
Asymmetry parameter computed @ WL1 |
1 |
?RRTMG? |
RadAsymWL2_<name|wc> |
Asymmetry parameter computed @ WL2 |
1 |
?RRTMG? |
RadAsymWL3_<name|wc> |
Asymmetry parameter computed @ WL3 |
1 |
?RRTMG? |
RadAllSkyLWSurf_<name|wc> |
All-sky longwave radiation at the surface |
W/m2 |
?RRTMG? |
RadAllSkyLWTOA_<name|wc> |
All-sky longwave radiation at top-of-atmosphere |
W/m2 |
?RRTMG? |
RadAllSkyLWTrop_<name|wc> |
All-sky longwave radiation at the tropopause |
W/m2 |
?RRTMG? |
RadAllSkySWSurf_<name|wc> |
All-sky shortwave radiation at the surface |
W/m2 |
?RRTMG? |
RadAllSkySWTOA_<name|wc> |
All-sky shortwave radiation at top-of-atmosphere |
W/m2 |
?RRTMG? |
RadAllSkySWTrop_<name|wc> |
All-sky shortwave radiation at the tropopause |
W/m2 |
?RRTMG? |
RadClrSkyLWSurf_<name|wc> |
Clear-sky longwave radiation at the surface |
W/m2 |
?RRTMG? |
RadClrSkyLWTOA_<name|wc> |
Clear-sky longwave radiation at top-of-atmosphere |
W/m2 |
?RRTMG? |
RadClrSkyLWTrop_<name|wc> |
Clear-sky longwave radiation at the tropopause |
W/m2 |
?RRTMG? |
RadClrSkySWSurf_<name|wc> |
Clear-sky shortwave radiation at the surface |
W/m2 |
?RRTMG? |
RadClrSkySWTOA_<name|wc> |
Clear-sky shortwave radiation at top-of-atmosphere |
W/m2 |
?RRTMG? |
RadClrSkySWTrop_<name|wc> |
Clear-sky shortwave radiation at the tropopause |
W/m2 |
?RRTMG? |
RadSSAWL1_<name|wc> |
Single-scattering albedo computed @ WL1 |
1 |
?RRTMG? |
RadSSAWL2_<name|wc> |
Single-scattering albedo computed @ WL2 |
1 |
?RRTMG? |
RadSSAWL3_<name|wc> |
Single-scattering albedo computed @ WL3 |
1 |
?RRTMG? |
RxnConst
The RxnConst collection contains reaction rate constants from the KPP solver.
# It is best to list individual reactions to avoid using too much memory.
# Reactions should be listed as "RxnConst_EQnnn", where nnn is the reaction
# index as listed in KPP/fullchem/gckpp_Monitor.F90 (pad zeroes as needed).
#
# The units of reaction rate constants vary according to the number of reactants
# in the reaction.
#
# Available for the fullchem simulations.
RxnConst.template: '%y4%m2%d2_%h2%n2z.nc4',
RxnConst.frequency: ${RUNDIR_HIST_TIME_AVG_FREQ}
RxnConst.duration: ${RUNDIR_HIST_TIME_AVG_DUR}
RxnConst.mode: 'time-averaged'
RxnConst.fields: 'RxnConst_EQ001 ',
'RxnConst_EQ002 ',
::
List of diagnostic fields in the RxnRate collection
Units are
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
RxnConst_EQnnn[17] |
Rate constant for KPP
reaction |
?RXN? |
Notes for the RxnRates collection
See the gckpp_Monitor.F90
file to get a numbered list of reactions.
Units are { (cm3/molec)**(nreactants-1) }/s
RxnRates
The RxnRates collection contains reaction rates (aka equation 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.
Sample definition section for HISTORY.rc
#
# 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,
# KPP/carbon/gckpp_Monitor.F90, and 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 ',
::
List of diagnostic fields in the RxnRate collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
RxnRate_EQnnn[19] |
Rate for KPP reaction |
molec/cm3/s |
?RXN? |
Notes for the RxnRates collection
See the gckpp_Monitor.F90
file to get a numbered list of reactions.
SatDiagn
The SatDiagn collection contains diagnostic quantities that will be sampled within a specified local time range. This is to mimic the overpass sampling times of sun-synchronus satellite instruments.
Tip
Set the the hours (local time) for the averaging interval with:
SatDiagn.hrrange: 11.98 15.02
Sample definition section for HISTORY.rc
SatDiagn.template: '%y4%m2%d2_%h2%n2z.nc4',
SatDiagn.frequency: 00000001 000000
SatDiagn.duration: 00000100 000000
SatDiagn.hrrange: 11.98 15.02
SatDiagn.mode: 'time-averaged'
SatDiagn.fields: 'SatDiagnConc_O3 ',
'SatDiagnOH ',
'SatDiagnRH ',
'SatDiagnAirDen ',
'SatDiagnBoxHeight ',
'SatDiagnPEdge ',
'SatDiagnTROPP ',
'SatDiagnPBLHeight ',
'SatDiagnPBLTop ',
'SatDiagnTAir ',
'SatDiagnGWETROOT ',
'SatDiagnGWETTOP ',
'SatDiagnPARDR ',
'SatDiagnPARDF ',
'SatDiagnPRECTOT ',
'SatDiagnSLP ',
'SatDiagnSPHU ',
'SatDiagnTS ',
'SatDiagnPBLTOPL ',
'SatDiagnMODISLAI ',
'SatDiagnWetLossLS_ ',
'SatDiagnWetLossConv_ ',
'SatDiagnJval_ ',
'SatDiagnJvalO3O1D ',
'SatDiagnJvalO3O3P ',
'SatDiagnDryDep_ ',
'SatDiagnDryDepVel_ ',
'SatDiagnOHreactivity ',
'SatDiagnColEmis_ ',
'SatDiagnSurfFlux_ ',
'SatDiagnProd_?PRD? ',
'SatDiagnLoss_?LOS? ',
'SatDiagnRxnRate_EQnnn',
::
List of diagnostic fields in the SatDiagn collection
Diagnostic field |
Description |
Units |
Wildcards |
---|---|---|---|
SatDiagnAirDen |
Air density |
molec/cm3 |
|
SatDiagnBoxHeight |
Grid box height |
m |
|
SatDiagnColEmis_<name|wc> |
Column emissions |
kg/m2/s |
?ADV? |
SatDiagnConc_<name|wc> |
Dry mixing ratio of species |
mol/mol |
?ADV? |
SatDiagnDryDep_<name|wc> |
Dry deposition flux of species |
molec/cm2/s |
?DRY? |
SatDiagnDryDepVel_<name|wc> |
Dry deposition velocity of species |
cm/s |
?DRY? |
SatDiagnGWETROOT |
Root zone soil moisture |
1 |
|
SatDiagnGWETTOP |
Topsoil moisture (or |
1 |
|
SatDiagnJVal_<name|wc> |
Photolysis rate |
1/s |
?PHO? |
SatDiagnJvalO3O1D |
Photolysis rate for O3 \(\rightarrow\) O1D |
1/s |
|
SatDiagnJvalO3O3P |
Photolysis rate for O3 \(\rightarrow\) O1D |
1/s |
|
SatDiagnLoss_<name|wc> |
Chemical loss of species or families |
molec/cm3/s |
?LOS? |
SatDiagnMODISLAI |
MODIS daily LAI |
m2/m2 |
|
SatDiagnOH |
OH number density |
molec/cm3 |
|
SatDiagnOHreactivity |
OH reactivity |
1/s |
|
SatDiagnPARDF |
Diffuse photosynthetically active radiation |
W/m2 |
|
SatDiagnPARDR |
Direct photosynthetically active radiation |
W/m2 |
|
SatDiagnPBLHeight |
PBL Height |
m |
|
SatDiagnPBLTop |
PBL Top |
m |
|
SatDiagnPBLTOPL |
PBL top height |
level |
|
SatDiagnPRECTOT |
Total precipitation at surface |
mm/day |
|
SatDiagnProd_<name|wc> |
Chemical production of species or families |
molec/cm3/s |
?PRD? |
SatDiagnRH |
Relative humidity |
% |
|
SatDiagnRxnRate_EQnnn |
Rate for chemical reaction
|
molec/cm3/s |
?RXN? |
SatDiagnSLP |
Sea level pressure |
hPa |
|
SatDiagnSPHU |
Specific humidity interpolated to current time |
g H2O/kg air |
|
SatDiagnSurfFlux_<name|wc> |
Total surface fluxes (emis - drydep) from surface to top of PBL |
kg/m2/s |
?ADV? |
SatDiagnTAir |
Air temperature |
K |
|
SatDiagnTROPP |
Tropopause pressure |
hPa |
|
SatDiagnTS |
Surface temperature at 2m |
K |
|
SatDiagnWetLossConv_<name|wc> |
Loss of soluble species in convective updrafts |
kg/s |
?WET? |
SatDiagnWetLossLS_<name|wc> |
Loss of soluble species in large-scale precipitation |
kg/s |
?WET? |
SatDiagnEdge
The SatDiagn collection contains diagnostic quantities (placed on level edges) that will be sampled within a specified local time range. This is to mimic the overpass sampling times of sun-synchronus satellite instruments.
Tip
Set the the hours (local time) for the averaging interval with:
SatDiagn.hrrange: 11.98 15.02
Sample definition section for HISTORY.rc
SatDiagnEdge.template: '%y4%m2%d2_%h2%n2z.nc4',
SatDiagnEdge.frequency: 00000001 000000
SatDiagnEdge.duration: 00000100 000000
SatDiagnEdge.hrrange: 11.98 15.02
SatDiagnEdge.mode: 'time-averaged'
SatDiagnEdge.fields: 'SatDiagnConc_PEDGE',
::
List of diagnostic fields in the SatDiagnEdge collection
Diagnostic field |
Description |
Units |
---|---|---|
SatDiagnPEDGE |
Pressure at grid box edges |
hPa |
SpeciesConc
The SpeciesConc diagnostic collection contains species concentrations.
Sample definition section for HISTORY.rc
SpeciesConc.template: '%y4%m2%d2_%h2%n2z.nc4',
SpeciesConc.format: 'CFIO',
SpeciesConc.frequency: 00000100 000000
SpeciesConc.duration: 00000100 000000
SpeciesConc.mode: 'time-averaged'
SpeciesConc.fields: 'SpeciesConcVV_?ALL?',
'SpeciesConcMND_?ALL?',
::
List of diagnostic fields in the SpeciesConc collection
Diagnostic field |
Description |
Units |
Wildcards |
---|---|---|---|
SpeciesConcMND_<name|wc> |
Species concentration |
molec/cm3 |
can be used with all wildcards |
SpeciesConcVV_<name|wc> |
Species concentration |
mol/mol dry air |
can be used with all wildcards |
StateChm
The ‘’’StateChm ‘’’ collection contains quantities from State_Chm, the Chemistry State object (other than the species concentrations, which are stored in the SpeciesConc collection).
Sample definition section for HISTORY.rc
StateChm.template: '%y4%m2%d2_%h2%n2z.nc4',
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',
::
List of diagnostic fields in the StateChm collection
Diagnostic field |
Description |
Units |
---|---|---|
Chm_AcidPurSav |
ISORROPIA acidpur concentration |
M |
Chm_BiSulSav |
ISORROPIA bisulfate general acid concentration |
M |
Chm_fupdateHOBr |
Correction factor for HOBr removal by SO2 grid box (wet air) |
mol/L |
Chm_HplusSav |
ISORROPIA H+ concentration |
M |
Chm_HSO3AQ |
Cloud bisulfite concentration |
mol/L |
Chm_NaRatSav |
ISORROPIA Na+ concentration |
M |
Chm_phCloud |
Cloud pH |
1 |
Chm_phSav |
ISORROPIA aerosol pH |
1 |
Chm_SO3AQ |
Cloud sulfite concentration |
mol/L |
Chm_SulRatSav |
ISORROPIA sulfate concentration |
M |
Chm_SSalk |
Sea salt alkalinity |
|
Chm_WaterSav |
ISORROPIA aerosol water |
\({\mu}g/m3\) |
StateMet
The StateMet collection contains met fields and other derived quantities that are carried in the State_Met object.
Sample definition section for HISTORY.rc
StateMet.template: '%y4%m2%d2_%h2%n2z.nc4',
StateMet.frequency: 00000100 000000
StateMet.duration: 00000100 000000
StateMet.mode: 'time-averaged'
StateMet.fields: 'Met_AD ',
'Met_AIRDEN ',
'Met_AIRVOL ',
'Met_ALBD ',
'Met_AREAM2 ',
'Met_AVGW ',
'Met_BXHEIGHT ',
'Met_ChemGridLev',
'Met_CLDF ',
'Met_CLDFRC ',
'Met_CLDTOPS ',
'Met_DELP ',
'Met_DQRCU ',
'Met_DQRLSAN ',
'Met_DTRAIN ',
'Met_EFLUX ',
'Met_FRCLND ',
'Met_FRLAKE ',
'Met_FRLAND ',
'Met_FRLANDIC ',
'Met_FROCEAN ',
'Met_FRSEAICE ',
'Met_FRSNO ',
'Met_GWETROOT ',
'Met_GWETTOP ',
'Met_HFLUX ',
'Met_LAI ',
'Met_LWI ',
'Met_PARDR ',
'Met_PARDF ',
'Met_PBLTOPL ',
'Met_PBLH ',
'Met_PHIS ',
'Met_PMID ',
'Met_PMIDDRY ',
'Met_PRECANV ',
'Met_PRECCON ',
'Met_PRECLSC ',
'Met_PRECTOT ',
'Met_PS1DRY ',
'Met_PS1WET ',
'Met_PS2DRY ',
'Met_PS2WET ',
'Met_PSC2WET ',
'Met_PSC2DRY ',
'Met_QI ',
'Met_QL ',
'Met_OMEGA ',
'Met_OPTD ',
'Met_REEVAPCN ',
'Met_REEVAPLS ',
'Met_SLP ',
'Met_SNODP ',
'Met_SNOMAS ',
'Met_SPHU ',
'Met_SPHU1 ',
'Met_SPHU2 ',
'Met_SUNCOS ',
'Met_SUNCOSmid ',
'Met_SWGDN ',
'Met_T ',
'Met_TAUCLI ',
'Met_TAUCLW ',
'Met_THETA ',
'Met_TMPU1 ',
'Met_TMPU2 ',
'Met_TO3 ',
'Met_TropHt ',
'Met_TropLev ',
'Met_TropP ',
'Met_TS ',
'Met_TSKIN ',
'Met_TV ',
'Met_U ',
'Met_U10M ',
'Met_USTAR ',
'Met_UVALBEDO ',
'Met_V ',
'Met_V10M ',
'Met_Z0 ',
::
List of diagnostic fields in the StateMet collection
Diagnostic field |
Description |
Units |
---|---|---|
Met_AD |
Dry air mass |
kg |
Met_AIRDEN |
Dry air density |
kg/m3 |
Met_AIRVOL |
Grid box volume, dry air |
m3 |
Met_ALBD |
Surface albedo |
1 |
Met_AREAM2 |
Grid box area |
m2 |
Met_AVGW |
Water vapor volume mixing ratio |
vol H2O/vol dry air |
Met_BXHEIGHT |
Grid box height |
m |
Met_ChemGridLev |
Chemistry grid level |
1 |
Met_CLDF |
3-D cloud fraction |
|
Met_CLDFRC |
Column cloud fraction |
1 |
Met_CLDTOPS |
Maximum cloud top height |
1 |
Met_DELP |
Delta-pressure between top and bottom edges of grid box (wet air) |
hPa |
Met_DQRCU |
Convective precipitation production rate (dry air) |
kg/kg/s |
Met_DTRAIN |
Detrainment flux |
kg/m2/s |
Met_EFLUX |
Latent heat flux |
W/m2 |
Met_FRCLND |
Olson land fraction |
1 |
Met_FRLAKE |
Fraction of grid box covered by lakes |
1 |
Met_FRLAND |
Fraction of grid box covered by land |
1 |
Met_FRLANDIC |
Fraction of grid box covered by land ice |
1 |
Met_FROCEAN |
Fraction of grid box covered by ocean |
1 |
Met_FRSEAICE |
Fraction of grid box covered by sea ice |
1 |
Met_FRSNO |
Fraction of grid box covered by snow |
1 |
Met_GWETROOT |
Root soil moisture |
1 |
Met_GWETTOP |
Topsoil moisture |
1 |
Met_HFLUX |
Sensible heat flux |
W/m2 |
Met_LAI |
Leaf area index from met field archive |
m2/m2 |
Met_LWI |
Land-water-ice indices |
1 |
Met_PARDF |
Diffuse photosynthetically active radiation |
W/m2 |
Met_PARDR |
Diffuse photosynthetically active radiation |
W/m2 |
Met_PBLTOPL |
PBL top layer |
1 |
Met_PBLH |
PBL height |
m |
Met_PHIS |
Surface geopotential height |
m |
Met_PMID |
Pressure at midpoint of model layers, defined as arithmetic average of edge pressures (wet air) |
hPa |
Met_PMIDDRY |
Pressure at midpoint of model layers, defined as arithmetic average of edge pressures (dry air) |
hPa |
Met_PRECANV |
Anvil precipitation (at surface) |
mm/day |
Met_PRECCON |
Convective precipitation (at surface) |
mm/day |
Met_PRECLSC |
Large-scale precipitation (at surface) |
mm/day |
Met_PRECTOT |
Total precipitation (at surface) |
mm/day |
Met_PS1DRY |
Instantaneous surface pressure at start of 3-hr met field interval (dry air) |
hPa |
Met_PS2DRY |
Instantaneous surface pressure at end of 3-hr met field interval (dry air) |
hPa |
Met_PSC2DRY |
Surface pressure interpolated to current time (dry air) |
hPa |
Met_PS1WET |
Instantaneous surface pressure at start of 3-hr met field interval (wet air) |
hPa |
Met_PS2WET |
Instantaneous surface pressure at end of 3-hr met field interval (wet air) |
hPa |
Met_PSC2WET |
Surface pressure interpolated to current time (wet air) |
hPa |
Met_QI |
Ice mixing ratio (dry air) |
kg/kg dry air |
Met_QL |
Liquid water mixing ratio (dry air) |
kg/kg dry air |
Met_OMEGA |
Updraft velocity |
Pa/s |
Met_OPTD |
Visible optical depth |
1 |
Met_REEVAPCN |
Evaporation of convective precipitation (dry air) |
kg/kg/s |
Met_REEVAPLS |
Evaporation of large-scale + anvil precipitation (dry air) |
kg/kg/s |
Met_SLP |
Sea level pressure |
hPa |
Met_SNODP |
Snow depth |
m |
Met_SNOMAS |
Snow mass |
kg/m2 |
Met_SPHU1 |
Instantaneous specific humidity at start of 3 hr met field interval (wet air) |
kg/kg |
Met_SPHU2 |
Instantaneous specific humidity at end of 3-hr met field interval (wet air) |
kg/kg |
Met_SPHU |
Specific humidity interpolated to current time (wet air) |
g H2O/kg air |
Met_SUNCOS |
Cosine of solar zenith angle at current time |
1 |
Met_SUNCOSMID |
Cosine of solar zenith angle at midpoint of chemistry timestep |
1 |
Met_SWGDN |
Incident shortwave radiation at ground |
W/m2 |
Met_TAUCLI |
Visible optical depth of ice clouds |
1 |
Met_TAUCLW |
Visible optical depth of water clouds |
1 |
Met_THETA |
Potential temperature |
K |
Met_TMPU1 |
Instantaneous temperature at start of 3-hr met field interval |
K |
Met_TMPU2 |
Instantaneous temperature at end of 3-hr met field interval |
K |
Met_T |
Temperature interpolated to current time |
K |
Met_TO3 |
Total overhead ozone column |
Dobsons |
Met_TropHt |
Tropopause height |
u |
Met_TropLev |
Tropopause level |
1 |
Met_TROPP |
Tropopause pressure |
hPa |
Met_TS |
Surface temperature |
K |
Met_TSKIN |
Surface skin temperature |
K |
Met_U |
East-west ccomponent of wind |
m/s |
Met_U10M |
East-west component of wind at 10 m height above surface |
m/s |
Met_USTAR |
Friction velocity |
m/s |
Met_UVALBEDO |
Ultraviolet surface albedo |
1 |
Met_V |
North-south ccomponent of wind |
m/s |
Met_V10M |
North-south component of wind at 10 m height above surface |
m/s |
Met_Z0 |
Surface roughness height |
m |
FracOfTimeInTrop |
Fraction of time spent in the troposphere |
1 |
StratBM
The StratBM collection contains diagnostic fields for GEOS-Chem 10-year stratospheric benchmark simulations. Unless you are involved with benchmarking GEOS-Chem, you may leave this collection deactivated.
Sample definition section for the StratBM collection
StratBM.template: '%y4%m2%d2_%h2%n2z.nc4',
StratBM.frequency: 00000000 010000
StratBM.duration: 00000001 000000
StratBM.mode: 'time-averaged'
StratBM.fields: 'SpeciesConcVV_NO2 ',
'SpeciesConcVV_O3 ',
'SpeciesConcVV_ClO ',
'Met_PSC2WET ',
'Met_BXHEIGHT ',
'Met_AIRDEN ',
'Met_AD ',
::
List of diagnostic fields in the StateMet collection
Diagnostic field |
Description |
Units |
---|---|---|
Met_AD |
Dry air mass |
kg |
Met_AIRDEN |
Dry air density |
kg/m3 |
Met_BXHEIGHT |
Grid box height |
m |
Met_PSC2WET |
Surface pressure interpolated to current time (wet air) |
hPa |
SpeciesConcVV_ClO |
ClO concentration |
mol/mol dry air |
SpeciesConcVV_NO2 |
NO2 concentration |
mol/mol dry air |
SpeciesConcVV_O3 |
O3 concentration |
mol/mol dry air |
Tomas
The TOMAS collection contains diagnostic fields for fullchem simulations with TOMAS aerosol microphysics.
Tomas.template: '%y4%m2%d2_%h2%n2z.nc4',
Tomas.format: 'CFIO',
Tomas.timestampStart: .true.
Tomas.monthly: 0
Tomas.frequency: 010000
Tomas.duration: 010000
Tomas.mode: 'time-averaged'
Tomas.fields: 'TomasH2SO4 ',
'TomasH2SO4mass_?TOMASBIN? ',
#----------------------------------------------
# NOTE: for GEOS-Chem Classic you can use the
# ?TOMASBIN? wildcard. For GCHP you will need
# to list each diagnostic field individually
# such as is shown below:
#'TomasH2SO4mass_bin01 ',
#'TomasH2SO4mass_bin02 ',
#'TomasH2SO4mass_bin03 ',
#'TomasH2SO4mass_bin04 ',
#'TomasH2SO4mass_bin05 ',
#'TomasH2SO4mass_bin06 ',
#'TomasH2SO4mass_bin07 ',
#'TomasH2SO4mass_bin08 ',
#'TomasH2SO4mass_bin09 ',
#'TomasH2SO4mass_bin10 ',
#'TomasH2SO4mass_bin11 ',
#'TomasH2SO4mass_bin12 ',
#'TomasH2SO4mass_bin13 ',
#'TomasH2SO4mass_bin14 ',
#'TomasH2SO4mass_bin15 ',
#----------------------------------------------
'TomasH2SO4number_?TOMASBIN? ',
'TomasCOAG ',
'TomasCOAGmass_?TOMASBIN ',
'TomasCOAGnumber_?TOMASBIN? ',
'TomasNUCL ',
'TomasNUCRATEFN ',
'TomasNUCLmass_?TOMASBIN? ',
'TomasNUCLnumber_?TOMASBIN? ',
'TomasNUCRATEnumber_?TOMASBIN? ',
'TomasAQOX ',
'TomasAQOXmass_?TOMASBIN? ',
'TomasAQOXnumber_?TOMASBIN? ',
'TomasMNFIX ',
'TomasMNFIXmass_?TOMASBIN? ',
'TomasMNFIXnumber_?TOMASBIN? ',
'TomasMNFIXh2so4mass_?TOMASBIN? ',
'TomasMNFIXh2so4number_?TOMASBIN? ',
'TomasMNFIXcoagmass_?TOMASBIN? ',
'TomasMNFIXcoagnumber_?TOMASBIN? ',
'TomasMNFIXaqoxmass_?TOMASBIN? ',
'TomasMNFIXaqoxnumber_?TOMASBIN? ',
'TomasMNFIXezwat1number_?TOMASBIN?',
'TomasMNFIXezwat2mass_?TOMASBIN? ',
'TomasMNFIXezwat2number_?TOMASBIN?',
'TomasMNFIXezwat3mass_?TOMASBIN? ',
'TomasMNFIXezwat3number_?TOMASBIN?',
'TomasMNFIXcheck1mass_?TOMASBIN? ',
'TomasMNFIXcheck1number_?TOMASBIN?',
'TomasMNFIXcheck2mass_?TOMASBIN? ',
'TomasMNFIXcheck2number_?TOMASBIN?',
'TomasMNFIXcheck3mass_?TOMASBIN? ',
'TomasMNFIXcheck3number_?TOMASBIN?',
'TomasSOA ',
'TomasSOAmass_?TOMASBIN? ',
'TomasSOAnumber_?TOMASBIN? ',
::
List of diagnostic fields for the Tomas collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
TomasAQOX |
Tomas aqueous oxidation rate |
1 |
|
TomasAQOXmass_<name|wc> |
TomasAQOX mass rate |
kg/kg/s |
|
TomasAQOXnumber_<name|wc> |
TomasAQOX number rate |
#/kg/s |
|
TomasCOAG |
Tomas coagulation rate |
1 |
|
TomasCOAGmass_<name|wc> |
TOMASCOAG mass rate |
kg/kg/s |
|
TomasCOAGnumber_<name|wc> |
TomasCOAG number rate |
#/kg/s |
|
TomasH2SO4 |
Tomas condensation rate |
1 |
|
TomasH2SO4mass_<name|wc> |
TomasH2SO4 mass rate |
kg/kg/s |
|
TomasH2SO4number_<name|wc> |
TomasH2SO4 number rate |
#/kg/s |
|
TomasMNFIX |
Tomas error rate |
1 |
|
TomasMNFIXmass_<name|wc> |
TomasMNFIX mass rate |
kg/kg/s |
|
TomasMNFIXnumber_<name|wc> |
TomasMNFIX number rate |
#/kg/s |
|
TomasMNFIXaqoxmass_<name|wc> |
TOMASMNFIXAQOX mass rate |
kg/kg/s |
|
TomasMNFIXaqoxnumber_ <name|wc> |
TOMASMNFIXAQOX number rate |
#/kg/s |
|
TomasMNFIXcoagmass_<name|wc> |
TomasMNFIXCOAG mass rate |
kg/kg/s |
|
TomasMNFIXcoagnumber_<name|wc> |
TomasMNFIXCOAG number rate |
#/kg/s |
|
TomasMNFIXcheck1mass_<name|wc> |
TOMASMNFIXCHECK1 mass rate |
kg/kg/s |
|
TomasMNFIXcheck1number_<name|wc> |
TOMASMNFIXCHECK1 number rate |
#/kg/s |
|
TomasMNFIXcheck2mass_<name|wc> |
TOMASMNFIXCHECK2 mass rate |
kg/kg/s |
|
TomasMNFIXcheck2number_<name|wc> |
TOMASMNFIXCHECK2 number rate |
#/kg/s |
|
TomasMNFIXcheck3mass_<name|wc> |
TOMASMNFIXCHECK3 mass rate |
kg/kg/s |
|
TomasMNFIXcheck3number_<name|wc> |
TOMASMNFIXCHECK3 number rate |
#/kg/s |
|
TomasMNFIXezwat1mass_<name|wc> |
TOMASMNFIXEZWAT1 mass rate |
kg/kg/s |
|
TomasMNFIXezwat1number_<name|wc> |
TOMASMNFIXEZWAT1 number rate |
#/kg/s |
|
TomasMNFIXezwat2mass_<name|wc> |
TOMASMNFIXEZWAT2 mass rate |
kg/kg/s |
|
TomasMNFIXezwat2number_<name|wc> |
TOMASMNFIXEZWAT2 number rate |
#/kg/s |
|
TomasMNFIXezwat3mass_<name|wc> |
TOMASMNFIXEZWAT3 mass rate |
kg/kg/s |
|
TomasMNFIXezwat3number_<name|wc> |
TOMASMNFIXEZWAT3 number rate |
#/kg/s |
|
TomasMNFIXh2so4mass_<name|wc> |
TomasMNFIXH2SO4 mass rate |
kg/kg/s |
|
TomasMNFIXh2so4number_<name|wc> |
TomasMNFIXH2SO4 number rate |
#/kg/s |
|
TomasNUCL |
Tomas nucleation rate |
1 |
|
TomasNUCLmass_<name|wc> |
TomasNUCL mass rate |
kg/kg/s |
|
TomasNUCLnumber_<name|wc> |
TomasNUCL number rate |
#/kg/s |
|
TomasSOA |
TomasSOA rate |
1 |
|
TomasSOAmass_<name|wc> |
TomasSOA mass rate |
kg/kg/s |
|
TomasSOAnumber_<name|wc> |
TomasSOA number rate |
#/kg/s |
Notes for the Tomas collection
This diagnostic field can use the ?TOMASBIN? wildcard (for GEOS-Chem Classic only).
UVFlux
The UVFlux diagnostic contains diffuse, direct, and net UV fluxes at each of the photolysis wavelength bins.
Sample definition section for HISTORY.rc
UVFlux.template: '%y4%m2%d2_%h2%n2z.nc4',
UVFlux.frequency: 00000100 000000
UVFlux.duration: 00000100 000000
UVFlux.mode: 'time-averaged'
UVFlux.template: 'UVFluxDiffuse_?UVFLX?',
'UVFluxDirect_?UVFLX? ',
'UVFluxNet_?UVFLX? ',
::
List of diagnostic fields in the UvFlux collection
Diagnostic field |
Description |
Units |
Wildcards |
---|---|---|---|
UVFluxDiffuse_<name|wc> |
Diffuse UV flux in wavelength bin |
W/m2 |
?UVFLX? |
UVFluxDirect_<name|wc> |
Direct UV flux in wavelength bin |
W/m2 |
?UVFLX? |
UVFluxNet_<name|wc> |
Net UV flux in wavelength bin |
W/m2 |
?UVFLX? |
WetLossConv
The WetLossConv collection contains diagnostics fluxes of soluble species lost to wet scavenging in convective updrafts.
Sample definition section for HISTORY.rc
WetLossConv.template: '%y4%m2%d2_%h2%n2z.nc4',
WetLossConv.frequency: 00000100 000000
WetLossConv.duration: 00000100 000000
WetLossConv.mode: 'time-averaged'
WetLossConv.fields: 'WetLossConv_?WET? ',
'WetLossConvFrac_?WET?',
::
List of diagnostic fields in the WetLossConv collection
Diagnostic field |
Description |
Units |
Wildcard |
---|---|---|---|
WetLossConv_<name|wc> |
Loss of soluble species scavenged by cloud updrafts in moist convection |
kg/s |
?WET? |
WetLossConvFrac_<name|wc> |
Fraction of species scavenged by cloud updrafts in moist convection |
1 |
?WET? |
WetLossLS
The WetLossLS collection contains diagnostics fluxes of soluble species lost to rainout and washout in large-scale wet deposition.
Sample definition section for HISTORY.rc
WetLossLS.template: '%y4%m2%d2_%h2%n2z.nc4',
WetLossLS.frequency: 00000100 000000
WetLossLS.duration: 00000100 000000
WetLossLS.mode: 'time-averaged'
WetLossLS.fields: 'WetLossLS_?WET?',
::
List of diagnostic fields in the WetLossLS collection
Diagnostic field |
Description |
Units | Wildcard |
||
---|---|---|---|---|
WetLossLS_<name|wc> |
Loss of soluble species in large-scale precipitation |
kg/s |
?WET? |
Adding new History diagnostics
To add your own diagnostics we recommend that you find a similar
existing diagnostic and use its implementation as a template. Most of
the work is done in Headers/state_diag_mod.F90
. Briefly, the
following updates to that file are essential for adding in your own
netCDF diagnostics:
Declare diagnostic array at top of module.
Set diagnostic array pointer to null in
Zero_State_Diag
subroutine.
Create a section in
Init_State_Diag
subroutine to allocate and register the array.
Deallocate the diagnostic array in subroutine
Cleanup_State_Diag
.
Add an
IF
block for the diagnostic within subroutineGet_Metadata_State_Diag
to define its metadata, making sure to list the diagnostic name with all capital letters.
Good diagnostics to use as templates are SpeciesConcVV
for
3-dimensional arrays that are for all species and DryDepVel
for
2-dimensional arrays that are for a subset of species. If your
diagnostic is not per species then AODDust
is a good
diagnostic to look at. Search the file
Headers/state_diag_mod.F90
for any of these strings to find
all instances of code related to their implementation.
Note that information about the dimensions and species collection the
diagnostic will include must be specified when allocating the array in
Init_State_Diag
and in Get_Metadata_State_Diag
. In the
latter subroutine the Rank
is the integer number of dimensions
of the diagnostic (not including species) and the TagID
string, if any, specifies the species collection to output the
diagnostic per. Each TagID
is defined in subroutine
Get_TagInfo. Each TagID string can also be used as a wildcard within
HISTORY.rc
to simplify diagnostic name specification (for
GEOS-Chem Classic only).
Once you have implemented your diagnostic in
Headers/state_diag_mod.F90
, try adding it to HISTORY.rc and
running. You should get your diagnostic output in the netCDF output
file as all zeros. The next step is to populate the array with
whatever value you want to output. You should do this by passing the
State_Diag
array to the location where you want to set the
values. Then write code to fill the array. A simple test of your
understanding is to initially set values to a constant other than zero
and see if the output matches what you set the arrays to.
For additional help implementing your own GEOS-Chem diagnostics please contact the GEOS-Chem Support Team.