Skip to content

Commit

Permalink
Apply vertical emission distributions to CEDS energy, industry, and p…
Browse files Browse the repository at this point in the history
…ower emissions

Barron Henderson provided scale factors for vertically distributing CEDS
anthropogenic emissions. To vertically allocate emissions via scale factors,
we copy emissions vertically. Scale factors can then be applied to determine
how much emissions to apply to each level.

In HEMCO_Config.rc the ratio of emissions to apply to each level for CEDS
energy, industry, and ship emissions is read from the gc_levels.nc file:

  #==============================================================================
  # --- CEDS vertical partitioning ---
  #==============================================================================
  (((CEDS
  315 ENERGY_LEVS   gc_layers.nc g_energy   2017/1/1/0 C xyz 1 1
  316 INDUSTRY_LEVS gc_layers.nc g_industry 2017/1/1/0 C xyz 1 1
  317 SHIP_LEVS     gc_layers.nc cmv_c3     2017/1/1/0 C xyz 1 1
  )))CEDS

These scale factors are applied to the CEDS_*_ENE, CEDS_*_IND, and
CEDS_*_SHIP fields. These fields are 2D in the CEDS files, but we now
can specify dimension xyL* instead of xy to tell HEMCO to copy the field
into each emissions level:

  0 CEDS_CO_ENE $ROOT/CEDS/v2020-08/$YYYY/CO-em-total-anthro_CEDS_$YYYY.nc  CO_ene  1970-2017/1-12/1/0 C xyL* kg/m2/s CO 26/37/35/315 1  5
  0 CEDS_CO_IND $ROOT/CEDS/v2020-08/$YYYY/CO-em-total-anthro_CEDS_$YYYY.nc  CO_ind  1970-2017/1-12/1/0 C xyL* kg/m2/s CO 26/316       1  5
  0 CEDS_CO_SHP $ROOT/CEDS/v2020-08/$YYYY/CO-em-total-anthro_CEDS_$YYYY.nc  CO_shp  1970-2017/1-12/1/0 C xyL* kg/m2/s CO 26/317       10 5

To properly save out these emissions, the dimensions in HEMCO_Diagn.rc
have be changed from dimension = 2 to 3. For example:

  EmisCO_Anthro      CO     0    1   -1   3   kg/m2/s  CO_emission_flux_from_anthropogenic
  EmisCO_Ship        CO     0    10  -1   3   kg/m2/s  CO_emission_flux_from_ships
  InvCEDS_CO         CO     0    1   5    3   kg/m2/s  CO_emission_flux_from_CEDS_inventory

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
  • Loading branch information
msulprizio committed May 30, 2021
1 parent aa76fb7 commit 87eb68e
Show file tree
Hide file tree
Showing 9 changed files with 1,082 additions and 1,077 deletions.
220 changes: 110 additions & 110 deletions run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.aerosol

Large diffs are not rendered by default.

708 changes: 354 additions & 354 deletions run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.fullchem

Large diffs are not rendered by default.

136 changes: 67 additions & 69 deletions run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.tagCO

Large diffs are not rendered by default.

32 changes: 16 additions & 16 deletions run/GCClassic/HEMCO_Diagn.rc.templates/HEMCO_Diagn.rc.aerosol
Expand Up @@ -36,11 +36,11 @@ EmisBCPI_Total BCPI -1 -1 -1 3 kg/m2/s BCPI_emission_flux_from_a
EmisBCPI_Aircraft BCPI 0 20 -1 3 kg/m2/s BCPI_emission_flux_from_aircraft
EmisBCPI_Anthro BCPI 0 1 -1 3 kg/m2/s BCPI_emission_flux_from_anthropogenic
EmisBCPI_BioBurn BCPI 111 -1 -1 2 kg/m2/s BCPI_emission_flux_from_biomass_burning
EmisBCPI_Ship BCPI 0 10 -1 2 kg/m2/s BCPI_emission_flux_from_ships
EmisBCPO_Total BCPO -1 -1 -1 2 kg/m2/s BCPO_emission_flux_from_all_sectors
EmisBCPO_Anthro BCPO 0 1 -1 2 kg/m2/s BCPO_emission_flux_from_anthropogenic
EmisBCPI_Ship BCPI 0 10 -1 3 kg/m2/s BCPI_emission_flux_from_ships
EmisBCPO_Total BCPO -1 -1 -1 3 kg/m2/s BCPO_emission_flux_from_all_sectors
EmisBCPO_Anthro BCPO 0 1 -1 3 kg/m2/s BCPO_emission_flux_from_anthropogenic
EmisBCPO_BioBurn BCPO 111 -1 -1 2 kg/m2/s BCPO_emission_flux_from_biomass_burning
EmisBCPO_Ship BCPO 0 10 -1 2 kg/m2/s BCPO_emission_flux_from_ships
EmisBCPO_Ship BCPO 0 10 -1 3 kg/m2/s BCPO_emission_flux_from_ships

###############################################################################
##### DMS emissions #####
Expand Down Expand Up @@ -77,33 +77,33 @@ EmisMTPO_Biogenic MTPO 0 3 -1 2 kg/m2/s MTPO_emission_flux_from_b
###############################################################################
##### NH3 emissions #####
###############################################################################
EmisNH3_Total NH3 -1 -1 -1 2 kg/m2/s NH3_emission_flux_from_all_sectors
EmisNH3_Anthro NH3 0 1 -1 2 kg/m2/s NH3_emission_flux_from_anthropogenic
EmisNH3_Total NH3 -1 -1 -1 3 kg/m2/s NH3_emission_flux_from_all_sectors
EmisNH3_Anthro NH3 0 1 -1 3 kg/m2/s NH3_emission_flux_from_anthropogenic
EmisNH3_BioBurn NH3 111 -1 -1 2 kg/m2/s NH3_emission_flux_from_biomass_burning
EmisNH3_Natural NH3 0 3 -1 2 kg/m2/s NH3_emission_flux_from_natural_sources
EmisNH3_Seabirds NH3 0 30 -1 2 kg/m2/s NH3_emission_flux_from_seabirds
EmisNH3_Ship NH3 0 10 -1 2 kg/m2/s NH3_emission_flux_from_ships
EmisNH3_Ship NH3 0 10 -1 3 kg/m2/s NH3_emission_flux_from_ships

###############################################################################
##### OCPI and OCPO emissions #####
###############################################################################
EmisOCPI_Total OCPI -1 -1 -1 3 kg/m2/s OCPI_emission_flux_from_all_sectors
EmisOCPI_Aircraft OCPI 0 20 -1 3 kg/m2/s OCPI_emission_flux_from_aircraft
EmisOCPI_Anthro OCPI 0 1 -1 2 kg/m2/s OCPI_emission_flux_from_anthropogenic
EmisOCPI_Anthro OCPI 0 1 -1 3 kg/m2/s OCPI_emission_flux_from_anthropogenic
EmisOCPI_BioBurn OCPI 111 -1 -1 2 kg/m2/s OCPI_emission_flux_from_biomass_burning
EmisOCPI_Ship OCPI 0 10 -1 2 kg/m2/s OCPI_emission_flux_from_ships
EmisOCPO_Total OCPO -1 -1 -1 2 kg/m2/s OCPO_emission_flux_from_all_sectors
EmisOCPO_Anthro OCPO 0 1 -1 2 kg/m2/s OCPO_emission_flux_from_anthropogenic
EmisOCPI_Ship OCPI 0 10 -1 3 kg/m2/s OCPI_emission_flux_from_ships
EmisOCPO_Total OCPO -1 -1 -1 3 kg/m2/s OCPO_emission_flux_from_all_sectors
EmisOCPO_Anthro OCPO 0 1 -1 3 kg/m2/s OCPO_emission_flux_from_anthropogenic
EmisOCPO_BioBurn OCPO 111 -1 -1 2 kg/m2/s OCPO_emission_flux_from_biomass_burning
EmisOCPO_Ship OCPO 0 10 -1 2 kg/m2/s OCPO_emission_flux_from_ships
EmisOCPO_Ship OCPO 0 10 -1 3 kg/m2/s OCPO_emission_flux_from_ships

###############################################################################
##### pFe emissions #####
###############################################################################
EmispFe_Total pFe -1 -1 -1 3 kg/m2/s pFe_emission_flux_from_all_sectors
EmispFe_Anthro pFe 0 1 -1 3 kg/m2/s pFe_emission_flux_from_anthropogenic
EmispFe_BioBurn pFe 111 -1 -1 2 kg/m2/s pFe_emission_flux_from_GFED_inventory
EmispFe_Ship pFe 0 10 -1 2 kg/m2/s pFe_emission_flux_from_ships
EmispFe_Ship pFe 0 10 -1 3 kg/m2/s pFe_emission_flux_from_ships

###############################################################################
##### SESQ emissions #####
Expand All @@ -130,15 +130,15 @@ EmisSO2_Anthro SO2 0 1 -1 3 kg/m2/s SO2_emission_flux_from_an
EmisSO2_BioBurn SO2 111 -1 -1 2 kg/m2/s SO2_emission_flux_from_biomass_burning
EmisSO2_VolcErupt SO2 117 51 -1 3 kg/m2/s SO2_emission_flux_from_eruptive_volcano
EmisSO2_VolcDegas SO2 117 52 -1 3 kg/m2/s SO2_emission_flux_from_noneruptive_volcano
EmisSO2_Ship SO2 0 10 -1 2 kg/m2/s SO2_emission_flux_from_ships
EmisSO2_Ship SO2 0 10 -1 3 kg/m2/s SO2_emission_flux_from_ships

###############################################################################
##### SO4 emissions #####
###############################################################################
EmisSO4_Total SO4 -1 -1 -1 3 kg/m2/s SO4_emission_flux_from_all_sectors
EmisSO4_Aircraft SO4 0 20 -1 3 kg/m2/s SO4_emission_flux_from_aircraft
EmisSO4_Anthro SO4 0 1 -1 3 kg/m2/s SO4_emission_flux_from_anthropogenic
EmisSO4_Ship SO4 0 10 -1 2 kg/m2/s SO4_emission_flux_from_ship
EmisSO4_Ship SO4 0 10 -1 3 kg/m2/s SO4_emission_flux_from_ship

###############################################################################
##### SOAP sources #####
Expand All @@ -148,7 +148,7 @@ EmisSOAP_Aircraft SOAP 0 20 -1 3 kg/m2/s SOAP_emission_flux_from_a
EmisSOAP_Anthro SOAP 0 1 -1 3 kg/m2/s SOAP_emission_flux_from_anthropogenic
EmisSOAP_Biogenic SOAP 0 3 -1 2 kg/m2/s SOAP_emission_flux_from_monoterpenes
EmisSOAP_BioBurn SOAP 111 -1 -1 2 kg/m2/s SOAP_emission_flux_from_biomass_burning
EmisSOAP_Ship SOAP 0 10 -1 2 kg/m2/s SOAP_emission_flux_from_ships
EmisSOAP_Ship SOAP 0 10 -1 3 kg/m2/s SOAP_emission_flux_from_ships

###############################################################################
##### SOAS sources #####
Expand Down

0 comments on commit 87eb68e

Please sign in to comment.