Skip to content

Commit

Permalink
Update HEMCO_Config.rc for carbon simulations read chemistry inputs b…
Browse files Browse the repository at this point in the history 
…ased on species defines

Logical switches have been added to HEMCO_Config.rc.carbon to only read
the data necessary based on the carbon species used. For example, GLOBAL_OH
is not needed for CO2-only simulations and will therefore not be read anymore.
This should speed up single-species carbon simulations.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
  • Loading branch information
@msulprizio
msulprizio committed Oct 16, 2024
1 parent 7043157 commit e29adfd
Show file tree
Hide file tree
Showing 3 changed files with 35 additions and 22 deletions.
9 changes: 9 additions & 0 deletions CHANGELOG.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4,6 +4,15 @@ This file documents all notable changes to the GEOS-Chem repository starting in

The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/), and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).

## [Unreleased] - TBD
### Changed
- Renamed `Emiss_Carbon_Gases` to `CO2_Production` in `carbon_gases_mod.F90`

### Fixed
- Added a fix to skip the call to KPP when only CO2 is defined in the carbon simulation
- Added fix to turn on ship emissions for CO2 in the carbon simulation
- Update `HEMCO_Config.rc` for carbon simulation to read data based on carbon species used

## [Unreleased] - TBD
### Added
- Added vectors `State_Chm%KPP_AbsTol` and `State_Chm%KPP_RelTol`
Expand Down
24 changes: 13 additions & 11 deletions run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1262,6 +1262,9 @@ Mask fractions: false
#------------------------------------------------------------------------------
# --- Global OH fields ---
#------------------------------------------------------------------------------
(((USE_CH4_DATA.or.USE_CO_DATA
# --- OH from GEOS-Chem v5-07 [kg/m3], needed for CH4/IMI ---
(((GLOBAL_OH_GCv5
Expand All @@ -1282,10 +1285,14 @@ ${RUNDIR_GLOBAL_OH}
))).not.GLOBAL_OH_GCv5
)))GLOBAL_OH_GC14
)))USE_CH4_DATA.or.USE_CO_DATA
#------------------------------------------------------------------------------
# --- Quantities needed for CH4 chemistry ---
#------------------------------------------------------------------------------
(((USE_CH4_DATA
# --- Global CH4 loss frequencies [1/s] -----------
(((CH4_LOSS_FREQ
${RUNDIR_CH4_LOSS}
Expand All @@ -1296,9 +1303,12 @@ ${RUNDIR_CH4_LOSS}
${RUNDIR_GLOBAL_Cl}
)))GLOBAL_CL
)))USE_CH4_DATA
#------------------------------------------------------------------------------
# --- Quantities needed for CO chemistry ---
#------------------------------------------------------------------------------
(((USE_CO_DATA
# -- P(CO) from CH4 and NMVOC from the last 10-yr benchmark [molec/cm3/s] ---
(((PROD_CO_CH4
Expand All @@ -1310,33 +1320,25 @@ ${RUNDIR_PCO_NMVOC}
# --- GMI chemistry: prod/loss rates (for strato-/mesosphere) ---
# --- Units: prod [v/v/s]; loss [1/s] ---
(((USE_CO_DATA
(((GMI_PROD_LOSS
* GMI_LOSS_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc loss 2005/1-12/1/0 C xyz s-1 CO - 1 1
* GMI_PROD_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc prod 2005/1-12/1/0 C xyz v/v/s CO - 1 1
)))GMI_PROD_LOSS
)))USE_CO_DATA
# If CO is not an advected species, then just read as a chemistry input
(((.not.USE_CO_DATA
(((GMI_PROD_LOSS
* GMI_LOSS_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc loss 2005/1-12/1/0 C xyz s-1 * - 1 1
* GMI_PROD_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc prod 2005/1-12/1/0 C xyz v/v/s * - 1 1
)))GMI_PROD_LOSS
))).not.USE_CO_DATA
)))USE_CO_DATA
#------------------------------------------------------------------------------
# --- Quantities needed for CO2 chemistry ---
#------------------------------------------------------------------------------
(((USE_CO2_DATA
# --- CHEMICAL PRODUCTION FROM CO OXIDATION ------
# --- Recommended for use in forward modelling ---
# --- Optional for inversion/assimilation --------
(((USE_CO2_DATA
(((CO2_COPROD
${RUNDIR_CO2_COPROD}
)))CO2_COPROD
)))USE_CO2_DATA
)))CHEMISTRY_INPUT
Expand Down
24 changes: 13 additions & 11 deletions run/GCHP/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1262,6 +1262,9 @@ Mask fractions: false
#------------------------------------------------------------------------------
# --- Global OH fields ---
#------------------------------------------------------------------------------
(((USE_CH4_DATA.or.USE_CO_DATA
# --- OH from GEOS-Chem v5-07 [kg/m3], needed for CH4/IMI ---
(((GLOBAL_OH_GCv5
Expand All @@ -1282,10 +1285,14 @@ ${RUNDIR_GLOBAL_OH}
))).not.GLOBAL_OH_GCv5
)))GLOBAL_OH_GC14
)))USE_CH4_DATA.or.USE_CO_DATA
#------------------------------------------------------------------------------
# --- Quantities needed for CH4 chemistry ---
#------------------------------------------------------------------------------
(((USE_CH4_DATA
# --- Global CH4 loss frequencies [1/s] -----------
(((CH4_LOSS_FREQ
${RUNDIR_CH4_LOSS}
Expand All @@ -1296,9 +1303,12 @@ ${RUNDIR_CH4_LOSS}
${RUNDIR_GLOBAL_Cl}
)))GLOBAL_CL
)))USE_CH4_DATA
#------------------------------------------------------------------------------
# --- Quantities needed for CO chemistry ---
#------------------------------------------------------------------------------
(((USE_CO_DATA
# -- P(CO) from CH4 and NMVOC from the last 10-yr benchmark [molec/cm3/s] ---
(((PROD_CO_CH4
Expand All @@ -1310,33 +1320,25 @@ ${RUNDIR_PCO_NMVOC}
# --- GMI chemistry: prod/loss rates (for strato-/mesosphere) ---
# --- Units: prod [v/v/s]; loss [1/s] ---
(((USE_CO_DATA
(((GMI_PROD_LOSS
* GMI_LOSS_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc loss 2005/1-12/1/0 C xyz s-1 CO - 1 1
* GMI_PROD_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc prod 2005/1-12/1/0 C xyz v/v/s CO - 1 1
)))GMI_PROD_LOSS
)))USE_CO_DATA
# If CO is not an advected species, then just read as a chemistry input
(((.not.USE_CO_DATA
(((GMI_PROD_LOSS
* GMI_LOSS_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc loss 2005/1-12/1/0 C xyz s-1 * - 1 1
* GMI_PROD_CO $ROOT/GMI/v2022-11/gmi.clim.CO.geos5.2x25.nc prod 2005/1-12/1/0 C xyz v/v/s * - 1 1
)))GMI_PROD_LOSS
))).not.USE_CO_DATA
)))USE_CO_DATA
#------------------------------------------------------------------------------
# --- Quantities needed for CO2 chemistry ---
#------------------------------------------------------------------------------
(((USE_CO2_DATA
# --- CHEMICAL PRODUCTION FROM CO OXIDATION ------
# --- Recommended for use in forward modelling ---
# --- Optional for inversion/assimilation --------
(((USE_CO2_DATA
(((CO2_COPROD
${RUNDIR_CO2_COPROD}
)))CO2_COPROD
)))USE_CO2_DATA
)))CHEMISTRY_INPUT
Expand Down

0 comments on commit e29adfd

Please sign in to comment.