diff --git a/CHANGELOG.md b/CHANGELOG.md index 7ae928fca..f61067d70 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -30,6 +30,7 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/), - Allocated `State_Diag%SatDiagnPEDGE` ffield with vertical dimension `State_Grid%NZ+1` - Modified `run/GCClassic/cleanRunDir.sh` to skip removing bpch files, as well as now removing `fort.*` and `OutputDir/*.txt` files - Edited `run/shared/kpp_standalone_interface.yml` to include additional entries under `active cells` and `locations` +- Modified tagCO simulation to use GFED4 biomass burning emissions and GEOS-Chem v5 OH fields for consistency with carbon simulation ### Fixed - Added a fix to skip the call to KPP when only CO2 is defined in the carbon simulation @@ -43,6 +44,12 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/), - Fixed errors in adjoint-only code preventing successful adjoint build - Fixed zero convective precipitation and high cloud base in runs using GEOS-FP (>=01Jun2020) or GEOS-IT - Updated GEOS-only code for compatibility with GEOS-Chem 14.5 +- Fixed typos in `HEMCO_Config.rc` for CH4 simulations causing mobile combustion emissions to be double counted +- Fixed handling of FIRST flag in carbon_gases_mod.F to limit log prints to first timestep only +- Removed extraneous pressure correction in GCHP carbon simulations by adding 'activate: true' to geoschem_config.yml +- Fixed bug in GC-Classic OCS emissions where unit conversion of km2 to m2 occurred twice +- Changed dimension of EmisOCS_Total from 3D to 2D since all emissions for all sectors are 2D +- Added fixes to only apply archived PCO_CH4 field for carbon simulations with CO only ### Removed - Removed duplicate `WD_RetFactor` tag for HgClHO2 in `species_database.yml` diff --git a/GeosCore/carbon_gases_mod.F90 b/GeosCore/carbon_gases_mod.F90 index 4a91aeb40..7a30db321 100644 --- a/GeosCore/carbon_gases_mod.F90 +++ b/GeosCore/carbon_gases_mod.F90 @@ -238,6 +238,9 @@ SUBROUTINE CO2_Production( Input_Opt, State_Chm, State_Diag, & Spc => NULL() ENDIF + ! Reset first-time flag + FIRST = .FALSE. + ! Free pointers for safety's sake Spc => NULL() Ptr2D => NULL() diff --git a/GeosCore/input_mod.F90 b/GeosCore/input_mod.F90 index 0812c9924..1e2c1c180 100644 --- a/GeosCore/input_mod.F90 +++ b/GeosCore/input_mod.F90 @@ -2020,7 +2020,7 @@ SUBROUTINE Config_CO( Config, Input_Opt, RC ) !------------------------------------------------------------------------ ! Use P(CO) from CH4 (archived from a fullchem simulation)? !------------------------------------------------------------------------ - key = "CO_simulation_options%use_fullchem_PCO_from_CH4" + key = "CO_simulation_options%use_archived_PCO_from_CH4" v_bool = MISSING_BOOL CALL QFYAML_Add_Get( Config, key, v_bool, "", RC ) IF ( RC /= GC_SUCCESS ) THEN @@ -2033,7 +2033,7 @@ SUBROUTINE Config_CO( Config, Input_Opt, RC ) !------------------------------------------------------------------------ ! Use P(CO) from NMVOC (archived from a fullchem simulation)? !------------------------------------------------------------------------ - key = "CO_simulation_options%use_fullchem_PCO_from_NMVOC" + key = "CO_simulation_options%use_archived_PCO_from_NMVOC" v_bool = MISSING_BOOL CALL QFYAML_Add_Get( Config, key, v_bool, "", RC ) IF ( RC /= GC_SUCCESS ) THEN @@ -2050,8 +2050,8 @@ SUBROUTINE Config_CO( Config, Input_Opt, RC ) WRITE(6,90 ) 'TAGGED CO SIMULATION SETTINGS' WRITE(6,95 ) '(overwrites any other settings related to CO)' WRITE(6,95 ) '---------------------------------------------' - WRITE(6,100) 'Use full chem. P(CO) from CH4? :', Input_Opt%LPCO_CH4 - WRITE(6,100) 'Use full chem. P(CO) from NMVOC? :', Input_Opt%LPCO_NMVOC + WRITE(6,100) 'Use archived P(CO) from CH4? :', Input_Opt%LPCO_CH4 + WRITE(6,100) 'Use archived P(CO) from NMVOC? :', Input_Opt%LPCO_NMVOC ENDIF ! FORMAT statements @@ -2113,9 +2113,9 @@ SUBROUTINE Config_CO2( Config, Input_Opt, RC ) thisLoc = ' -> at Config_CO2 (in module GeosCore/input_mod.F90)' !------------------------------------------------------------------------ - ! Turn on CO2 3D chemical source and surface correction? + ! Use archived fields of CO2 production from CO oxidation? !------------------------------------------------------------------------ - key = "CO2_simulation_options%sources%3D_chemical_oxidation_source" + key = "CO2_simulation_options%sources%use_archived_PCO2_from_CO" v_bool = MISSING_BOOL CALL QFYAML_Add_Get( Config, key, v_bool, "", RC ) IF ( RC /= GC_SUCCESS ) THEN @@ -2158,7 +2158,7 @@ SUBROUTINE Config_CO2( Config, Input_Opt, RC ) WRITE( 6,90 ) 'CO2 SIMULATION SETTINGS' WRITE( 6,95 ) '(overwrites any other settings related to CO2)' WRITE( 6,95 ) '----------------------------------------------' - WRITE( 6,100 ) 'CO2 from oxidation (CO,CH4,..):', Input_Opt%LCHEMCO2 + WRITE( 6,100 ) 'Use archived P(CO2) from CO? :', Input_Opt%LCHEMCO2 WRITE( 6, 95 ) 'Tagged CO2 settings' WRITE( 6,100 ) ' Tag Biosphere/Ocean CO2 :', Input_Opt%LBIOSPHTAG WRITE( 6,100 ) ' Tag Fossil Fuel CO2 :', Input_Opt%LFOSSILTAG diff --git a/KPP/carbon/carbon_Funcs.F90 b/KPP/carbon/carbon_Funcs.F90 index cfe12a7d2..734604aa2 100644 --- a/KPP/carbon/carbon_Funcs.F90 +++ b/KPP/carbon/carbon_Funcs.F90 @@ -193,7 +193,7 @@ SUBROUTINE carbon_ComputeRateConstants( & C(ind_FixedCl) = ConcClMnd ! CH4 + offline OH reaction rate [1/s] - ! This is a pseudo-2nd order rate appropriate for CH4 + OH + ! Use rates saved from full-chemistry run (if CH4 is not advected) IF ( PCO_fr_CH4_use ) THEN k_Trop(1) = PCO_fr_CH4 * OHdiurnalFac k_Trop(1) = SafeDiv( k_Trop(1), C(ind_CH4)*C(ind_FixedOH), 0.0_dp ) diff --git a/KPP/carbon/gckpp_Function.F90 b/KPP/carbon/gckpp_Function.F90 index 31145e8f7..05a1fe789 100644 --- a/KPP/carbon/gckpp_Function.F90 +++ b/KPP/carbon/gckpp_Function.F90 @@ -2,7 +2,7 @@ ! ! The ODE Function of Chemical Model File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Global.F90 b/KPP/carbon/gckpp_Global.F90 index ab0d870d2..35126180f 100644 --- a/KPP/carbon/gckpp_Global.F90 +++ b/KPP/carbon/gckpp_Global.F90 @@ -2,7 +2,7 @@ ! ! Global Data Module File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Initialize.F90 b/KPP/carbon/gckpp_Initialize.F90 index b04d661ad..fc2a3904c 100644 --- a/KPP/carbon/gckpp_Initialize.F90 +++ b/KPP/carbon/gckpp_Initialize.F90 @@ -2,7 +2,7 @@ ! ! Initialization File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Integrator.F90 b/KPP/carbon/gckpp_Integrator.F90 index 3a267e9fa..15e1ee7dd 100644 --- a/KPP/carbon/gckpp_Integrator.F90 +++ b/KPP/carbon/gckpp_Integrator.F90 @@ -2,7 +2,7 @@ ! ! Numerical Integrator (Time-Stepping) File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Jacobian.F90 b/KPP/carbon/gckpp_Jacobian.F90 index a79264e6f..65cfcc29d 100644 --- a/KPP/carbon/gckpp_Jacobian.F90 +++ b/KPP/carbon/gckpp_Jacobian.F90 @@ -2,7 +2,7 @@ ! ! The ODE Jacobian of Chemical Model File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_LinearAlgebra.F90 b/KPP/carbon/gckpp_LinearAlgebra.F90 index 655002792..afd23849e 100644 --- a/KPP/carbon/gckpp_LinearAlgebra.F90 +++ b/KPP/carbon/gckpp_LinearAlgebra.F90 @@ -2,7 +2,7 @@ ! ! Linear Algebra Data and Routines File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Monitor.F90 b/KPP/carbon/gckpp_Monitor.F90 index 7063190ad..126258534 100644 --- a/KPP/carbon/gckpp_Monitor.F90 +++ b/KPP/carbon/gckpp_Monitor.F90 @@ -2,7 +2,7 @@ ! ! Utility Data Module File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Parameters.F90 b/KPP/carbon/gckpp_Parameters.F90 index 70d35690b..874e5dfeb 100644 --- a/KPP/carbon/gckpp_Parameters.F90 +++ b/KPP/carbon/gckpp_Parameters.F90 @@ -2,7 +2,7 @@ ! ! Parameter Module File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Rates.F90 b/KPP/carbon/gckpp_Rates.F90 index 8d0b302e7..0195411f7 100644 --- a/KPP/carbon/gckpp_Rates.F90 +++ b/KPP/carbon/gckpp_Rates.F90 @@ -2,7 +2,7 @@ ! ! The Reaction Rates File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/KPP/carbon/gckpp_Util.F90 b/KPP/carbon/gckpp_Util.F90 index 7fd3d92fb..4fd4fee48 100644 --- a/KPP/carbon/gckpp_Util.F90 +++ b/KPP/carbon/gckpp_Util.F90 @@ -2,7 +2,7 @@ ! ! Auxiliary Routines File ! -! Generated by KPP-3.0.0 symbolic chemistry Kinetics PreProcessor +! Generated by KPP-3.1.1 symbolic chemistry Kinetics PreProcessor ! (https:/github.com/KineticPreProcessor/KPP ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) diff --git a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.CH4 b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.CH4 index b914853c3..840fc4780 100644 --- a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.CH4 +++ b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.CH4 @@ -135,6 +135,7 @@ VerboseOnCores: root # Accepted values: root all )))UseTotalPriorEmis (((.not.UseTotalPriorEmis + #============================================================================== # --- Gridded GHGI v2 (Maasakkers et al., submitted to ES&T, 2023) --- # @@ -251,7 +252,7 @@ VerboseOnCores: root # Accepted values: root all 0 GHGI_EE_COAST_GAS_TRANSMISSION - - - - - - CH4 1009 2 5 0 GHGI_EE_GAS_POSTMETER $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_Supp_1B2b_PostMeter 2012-2020/1/1/0 C xy molec/cm2/s CH4 1008 2 100 -### Coal ### +### Coal ### 0 GHGI_EE_COAL_UNDERGROUND $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1B1a_Underground_Coal 2012-2020/1/1/0 C xy molec/cm2/s CH4 1008 3 100 0 GHGI_EE_COAST_COAL_UNDERGROUND - - - - - - CH4 1009 3 5 0 GHGI_EE_COAL_SURFACE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1B1a_Surface_Coal 2012-2020/1/1/0 C xy molec/cm2/s CH4 1008 3 100 @@ -281,11 +282,11 @@ VerboseOnCores: root # Accepted values: root all ### Rice ### 0 GHGI_EE_RICE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_3C_Rice_Cultivation 2012-2020/1-12/1/0 C xy molec/cm2/s CH4 58/1008 7 100 -0 GHGI_EE_COAST_RICE - - - - -y - CH4 58/1009 7 1 +0 GHGI_EE_COAST_RICE - - - - - - CH4 58/1009 7 1 ### Other Anthro ### 0 GHGI_EE_OTHER_MCOMB $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1A_Combustion_Mobile 2012-2020/1/1/0 C xy molec/cm2/s CH4 1008 8 100 -0 GHGI_EE_COAST_OTHER_MCOMB - - - - - - CH4 1008 8 100 +0 GHGI_EE_COAST_OTHER_MCOMB - - - - - - CH4 1009 8 1 0 GHGI_EE_OTHER_SCOMB $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1A_Combustion_Stationary 2012-2020/1-12/1/0 C xy molec/cm2/s CH4 50/1008 8 100 0 GHGI_EE_COAST_OTHER_SCOMB - - - - - - CH4 50/1009 8 1 0 GHGI_EE_OTHER_PIND $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_2B8_Industry_Petrochemical 2012-2020/1/1/0 C xy molec/cm2/s CH4 1008 8 100 @@ -415,7 +416,7 @@ VerboseOnCores: root # Accepted values: root all )))EDGARv8 #============================================================================== -# CEDS (historical) or Shared Socioeconomic Pathways (future) +# --- CEDS (historical) or Shared Socioeconomic Pathways (future) --- #============================================================================== (((CMIP6_SFC_LAND_ANTHRO 0 CMIP6_CH4_AGR $ROOT/CMIP6/v2021-01/$GCAPSCENARIO/$GCAPSCENARIO_$YYYY.nc4 CH4_agr 1970-2014/1-12/1/0 C xy kg/m2/s CH4 - 4 1 @@ -431,7 +432,7 @@ VerboseOnCores: root # Accepted values: root all )))CMIP6_SHIP #============================================================================== -# BB4MIPs historical / SSP future biomass burning inventories +# --- BB4MIPs historical / SSP future biomass burning inventories --- #============================================================================== (((BB4MIPS 0 CMIP6_BB_CH4 $ROOT/CMIP6/v2021-01/$GCAPSCENARIO/$GCAPSCENARIO_$YYYY.nc4 CH4_bbn 1750-2015/1-12/1/0 C xyL=1:PBL kg/m2/s CH4 75 9 1 @@ -450,7 +451,7 @@ VerboseOnCores: root # Accepted values: root all #============================================================================== # --- CH4: JPL WetCHARTs v1.0 (Bloom et al., https://doi.org/10.3334/ORNLDAAC/1502) --- # -# # Use updated files (v2024-01); these are COARDS-compliant. +# Use updated files (v2024-01); these are COARDS-compliant. #============================================================================== (((JPL_WETCHARTS 0 JPLW_CH4 $ROOT/CH4/v2024-01/JPL_WetCharts/HEensemble/JPL_WetCharts_2010-2019.Ensemble_Mean.0.5x0.5.nc emi_ch4 2010-2019/1-12/1/0 C xy molec/cm2/s CH4 - 10 1 @@ -496,7 +497,7 @@ VerboseOnCores: root # Accepted values: root all )))FUNG_SOIL_ABSORPTION #============================================================================== -# --- # Soil absorption from MeMo model (Murguia-Flores et al. 2018, GMD) --- +# --- Soil absorption from MeMo model (Murguia-Flores et al. 2018, GMD) --- # # - Multiply soil absorption by -1 to get a "negative" flux #============================================================================== @@ -623,7 +624,7 @@ VerboseOnCores: root # Accepted values: root all # --- GEOS-Chem boundary condition file --- #============================================================================== (((GC_BCs -* BC_ $ROOT/SAMPLE_BCs/v2021-07/CH4/GEOSChem.BoundaryConditions.$YYYY$MM$DD_0000z.nc4 SpeciesBC_?ADV? 1900-2100/1-12/1-31/* EFY xyz 1 * - 1 1 +* BC_ $ROOT/SAMPLE_BCs/v2021-07/CH4/GEOSChem.BoundaryConditions.$YYYY$MM$DD_0000z.nc4 SpeciesBC_?ADV? 1900-2100/1-12/1-31/* EFY xyz 1 * - 1 1 )))GC_BCs #============================================================================== @@ -945,7 +946,6 @@ ${RUNDIR_GLOBAL_Cl} # Country/region masks #============================================================================== -# Use updated files (v2024-01); these are COARDS-compliant. (((Scarpelli_Mexico 1001 MEX_MASK $ROOT/MASKS/v2024-04/Mexico_Mask.01x01.nc MASK 2000/1/1/0 C xy 1 1 -118/17/-95/33 1010 MEX_MASK_MIRROR $ROOT/MASKS/v2024-04/Mexico_Mask_Mirror.01x01.nc MASK 2000/1/1/0 C xy 1 1 -118/17/-95/33 diff --git a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon index 95c72aca8..445ffe4c3 100644 --- a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon +++ b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon @@ -240,8 +240,8 @@ Mask fractions: false 0 GHGI_COAST_WASTEWATER_IND - - - - - - CH4 1009 6 1 ### Rice ### -0 GHGI_RICE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_3C_Rice_Cultivation 2012-2018/1-12/1/0 EFY xy molec/cm2/s CH4 58/1008 7 100 -0 GHGI_COAST_RICE - - - - - - CH4 58/1009 7 1 +0 GHGI_RICE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_3C_Rice_Cultivation 2012-2018/1-12/1/0 EFY xy molec/cm2/s CH4 58/1008 7 100 +0 GHGI_COAST_RICE - - - - - - CH4 58/1009 7 1 ### Other Anthro ### 0 GHGI_OTHER_MCOMB $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1A_Combustion_Mobile 2012-2018/1/1/0 EFY xy molec/cm2/s CH4 1008 8 100 @@ -324,7 +324,7 @@ Mask fractions: false ### Rice ### 0 GHGI_EE_RICE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_3C_Rice_Cultivation 2012-2020/1-12/1/0 C xy molec/cm2/s CH4 58/1008 7 100 -0 GHGI_EE_COAST_RICE - - - - - - CH4 58/1009 7 1 +0 GHGI_EE_COAST_RICE - - - - - - CH4 58/1009 7 1 ### Other Anthro ### 0 GHGI_EE_OTHER_MCOMB $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1A_Combustion_Mobile 2012-2020/1/1/0 C xy molec/cm2/s CH4 1008 8 100 @@ -399,7 +399,7 @@ Mask fractions: false )))Scarpelli_Canada #============================================================================== -# --- CH4: Global Fuel Exploitation Inventory (GFEI v2, Scarpelli et al., 2021) +# --- CH4: Global Fuel Exploitation Inventory (GFEI v2, Scarpelli et al., 2021) --- # # This inventory will replace EDGAR (oil, gas, & coal) #============================================================================== @@ -481,11 +481,11 @@ Mask fractions: false )))BB4MIPS #============================================================================== -# --- CH4: biomass burning --- +# --- CH4: QFED2 biomass burning --- +# +# Use QFED v2.4r8 for CH4. The data in QFED/v2014-09 are for v2.4r5 and have +# been found to be buggy for CH4 (same emissions on each day). #============================================================================== - -# --- QFED v2.4r8 --- -# NOTE: Use this instead of data in QFED/v2014-09 (these are buggy) (((QFED2 0 QFED_CH4 $ROOT/CH4/v2017-10/QFED/QFEDv2.4r8.emis_ch4.$YYYY.nc emi_ch4 2009-2015/1-12/1-31/0 C xy molec/cm2/s CH4 - 9 2 )))QFED2 @@ -986,19 +986,19 @@ Mask fractions: false (((USE_OCS_DATA (((OCS_ANTHRO_FLUX -0 OCS_ANTHRO $ROOT/OCS_FLUX/v2022-11/anthro_v3/v3_anthro/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 1 1 +0 OCS_ANTHRO $ROOT/OCS_FLUX/v2022-11/anthro_v3/v3_anthro/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 1 1 )))OCS_ANTHRO_FLUX (((OCS_BIOMASS_FLUX -0 OCS_BIOMASS $ROOT/OCS_FLUX/v2022-11/Biomassburn/CMSV2-CO2-GFED3-ZC/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 2 1 +0 OCS_BIOMASS $ROOT/OCS_FLUX/v2022-11/Biomassburn/CMSV2-CO2-GFED3-ZC/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 2 1 )))OCS_BIOMASS_FLUX (((OCS_MISSING_OCEAN_FLUX -0 OCS_MISSING_OCEAN $ROOT/OCS_FLUX/v2022-11/MissingOcean/LUKAI/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 3 1 +0 OCS_MISSING_OCEAN $ROOT/OCS_FLUX/v2022-11/MissingOcean/LUKAI/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 3 1 )))OCS_MISSING_OCEAN_FLUX (((OCS_OCEAN_FLUX -0 OCS_OCEAN $ROOT/OCS_FLUX/v2022-11/OceanCOS/Kettle/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 4 1 +0 OCS_OCEAN $ROOT/OCS_FLUX/v2022-11/OceanCOS/Kettle/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 4 1 )))OCS_OCEAN_FLUX )))USE_OCS_DATA @@ -1398,7 +1398,7 @@ ${RUNDIR_CO2_COPROD} (((EMISSIONS #============================================================================== -# --- Scale factors for CH4 emissions ---- +# --- CH4: Seasonal scaling factors ---- #============================================================================== (((USE_CH4_DATA @@ -1419,7 +1419,7 @@ ${RUNDIR_CO2_COPROD} )))GHGI_v2.or.GHGI_v2_Express_Ext #------------------------------------------------------------------------------ -# --- Scarpelli et al Mexico manure & rice scale factors --- +# --- Scarpelli_Mexico manure & rice scale factors --- #------------------------------------------------------------------------------ (((Scarpelli_Mexico 10 MANURE_SF $ROOT/CH4/v2017-10/Seasonal_SF/EMICH4_Manure_ScalingFactors.WithClimatology.nc sf_ch4 2008-2016/1-12/1/0 C xy 1 1 @@ -1555,11 +1555,6 @@ ${RUNDIR_CO2_COPROD} #------------------------------------------------------------------------------ 70 StoOCS MATH:60.07/32.06 - - - xy unitless 1 -#------------------------------------------------------------------------------ -# --- km-2 to m-2 --- -#------------------------------------------------------------------------------ -71 Km2ToM2 1.0e-6 - - - xy unitless 1 - )))USE_OCS_DATA #============================================================================== diff --git a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.tagCO b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.tagCO index a9e02c190..49c972ac6 100644 --- a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.tagCO +++ b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.tagCO @@ -116,7 +116,7 @@ Mask fractions: false --> CMIP6_AIRCRAFT : false # 1850-2100 --> CMIP6_SHIP : false # 1850-2100 # ----- BIOMASS BURNING EMISSIONS ----------------------------- - --> QFED2 : true # 2000-2020 + --> QFED2 : false # 2000-2020 --> GFAS : false # 2003-2021 --> BB4MIPS : false # 1850-2100 # ----- NON-EMISSIONS DATA ------------------------------------ @@ -139,6 +139,11 @@ Mask fractions: false --> Monoterp to SOAS : 0.050 --> Othrterp to SOAP : 0.050 --> Othrterp to SOAS : 0.050 +111 GFED : on CH4/CO/CO2 + --> GFED4 : true + --> GFED_daily : true + --> GFED_3hourly : false + --> Scaling_CO : 1.05 ### END SECTION EXTENSION SWITCHES ### @@ -838,6 +843,28 @@ Mask fractions: false 108 CLM4_PFT_CROP $ROOT/MEGAN/v2018-05/CLM4_PFT.geos.025x03125.v201805.nc PFT_CROP 2000/1/1/0 C xy 1 * - 1 1 )))MEGAN +#============================================================================== +# --- GFED biomass burning emissions (Extension 111) +# NOTE: These are the base emissions in kgDM/m2/s. +#============================================================================== + +(((GFED4 +111 GFED_TEMP $ROOT/GFED4/v2023-03/$YYYY/GFED4_gen.025x025.$YYYY$MM.nc DM_TEMP 2010-2023/1-12/01/0 RF xy kgDM/m2/s * - 1 1 +111 GFED_AGRI $ROOT/GFED4/v2023-03/$YYYY/GFED4_gen.025x025.$YYYY$MM.nc DM_AGRI 2010-2023/1-12/01/0 RF xy kgDM/m2/s * - 1 1 +111 GFED_DEFO $ROOT/GFED4/v2023-03/$YYYY/GFED4_gen.025x025.$YYYY$MM.nc DM_DEFO 2010-2023/1-12/01/0 RF xy kgDM/m2/s * - 1 1 +111 GFED_BORF $ROOT/GFED4/v2023-03/$YYYY/GFED4_gen.025x025.$YYYY$MM.nc DM_BORF 2010-2023/1-12/01/0 RF xy kgDM/m2/s * - 1 1 +111 GFED_PEAT $ROOT/GFED4/v2023-03/$YYYY/GFED4_gen.025x025.$YYYY$MM.nc DM_PEAT 2010-2023/1-12/01/0 RF xy kgDM/m2/s * - 1 1 +111 GFED_SAVA $ROOT/GFED4/v2023-03/$YYYY/GFED4_gen.025x025.$YYYY$MM.nc DM_SAVA 2010-2023/1-12/01/0 RF xy kgDM/m2/s * - 1 1 + +(((GFED_daily +111 GFED_FRAC_DAY $ROOT/GFED4/v2023-03/$YYYY/GFED4_dailyfrac_gen.025x025.$YYYY$MM.nc GFED_FRACDAY 2010-2023/1-12/1-31/0 RF xy 1 * - 1 1 +)))GFED_daily + +(((GFED_3hourly +111 GFED_FRAC_3HOUR $ROOT/GFED4/v2023-03/$YYYY/GFED4_3hrfrac_gen.025x025.$YYYY$MM.nc GFED_FRAC3HR 2010-2023/1-12/1/0-23 RF xy 1 * - 1 1 +)))GFED_3hourly +)))GFED4 + )))EMISSIONS ############################################################################### @@ -878,7 +905,8 @@ Mask fractions: false # --- Quantities needed for CO chemistry --- #============================================================================== (((GLOBAL_OH -${RUNDIR_GLOBAL_OH} +# Revert to GEOS-Chem v5 OH for validation with carbon simulation +* GLOBAL_OH $ROOT/OH/v2022-11/v5-07-08/OH_3Dglobal.geos5.72L.4x5.nc OH 1985/1-12/1/0 C xyz kg/m3 * - 1 1 )))GLOBAL_OH (((PROD_CO_CH4 ${RUNDIR_PCO_CH4} diff --git a/run/GCClassic/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon b/run/GCClassic/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon index 47504a22c..6d1f0b2e9 100644 --- a/run/GCClassic/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon +++ b/run/GCClassic/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon @@ -60,8 +60,8 @@ EmisCO2_BiomassBurn CO2 111 -1 -1 2 kg/m2/s CO2_biomass_burning_emi ############################################################################### ##### OCS sources ##### ############################################################################### -EmisOCS_Total OCS 0 -1 -1 3 kg/m2/s OCS_emission_flux_from_all_sectors +EmisOCS_Total OCS 0 -1 -1 2 kg/m2/s OCS_emission_flux_from_all_sectors EmisOCS_Anthro OCS 0 1 -1 2 kg/m2/s OCS_emission_flux_from_anthropogenic -EmisOCS_Bioburn OCS 0 2 -1 2 kg/m2/s OCS_emission_flux_from_biomass_burning +EmisOCS_BiomassBurn OCS 0 2 -1 2 kg/m2/s OCS_emission_flux_from_biomass_burning EmisOCS_MissingOcean OCS 0 3 -1 2 kg/m2/s OCS_emission_flux_from_missing_ocean EmisOCS_Ocean OCS 0 4 -1 2 kg/m2/s OCS_emission_flux_from_ocean diff --git a/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.CO2 b/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.CO2 index 54e8554a8..24fd480d1 100644 --- a/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.CO2 +++ b/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.CO2 @@ -73,7 +73,7 @@ operations: CO2_simulation_options: sources: - 3D_chemical_oxidation_source: true + use_archived_PCO2_from_CO: true tagged_species: tag_bio_and_ocean_CO2: false diff --git a/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.carbon b/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.carbon index c35423f8b..fcf79603d 100644 --- a/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.carbon +++ b/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.carbon @@ -86,8 +86,8 @@ CH4_simulation_options: #============================================================================ CO_simulation_options: - use_fullchem_PCO_from_CH4: true - use_fullchem_PCO_from_NMVOC: true + use_archived_PCO_from_CH4: false + use_archived_PCO_from_NMVOC: true #============================================================================ # Options for CO2 @@ -96,7 +96,7 @@ CO_simulation_options: CO2_simulation_options: sources: - 3D_chemical_oxidation_source: true + use_archived_PCO2_from_CO: false tagged_species: tag_bio_and_ocean_CO2: false diff --git a/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.tagCO b/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.tagCO index 6ac1a9500..e47984860 100644 --- a/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.tagCO +++ b/run/GCClassic/geoschem_config.yml.templates/geoschem_config.yml.tagCO @@ -84,8 +84,8 @@ operations: #============================================================================ CO_simulation_options: - use_fullchem_PCO_from_CH4: true - use_fullchem_PCO_from_NMVOC: true + use_archived_PCO_from_CH4: true + use_archived_PCO_from_NMVOC: true #============================================================================ # Settings for diagnostics (other than HISTORY and HEMCO) diff --git a/run/GCHP/ExtData.rc.templates/ExtData.rc.carbon b/run/GCHP/ExtData.rc.templates/ExtData.rc.carbon index c499b9161..a815fd6c1 100644 --- a/run/GCHP/ExtData.rc.templates/ExtData.rc.carbon +++ b/run/GCHP/ExtData.rc.templates/ExtData.rc.carbon @@ -292,7 +292,7 @@ FOSSILCO2_ODIAC kg/m2/s N Y F%y4-%m2-01T00:00:00 none none CO2 ./HcoDir/CO2/v2 OCEANCO2_SCALED_MONTHLY kg/m2/s N Y F%y4-%m2-01T00:00:00 none none CO2 ./HcoDir/CO2/v2022-11/OCEAN/Scaled_Ocean_CO2_monthly.nc # --- Balanced biosphere exchange --- -SIB_BBIO_CO2 kg/m2/s N Y F1985-%m2-%d2T%h2-01-01 none none CO2 ./HcoDir/CO2/v2022-11/BIO/SiB3_3hr_NEP.nc +SIB_BBIO_CO2 kg/m2/s N Y F%y4-%m2-%d2T%h2:00:00 none none CO2 ./HcoDir/CO2/v2022-11/BIO/SiB3_3hr_NEP.nc # --- Net Terrestrial Exchange --- CO2_NET_TERRESTRIAL kg/m2/s N Y - none none CO2 ./HcoDir/CO2/v2022-11/BIO/Net_terrestrial_exch_5.29Pg.generic.1x1.nc @@ -319,10 +319,10 @@ CO2_COPROD kgC/m3/s N Y F%y4-%m2-01T00:00:00 none none LCO ./HcoDir/CO2/ #============================================================================== # --- OCS emission fluxes --- #============================================================================== -OCS_ANTHRO kgS/km2/s N Y F%y4-%m2-01T00:00:00 none none COS_Flux ./HcoDir/OCS_FLUX/v2022-11/anthro_v3/v3_anthro/%y4/%m2.nc -OCS_BIOMASS kgS/km2/s N Y F%y4-%m2-01T00:00:00 none none COS_Flux ./HcoDir/OCS_FLUX/v2022-11/Biomassburn/CMSV2-CO2-GFED3-ZC/%y4/%m2.nc -OCS_MISSING_OCEAN kgS/km2/s N Y F%y4-%m2-01T00:00:00 none none COS_Flux ./HcoDir/OCS_FLUX/v2022-11/MissingOcean/LUKAI/%y4/%m2.nc -OCS_OCEAN kgS/km2/s N Y F%y4-%m2-01T00:00:00 none none COS_Flux ./HcoDir/OCS_FLUX/v2022-11/OceanCOS/Kettle/%y4/%m2.nc +OCS_ANTHRO kgS/km2/s N Y F%y4-%m2-01T00:00:00 none 1.0e-6 COS_Flux ./HcoDir/OCS_FLUX/v2022-11/anthro_v3/v3_anthro/%y4/%m2.nc +OCS_BIOMASS kgS/km2/s N Y F%y4-%m2-01T00:00:00 none 1.0e-6 COS_Flux ./HcoDir/OCS_FLUX/v2022-11/Biomassburn/CMSV2-CO2-GFED3-ZC/%y4/%m2.nc +OCS_MISSING_OCEAN kgS/km2/s N Y F%y4-%m2-01T00:00:00 none 1.0e-6 COS_Flux ./HcoDir/OCS_FLUX/v2022-11/MissingOcean/LUKAI/%y4/%m2.nc +OCS_OCEAN kgS/km2/s N Y F%y4-%m2-01T00:00:00 none 1.0e-6 COS_Flux ./HcoDir/OCS_FLUX/v2022-11/OceanCOS/Kettle/%y4/%m2.nc ############################################################################### ### diff --git a/run/GCHP/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon b/run/GCHP/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon index a7e4c97ce..89c8c66e8 100644 --- a/run/GCHP/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon +++ b/run/GCHP/HEMCO_Config.rc.templates/HEMCO_Config.rc.carbon @@ -240,8 +240,8 @@ Mask fractions: false 0 GHGI_COAST_WASTEWATER_IND - - - - - - CH4 1009 6 1 ### Rice ### -0 GHGI_RICE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_3C_Rice_Cultivation 2012-2018/1-12/1/0 EFY xy molec/cm2/s CH4 58/1008 7 100 -0 GHGI_COAST_RICE - - - - - - CH4 58/1009 7 1 +0 GHGI_RICE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_3C_Rice_Cultivation 2012-2018/1-12/1/0 EFY xy molec/cm2/s CH4 58/1008 7 100 +0 GHGI_COAST_RICE - - - - - - CH4 58/1009 7 1 ### Other Anthro ### 0 GHGI_OTHER_MCOMB $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1A_Combustion_Mobile 2012-2018/1/1/0 EFY xy molec/cm2/s CH4 1008 8 100 @@ -324,7 +324,7 @@ Mask fractions: false ### Rice ### 0 GHGI_EE_RICE $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_3C_Rice_Cultivation 2012-2020/1-12/1/0 C xy molec/cm2/s CH4 58/1008 7 100 -0 GHGI_EE_COAST_RICE - - - - - - CH4 58/1009 7 1 +0 GHGI_EE_COAST_RICE - - - - - - CH4 58/1009 7 1 ### Other Anthro ### 0 GHGI_EE_OTHER_MCOMB $ROOT/CH4/v2023-07/Gridded_GHGI_v2/Express_Extension_Gridded_GHGI_Methane_v2_$YYYY.nc emi_ch4_1A_Combustion_Mobile 2012-2020/1/1/0 C xy molec/cm2/s CH4 1008 8 100 @@ -399,7 +399,7 @@ Mask fractions: false )))Scarpelli_Canada #============================================================================== -# --- CH4: Global Fuel Exploitation Inventory (GFEI v2, Scarpelli et al., 2021) +# --- CH4: Global Fuel Exploitation Inventory (GFEI v2, Scarpelli et al., 2021) --- # # This inventory will replace EDGAR (oil, gas, & coal) #============================================================================== @@ -481,11 +481,11 @@ Mask fractions: false )))BB4MIPS #============================================================================== -# --- CH4: biomass burning --- +# --- CH4: QFED2 biomass burning --- +# +# Use QFED v2.4r8 for CH4. The data in QFED/v2014-09 are for v2.4r5 and have +# been found to be buggy for CH4 (same emissions on each day). #============================================================================== - -# --- QFED v2.4r8 --- -# NOTE: Use this instead of data in QFED/v2014-09 (these are buggy) (((QFED2 0 QFED_CH4 $ROOT/CH4/v2017-10/QFED/QFEDv2.4r8.emis_ch4.$YYYY.nc emi_ch4 2009-2015/1-12/1-31/0 C xy molec/cm2/s CH4 - 9 2 )))QFED2 @@ -986,19 +986,19 @@ Mask fractions: false (((USE_OCS_DATA (((OCS_ANTHRO_FLUX -0 OCS_ANTHRO $ROOT/OCS_FLUX/v2022-11/anthro_v3/v3_anthro/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 1 1 +0 OCS_ANTHRO $ROOT/OCS_FLUX/v2022-11/anthro_v3/v3_anthro/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 1 1 )))OCS_ANTHRO_FLUX (((OCS_BIOMASS_FLUX -0 OCS_BIOMASS $ROOT/OCS_FLUX/v2022-11/Biomassburn/CMSV2-CO2-GFED3-ZC/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 2 1 +0 OCS_BIOMASS $ROOT/OCS_FLUX/v2022-11/Biomassburn/CMSV2-CO2-GFED3-ZC/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 2 1 )))OCS_BIOMASS_FLUX (((OCS_MISSING_OCEAN_FLUX -0 OCS_MISSING_OCEAN $ROOT/OCS_FLUX/v2022-11/MissingOcean/LUKAI/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 3 1 +0 OCS_MISSING_OCEAN $ROOT/OCS_FLUX/v2022-11/MissingOcean/LUKAI/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 3 1 )))OCS_MISSING_OCEAN_FLUX (((OCS_OCEAN_FLUX -0 OCS_OCEAN $ROOT/OCS_FLUX/v2022-11/OceanCOS/Kettle/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70/71 4 1 +0 OCS_OCEAN $ROOT/OCS_FLUX/v2022-11/OceanCOS/Kettle/$YYYY/$MM.nc COS_Flux 2012-2016/1-12/1/0 C xy kgS/km2/s OCS 70 4 1 )))OCS_OCEAN_FLUX )))USE_OCS_DATA @@ -1398,7 +1398,7 @@ ${RUNDIR_CO2_COPROD} (((EMISSIONS #============================================================================== -# --- Scale factors for CH4 emissions ---- +# --- CH4: Seasonal scaling factors ---- #============================================================================== (((USE_CH4_DATA @@ -1419,7 +1419,7 @@ ${RUNDIR_CO2_COPROD} )))GHGI_v2.or.GHGI_v2_Express_Ext #------------------------------------------------------------------------------ -# --- Scarpelli et al Mexico manure & rice scale factors --- +# --- Scarpelli_Mexico manure & rice scale factors --- #------------------------------------------------------------------------------ (((Scarpelli_Mexico 10 MANURE_SF $ROOT/CH4/v2017-10/Seasonal_SF/EMICH4_Manure_ScalingFactors.WithClimatology.nc sf_ch4 2008-2016/1-12/1/0 C xy 1 1 @@ -1555,11 +1555,6 @@ ${RUNDIR_CO2_COPROD} #------------------------------------------------------------------------------ 70 StoOCS MATH:60.07/32.06 - - - xy unitless 1 -#------------------------------------------------------------------------------ -# --- km-2 to m-2 --- -#------------------------------------------------------------------------------ -71 Km2ToM2 1.0e-6 - - - xy unitless 1 - )))USE_OCS_DATA #============================================================================== diff --git a/run/GCHP/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon b/run/GCHP/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon index 73c1d3129..cb3cbfac7 100644 --- a/run/GCHP/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon +++ b/run/GCHP/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.carbon @@ -60,7 +60,7 @@ EmisCO2_BiomassBurn CO2 111 -1 -1 2 kg/m2/s CO2_biomass_burning_emi ############################################################################### ##### OCS sources ##### ############################################################################### -EmisOCS_Total OCS 0 -1 -1 3 kg/m2/s OCS_emission_flux_from_all_sectors +EmisOCS_Total OCS 0 -1 -1 2 kg/m2/s OCS_emission_flux_from_all_sectors EmisOCS_Anthro OCS 0 1 -1 2 kg/m2/s OCS_emission_flux_from_anthropogenic EmisOCS_Bioburn OCS 0 2 -1 2 kg/m2/s OCS_emission_flux_from_biomass_burning EmisOCS_MissingOcean OCS 0 3 -1 2 kg/m2/s OCS_emission_flux_from_missing_ocean diff --git a/run/GCHP/geoschem_config.yml.templates/geoschem_config.yml.carbon b/run/GCHP/geoschem_config.yml.templates/geoschem_config.yml.carbon index 7c68d6a07..11e4ed0b2 100644 --- a/run/GCHP/geoschem_config.yml.templates/geoschem_config.yml.carbon +++ b/run/GCHP/geoschem_config.yml.templates/geoschem_config.yml.carbon @@ -31,6 +31,7 @@ operations: use_non_local_pbl: ${RUNDIR_USE_NLPBL} transport: + activate: true transported_species: - CH4 - CO diff --git a/run/shared/singleCarbonSpecies.sh b/run/shared/singleCarbonSpecies.sh index a727dc2e8..ac9a166ff 100755 --- a/run/shared/singleCarbonSpecies.sh +++ b/run/shared/singleCarbonSpecies.sh @@ -31,7 +31,7 @@ function isItemInList() { # # Arguments: # ${1}: The item - # ${2}: Thie list + # ${2}: The list # # Returns (via $?) # 0 if item is in the list @@ -110,17 +110,22 @@ function updateGeosChemConfig() { sed -i -e "${cmd}" "${file}" done - # NOTE: CH4 options are already deactivated - # in the out-of-the-box geoschem_config.yml - - # If CO2 is in the exclude list, turn off CO2 options + # If CO2 is in the include list, turn on CO2 production options isItemInList "CO2" "${1}" - if [[ $? == 0 ]]; then - keys=("3D_chemical_oxidation_source" \ - "tag_bio_and_ocean_CO2" \ - "tag_land_fossil_fuel_CO2" ) + if [[ $? == 1 ]]; then + keys=("use_archived_PCO2_from_CO" ) + for key in ${keys[@]}; do + keyValueUpdate "${key}" "false" "true" "${file}" + done + fi + + # If CO is in the include list, turn on CO production options + isItemInList "CO" "${1}" + if [[ $? == 1 ]]; then + keys=("use_archived_PCO_from_CH4" \ + "use_archived_PCO_from_NMVOC" ) for key in ${keys[@]}; do - keyValueUpdate "${key}" "true" "false" "${file}" + keyValueUpdate "${key}" "false" "true" "${file}" done fi }