Difference between revisions of "GEOS-Chem v11-02"

From Geos-chem
Jump to: navigation, search
(What's new in this version)
 
(226 intermediate revisions by 4 users not shown)
Line 1: Line 1:
== Overview ==
 
  
=== History ===
+
----
 +
----
 +
<big><strong>GEOS-Chem v11-02-final</strong> '''will also carry the designation''' <strong>GEOS-Chem 12.0.0</strong>'''.'''  We are migrating to a purely numeric versioning system in order to adhere more closely to software development best practices. For a complete description of the new versioning system, please see [[GEOS-Chem version numbering system|our ''GEOS-Chem version numbering system'' wiki page]].</big>
 +
----
 +
----
 +
 
 +
 
 +
== Version history ==
  
 
The table below shows the previous, current, and successive versions of GEOS-Chem:
 
The table below shows the previous, current, and successive versions of GEOS-Chem:
Line 7: Line 13:
 
{| border=1 cellspacing=0 cellpadding=5  
 
{| border=1 cellspacing=0 cellpadding=5  
 
|-align="center"
 
|-align="center"
!width="250px"|Previous version
+
!width="220px"|Previous version
!width="250px" bgcolor="#CCFFFF"|This version
+
!width="220px" bgcolor="#CCFFFF"|This version
!width="250px"|Next version
+
!width="220px"|Next version
  
|-align="center"
+
|-align="center" valign="top"
|[[GEOS-Chem v11-01]]
+
|[[GEOS-Chem_v11-01#v11-01_public_release|GEOS-Chem v11-01-public]]
|bgcolor="#CCFFFF"|'''GEOS-Chem v11-02'''
+
|bgcolor="#CCFFFF"|[[GEOS-Chem_v11-02#GEOS-Chem_v11-02-final_.28aka_12.0.0.29|'''GEOS-Chem 12.0.0<br>(aka v11-02-final)''']]<br>
|TBD
+
|[[GEOS-Chem 12#12.1.0|'''GEOS-Chem 12.1.0''']]
  
|-align="center"
+
|-align="center" valign="top"
|PUBLIC RELEASE 01 Feb 2017
+
|Released 01 Feb 2017
|bgcolor="#CCFFFF"|'''TBD'''
+
|bgcolor="#CCFFFF"|'''Released 10 Aug 2018'''
|TBD
+
|Release TBD
 +
 
 +
|-align="center" valign="top"
 +
|[https://geoschem.github.io/gcclassic-manpage-archive/man.v11-01/index.html GC v11-01 online manual]
 +
|bgcolor="#CCFFFF"|'''[https://geoschem.github.io/gcclassic-manpage-archive/man.GC_12/index.html GC 12 online manual]'''
 +
|[[GEOS-Chem 12 benchmark history|GC 12 benchmark history]]
  
|-align="center"
 
|[[GEOS-Chem v11-01 benchmark history|View v11-01 benchmark history]]
 
|bgcolor="#CCFFFF"|'''[[GEOS-Chem v11-02 benchmark history|View v11-02 benchmark history]]'''
 
|TBD
 
 
|}
 
|}
  
=== What's new in this version ===
+
The sections below describe each GEOS-Chem version in the v11-02 series. Versions v11-02a, v11-02b, etc. were validated with 1-month benchmark simulations (and 1-year benchmark simulations if necessary).
  
''NOTE: v11-02a, v11-02b, etc. denote 1-month benchmark simulations, which are designed to evaluate GEOS-Chem's performance at intermediate stages of development.''
+
=== GEOS-Chem v11-02-final (aka 12.0.0) ===
  
==== v11-02a ====
+
<span style="color:green">'''''This version was released on 10 Aug 2018'''''</span>
  
Please see the following link for complete information about the validation of GEOS-Chem v11-02a:
+
For a list of the features that were added to GEOS-Chem during the v11-02 "release candidate" testing period, please see visit [[GEOS-Chem_12#12.0.0|this section of our ''GEOS-Chem 12'' wiki page]].
#[[GEOS-Chem_v11-02_benchmark_history#v11-02a|''Approval form for 1-month benchmark simulation v11-02a'']]
+
 
#[[GEOS-Chem_v11-02_benchmark_history#v11-02a-Run0|''Results for 1-year benchmark simulation v11-02a-Run0'']]
+
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 14:33, 21 June 2018 (UTC)
#[[GEOS-Chem_v11-02_benchmark_history#v11-02a-Run0|''Results for 1-year benchmark simulation v11-02a-Run1'']] (final version)
+
 
 +
=== GEOS-Chem v11-02 release candidate ===
 +
 
 +
NOTE: These updates were not received in time for the [[#v11-02f|v11-02f]] benchmarks, but were included in the '''v11-02-release-candidate''' version (aka '''v11-02-rc''').
  
 
{| border=1 cellspacing=0 cellpadding=5  
 
{| border=1 cellspacing=0 cellpadding=5  
Line 43: Line 53:
 
!width="200px"|Submitted by
 
!width="200px"|Submitted by
 
!width="100px"|Type
 
!width="100px"|Type
!width="150px"|Status
+
!width="300px"|Status
  
 
|-
 
|-
Line 49: Line 59:
  
 
|-valign="top"
 
|-valign="top"
|[[Updates_in_JPL_Publication_15-10|Update chemistry rate constants based on ''JPL Publication 15-10'']]
+
|[[Secondary_organic_aerosols#Only_add_ISOAAQ_species_to_PM2.5_diagnostics_for_simulations_using_the_complex_SOA_option|Only add ISOAAQ species to the PM2.5 diagnostics when using the Complex SOA option]]
|Barron Henderson (US EPA),<br>Mat Evans (U. York), &<br>[[Chemistry Issues|Oxidants and Chemistry WG]]
+
|Aerosols WG
|Science
+
|Bug fix
|Approved 12 May 2017
+
|''Released 22 Jun 2018''
  
 
|-valign="top"
 
|-valign="top"
|[[Chemistry_Issues#Fixes_to_correct_ALK4_lumping_issue|Fixes to correct ALK4 lumping issue]]
+
|[[Secondary_organic_aerosols#Only add ISOAAQ species to OCPISOA when using the complex SOA option|Only add ISOAAQ species to OCPISOA when using the complex SOA option]]
|Barron Henderson (US EPA)
+
|Aerosols WG
 +
|Bug fix
 +
|" "
 +
 
 +
|-valign="top"
 +
|[[#Disable Xiao propane emissions to avoid double-counting with CEDS|Disable Xiao propane emissions to avoid double-counting with CEDS]]
 +
|Christoph Keller (NASA/GMAO)<br>[[GCST]]
 +
|Bug fix
 +
|" "
 +
 
 +
|-valign="top"
 +
|[[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions#Grid_offset_bug_fix_in_high-resolution_MEGAN_files_at_southern_mid-latitudes|Grid offset bug fix in high-resolution MEGAN files at southern mid-latitudes]]
 +
|Jenny Fisher (Wollongong) <br>Dylan Millet (UMN)
 +
|Bug fix
 +
|" "
 +
 
 +
|-
 +
!colspan="4" bgcolor="#CCFFFF"|Features NOT affecting the full-chemistry simulation:
 +
 
 +
|-valign="top"
 +
|Minor fixes (including updated parameters) for the [[Mercury]] simulation.
 +
|Colin Thackray (Harvard)
 +
|Bug Fix<br>Science
 +
|" "
 +
 
 +
|}
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 19:16, 21 June 2018 (UTC)
 +
 
 +
=== v11-02f ===
 +
 
 +
Please see the following link for complete information about the validation of GEOS-Chem v11-02f:
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02f|''Approval form for 1-month benchmark simulation v11-02f'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02f-Run0|''Results for 1-year benchmark simulation v11-02f-Run0'']]
 +
 
 +
{| border=1 cellspacing=0 cellpadding=5
 +
|-bgcolor="#CCCCCC"
 +
!width="500px"|Feature
 +
!width="200px"|Submitted by
 +
!width="100px"|Type
 +
!width="300px"|Status
 +
 
 +
|-
 +
!colspan="4" bgcolor="#CCFFFF"|Features affecting the full-chemistry simulation:
 +
 
 +
|-valign="top"
 +
|[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02f|Update HEMCO from v2.1.001 to v2.1.005]]
 +
*HEMCO v2.1.005 requires [[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions#Global_0.25x0.25_MEGAN_input_and_MODIS_LAI|updated MODIS LAI data]]
 +
|Christoph Keller (NASA GMAO)
 +
|Structural
 +
|''1-month benchmark:''<br>Approved on 10 May 2018
 +
 
 +
''1-year benchmark:''<br>Approved on 17 May 2018
 +
 
 +
|-valign="top"
 +
|[[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions#Global_0.25x0.25_MEGAN_input_and_MODIS_LAI|Use global 0.25x0.25 files for MEGAN input and MODIS LAI]]
 +
|Jenny Fisher (U. Wollongong)<br>Eloise Marais (U. Birmingham)<br>Dylan Millet (UMN)
 +
|Bug fix/<br>Science
 +
|" "
 +
 
 +
|-valign="top"
 +
|[[Lightning_NOx_emissions|Update lightning BETA values and OTD/LIS files for consistency with latest met fields]]
 +
|Lee Murray (Rochester)
 
|Science
 
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[PAN|PAN updates]], including:
+
|[[CEDS_anthropogenic_emissions|CEDS global anthropogenic emissions (1950-2014)]]
*Updates to emissions, deposition, species, chemistry
+
*This inventory will be the default global anthropogenic emissions inventory
*Ocean+biogenic ALD2 emissions from Millet al. [2010]
+
*Includes a fix for NOx unit conversion
|Emily Fischer (CSU)
+
|Lu Shen (Harvard)<br>Melissa Sulprizio (Harvard)
 
|Science
 
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[EPA/NEI11_North_American_emissions#Monthly_mean_emissions|Monthly mean NEI2011 emissions]]
+
|[[EDGAR_v4.3_anthropogenic_emissions|EDGAR v4.3.1 emissions]]
|[[GCST]] &<br>Katie Travis (Harvard)
+
*This inventory will be included as an option
 +
|Chi Li (Dalhousie)
 
|Science
 
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Bug fixes in the GEOS-Chem sulfate module:
+
|[[GEOS-Chem_v11-02#Use_1x1_GEIA_NH3_emissions|Use 1x1 GEIA NH3 emissions]]
*[[Sulfate_aerosols#Fix_error_in_production_of_SO4s.2C_NITs_in_SEASALT_CHEM_routine|Fix error in production of SO4s and NITs in <tt>SEASALT_CHEM</tt>]]
+
|Eleanor Morris (York)
*[[Sulfate_aerosols#Fix_bug_in_CHEM_NIT_routine|Fix bug in <tt>CHEM_NIT</tt>]]
+
*[[Sulfate_aerosols#Fix_for_sulfate_production_in_HET_DROP_CHEM|Fix for sulfate production in <tt>HET_DROP_CHEM</tt>]]
+
*[[Sulfate_aerosols#Fix_bugs_in_sulfate_chemistry_routines|Fix bugs in sulfate chemistry routines]]
+
|<br>Prasad Kasibhatla (Duke)<br>Prasad Kasibhatla (Duke)<br>Qianjie Chen (UW)<br>Viral Shah (UW)
+
 
|Bug fix
 
|Bug fix
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Dry_deposition#Bug_in_computation_of_aerodynamic_resistance_RA|Fix bug in dry deposition aerodynamic resistance]]
+
|[[CAC_anthropogenic_emissions#Historical_Canadian_emissions|Historical Canadian emissions]]
|Brian Boys (Dalhousie)
+
|Chi Li (Dalhousie)
|Bug fix
+
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Acetone#Fix_acetone_parameterization_in_hcox_seaflux_mod.F90|Fix acetone parameterization in <tt>hcox_seaflux_mod.F90</tt>]]
+
|[[Anthropogenic_emissions#C2H6|Global ethane emissions]]
|[[GCST]]
+
|Zitely Tzompa-Sosa (CSU)
|Bug fix
+
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_v11-01#Bug_in_cos.28SZA.29_for_start_of_timestep|Bug fix in cos(SZA) for start of timestep]]
+
|[[DICE-Africa anthropogenic emissions inventory|Inefficient combustion emissions in Africa (DICE-Africa inventory)]]
|Lizzie Lundgren ([[GCST]])
+
|Eloise Marais (U. Birmingham)<br>Christine Wiedinmyer (NCAR)
|Bug fix
+
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_v11-02#Include_TOA_pressure_when_calculating_dry_pressure_edges|Include TOA pressure when calculating dry pressure edges]]
+
|[[Trash burning emissions]]
|Seb Eastham (Harvard)
+
|Eloise Marais (U. Birmingham)<br>Christine Wiedinmyer (NCAR)
 
|Science
 
|Science
 
|" "
 
|" "
  
|-
+
|-valign="top"
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
+
|[[GFED4_biomass_burning_emissions#Update_GFED4_emission_factors|Update GFED4 emission factors]]
 +
|Katie Travis (MIT)
 +
|Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[ISORROPIA_II#Implement_ISORROPIA_v2.0_as_a_Fortran_module|Implement ISORROPIA v2.0 as a Fortran module]]
+
|[[GFED3_biomass_burning_emissions|Retire GFED3 biomass burning emissions]]
|Seb Eastham (Harvard) &<br>[[GCST]]
+
|Katie Travis (MIT)
 
|Structural
 
|Structural
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Updates to the HEMCO emissions component:
+
|[[NH3_emissions#NH3_emissions_from_arctic_seabirds|Ammonia emissions from Arctic seabirds]]
*[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02a|Update to HEMCO v2.0.004]]
+
|Betty Croft (Dalhousie)
*[[Implementation_of_HEMCO_in_GEOS-Chem#Fixed_bug_in_computation_of_local_time_in_HCO_GetSuncos|Fixed bug in computation of local time in <tt>HCO_GetSuncos</tt>]]
+
|Science
*[[EPA/NEI11_North_American_emissions#Default_US_emissions_to_NEI2011_after_2013|Default US emissions to NEI2011 after 2013]]
+
*[[Implementation_of_HEMCO_in_GEOS-Chem#Read_default_DEP_RESERVOIR_fields_from_file_when_not_found_in_HEMCO_restart_file|Read default <tt>DEP_RESERVOIR</tt> fields from file]]
+
|<br>Christoph Keller (NASA GMAO)<br>Seb Eastham (Harvard)<br>Jessica Morena (Dalhousie)<br>Brian Boys (Dalhousie) & [[GCST]]
+
|Structural<br>& Bug fix
+
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Updates to gain computational speedup:
+
|Fixes for prod/loss rates used in linearized chemistry routines:
*[[FlexChem#In_GEOS-Chem_v11-02_and_later_versions|Add simplified prod/loss families to KPP]]
+
*[[GEOS-Chem_v11-02#Restore_use_of_prod.2Floss_rates_from_GMI_in_the_mesosphere_for_UCX-based_simulations|Restore use of prod/loss rates from GMI in the mesosphere for UCX-based simulations]]
*[[Cloud_convection#Computational_bottleneck_in_the_v11-01_convection_module|Remove computational bottleneck in the convection module]]
+
*[[GEOS-Chem_v11-02#Fix_bug_in_units_of_UCX_P.2FL_rates|Fix bug in units of UCX prod/loss rates]]
|<br>Mike Long (Harvard)<br>Bob Yantosca ([[GCST]])
+
|Chris Holmes (FSU)<br>Christoph Keller (NASA GMAO)
|Structural<br>& Bug fix
+
|Bug fix
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Mercury#Fixed_typo_in_nested_grid_transport_routine_INIT_WINDOW|Fixed typo in <tt>INIT_WINDOW</tt>]]
+
|[[UCX_chemistry_mechanism#Fix_for_UCX_mesospheric_NOy_routine|Add fix for UCX mesospheric NOy routine]]
|Bob Yantosca ([[GCST]])
+
|Sebastian Eastham (MIT)
 
|Bug fix
 
|Bug fix
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|netCDF file I/O updates:
+
|[[Physical_properties_of_GEOS-Chem_species#Fix_Henry.27s_law_parameters_for_HCl|Fix Henry's Law parameters for HCl]]
*[[The NcdfUtilities package#Enable_compression_in_netCDF-4_output_files|Enable data compression in netCDF-4 output files]]
+
|Xuan Wang (Harvard)
*[[The NcdfUtilities package#Routine DO_ERR_OUT now returns a non-zero error code|Routine <tt>DO_ERR_OUT</tt> now returns a non-zero error code]]
+
|Bug fix
*[[Implementation of HEMCO in GEOS-Chem#HEMCO diagnostic and restart files now have an unlimited time dimension|HEMCO files now have an unlimited time dimension]]
+
|<br>Chris Holmes (Florida State)<br>Andy Jacobson (NOAA)<br>[[GCST]]
+
|Bug fix &<br>Structural
+
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Makefile and build sequence updates:
+
|[[Halogen_chemistry_mechanism#Fix_for_N2O5_uptake|Fix for N2O5 uptake]]
*[[GEOS-Chem_Makefile_Structure#Bug_fix:_Specifying_NO_REDUCED.3Dno_now_compiles_GEOS-Chem_for_reduced_grids|Specifying <tt>NO_REDUCED=no</tt> now compiles GEOS-Chem for reduced grids]]
+
|Xuan Wang (Harvard)
*[[GEOS-Chem_Makefile_Structure#Removed_the_COMPILER_variable_from_Makefile_header.mk_for_a_cleaner_build_sequence|Removed the <tt>COMPILER</tt> variable from <tt>Makefile_header.mk</tt>]]
+
|Bug fix
|<br>Jiawei Zhang (Harvard)<br>[[GCST]]
+
|Bug fix &<br>Structural
+
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[UCX_chemistry_mechanism#Bug_fixes_for_running_UCX_in_ESMF_environment|Bug fixes for running UCX in ESMF environment]]
+
|[[GEOS-Chem_v11-02#Fix_for_organic_nitrate_chemistry|Add fix for organic nitrate chemistry]]
|Christoph Keller (NASA GMAO)
+
|Xuan Wang (Harvard)
 
|Bug fix
 
|Bug fix
 
|" "
 
|" "
 +
 +
|-
 +
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
  
 
|-valign="top"
 
|-valign="top"
|Bug fixes for diagnostics:
+
|[[CH4_simulation#Updated_CH4_simulation_in_v11-02_and_later|Updated CH4 simulation]]
*[[Particulate_matter_in_GEOS-Chem#Save_out_PM2.5_diagnostic_at_STP_conditions|Save out PM2.5 diagnostic at STP conditions]]
+
*Includes updates to the [[planeflight diagnostic]]
*[[GEOS-Chem_v11-01#Bugs_in_convective_flux_diagnostics_.28ND14.29|Bug fixes in convective mass flux diagnostic]]
+
|Bram Maasakkers (Harvard)<br>Melissa Sulprizio (Harvard)
*[[GEOS-Chem_v11-01#P.28OH.29_not_saved_out_in_ND22_diagnostic|Restore P(OH) in ND22 diagnostic]]
+
|Science
*[[Implementation_of_HEMCO_in_GEOS-Chem#Make_anthropogenic_emissions_diagnostics_3D|Make anthropogenic emissions diagnostics 3D]]
+
*[[Tagged_CO_simulation#Fix_ND65_bugs_in_tagged_CO_simulation|Fix ND65 bugs in tagged CO simulation]]
+
*[[Mineral_dust_aerosols#Bug_in_ND21_diagnostic_indexing_for_size-resolved_dust_species|Fix bug in ND21 diagnostic indexing for dust species]]
+
|<br>Aaron van Donkelaar (Dalhousie)<br>Jenny Fisher (U. Wollongong)<br>[[GCST]]<br>Jenny Fisher (U. Wollongong)<br>Chris Holmes (Florida State)
+
|Bug fix
+
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Removal of obsolete variables:
+
|[[Tagged_CO_simulation#New_tagged_CO_simulation_in_v11-02_and_later|Updated Tagged CO simulation]]
*[[#Removed_the_NNPAR_parameter_from_CMN_SIZE_mod.F|Removed the <tt>NNPAR</tt> parameter from <tt>CMN_SIZE_mod.F</tt>]]
+
|Jenny Fisher (Wollongong)
*[[Wet_deposition#Removal_of_obsolete_variables_NSOL_and_IDWETD|Removed obsolete variables <tt>NSOL</tt> and <tt>IDWETD</tt>]]
+
|Science
*[[Derived type objects used by GEOS-Chem#Summary of fields removed from Input_Opt|Removed obsolete fields of the <tt>Input_Opt</tt> object]]
+
|[[GCST]]
+
|Structural
+
 
|" "
 
|" "
  
|}
+
|-valign="top"
 +
|[http://www.atmos-chem-phys.net/17/6353/2017/acp-17-6353-2017.pdf Updates to Hg chemistry from Horowitz et al. (2017)]
 +
|Hannah Horowitz (U. Washington)<br>Colin Thackray (Harvard)<br>Rebecca Stern (Harvard)
 +
|Science
 +
|" "
  
==== v11-02b ====
+
|-valign="top"
 
+
|[[MERRA-2_implementation_details#Default_timestep_for_MERRA-2_nested_grid_simulations|Change default timestep for MERRA-2 nested grid simulations]]
Please see the following link for complete information about the validation of GEOS-Chem v11-02b:
+
|Bram Maasakkers (Harvard)
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b|''Approval form for 1-month benchmark simulation v11-02b'']]
+
|Science
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b_HP_1-month|''Approval form for 1-month benchmark simulation v11-02b with high performance option (v11-02b-HP)'']]
+
|" "
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b-HP-Run0|''Results for 1-year standard benchmark simulation of v11-02b-HP Run0'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b-HP-Run1|''Results for 1-year standard benchmark simulation of v11-02b-HP Run1'']] (uses offline archived lightning and dust emissions)
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b_HP_1-year_Rn-Pb-Be|''Results for 1-year RnPbBe benchmark simulation of v11-02b-HP'']]
+
 
+
{| border=1 cellspacing=0 cellpadding=5
+
|-bgcolor="#CCCCCC"
+
!width="500px"|Feature
+
!width="200px"|Submitted by
+
!width="100px"|Type
+
!width="150px"|Status
+
 
+
|-
+
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
+
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem HP|Source code updates for high performance option (GCHP)]]
+
|[[Biomass_burning_emissions#Fix_bug_in_biomass_burning_emissions_when_GFED_is_turned_off|Fix bug in biomass burning emissions when GFED is turned off]]
|Seb Eastham (Harvard)<br>Lizzie Lundgren (Harvard)<br>Mike Long (Harvard)<br>Jiawei Zhuang (Harvard)<br>Bob Yantosca (Harvard)
+
|Corey Trujillo (U. Colorado)<br>Therese Carter (MIT)
|Structural
+
|Bug fix
|Approved 16 Jun 2017
+
|" "
  
 
|-valign="top"
 
|-valign="top"
|Bug fixes for diagnostics:
+
|[[RCP_future_emissions_scenarios#Fix_missing_HCOOH_factors_for_RCP_in_HEMCO_Config.rc|Fix missing HCOOH factors for RCP in HEMCO_Config.rc]]
*[[GEOS-Chem_v11-02#ND26_vertical_flux_diagnostic_does_not_yield_mass_balance_in_advection|Use FZ array to compute ND26 vertical flux diagnostic]]
+
|Seb Eastham (MIT)
*[[GEOS-Chem_v11-02#ND40_planeflight_diagnostic_not_always_written_out_when_turned_on|Always write ND40 planeflight diagnostic when enabled]]
+
|<br>Ilya Stanevic (Toronto)<br>Lee Murray (Rochester)
+
 
|Bug fix
 
|Bug fix
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Tagged_O3_simulation#Write_initial_and_final_Ox_mass_to_file|Write initial and final Ox mass to file when using the tagged O3 simulation]]
+
|[[Tagged_O3_simulation#Fix_bug_in_calculation_of_L.28Ox.29|Fix bug in calculation of L(Ox) for tagged Ox simulations]]
|Bob Yantosca ([[GCST]])
+
|Lu Shen (Harvard)
|Structural
+
|Bug fix
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_v11-02#CO2_emissions_unit_conversion_not_consistent_with_new_advection_pressure|Convert CO2 emissions units using dry pressure used in advection]]
+
|[[GEOS-Chem_v11-02#KPP_errors_with_halogen_chemistry_and_tropchem-based_mechanisms|Change tolerance value as a temporary fix for convergence errors in 2x2.5 tropchem simulations]]
|Meemong Lee (JPL)
+
|Katie Travis (MIT)<br>[[Chemistry_Issues|Chemistry WG]]
|Science
+
|Bug fix
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[QFED_biomass_burning_emissions|Add QFED emissions for 2014-2016]]
+
|[[GEOS-Chem_v11-02#Fix_compatibility_issue_with_ifort17|Fix compatibility issue with ifort17]]
|Christoph Keller (NASA GMAO)
+
|Bob Yantosca ([[GEOS-Chem_Support_Team|GCST]])
|Science
+
|Bug fix
 
|" "
 
|" "
  
 
|}
 
|}
  
==== v11-02c ====
+
=== v11-02e ===
  
Please see the following link for complete information about the validation of GEOS-Chem v11-02c:
+
Please see the following link for complete information about the validation of GEOS-Chem v11-02e:
#[[GEOS-Chem_v11-02_benchmark_history#v11-02c|''Approval form for 1-month benchmark simulation v11-02c'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e|''Approval form for 1-month benchmark simulation v11-02e'']]
#[[GEOS-Chem_v11-02_benchmark_history#v11-02c-Run0|''Results for 1-year benchmark simulation v11-02c-Run0'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e_using_spun-up_restart_file|''Approval form for 1-month benchmark simulation v11-02e using a spun-up restart file'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e-Run0|''Results for 1-year benchmark simulation v11-02e-Run0'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e-RnPbBePasv--Run0|''Results for 1-year Rn-Pb-Be benchmark simulation with fixed GEOS-FP met fields'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e-RnPbBePasv-Run1|''Results for 1-year Rn-Pb-Be benchmark simulation using 72 levels'']]
 +
 
 +
We also completed a second set of benchmark simulations for GEOS-Chem v11-02e using 2016 GEOS-FP met fields. Please see the following links:
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e_using_2016_GEOS-FP_met_fields|''Approval form for 1-month benchmark simulation v11-02e using 2016 GEOS-FP met fields'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e-Run1|''Results for 1-year benchmark simulation v11-02e-Run1'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e-RnPbBePasv-Run2|''Results for 1-year Rn-Pb-Be benchmark simulation using 2016 GEOS-FP met fields'']]
  
 
{| border=1 cellspacing=0 cellpadding=5  
 
{| border=1 cellspacing=0 cellpadding=5  
Line 251: Line 306:
 
!width="200px"|Submitted by
 
!width="200px"|Submitted by
 
!width="100px"|Type
 
!width="100px"|Type
!width="150px"|Status
+
!width="300px"|Status
  
 
|-
 
|-
Line 257: Line 312:
  
 
|-valign="top"
 
|-valign="top"
|[[Secondary_organic_aerosols#Simple_SOA_scheme|Enhance default GEOS-Chem simple SOA]]
+
|[[NOx-Ox-HC-aerosol#Monthly mean surface distributions|Monthly mean surface methane distributions]]
|Sal Farina (Colorado State)<br>Aerosols Working Group Chairs
+
|Lee Murray (Rochester)
 
|Science
 
|Science
|''1-month benchmark:''<br> Approved on 07 Sep 2017
+
|''1-month benchmark:''<br>Approved on 15 Mar 2018
  
''1-year benchmark:''<br> Approved 21 Sep 2017
+
''1-year benchmark:''<br>Approved 24 Mar 2018
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_chemistry_mechanisms#Updated_isoprene_and_monoterpene_chemistry|Updates to isoprene and monoterpene chemistry]]
+
|UCX updates, including:
|Katie Travis (MIT)<br>Jenny Fisher (U. Wollongong)<br>Elo&iuml;se Marais (U. Birmingham)<br>Christopher Chan Miller (Harvard)<br>Kelvin Bates (Caltech)<br>Rebecca Schwantes(Caltech)
+
*[[UCX_chemistry_mechanism#v11-02|Remove initial stratospheric 2D mixing ratio option]]
 +
*[[UCX_chemistry_mechanism#UCX_stratospheric_water_boundary_condition_update|UCX stratospheric water boundary condition update]]
 +
*[[UCX_chemistry_mechanism#Remove_UCX_.23ifdefs_and_replace_with_Input_Opt.25LUCX|Remove <tt>UCX=y</tt> as compile switch and use <tt>Input_Opt%LUCX</tt> instead]]
 +
|<br>Seb Eastham (Harvard)<br>Chris Holmes (Florida State)<br>Melissa Sulprizio ([[GCST]])
 
|Science
 
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Secondary_organic_aerosols#SOA_formation_from_aqueous_isoprene_uptake|Add aqueous isoprene uptake to SOA scheme]]
+
|[[Stratospheric_chemistry#Prod.2Floss_rates_from_UCX|Use monthly mean stratospheric prod/loss rates from UCX simulation for all species]]
|Eloise Marais (U. Birmingham)
+
|Melissa Sulprizio ([[GCST]])
 
|Science
 
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Chemistry_Issues#Fixes_for_carbon_creating_reactions|Carbon balance (fix C creation)]]
+
|[[Aerosol_optical_properties#Add_BC_absorption_enhancement_factor|Update density of BC to 1.8 and add absorption enhancement factor in input.geos]]
|Sarah Safieddine (MIT)
+
|Xuan Wang (MIT)
|Bug fix
+
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[PAN#Bug_fixes_for_EOH_and_MGLY|Fix bugs for EOH and MGLY following implementation of PAN updates in v11-02a]]
+
|[[Sulfate_aerosols#Metal_catalyzed_oxidation_of_SO2|Metal catalyzed oxidation of SO2 as a switch]] (on by default)
|Melissa Sulprizio ([[GCST]])
+
|Viral Shah (UW)<br>Becky Alexander (UW)
|Bug fix
+
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02c|Update HEMCO from v2.0.004 to v2.1.001]]
+
|Add SOA updates and bug fixes, including:
|Christoph Keller (NASA GMAO)
+
*[[GEOS-Chem_v11-02#Remove_MONX_from_chemical_mechanisms|Remove MONX from chemical mechanisms]]
|Structural
+
*[[Secondary_organic_aerosols#Add_MTPO_as_an_advected_and_chemical_species_to_all_full-chemistry_simulations|Add MTPO as a species in full-chemistry simulations]]
 +
*[[Secondary_organic_aerosols#Update_simple_SOA_entries_in_HEMCO_to_follow_MTPA.2BLIMO.2BMTPO|Updating SOAP and SOAS entries in HEMCO to follow MTPA+LIMO+MTPO rather than MONX]]
 +
*[[Secondary_organic_aerosols#Calculation_of_PM2.5.2C_AOD.2C_and_aerosol_mass|Use simple SOA in calculation of PM2.5, AOD, and aerosol mass by default]]
 +
*[[Secondary_organic_aerosols#SOA_formation_from_aqueous_isoprene_uptake|Always use SOA formation from aqueous isoprene uptake (remove switch from input.geos)]]
 +
*[[Secondary_organic_aerosols#Fixes_for_isoprene_SOA_for_consistency_with_Marais_et_al._.282016.29|Fixes for isoprene SOA for consistency with Marais et al. (2016)]]
 +
|Sam Silva (MIT)<br>[[Chemistry_Issues|Chemistry Working Group]]<br>[[Aerosols Working Group]]<br>[[GCST]]<br>Eloise Marais (U. Birmingham)
 +
|Science & Bug fix
 
|" "
 
|" "
  
|-
+
|-valign="top"
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
+
|[[GEOS-Chem_v11-02#Fix_dust_rainout_settings_to_match_BCPO|Fix dust rainout settings to match BCPO]]
 +
|Colette Heald (MIT)
 +
|Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|HEMCO updates:
+
|[[Halogen_chemistry_mechanism#Fixes_in_gckpp_HetRates.F90_following_implementation_in_v11-02d|Fixes in gckpp_HetRates.F90 for halogen chemistry]]
*[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02c|Add a HEMCO-standalone run directory for benchmarking purposes]]
+
|Lei Zhu (Harvard)<br>Seb Eastham (MIT)
*[[Implementation_of_HEMCO_in_GEOS-Chem#Now_use_YYYYMMDDhhmm_format_for_time_stamp_values|Now use <tt>YYYYMMDDhhmm</tt> for time stamp values]]
+
|Bug fix
*[[Implementation_of_HEMCO_in_GEOS-Chem#Avoid_segmentation_fault_in_DustGinoux_extension|Add error trap in to avoid a segmentation fault when the DustGinoux extension is turned off]]
+
|<br>Christoph Keller (NASA GMAO)<br>Andy Jacobson (NOAA)<br>Paulo Tuccella (L'Aquila)
+
|Structural
+
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Fixes for several minor issues:
+
|[[Sulfate_aerosols#Add_bug_fix_in_sulfate_mod.F_involving_reorder_IF_statements_for_Fahey_and_Seinfeld_decision_algorithm|Add bug fix in sulfate_mod.F involving reorder IF statements for Fahey and Seinfeld decision algorithm]]
*[[#Bug fixes for the ND50 timeseries diagnostic|Bug fixes for the ND50 timeseries diagnostic]]
+
|Chris Holmes (FSU)<br>Becky Alexander (UW)
*[[#Fixed_an_incorrect_format_statement_in_input_mod.F|Fixed an incorrect format statement in <tt>input_mod.F</tt>]]
+
*[[#Add MERRA2 to #if and #elif statements where it had been omitted|Add MERRA2 to <tt>#if</tt> and <tt>#elif</tt> statements where it had been omitted]]
+
*[[Mercury#Bug_fixes_for_Hg_emissions|Bug fixes for Hg emissions]]
+
*[[Aerosol_emissions#Bug_fix_for_MOPO_and_MOPI_definitions_in_species_database|Fix bug in species definitions for marine POA simulation]]
+
*[[Secondary_organic_aerosols#Fix_diagnostics_bugs_in_the_SOA-SVPOA_simulation|Fix diagnostic bugs in the SOA-SVPOA simulation]]
+
|<br>Chris Holmes (Florida State)<br>Chris Holmes (Florida State)<br>Jiawei Zhuang (Harvard)<br>Amanda Giang (MIT)<br>Katie Travis (MIT)<br>Katie Travis(MIT)
+
 
|Bug fix
 
|Bug fix
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Fixes for the TOMAS simulation:
+
|[[Downloading_GEOS-Chem_source_code_and_data#Inconsistency_in_GEOS-FP_files_at_Harvard_and_Dalhousie_for_July_2013|Sync GEOS-FP files on Harvard ftp with files at Dalhousie]]
*[[TOMAS_aerosol_microphysics#Fixes_for_TOMAS_simulation_in_v11-02c|Fixes for TOMAS simulation in v11-02c]]
+
|[[GCST]]
*[[TOMAS_aerosol_microphysics#Typo_in_wetscav_mod.F_for_TOMAS30_simulation|Fix typo in wetscav_mod.F for TOMAS30]]
+
*Remove a couple of array temporaries and an out-of-bounds error
+
|<br>Sal Farina (Colorado State)<br>Jack Kodros (Colorado State)<br>Bob Yantosca ([[GCST]])
+
 
|Bug fix
 
|Bug fix
 
|" "
 
|" "
 +
 +
|-
 +
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
  
 
|-valign="top"
 
|-valign="top"
|[[Stratospheric_chemistry#STE_fluxes|Fix STE flux diagnostics and add to benchmark procedure]]
+
|[[Particulate_matter_in_GEOS-Chem#Option_to_include_spatially_and_seasonally_varying_OM.2FOC|Spatially varying OM/OC as a switch]] (off by default)
 +
|Sajeev Philip (Dalhousie/NASA Ames)
 +
|Science
 +
|" "
 +
 
 +
|-valign="top"
 +
|[[Rn-Pb-Be_simulation#Vertical_Grid|Update Rn-Pb-Be-Pasv simulations to use 72 levels by default]]
 
|Melissa Sulprizio ([[GCST]])
 
|Melissa Sulprizio ([[GCST]])
|Bug fix &<br>Benchmarking
+
|Science
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|Initial structural modifications for netCDF diagnostics:
+
|[[ISORROPIA_II#Bug_fixes_for_ISORROPIA_II_stable_mode|Bug fixes for ISORROPIA II stable mode]]
*Introduce <tt>Headers/State_Diag</tt> as a stub module (for now)
+
|Jingyuan Shao (UW)<br>Shaojie Song (Harvard)
*Add a registry object into <tt>State_Met</tt>, <tt>State_Diag</tt>, and <tt>State_Chm</tt> in order to obtain a pointer to any module variable (or slice) by looking up its name
+
|Bug fix
*Add new module <tt>Headers/registry_mod.F90</tt> which contains derived types and routines for registering module variables.
+
|Bob Yantosca ([[GCST]])
+
|Structural
+
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_v11-02#Code_cleanup_suggestions|Combine timestep settings in input.geos in a Timesteps menu]]
+
|[[GEOS-Chem_v11-02#Change_timestep_units_from_minutes_to_seconds|Change input.geos timestep units from minutes to seconds]]
|Melissa Sulprizio ([[GCST]])
+
|Lizzie Lundgren ([[GCST]])
 
|Structural
 
|Structural
 
|" "
 
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_benchmarking#Update_2005.2F2009_CMDL_CO_data_to_2013.2F2014_GMD_CO_data|Update CO data used in 1-year benchmark plots]]
+
|[[GEOS-Chem_v11-02#Fixes_for_two_new_ND21_binary_diagnostics_introduced_in_v11-02|Fixes for ND21 binary diagnostics]]
|Jenny Fisher (U. Wollongong)
+
|Lizzie Lundgren ([[GCST]])
|Benchmarking
+
|Bug fix
 +
|" "
 +
 
 +
|-valign="top"
 +
|[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02e|Fix bug in HEMCO reference time code]]
 +
|Barron Henderson (US EPA)
 +
|Bug fix
 
|" "
 
|" "
  
 
|}
 
|}
  
==== v11-02d ====
+
=== v11-02d ===
  
 
Please see the following link for complete information about the validation of GEOS-Chem v11-02d:
 
Please see the following link for complete information about the validation of GEOS-Chem v11-02d:
 
#[[GEOS-Chem_v11-02_benchmark_history#v11-02d|''Approval form for 1-month benchmark simulation v11-02d'']]
 
#[[GEOS-Chem_v11-02_benchmark_history#v11-02d|''Approval form for 1-month benchmark simulation v11-02d'']]
#[[GEOS-Chem_v11-02_benchmark_history#v11-02d-Run0|''Results for 1-year benchmark simulation v11-02d-Run0'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02d-Run2|''Results for 1-year benchmark simulation v11-02d-Run2'']]
  
 
{| border=1 cellspacing=0 cellpadding=5  
 
{| border=1 cellspacing=0 cellpadding=5  
Line 367: Line 434:
 
!width="200px"|Submitted by
 
!width="200px"|Submitted by
 
!width="100px"|Type
 
!width="100px"|Type
!width="150px"|Status
+
!width="300px"|Status
  
 
|-
 
|-
Line 376: Line 443:
 
|Tom&aacute;s Sherwen (York)<br>Johan Schmidt (Harvard)<br>Seb Eastham (MIT)<br>Lei Zhu (Harvard)<br>[[Chemistry Issues|Oxidants and Chemistry WG]]
 
|Tom&aacute;s Sherwen (York)<br>Johan Schmidt (Harvard)<br>Seb Eastham (MIT)<br>Lei Zhu (Harvard)<br>[[Chemistry Issues|Oxidants and Chemistry WG]]
 
|Science
 
|Science
|Pending approval
+
|Approved 12 Feb 2018
  
 
|-valign="top"
 
|-valign="top"
Line 424: Line 491:
  
 
|-valign="top"
 
|-valign="top"
|[[List of diagnostics for v11-02|Structural updates for netCDF diagnostics]]
+
|[[Guide to GEOS-Chem History diagnostics|Structural updates for netCDF diagnostics]]
 
|[[GCST]]
 
|[[GCST]]
 
|Structural
 
|Structural
Line 431: Line 498:
 
|}
 
|}
  
==== v11-02e ====
+
=== v11-02c ===
  
Please see the following link for complete information about the validation of GEOS-Chem v11-02e:
+
Please see the following link for complete information about the validation of GEOS-Chem v11-02c:
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e|''Approval form for 1-month benchmark simulation v11-02e'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02c|''Approval form for 1-month benchmark simulation v11-02c'']]
#[[GEOS-Chem_v11-02_benchmark_history#v11-02e-Run0|''Results for 1-year benchmark simulation v11-02e-Run0'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02c-Run0|''Results for 1-year benchmark simulation v11-02c-Run0'']]
  
 
{| border=1 cellspacing=0 cellpadding=5  
 
{| border=1 cellspacing=0 cellpadding=5  
Line 442: Line 509:
 
!width="200px"|Submitted by
 
!width="200px"|Submitted by
 
!width="100px"|Type
 
!width="100px"|Type
!width="150px"|Status
+
!width="300px"|Status
  
 
|-
 
|-
Line 448: Line 515:
  
 
|-valign="top"
 
|-valign="top"
|[[Downloading_GEOS-Chem_source_code_and_data#Inconsistency_in_GEOS-FP_files_at_Harvard_and_Dalhousie_for_July_2013|Sync GEOS-FP files on Harvard ftp with files at Dalhousie]]
+
|[[Secondary_organic_aerosols#Simple_SOA_scheme|Enhance default GEOS-Chem simple SOA]]
|[[GCST]]
+
|Sal Farina (Colorado State)<br>Aerosols Working Group Chairs
|Bug fix
+
|
+
 
+
|-valign="top"
+
|[[NOx-Ox-HC-aerosol#Monthly mean surface distributions|Monthly mean surface methane distributions]]
+
|Lee Murray (NASA GISS/LDEO)
+
 
|Science
 
|Science
|
+
|''1-month benchmark:''<br> Approved on 07 Sep 2017
*Delivered to [[GEOS-Chem Support Team|GCST]] (Nov 2017)
+
 
 +
''1-year benchmark:''<br> Approved 21 Sep 2017
  
 
|-valign="top"
 
|-valign="top"
|UCX updates, including:
+
|[[GEOS-Chem_chemistry_mechanisms#Updated_isoprene_and_monoterpene_chemistry|Updates to isoprene and monoterpene chemistry]]
*[[UCX_chemistry_mechanism#v11-02|Remove initial stratospheric 2D mixing ratio option]]
+
|Katie Travis (MIT)<br>Jenny Fisher (U. Wollongong)<br>Elo&iuml;se Marais (U. Birmingham)<br>Christopher Chan Miller (Harvard)<br>Kelvin Bates (Caltech)<br>Rebecca Schwantes(Caltech)
*[[UCX_chemistry_mechanism#UCX_stratospheric_water_boundary_condition_update|UCX stratospheric water boundary condition update]]
+
*Remove <tt>UCX=y</tt> as compile switch and use <tt>Input_Opt%LUCX</tt> instead
+
|Seb Eastham (Harvard)<br>Chris Holmes (Florida State)<br>[[GCST]]
+
 
|Science
 
|Science
|
+
|" "
*Delivered to GCST (Jan 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[[Stratospheric_chemistry#Prod.2Floss_rates_from_UCX|Use monthly mean stratospheric prod/loss rates from UCX simulation for all species]]
+
|[[Secondary_organic_aerosols#SOA_formation_from_aqueous_isoprene_uptake|Add aqueous isoprene uptake to SOA scheme]]
|[[GCST]]
+
|Eloise Marais (U. Birmingham)
 
|Science
 
|Science
|
+
|" "
*Delivered to GCST (Nov 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[[Sulfate_aerosols#Metal_catalyzed_oxidation_of_SO2|Metal catalyzed oxidation of SO2]] (as a switch)
+
|[[Chemistry_Issues#Fixes_for_carbon_creating_reactions|Carbon balance (fix C creation)]]
|Viral Shah (UW)<br>Becky Alexander (UW)
+
|Sarah Safieddine (MIT)
|Science
+
|Bug fix
|
+
|" "
*Delivered to GCST (Dec 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[[Aerosol_optical_properties#Add_BC_absorption_enhancement_factor|Update density of BC to 1.8 and add absorption enhancement factor in input.geos]]
+
|[[PAN#Bug_fixes_for_EOH_and_MGLY|Fix bugs for EOH and MGLY following implementation of PAN updates in v11-02a]]
|Xuan Wang (MIT)
+
|Melissa Sulprizio ([[GCST]])
|Science
+
|Bug fix
|
+
|" "
*Delivered to GCST (Jan 2016)
+
  
 
|-valign="top"
 
|-valign="top"
|[[Particulate_matter_in_GEOS-Chem#Option_to_include_spatially_and_seasonally_varying_OM.2FOC|Spatially varying OM/OC]]
+
|[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02c|Update HEMCO from v2.0.004 to v2.1.001]]
|Sajeev Philip (Dalhousie/NASA Ames)
+
|Christoph Keller (NASA GMAO)
|Science
+
|Structural
|
+
|" "
*Delivered to GCST (May 2017)
+
 
 +
|-
 +
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
  
 
|-valign="top"
 
|-valign="top"
|Add SOA updates and bug fixes, including:
+
|HEMCO updates:
*[[GEOS-Chem_v11-02#Remove_MONX_from_chemical_mechanisms|Remove MONX from chemical mechanisms]]
+
*[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02c|Add a HEMCO-standalone run directory for benchmarking purposes]]
*[[Secondary_organic_aerosols#Add_MTPO_as_an_advected_species_to_all_full-chemistry_simulations|Add MTPO as an advected species in full-chemistry simulations]]
+
*[[Implementation_of_HEMCO_in_GEOS-Chem#Now_use_YYYYMMDDhhmm_format_for_time_stamp_values|Now use <tt>YYYYMMDDhhmm</tt> for time stamp values]]
*[[Secondary_organic_aerosols#Update_simple_SOA_entries_in_HEMCO_to_follow_MTPA.2BLIMO.2BMTPO|Updating SOAP and SOAS entries in HEMCO to follow MTPA+LIMO+MTPO rather than MONX]]
+
*[[Implementation_of_HEMCO_in_GEOS-Chem#Avoid_segmentation_fault_in_DustGinoux_extension|Add error trap in to avoid a segmentation fault when the DustGinoux extension is turned off]]
*[[Secondary_organic_aerosols#Calculation_of_PM2.5.2C_AOD.2C_and_aerosol_mass|Use simple SOA in calculation of PM2.5, AOD, and aerosol mass by default]]
+
|<br>Christoph Keller (NASA GMAO)<br>Andy Jacobson (NOAA)<br>Paulo Tuccella (L'Aquila)
*[[GEOS-Chem_v11-02#Fixes_for_isoprene_SOA_for_consistency_with_Marais_et_al._.282016.29|Fixes for isoprene SOA for consistency with Marais et al. (2016)]]
+
|Structural
|Sam Silva (MIT)<br>[[Chemistry_Issues|Chemistry Working Group]]<br>[[Aerosols Working Group]]<br>[[GCST]]<br>Eloise Marais (U. Birmingham)
+
|" "
|Science & Bug fix
+
|
+
*Delivered to GCST (Nov 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_v11-02#Bug_fixes_for_aerosols_in_v11-01|Add bug fixes for aerosols]]
+
|Fixes for several minor issues:
|Jingyuan Shao (UW)
+
*[[#Bug fixes for the ND50 timeseries diagnostic|Bug fixes for the ND50 timeseries diagnostic]]
 +
*[[#Fixed_an_incorrect_format_statement_in_input_mod.F|Fixed an incorrect format statement in <tt>input_mod.F</tt>]]
 +
*[[#Add MERRA2 to #if and #elif statements where it had been omitted|Add MERRA2 to <tt>#if</tt> and <tt>#elif</tt> statements where it had been omitted]]
 +
*[[Mercury#Bug_fixes_for_Hg_emissions|Bug fixes for Hg emissions]]
 +
*[[Aerosol_emissions#Bug_fix_for_MOPO_and_MOPI_definitions_in_species_database|Fix bug in species definitions for marine POA simulation]]
 +
*[[Secondary_organic_aerosols#Fix_diagnostics_bugs_in_the_SOA-SVPOA_simulation|Fix diagnostic bugs in the SOA-SVPOA simulation]]
 +
|<br>Chris Holmes (Florida State)<br>Chris Holmes (Florida State)<br>Jiawei Zhuang (Harvard)<br>Amanda Giang (MIT)<br>Katie Travis (MIT)<br>Katie Travis(MIT)
 
|Bug fix
 
|Bug fix
|
+
|" "
*Delivered to GCST (Nov 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_v11-02#Fix_dust_rainout_settings_to_match_BCPO|Fix dust rainout settings to match BCPO]]
+
|Fixes for the TOMAS simulation:
|Colette Heald (MIT)
+
*[[TOMAS_aerosol_microphysics#Fixes_for_TOMAS_simulation_in_v11-02c|Fixes for TOMAS simulation in v11-02c]]
 +
*[[TOMAS_aerosol_microphysics#Typo_in_wetscav_mod.F_for_TOMAS30_simulation|Fix typo in wetscav_mod.F for TOMAS30]]
 +
*Remove a couple of array temporaries and an out-of-bounds error
 +
|<br>Sal Farina (Colorado State)<br>Jack Kodros (Colorado State)<br>Bob Yantosca ([[GCST]])
 
|Bug fix
 
|Bug fix
|
+
|" "
*Delivered to GCST (Aug 2017)
+
  
|-
+
|-valign="top"
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
+
|[[Stratospheric_chemistry#STE_fluxes|Fix STE flux diagnostics and add to benchmark procedure]]
 +
|Melissa Sulprizio ([[GCST]])
 +
|Bug fix &<br>Benchmarking
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Chemistry_Issues#Correcting_ozone_from_the_height_of_the_lowest_model_level_to_10m|Add ND49 timeseries diagnostic for correcting ozone from the height of the lowest model level to 10m]]
+
|Initial structural modifications for netCDF diagnostics:
 +
*Introduce <tt>Headers/State_Diag</tt> as a stub module (for now)
 +
*Add a registry object into <tt>State_Met</tt>, <tt>State_Diag</tt>, and <tt>State_Chm</tt> in order to obtain a pointer to any module variable (or slice) by looking up its name
 +
*Add new module <tt>Headers/registry_mod.F90</tt> which contains derived types and routines for registering module variables.
 +
|Bob Yantosca ([[GCST]])
 
|Structural
 
|Structural
|Katie Travis (MIT)
+
|" "
|
+
*Delivered to GCST (Nov 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|Change input.geos timestep units from minutes to seconds
+
|[[GEOS-Chem_v11-02#Code_cleanup_suggestions|Combine timestep settings in input.geos in a Timesteps menu]]
 +
|Melissa Sulprizio ([[GCST]])
 
|Structural
 
|Structural
|GCST
+
|" "
|
+
 
*Delivered to GCST (Nov 2017)
+
|-valign="top"
 +
|[[GEOS-Chem_benchmarking#Update_2005.2F2009_CMDL_CO_data_to_2013.2F2014_GMD_CO_data|Update CO data used in 1-year benchmark plots]]
 +
|Jenny Fisher (U. Wollongong)
 +
|Benchmarking
 +
|" "
  
 
|}
 
|}
  
==== v11-02f ====
+
=== v11-02b ===
  
Please see the following link for complete information about the validation of GEOS-Chem v11-02f:
+
Please see the following link for complete information about the validation of GEOS-Chem v11-02b:
#[[GEOS-Chem_v11-02_benchmark_history#v11-02f|''Approval form for 1-month benchmark simulation v11-02f'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b|''Approval form for 1-month benchmark simulation v11-02b'']]
#[[GEOS-Chem_v11-02_benchmark_history#v11-02f-Run0|''Results for 1-year benchmark simulation v11-02f-Run0'']]
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b_HP_1-month|''Approval form for 1-month benchmark simulation v11-02b with high performance option (v11-02b-HP)'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b-HP-Run0|''Results for 1-year standard benchmark simulation of v11-02b-HP Run0'']]
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b-HP-Run1|''Results for 1-year standard benchmark simulation of v11-02b-HP Run1'']] (uses offline archived lightning and dust emissions)
 +
#[[GEOS-Chem_v11-02_benchmark_history#v11-02b_HP_1-year_Rn-Pb-Be|''Results for 1-year RnPbBe benchmark simulation of v11-02b-HP'']]
  
 
{| border=1 cellspacing=0 cellpadding=5  
 
{| border=1 cellspacing=0 cellpadding=5  
Line 554: Line 628:
 
!width="200px"|Submitted by
 
!width="200px"|Submitted by
 
!width="100px"|Type
 
!width="100px"|Type
!width="150px"|Status
+
!width="300px"|Status
  
 
|-
 
|-
!colspan="4" bgcolor="#CCFFFF"|Features affecting the full-chemistry simulation:
+
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
  
 
|-valign="top"
 
|-valign="top"
|[[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions#Global_0.25x0.25_MEGAN_input_and_MODIS_LAI|Use global 0.25x0.25 files for MEGAN input and MODIS LAI]]
+
|[[GEOS-Chem HP|Source code updates for high performance option (GCHP)]]
|Jenny Fisher (U. Wollongong)<br>Eloise Marais (U. Birmingham)<br>Dylan Millet (UMN)
+
|Seb Eastham (Harvard)<br>Lizzie Lundgren (Harvard)<br>Mike Long (Harvard)<br>Jiawei Zhuang (Harvard)<br>Bob Yantosca (Harvard)
|Bug fix/<br>Science
+
|
+
*Delivered to GCST (Nov 2017)
+
 
+
|-valign="top"
+
|[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02|Update to HEMCO v2.1.003]] (requires updated MODIS LAI data)
+
|Christoph Keller (NASA GMAO)
+
 
|Structural
 
|Structural
|
+
|Approved 16 Jun 2017
*Delivered to GCST (Jul 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[http://onlinelibrary.wiley.com/doi/10.1002/2015JD023134/abstract OMI-based Volcanic emissions 2005-2012]
+
|Bug fixes for diagnostics:
|Cui Ge (U. Iowa)<br>Jun Wang (U. Iowa)
+
*[[GEOS-Chem_v11-02#ND26_vertical_flux_diagnostic_does_not_yield_mass_balance_in_advection|Use FZ array to compute ND26 vertical flux diagnostic]]
|Science
+
*[[GEOS-Chem_v11-02#ND40_planeflight_diagnostic_not_always_written_out_when_turned_on|Always write ND40 planeflight diagnostic when enabled]]
|
+
|<br>Ilya Stanevic (Toronto)<br>Lee Murray (Rochester)
*Delivered to GCST (Nov 2017)
+
|Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[EDGAR_v4.3_anthropogenic_emissions|EDGAR v4.3.1 emissions]]
+
|[[Tagged_O3_simulation#Write_initial_and_final_Ox_mass_to_file|Write initial and final Ox mass to file when using the tagged O3 simulation]]
|Chi Li (Dalhousie)
+
|Bob Yantosca ([[GCST]])
|Science
+
|
+
*Delivered to GCST (Nov 2016)
+
 
+
|-valign="top"
+
|EPA-derived BC/OA emissions over US (1990-2012)
+
|David Ridley (MIT)<br>Chi Li (Dalhousie)
+
|Science
+
|
+
*Need to coordinate between Ridley and Li versions
+
 
+
|-valign="top"
+
|[[CAC_anthropogenic_emissions#Historical_emissions|Historical CAC emissions]]
+
|Chi Li (Dalhousie)
+
|Science
+
|
+
*Delivered to GCST (Nov 2016)
+
 
+
|-valign="top"
+
|[http://onlinelibrary.wiley.com/doi/10.1002/2016JD025767/full Global ethane emissions]
+
|Zitely Tzompa-Sosa (CSU)
+
|Science
+
|
+
*Delivered to GCST (May 2017)
+
 
+
|-valign="top"
+
|[https://www2.acom.ucar.edu/modeling/dice-africa Inefficient combustion emissions in Africa (DICE-Africa inventory)]
+
|Eloise Marais (U. Birmingham)<br>Christine Wiedinmyer (NCAR)
+
|Science
+
|
+
*Delivered to GCST (Nov 2017)
+
 
+
|-valign="top"
+
|Trash burning emissions
+
|Eloise Marais (U. Birmingham)<br>Christine Wiedinmyer (NCAR)
+
|Science
+
|
+
*Delivered to GCST (Nov 2017)
+
 
+
|-valign="top"
+
|[[GFED3_biomass_burning_emissions|Retire GFED3 biomass burning emissions]]
+
|[[GCST]]
+
 
|Structural
 
|Structural
|
+
|" "
*Ready to go
+
  
 
|-valign="top"
 
|-valign="top"
|MEIC emissions over China (2000-2015)
+
|[[GEOS-Chem_v11-02#CO2_emissions_unit_conversion_not_consistent_with_new_advection_pressure|Convert CO2 emissions units using dry pressure used in advection]]
|Lu Shen (Harvard)<br>Melissa Sulprizio (Harvard)
+
|Meemong Lee (JPL)
 
|Science
 
|Science
|
+
|" "
*Delivered to GCST (Nov 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[http://www.globalchange.umd.edu/ceds/ CEDS global anthropogenic emissions (1950-2014)]
+
|[[QFED_biomass_burning_emissions|Add QFED emissions for 2014-2016]]
|Lu Shen (Harvard)<br>Melissa Sulprizio (Harvard)
+
|Christoph Keller (NASA GMAO)
 
|Science
 
|Science
|
+
|" "
*Delivered to GCST (Nov 2017)
+
  
|-valign="top"
+
|}
|[[Mineral_dust_aerosols#Anthropogenic_PM2.5_dust_source_in_GEOS-Chem|Anthropogenic PM2.5 dust source]]
+
|Sajeev Philip (Dalhousie)
+
|Science
+
|
+
*Delivered to GCST (Mar 2017)
+
  
|-valign="top"
+
=== v11-02a ===
|[[FINNv1_biomass_burning_emissions|FINN biomass burning emissions for 2014-2016]]
+
|Barron Henderson (US EPA)
+
|Science
+
|
+
*Delivered to GCST (Nov 2017)
+
  
 
+
Please see the following link for complete information about the validation of GEOS-Chem v11-02a:
|-
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02a|''Approval form for 1-month benchmark simulation v11-02a'']]
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02a-Run0|''Results for 1-year benchmark simulation v11-02a-Run0'']]
 
+
#[[GEOS-Chem_v11-02_benchmark_history#v11-02a-Run0|''Results for 1-year benchmark simulation v11-02a-Run1'']] (final version)
|-valign="top"
+
|[[Tagged_CO_simulation#New_tagged_CO_simulation_in_v11-02_and_later|Updated Tagged CO simulation]]
+
|Jenny Fisher (Wollongong)
+
|Science
+
|
+
*Delivered to GCST (Nov 2017)
+
 
+
|-valign="top"
+
|[[CH4_simulation#Updated_CH4_simulation_in_v11-02_and_later|Updated CH4 simulation]]
+
|Bram Maasakkers (Harvard)<br>Melissa Sulprizio (Harvard)
+
|Science
+
|
+
*Delivered to GCST (Oct 2017)
+
 
+
|-valign="top"
+
|Updates to Hg chemistry from Horowitz et al. (2017)
+
|Hannah Horowitz (U. Washington)<br>Colin Thackray (Harvard)<br>Rebecca Stern (Harvard)
+
|Science
+
|
+
 
+
|}
+
 
+
==== Defer to v11-03 ====
+
  
 
{| border=1 cellspacing=0 cellpadding=5  
 
{| border=1 cellspacing=0 cellpadding=5  
Line 690: Line 679:
 
!width="200px"|Submitted by
 
!width="200px"|Submitted by
 
!width="100px"|Type
 
!width="100px"|Type
!width="150px"|Status
+
!width="300px"|Status
  
 
|-
 
|-
!colspan="4" bgcolor="#CCFFFF"|Transport updates
+
!colspan="4" bgcolor="#CCFFFF"|Features affecting the full-chemistry simulation:
  
 
|-valign="top"
 
|-valign="top"
|[https://www.geosci-model-dev-discuss.net/gmd-2017-125/ GEOS resolution & subgrid convection]
+
|[[Updates_in_JPL_Publication_15-10|Update chemistry rate constants based on ''JPL Publication 15-10'']]
|Karen Yu (Harvard)
+
|Barron Henderson (US EPA),<br>Mat Evans (U. York), &<br>[[Chemistry Issues|Oxidants and Chemistry WG]]
 
|Science
 
|Science
|
+
|Approved 12 May 2017
*Delivered to GCST (Oct 2017)
+
 
+
|-
+
!colspan="4" bgcolor="#CCFFFF"|Emissions updates
+
  
 
|-valign="top"
 
|-valign="top"
|Ammonia emissions from Arctic seabirds
+
|[[Chemistry_Issues#Fixes_to_correct_ALK4_lumping_issue|Fixes to correct ALK4 lumping issue]]
|Betty Croft (Dalhousie)
+
|Barron Henderson (US EPA)
 
|Science
 
|Science
|
+
|" "
*Delivered to GCST (Nov 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[[Sulfate_aerosols#Ocean_ammonia_emission_inventory|Ocean ammonia emission inventory]]
+
|[[PAN|PAN updates]], including:
|Fabien Paulot (NOAA/GFDL)
+
*Updates to emissions, deposition, species, chemistry
 +
*Ocean+biogenic ALD2 emissions from Millet al. [2010]
 +
|Emily Fischer (CSU)
 
|Science
 
|Science
|
+
|" "
*Data files delivered to GCST (Jan 2016)
+
*Needs to be implemented via HEMCO
+
  
 
|-valign="top"
 
|-valign="top"
|[http://www.atmos-chem-phys.net/16/2819/2016/ A snow NOx source from deep snowpack]
+
|[[EPA/NEI11_North_American_emissions#Monthly_mean_emissions|Monthly mean NEI2011 emissions]]
|Maria Zatko (UW)<br>Becky Alexander (UW)
+
|[[GCST]] &<br>Katie Travis (Harvard)
 
|Science
 
|Science
|
+
|" "
*Delivered to GCST (Nov 2017)
+
*Implemented in v9-01-01
+
*Needs to be added as a HEMCO extension
+
  
 
|-valign="top"
 
|-valign="top"
|[[MASAGE_NH3_inventory|MASAGE inventory in HEMCO as default]]
+
|Bug fixes in the GEOS-Chem sulfate module:
|Fabien Paulot (NOAA/GFDL)<br>Amos Tai (CUHK)
+
*[[Sulfate_aerosols#Fix_error_in_production_of_SO4s.2C_NITs_in_SEASALT_CHEM_routine|Fix error in production of SO4s and NITs in <tt>SEASALT_CHEM</tt>]]
|Science
+
*[[Sulfate_aerosols#Fix_bug_in_CHEM_NIT_routine|Fix bug in <tt>CHEM_NIT</tt>]]
|
+
*[[Sulfate_aerosols#Fix_for_sulfate_production_in_HET_DROP_CHEM|Fix for sulfate production in <tt>HET_DROP_CHEM</tt>]]
 +
*[[Sulfate_aerosols#Fix_bugs_in_sulfate_chemistry_routines|Fix bugs in sulfate chemistry routines]]
 +
|<br>Prasad Kasibhatla (Duke)<br>Prasad Kasibhatla (Duke)<br>Qianjie Chen (UW)<br>Viral Shah (UW)
 +
|Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[MASAGE_NH3_inventory#Non-agricultural_NH3_for_use_with_MASAGE|Non-agricultural NH3 for use with MASAGE]]
+
|[[Dry_deposition#Bug_in_computation_of_aerodynamic_resistance_RA|Fix bug in dry deposition aerodynamic resistance]]
|Amos Tai (CUHK)
+
|Brian Boys (Dalhousie)
|Science
+
|Bug fix
|
+
|" "
  
 
|-valign="top"
 
|-valign="top"
|Radon emission update
+
|[[Acetone#Fix_acetone_parameterization_in_hcox_seaflux_mod.F90|Fix acetone parameterization in <tt>hcox_seaflux_mod.F90</tt>]]
|Bo Zhang (NIA/NASA Langley)
+
|[[GCST]]
|Science
+
|Bug fix
|
+
|" "
*Delivered to GCST (Nov 2017)
+
*Implementation pending paper submission
+
  
 
|-valign="top"
 
|-valign="top"
|EPA NEI2011 Version EK
+
|[[GEOS-Chem_v11-01#Bug_in_cos.28SZA.29_for_start_of_timestep|Bug fix in cos(SZA) for start of timestep]]
|Zitely Tzompa-Sosa (CSU)<br>Emily Fischer (CSU)
+
|Lizzie Lundgren ([[GCST]])
|Science
+
|Bug fix
|
+
|" "
*Pending paper submission
+
  
 
|-valign="top"
 
|-valign="top"
|GFAS biomass burning emissions
+
|[[GEOS-Chem_v11-02#Include_TOA_pressure_when_calculating_dry_pressure_edges|Include TOA pressure when calculating dry pressure edges]]
|Mat Evans (York)
+
|Seb Eastham (Harvard)
 
|Science
 
|Science
|
+
|" "
  
 
|-
 
|-
!colspan="4" bgcolor="#CCFFFF"|Deposition/surface updates
+
!colspan="4" bgcolor="#CCFFFF"|Features '''not''' affecting the full-chemistry simulation:
  
 
|-valign="top"
 
|-valign="top"
|Wet scavenging in mixed-phase clouds
+
|[[ISORROPIA_II#Implement_ISORROPIA_v2.0_as_a_Fortran_module|Implement ISORROPIA v2.0 as a Fortran module]]
|Bo Zhang (NIA/NASA Langley)
+
|Seb Eastham (Harvard) &<br>[[GCST]]
|Science
+
|Structural
|
+
|" "
*Delivered to GCST (Nov 2017)
+
*Only applied/tested for Pb
+
*Implementation pending paper submission
+
  
 
|-valign="top"
 
|-valign="top"
|Harmonizing Henry's Law coefficients across wet and dry deposition
+
|Updates to the HEMCO emissions component:
|[[GCST]]
+
*[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-02a|Update to HEMCO v2.0.004]]
|Science
+
*[[Implementation_of_HEMCO_in_GEOS-Chem#Fixed_bug_in_computation_of_local_time_in_HCO_GetSuncos|Fixed bug in computation of local time in <tt>HCO_GetSuncos</tt>]]
|
+
*[[EPA/NEI11_North_American_emissions#Default_US_emissions_to_NEI2011_after_2013|Default US emissions to NEI2011 after 2013]]
 +
*[[Implementation_of_HEMCO_in_GEOS-Chem#Read_default_DEP_RESERVOIR_fields_from_file_when_not_found_in_HEMCO_restart_file|Read default <tt>DEP_RESERVOIR</tt> fields from file]]
 +
|<br>Christoph Keller (NASA GMAO)<br>Seb Eastham (Harvard)<br>Jessica Morena (Dalhousie)<br>Brian Boys (Dalhousie) & [[GCST]]
 +
|Structural<br>& Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|[[Aerosols_Working_Group#Sedimentation_of_stratospheric_aerosols|Sedimentation of stratospheric aerosols]]
+
|Updates to gain computational speedup:
|Sebastian Eastham (Harvard)
+
*[[FlexChem#In_GEOS-Chem_v11-02_and_later_versions|Add simplified prod/loss families to KPP]]
|Science
+
*[[Cloud_convection#Computational_bottleneck_in_the_v11-01_convection_module|Remove computational bottleneck in the convection module]]
|
+
|<br>Mike Long (Harvard)<br>Bob Yantosca ([[GCST]])
 +
|Structural<br>& Bug fix
 +
|" "
  
|-
+
|-valign="top"
!colspan="4" bgcolor="#CCFFFF"|Aerosol updates
+
|[[Mercury#Fixed_typo_in_nested_grid_transport_routine_INIT_WINDOW|Fixed typo in <tt>INIT_WINDOW</tt>]]
 +
|Bob Yantosca ([[GCST]])
 +
|Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|Sea salt aerosol alkalinity
+
|netCDF file I/O updates:
|Qianjie Chen (UW)<br>Becky Alexander (UW)<br>Xuan Wang (Harvard)
+
*[[The NcdfUtilities package#Enable_compression_in_netCDF-4_output_files|Enable data compression in netCDF-4 output files]]
|Science
+
*[[The NcdfUtilities package#Routine DO_ERR_OUT now returns a non-zero error code|Routine <tt>DO_ERR_OUT</tt> now returns a non-zero error code]]
|
+
*[[Implementation of HEMCO in GEOS-Chem#HEMCO diagnostic and restart files now have an unlimited time dimension|HEMCO files now have an unlimited time dimension]]
*UW version delivered to GCST (Nov 2017)
+
|<br>Chris Holmes (Florida State)<br>Andy Jacobson (NOAA)<br>[[GCST]]
 +
|Bug fix &<br>Structural
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|VBS SOA option
+
|Makefile and build sequence updates:
|Rokjin Park (SNU)
+
*[[GEOS-Chem_Makefile_Structure#Bug_fix:_Specifying_NO_REDUCED.3Dno_now_compiles_GEOS-Chem_for_reduced_grids|Specifying <tt>NO_REDUCED=no</tt> now compiles GEOS-Chem for reduced grids]]
|Science
+
*[[GEOS-Chem_Makefile_Structure#Removed_the_COMPILER_variable_from_Makefile_header.mk_for_a_cleaner_build_sequence|Removed the <tt>COMPILER</tt> variable from <tt>Makefile_header.mk</tt>]]
|
+
|<br>Jiawei Zhang (Harvard)<br>[[GCST]]
 +
|Bug fix &<br>Structural
 +
|" "
  
|-
+
|-valign="top"
!colspan="4" bgcolor="#CCFFFF"|Specialty simulation updates
+
|[[UCX_chemistry_mechanism#Bug_fixes_for_running_UCX_in_ESMF_environment|Bug fixes for running UCX in ESMF environment]]
 +
|Christoph Keller (NASA GMAO)
 +
|Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|Oxidant fields from v10 for 2006-2012
+
|Bug fixes for diagnostics:
*3D chemical source for CO2, CO
+
*[[Particulate_matter_in_GEOS-Chem#Save_out_PM2.5_diagnostic_at_STP_conditions|Save out PM2.5 diagnostic at STP conditions]]
*OH for tagged CO & CH4
+
*[[GEOS-Chem_v11-01#Bugs_in_convective_flux_diagnostics_.28ND14.29|Bug fixes in convective mass flux diagnostic]]
|Dylan Jones (UofT)
+
*[[GEOS-Chem_v11-01#P.28OH.29_not_saved_out_in_ND22_diagnostic|Restore P(OH) in ND22 diagnostic]]
|Science
+
*[[Implementation_of_HEMCO_in_GEOS-Chem#Make_anthropogenic_emissions_diagnostics_3D|Make anthropogenic emissions diagnostics 3D]]
|
+
*[[Tagged_CO_simulation#Fix_ND65_bugs_in_tagged_CO_simulation|Fix ND65 bugs in tagged CO simulation]]
*Will make data files available on ftp site
+
*[[Mineral_dust_aerosols#Bug_in_ND21_diagnostic_indexing_for_size-resolved_dust_species|Fix bug in ND21 diagnostic indexing for dust species]]
 +
|<br>Aaron van Donkelaar (Dalhousie)<br>Jenny Fisher (U. Wollongong)<br>[[GCST]]<br>Jenny Fisher (U. Wollongong)<br>Chris Holmes (Florida State)
 +
|Bug fix
 +
|" "
  
 
|-valign="top"
 
|-valign="top"
|Year 2010 ocean Hg boundary conditions from ocean model
+
|Removal of obsolete variables:
|Hannah Horowitz (Harvard)
+
*[[#Removed_the_NNPAR_parameter_from_CMN_SIZE_mod.F|Removed the <tt>NNPAR</tt> parameter from <tt>CMN_SIZE_mod.F</tt>]]
|Science
+
*[[Wet_deposition#Removal_of_obsolete_variables_NSOL_and_IDWETD|Removed obsolete variables <tt>NSOL</tt> and <tt>IDWETD</tt>]]
|
+
*[[Derived type objects used by GEOS-Chem#Summary of fields removed from Input_Opt|Removed obsolete fields of the <tt>Input_Opt</tt> object]]
 +
|[[GCST]]
 +
|Structural
 +
|" "
 +
 
 +
|}
 +
 
 +
== New data directories ==
 +
 
 +
The following [[HEMCO data directories]] have been added or updated for GEOS-Chem v11-02.  You will have to download the directories relevant to your simulation.
 +
 
 +
{| border=1 cellspacing=0 cellpadding=5
 +
|-bgcolor="#CCCCCC"
 +
!width="200px"|Inventory
 +
!width="400px"|What was added?
 +
!width="100px"|When added?
 +
!width="380px"|Directory
  
 
|-valign="top"
 
|-valign="top"
|Historical products emissions from Zhang et al. (2016), should become default
+
|[[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions|MEGAN]]
|Colin Thackray (Harvard)
+
|[[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions#Global_0.25x0.25_MEGAN_input_and_MODIS_LAI|Global 0.25x0.3125 CLM4 PFT and emission factor input files for MEGAN]]
|Science
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/MEGAN/v2017-07/</tt>
  
 
|-valign="top"
 
|-valign="top"
|AMAP/UNEP 2035 for Hg simulation
+
|[[MODIS_leaf_area_indices|MODIS LAI]]
|Helene Angot (MIT)
+
|[[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions#Global_0.25x0.25_MEGAN_input_and_MODIS_LAI|Updated 0.25 x 0.25 MODIS LAI files (2005-2011)]]
|Science
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>CHEM_INPUTS/MODIS_LAI_201707/</tt>
  
 
|-valign="top"
 
|-valign="top"
|Africa nested-grid for GEOS-FP and MERRA-2
+
|[[Lightning_NOx_emissions|Lightning NOx]]
|Eloise Marais (U. Birmingham)
+
|GEOS-FP OTD-LIS redistribution factors for Apr 2012 - Jul 2017 and MERRA-2 OTD-LIS redistribution factors for any date
|Science
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/LIGHTNOX/v2017-09/</tt>
*Delivered to GCST (Dec 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|West Africa nested-grid for GEOS-FP
+
|[[Volcanic_SO2_emissions_from_Aerocom|Volcanic SO2 emissions]]
|Mat Evans (York)<br>Eleanor Morris (York)
+
|[[Volcanic_SO2_emissions_from_Aerocom#OMI-based_volcanic_emissions|OMI-based volcanic emissions for 2005-2012]]
|Science
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/VOLCANO/v2018-03/</tt>
  
 
|-valign="top"
 
|-valign="top"
|[[Sulfate_aerosols#Tagged_sulfate_and_nitrate_simulation|Tagged sulfate simulation]]
+
|APEI
|Becky Alexander (UW)
+
|[[CAC_anthropogenic_emissions#Historical_Canadian_emissions|Historical Canadian emissions]]
|Science
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/APEI/v2016-11/</tt>
  
 
|-valign="top"
 
|-valign="top"
|Arsenic simulation
+
|[[EDGAR_v4.3_anthropogenic_emissions|EDGAR v4.3]]
|Shiliang Wu (MT)
+
|Annual global anthropogenic emissions for 1970-2010 at 0.1x0.1 resolution
|Science
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/EDGARv43/v2016-11/</tt>
 
+
|-
+
!colspan="4" bgcolor="#CCFFFF"|Structural updates
+
  
 
|-valign="top"
 
|-valign="top"
|Updated flux diagnostics
+
|C2H6_2010
|Chris Holmes (Florida State)<br>[[GCST]]
+
|[[Anthropogenic_emissions#C2H6|Global fossil fuel and biofuel emissions of C2H6 for 2010]]
|Structural/<br>Benchmarking
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/C2H6_2010/v2017-05/</tt>
  
 
|-valign="top"
 
|-valign="top"
|Additional modifications for [[FlexChem]], including:
+
|[[DICE-Africa_anthropogenic_emissions_inventory|DICE-Africa]]
*Rebuild the pre-built chemistry mechanisms (Standard, Tropchem, UCX, SOA, SOA-SVPOA) with [http://www.paratools.com/Kppa Kppa]
+
|Biofuel and diffuse anthropogenic emissions for Africa
*Add capability for custom chemistry mechanism by building KPP on-the-fly in GEOS-Chem
+
|[[#v11-02f|v11-02f]]
|Mike Long (Harvard)<br>[[GCST]]
+
|<tt>HEMCO/DICE_Africa/v2016-10/</tt>
|Structural
+
|
+
  
 
|-valign="top"
 
|-valign="top"
|Benchmark of stratospheric simulation
+
|[[Trash burning emissions]]
|Dylan Jones (U. Toronto)
+
|Global trash burning emissions from Wiedinmyer et al. (2014)
|Benchmarking
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/TrashEmis/v2015-03/</tt>
*1-yr benchmarks for GEOS-Chem v11-01d and later include comparisons with OSIRIS O3
+
*Waiting on OSIRIS NO2 data
+
*Waiting on comparisons with ACE-FTS data
+
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_species_units|Uniformity of units]]
+
|SEABIRD NH3
|[[GCST]]
+
|[[NH3_emissions#NH3_emissions_from_arctic_seabirds|NH3 emissions from arctic sea birds]]
|Structural
+
|[[#v11-02f|v11-02f]]
|
+
|<tt>HEMCO/NH3/v2018-04</tt>
*Phase 2 (unit uniformity throughout chemistry) to be implemented after Flexchem
+
*Phase 3 (unit uniformity in some specialty simulations) will be implemented in pieces with input from 3rd party developers
+
*Phase 4 (removal of area and area-dependent State_Met variables) will be ongoing
+
  
 
|-valign="top"
 
|-valign="top"
|ObsPack diagnostic as netCDF replacement for planeflight_mod.F
+
|OMOC
|Andy Jacobson (NOAA/ESRL)
+
|[[Particulate_matter_in_GEOS-Chem#Option_to_include_spatially_and_seasonally_varying_OM.2FOC|Spatially varying OM/OC ratio]]
|Structural
+
|[[#v11-02e|v11-02e]]
|
+
|<tt>HEMCO/OMOC/v2018-01</tt>
*Delivered to GCST (May 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|[[GEOS-Chem_Output_Files#NetCDF_File_Format|NetCDF Output Phase 3b: Updates to eliminate differences between single and multi-segmented GEOS-Chem runs]]  
+
|NOAA GMD CH4
|[[GCST]]<br>Christoph Keller
+
|[[NOx-Ox-HC-aerosol#Monthly_mean_surface_distributions|Monthly mean surface methane distributions]]
|Bug fix/Structural
+
|[[#v11-02e|v11-02e]]
|
+
|<tt>HEMCO/NOAA_GMD/v2018-01</tt>
 
+
|-
+
!colspan="4" bgcolor="#CCFFFF"|Expansion of model capabilities
+
  
 
|-valign="top"
 
|-valign="top"
|[[FlexGrid|FlexGrid: Use HEMCO to read met fields and enable nested simulations over any custom domains]]
+
|UCX
|Jiawei Zhuang (Harvard)<br>Jintai Lin (PKU)
+
|[[Stratospheric_chemistry#Prod.2Floss_rates_from_UCX|Input files for applying linearized stratospheric chemistry production and loss rates saved out from a UCX simulation using GEOS-Chem v11-02d]]
|Structural
+
|[[#v11-02e|v11-02e]]
|
+
|<tt>HEMCO/UCX/v2018-02/</tt>
*Delivered to GCST (Oct 2016)
+
  
 
|-valign="top"
 
|-valign="top"
|[[ISORROPIA_II#Investigating_persistent_noise_observed_in_ISORROPIA_output|Update ISORROPIA from v2.0 to v2.2, in order to remove persistent noise in output]]
+
|IODINE
|Seb Eastham (Harvard)
+
|[[Halogen_chemistry_mechanism#Source_code_and_data_files|Monthly emission fluxes of organic iodine compounds (CH3I, CH2I2, CH2ICl, CH2IBr) from Carlos Ordonez.]]
|Science
+
|[[#v11-02d|v11-02d]]
|
+
|<tt>HEMCO/IODINE/v2017-03/</tt>
*Delivered to GCST (May 2017)
+
  
 
|-valign="top"
 
|-valign="top"
|Archival of roughness length for temperature z0T
+
|GMI
|[[GCST]]
+
|New symbolic links:
|Met processing
+
*<tt>gmi.clim.RIPA.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25</tt>
|
+
*<tt>gmi.clim.RIPB.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25</tt>
 +
*<tt>gmi.clim.RIPD.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25</tt>
 +
*<tt>gmi.clim.IPMN.geos5.2x25.nc -> gmi.clim.PMN.geos5.2x25</tt>
 +
*<tt>gmi.clim.NPMN.geos5.2x25.nc -> gmi.clim.PMN.geos5.2x25</tt>
 +
|[[#v11-02c|v11-02c]]
 +
|<tt>HEMCO/GMI/v2015-02/</tt>
  
|}
+
|-valign="top"
 
+
|[[QFED_biomass_burning_emissions|QFED]]
=== New data directories ===
+
|Data for 2014-2016
 
+
|[[#v11-02b|v11-02b]]
The following [[HEMCO data directories]] have been added or updated for GEOS-Chem v11-02.  You will have to download the directories relevant to your simulation.
+
|<tt>HEMCO/QFED/v2014-09/</tt>
 
+
{| border=1 cellspacing=0 cellpadding=5
+
|-bgcolor="#CCCCCC"
+
!width="200px"|Inventory
+
!width="400px"|What was added?
+
!width="100px"|When added?
+
!width="380px"|Directory
+
  
 
|-valign="top"
 
|-valign="top"
Line 967: Line 961:
 
|[[#v11-02a|v11-02a]]
 
|[[#v11-02a|v11-02a]]
 
|<tt>HEMCO/SOILNOX/v2014-07/</tt>
 
|<tt>HEMCO/SOILNOX/v2014-07/</tt>
 
|-valign="top"
 
|[[QFED_biomass_burning_emissions|QFED]]
 
|Data for 2014-2016
 
|[[#v11-02b|v11-02b]]
 
|<tt>HEMCO/QFED/v2014-09/</tt>
 
 
|-valign="top"
 
|GMI
 
|New symbolic links:
 
*<tt>gmi.clim.RIPA.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25</tt>
 
*<tt>gmi.clim.RIPB.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25</tt>
 
*<tt>gmi.clim.RIPD.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25</tt>
 
*<tt>gmi.clim.IPMN.geos5.2x25.nc -> gmi.clim.PMN.geos5.2x25</tt>
 
*<tt>gmi.clim.NPMN.geos5.2x25.nc -> gmi.clim.PMN.geos5.2x25</tt>
 
|[[#v11-02c|v11-02c]]
 
|<tt>HEMCO/GMI/v2015-02/</tt>
 
 
|-valign="top"
 
|IODINE
 
|[[Halogen_chemistry_mechanism#Source_code_and_data_files|Monthly emission fluxes of organic iodine compounds (CH3I, CH2I2, CH2ICl, CH2IBr) from Carlos Ordonez.]]
 
|[[#v11-02d|v11-02d]]
 
|<tt>HEMCO/IODINE/v2017-03/</tt>
 
 
|-valign="top"
 
|UCX
 
|[[Stratospheric_chemistry#Prod.2Floss_rates_from_UCX|Input files for applying linearized stratospheric chemistry production and loss rates saved out from a UCX simulation using GEOS-Chem v11-02d.]]
 
|[[#v11-02d|v11-02d]]
 
|<tt>HEMCO/UCX/v2017-12/</tt>
 
  
 
|}
 
|}
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:31, 18 January 2018 (UTC)
 
  
 
== Validation ==
 
== Validation ==
  
In this section we provide information about the benchmarks and tests that we have done to validate GEOS-Chem v11-01.
+
In this section we provide information about the benchmarks and tests that we have done to validate GEOS-Chem v11-02.
  
 
=== 1-month and 1-year benchmarks ===
 
=== 1-month and 1-year benchmarks ===
Line 1,027: Line 990:
 
|[http://ftp.as.harvard.edu/gcgrid/geos-chem/unit_tests/v11-02/v11-02b.results.html v11-02b]
 
|[http://ftp.as.harvard.edu/gcgrid/geos-chem/unit_tests/v11-02/v11-02b.results.html v11-02b]
 
|18 May 2017
 
|18 May 2017
 +
 +
|-valign="top"
 +
|[http://ftp.as.harvard.edu/gcgrid/geos-chem/unit_tests/v11-02/v11-02c.results.html v11-02c]
 +
|5 Sep 2017
 +
 +
|-valign="top"
 +
|[http://ftp.as.harvard.edu/gcgrid/geos-chem/unit_tests/v11-02/v11-02d.results.html v11-02d]
 +
|5 Jan 2018
 +
 +
|-valign="top"
 +
|[http://ftp.as.harvard.edu/gcgrid/geos-chem/unit_tests/v11-02/v11-02e.results.html v11-02e]
 +
|07 Mar 2018
 +
 +
|-valign="top"
 +
|[http://ftp.as.harvard.edu/gcgrid/geos-chem/unit_tests/v11-02/v11-02f.results.html v11-02f]
 +
|02 May 2018
  
 
|}
 
|}
Line 1,032: Line 1,011:
 
== Previous issues now resolved in GEOS-Chem v11-02 ==
 
== Previous issues now resolved in GEOS-Chem v11-02 ==
  
The following bugs and/or technical issues have been resolved in v11-02:
+
The following bugs and/or technical issues have been resolved in v11-02.
  
 
=== Removed the NNPAR parameter from CMN_SIZE_mod.F ===
 
=== Removed the NNPAR parameter from CMN_SIZE_mod.F ===
Line 1,106: Line 1,085:
 
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02b|v11-02b]] and approved on 16 Jun 2017.'''''</span>
 
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02b|v11-02b]] and approved on 16 Jun 2017.'''''</span>
  
'''''Lee Murray (University of Rochester) wrote:'''''
+
[[Planeflight_diagnostic#ND40_planeflight_diagnostic_not_always_written_out_when_turned_on|We have moved the discussion of this issue]] to our ''Planeflight diagnostics'' wiki page.
 
+
<blockquote>I think the FlexChem implementation accidentally disabled the planeflight diagnostic for the full chemistry simulation. We used to initialize the planeflight diagnostic in chemdr.F for full chemistry simulations (now no more) and all the other versions in main.F. We should change main.F to always call <tt>SETUP_PLANEFLIGHT</tt> whenever the diagnostic is turned on (remove code in <span style="color:red">RED</span>, add code in <span style="color:green">GREEN</span>).</blockquote>
+
 
+
        ! Plane following diagnostic
+
        IF ( ND40 > 0 .and. DO_DIAG_WRITE ) THEN
+
       
+
            ! Call SETUP_PLANEFLIGHT routine if necessary
+
            IF ( ITS_A_NEW_DAY() ) THEN
+
             
+
<span style="color:red">!              ! If it's a full-chemistry simulation but LCHEM=F,
+
!              ! or if it's an offline simulation, call setup routine
+
!              IF ( ITS_A_FULLCHEM_SIM ) THEN
+
!                  IF ( .not. LCHEM ) THEN
+
!                    CALL SETUP_PLANEFLIGHT( am_I_Root, Input_Opt,
+
!    &                                      State_Chm, RC )
+
!                  ENDIF
+
!              ELSE
+
!                  CALL SETUP_PLANEFLIGHT( am_I_Root, Input_Opt,
+
!    &                                    State_Chm, RC )
+
!              ENDIF</span>
+
+
<span style="color:green">              CALL SETUP_PLANEFLIGHT( am_I_Root, Input_Opt, State_Chm, RC )</span>
+
+
            ENDIF
+
        ENDIF
+
 
+
--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 19:44, 10 May 2017 (UTC)
+
  
 
=== CO2 emissions unit conversion not consistent with new advection pressure ===
 
=== CO2 emissions unit conversion not consistent with new advection pressure ===
Line 1,274: Line 1,226:
 
<span style="color:green">'''''Suggestion 1 was included in [[GEOS-Chem v11-02#v11-02c|v11-02c]] and approved on 21 Sep 2017.'''''</span>
 
<span style="color:green">'''''Suggestion 1 was included in [[GEOS-Chem v11-02#v11-02c|v11-02c]] and approved on 21 Sep 2017.'''''</span>
  
<span style="color:darkorange">'''''Suggestion 2 is slated for inclusion later in [[GEOS-Chem v11-02]].'''''</span>
+
<span style="color:green">'''''Suggestion 2 (Git ID: [https://github.com/geoschem/geos-chem/commit/d06ae04ad44defd15854a9c4dadd36923dee0498 d06ae04a]) was included in [[GEOS-Chem 12#12.1.0|GEOS-Chem 12.1.0]], which was released on 26 Nov 2018.'''''</span>
  
 
'''''[[User:Chris Holmes|Chris Holmes]] wrote:'''''
 
'''''[[User:Chris Holmes|Chris Holmes]] wrote:'''''
Line 1,294: Line 1,246:
 
       ssa.dat
 
       ssa.dat
 
       ssc.dat
 
       ssc.dat
 +
 +
The .dat files required for FAST-JX can now be found in [http://ftp.as.harvard.edu/gcgrid/data/ExtData/CHEM_INPUTS/FAST_JX/v2018-09/ <tt>CHEM_INPUTS/FAST_JX/v2018-09/</tt>].
  
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:37, 19 July 2017 (UTC)
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:37, 19 July 2017 (UTC)
Line 1,299: Line 1,253:
 
=== Update CH4 latitude bands for 2014-2016 ===
 
=== Update CH4 latitude bands for 2014-2016 ===
  
<span style="color:darkorange">'''''This update will be included in [[GEOS-Chem v11-02#v11-02d|v11-02d]].'''''</span>
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02d|v11-02d]] (approved 12 Feb 2018).'''''</span>
  
 
The last year of latitudinal CH4 concentration data (based on CMDL flask observations) was 2013.  (Data is provided in 4 [[NOx-Ox-HC-Aer-Br_chemistry_mechanism#Latitude_bands|latitude bands]]: 90&deg;S-30&deg;S, 30&deg;S-0&deg;N:, 0&deg;N-30&deg;N, 30&deg;N-90&deg;N.) [[User:Katherine Travis|Katherine Travis]] has updated the subroutine <code>get_global_ch4.F</code> with additional years of latitudinal CH4 data from 2014 thru 2016.  
 
The last year of latitudinal CH4 concentration data (based on CMDL flask observations) was 2013.  (Data is provided in 4 [[NOx-Ox-HC-Aer-Br_chemistry_mechanism#Latitude_bands|latitude bands]]: 90&deg;S-30&deg;S, 30&deg;S-0&deg;N:, 0&deg;N-30&deg;N, 30&deg;N-90&deg;N.) [[User:Katherine Travis|Katherine Travis]] has updated the subroutine <code>get_global_ch4.F</code> with additional years of latitudinal CH4 data from 2014 thru 2016.  
Line 1,305: Line 1,259:
 
[http://www.ldeo.columbia.edu/~ltmurray/ Lee Murray] will provide additional latitudinal CH4 concentration data through present day, once it is available. See [[NOx-Ox-HC-aerosol#Latitude bands|this post on our ''NOx-Ox-HC-aerosol'' wiki page]] for more information.
 
[http://www.ldeo.columbia.edu/~ltmurray/ Lee Murray] will provide additional latitudinal CH4 concentration data through present day, once it is available. See [[NOx-Ox-HC-aerosol#Latitude bands|this post on our ''NOx-Ox-HC-aerosol'' wiki page]] for more information.
  
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 15:16, 5 September 2017 (UTC)
+
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 15:16, 5 September 2017 (UTC)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 21:59, 20 February 2018 (UTC)
  
 
=== Add bug fixes for GLYX, MGLY, IEPOX, and IMAE heterogeneous chemistry added in v11-02c ===
 
=== Add bug fixes for GLYX, MGLY, IEPOX, and IMAE heterogeneous chemistry added in v11-02c ===
  
<span style="color:darkorange">'''''This update will be included in [[GEOS-Chem v11-02#v11-02d|v11-02d]].'''''</span>
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02d|v11-02d]] (approved 12 Feb 2018).'''''</span>
  
 
'''''Sebastian Eastham wrote:'''''
 
'''''Sebastian Eastham wrote:'''''
Line 1,342: Line 1,296:
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 17:57, 20 September 2017 (UTC)
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 17:57, 20 September 2017 (UTC)
  
=== Add fix for IPMN reactions contributing to high P(Ox) and L(NOx) ===
+
=== Add fix for IPMN reactions contributing to high P(Ox) and L(Ox) ===
  
<span style="color:darkorange">'''''This update will be included in [[GEOS-Chem v11-02#v11-02d|v11-02d]].'''''</span>
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02d|v11-02d]] (approved 12 Feb 2018).'''''</span>
  
 
During the 1-year benchmark for v11-02d, Lu Hu pointed out an issue in the O3 budget where P(Ox) and L(Ox) were too high (~40000 Tg/yr compared to the expected 4000-5000 Tg/yr). This issue was tracked to the following reactions in the chemical mechanism:
 
During the 1-year benchmark for v11-02d, Lu Hu pointed out an issue in the O3 budget where P(Ox) and L(Ox) were too high (~40000 Tg/yr compared to the expected 4000-5000 Tg/yr). This issue was tracked to the following reactions in the chemical mechanism:
Line 1,355: Line 1,309:
 
:I'd suggest just using the same reaction rate values as for NPMN: non-isoprene PMN, as the reaction rates I used come from a rather speculative chemical mechanism anyway proposed by Lin et al., 2013 (http://www.pnas.org/content/suppl/2013/04/03/1221150110.DCSupplemental/pnas.201221150SI.pdf) and the isoprene SOA that these reactions form is relatively minor.
 
:I'd suggest just using the same reaction rate values as for NPMN: non-isoprene PMN, as the reaction rates I used come from a rather speculative chemical mechanism anyway proposed by Lin et al., 2013 (http://www.pnas.org/content/suppl/2013/04/03/1221150110.DCSupplemental/pnas.201221150SI.pdf) and the isoprene SOA that these reactions form is relatively minor.
  
This fix involves changing the reactions rates for the following IPMN reactions from:
+
This fix involves changing the reactions rates for the following IPMN reactions in the <tt>.eqn</tt> files from:
  
 
  IMAO3 + NO2 {+M} = IPMN : GCJPLPR(<span style="color:red">1.21E-11, 1.1E+00, 0.0, 7.7E-12, 0.2, 0.0, 0.6, 0.0, 0.0</span>);
 
  IMAO3 + NO2 {+M} = IPMN : GCJPLPR(<span style="color:red">1.21E-11, 1.1E+00, 0.0, 7.7E-12, 0.2, 0.0, 0.6, 0.0, 0.0</span>);
Line 1,369: Line 1,323:
 
  NPMN = NMAO3 + NO2 : GCJPLEQ(9.00E-29, 0.0E+00, 14000.0, 9.00E-28, 8.9E0, 0.0, 7.7E-12, 0.2, 0., 0.6, 0., 0.);
 
  NPMN = NMAO3 + NO2 : GCJPLEQ(9.00E-29, 0.0E+00, 14000.0, 9.00E-28, 8.9E0, 0.0, 7.7E-12, 0.2, 0., 0.6, 0., 0.);
  
and rebuilding the chemical mechanisms with KPP. In addition IPMN should be added to the Ox family in <tt>gckpp.kpp</tt>:
+
In addition IPMN should be added to the Ox family in <tt>gckpp.kpp</tt>:
  
 
  POx : O3 + NO2 + 2NO3 + PAN + NPMN + <span style="color:green">IPMN</span> + PPN + HNO4 + 3N2O5 + HNO3 + BrO + HOBr + BrNO2 + 2BrNO3 + MPN + ETHLN + ISN1 + ISOPNB + ISOPND + MACRN + MVKN + PROPNN + R4N2 + INPN + ISNP + INO2 + ISNOOA + ISNOOB + ISNOHOO + MAN2 + PRN1 + PRPN + R4N1 + PMNN + MACRNO2 + ClO + HOCl + ClNO2 + 2ClNO3 + 2Cl2O2 + 2OClO + O + O1D + IO + HOI + IONO + 2IONO2 + 2OIO + 2I2O2 + 3I2O3 + 4I2O4;
 
  POx : O3 + NO2 + 2NO3 + PAN + NPMN + <span style="color:green">IPMN</span> + PPN + HNO4 + 3N2O5 + HNO3 + BrO + HOBr + BrNO2 + 2BrNO3 + MPN + ETHLN + ISN1 + ISOPNB + ISOPND + MACRN + MVKN + PROPNN + R4N2 + INPN + ISNP + INO2 + ISNOOA + ISNOOB + ISNOHOO + MAN2 + PRN1 + PRPN + R4N1 + PMNN + MACRNO2 + ClO + HOCl + ClNO2 + 2ClNO3 + 2Cl2O2 + 2OClO + O + O1D + IO + HOI + IONO + 2IONO2 + 2OIO + 2I2O2 + 3I2O3 + 4I2O4;
 
  LOx : O3 + NO2 + 2NO3 + PAN + NPMN + <span style="color:green">IPMN</span> + PPN + HNO4 + 3N2O5 + HNO3 + BrO + HOBr + BrNO2 + 2BrNO3 + MPN + ETHLN + ISN1 + ISOPNB + ISOPND + MACRN + MVKN + PROPNN + R4N2 + INPN + ISNP + INO2 + ISNOOA + ISNOOB + ISNOHOO + MAN2 + PRN1 + PRPN + R4N1 + PMNN + MACRNO2 + ClO + HOCl + ClNO2 + 2ClNO3 + 2Cl2O2 + 2OClO + O + O1D + IO + HOI + IONO + 2IONO2 + 2OIO + 2I2O2 + 3I2O3 + 4I2O4;
 
  LOx : O3 + NO2 + 2NO3 + PAN + NPMN + <span style="color:green">IPMN</span> + PPN + HNO4 + 3N2O5 + HNO3 + BrO + HOBr + BrNO2 + 2BrNO3 + MPN + ETHLN + ISN1 + ISOPNB + ISOPND + MACRN + MVKN + PROPNN + R4N2 + INPN + ISNP + INO2 + ISNOOA + ISNOOB + ISNOHOO + MAN2 + PRN1 + PRPN + R4N1 + PMNN + MACRNO2 + ClO + HOCl + ClNO2 + 2ClNO3 + 2Cl2O2 + 2OClO + O + O1D + IO + HOI + IONO + 2IONO2 + 2OIO + 2I2O2 + 3I2O3 + 4I2O4;
 +
 +
After implementing those fixes, the chemical mechanisms will need to be [[FlexChem#Building_a_custom_chemical_mechanism|rebuilt with KPP]].
  
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:42, 25 January 2018 (UTC)
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:42, 25 January 2018 (UTC)
Line 1,378: Line 1,334:
 
=== Remove support for GCAP, GEOS-4, MERRA, GEOS-5  ===
 
=== Remove support for GCAP, GEOS-4, MERRA, GEOS-5  ===
  
<span style="color:darkorange">'''''This update will be included in [[GEOS-Chem v11-02#v11-02d|v11-02d]].'''''</span>
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02d|v11-02d]] (approved 12 Feb 2018).'''''</span>
  
 
Daniel Jacob wrote to the GEOS-Chem Steering Committee:
 
Daniel Jacob wrote to the GEOS-Chem Steering Committee:
Line 1,392: Line 1,348:
 
:I believe we currently have GCAP2 met fields stored at Harvard, and these should be maintained. GCAP2 consists of met fields from the GISS ModelE2.1 present-day and future. Lee Murray and Eric Leibensperger created these files, after a lot of work!
 
:I believe we currently have GCAP2 met fields stored at Harvard, and these should be maintained. GCAP2 consists of met fields from the GISS ModelE2.1 present-day and future. Lee Murray and Eric Leibensperger created these files, after a lot of work!
  
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:11, 2 November 2017 (UTC)
+
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:11, 2 November 2017 (UTC)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 21:59, 20 February 2018 (UTC)
  
 
=== Bug fixes for aerosols in v11-01 ===
 
=== Bug fixes for aerosols in v11-01 ===
  
<span style="color:darkorange">'''''These fixes will be included in [[GEOS-Chem v11-02#v11-02e|v11-02e]].'''''</span>
+
<span style="color:green">'''''This fix was included in [[GEOS-Chem v11-02#v11-02c|v11-02c]].'''''</span>
  
 
'''''Jingyuan Shao (UW) wrote:'''''
 
'''''Jingyuan Shao (UW) wrote:'''''
Line 1,416: Line 1,372:
 
       IF ( AS /= 0) CALL ALLOC_ERR( 'WAERSL' )
 
       IF ( AS /= 0) CALL ALLOC_ERR( 'WAERSL' )
 
       <span style="color:green">WAERSL = 0d0</span>
 
       <span style="color:green">WAERSL = 0d0</span>
 +
 +
The fix in aerosol_mod.F was included in [[GEOS-Chem v11-02#v11-02c|v11-02c]] with the implementation of [[Secondary_organic_aerosols#SOA_formation_from_aqueous_isoprene_uptake|aqueous isoprene uptake]] updates from Eloise Marais.
  
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:28, 15 November 2017 (UTC)
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:28, 15 November 2017 (UTC)
Line 1,421: Line 1,379:
 
=== Remove MONX from chemical mechanisms ===
 
=== Remove MONX from chemical mechanisms ===
  
<span style="color:darkorange">'''''This fix will be included in [[GEOS-Chem v11-02#v11-02e|v11-02e]].'''''</span>
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02e|v11-02e]] (approved 24 Mar 2018).'''''</span>
  
 
'''''Sam Silva (MIT) wrote:'''''
 
'''''Sam Silva (MIT) wrote:'''''
Line 1,442: Line 1,400:
 
:Looking back at Emily’s original PAN updates that she sent to us (based on v9-02h), it looks like MONX is a dead species in globchem.dat and the MONX -> 3CH2O reaction is not included. That may have been an oversight during the PAN updates implementation. Since Emily’s globchem.dat did not include the SOA chemistry by default and our standard/benchmark mechanism did, it wasn’t clear that we should remove the MONX species and reaction.  
 
:Looking back at Emily’s original PAN updates that she sent to us (based on v9-02h), it looks like MONX is a dead species in globchem.dat and the MONX -> 3CH2O reaction is not included. That may have been an oversight during the PAN updates implementation. Since Emily’s globchem.dat did not include the SOA chemistry by default and our standard/benchmark mechanism did, it wasn’t clear that we should remove the MONX species and reaction.  
  
The MONX species will be removed from GEOS-Chem mechanisms in v11-02e, upon the recommendation given to us by Emily Fischer.
+
The MONX species will be removed from GEOS-Chem mechanisms in v11-02e, upon the recommendation given to us by Emily Fischer. This update also requires the addition of the MTPO species to represents all the monoterpenes not covered by species MTPA. See [[Secondary_organic_aerosols#Add_MTPO_as_an_advected_and_chemical_species_to_all_full-chemistry_simulations|this post on our ''Secondary organic aerosols'' wiki page]] for more details.
  
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:28, 15 November 2017 (UTC)
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:28, 15 November 2017 (UTC)
Line 1,448: Line 1,406:
 
=== Fix dust rainout settings to match BCPO ===
 
=== Fix dust rainout settings to match BCPO ===
  
<span style="color:darkorange">'''''These fixes will be included in [[GEOS-Chem v11-02#v11-02e|v11-02e]].'''''</span>
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02e|v11-02e]] (approved 24 Mar 2018).'''''</span>
  
 
'''''Colette Heald (MIT) wrote:'''''
 
'''''Colette Heald (MIT) wrote:'''''
Line 1,456: Line 1,414:
 
--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 15:44, 30 November 2017 (UTC)
 
--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 15:44, 30 November 2017 (UTC)
  
===  Fixes for isoprene SOA for consistency with Marais et al. (2016) ===
+
=== Change timestep units from minutes to seconds ===
  
<span style="color:darkorange">'''''These fixes will be included in [[GEOS-Chem v11-02#v11-02e|v11-02e]].'''''</span>
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02e|v11-02e]] (approved 24 Mar 2018).'''''</span>
 +
 
 +
'''''Andrea Molod (GMAO) wrote:'''''
 +
 
 +
:As part of the side by side testing of GMI and GEOSChem chemistry in our coupled framework, we have run across a small but annoying issue. The "annoyance" is that the geoschem code does not accomodate a call at a time step that is not a multiple of a minute. our current model "heartbeat" is 450 seconds (for many reasons) at what is now a wide range of resolutions, and we end up having to run with geoschem at 300 seconds (for other reasons internal to GEOS).
 +
 
 +
To address this issue, the timestep duration specifications in the configuration file <tt>input.geos</tt> will be changed from minutes to seconds. The data type will stay as integer. This change has implications for the GEOS-Chem source code. Subroutines of form <tt>GET_TS_{component}</tt> now return seconds instead of minutes and usages of time units since a reference time must incorporate seconds as well as minutes in the calculation.
 +
 
 +
--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 23:19, 29 January 2018 (UTC)
 +
 
 +
===  Fixes for two new ND21 binary diagnostics introduced in v11-02 ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02e|v11-02e]] (approved 24 Mar 2018).'''''</span>
 +
 
 +
Two new ND21 diagnostics introduced in [[GEOS-Chem v11-02#v11-02c|v11-02c]] are assigned to the wrong tracer numbers in the AD21 array written to binary output. The optical depth of SOA from aqueous isoprene is written to tracer 34 but should be written to tracer 58. The aqueous aerosol volume is written to tracer 35 but should be written to tracer 59.
 +
 
 +
These tracer number errors resulted in an inconsistency with metadata supplied by GAMAP and prevented the diagnostics from being written to file during a full chemistry simulation. The new diagnostics also overwrite the dust AOD wavelength 2 diagnostics for bins 2 and 3 if run with RRTMG. A parallelization issue in setting one of the diagnostics is also present.
 +
 
 +
To get around these issues in v11-02c and v11-02d, replace AD21 index 34 with 58 and 35 with 59 where the diagnostics are set in <tt>aerosol_mod.F</tt> and declare local variable <tt>TAERVOL</tt> as private in the OMP loop in subroutine <tt>RDAER</tt>.
 +
 
 +
--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 23:19, 29 January 2018 (UTC)
 +
 
 +
=== Use spun-up restart file to initialize benchmark simulations ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02e|v11-02e]] (approved 24 Mar 2018).'''''</span>
 +
 
 +
Starting in [[GEOS-Chem v11-02#v11-02e|GEOS-Chem v11-02e]], we have replaced the default restart file ([http://ftp.as.harvard.edu/gcgrid/data/ExtData/SPC_RESTARTS/initial_GEOSChem_rst.4x5_standard.nc initial_GEOSChem_rst.4x5_standard.nc]) used to initialize 1-month benchmark simulations and UCX-based simulations in general. The previous default restart file had not been spun up in many years, so several species concentrations did not reflect updates to the model. The new restart file was obtained from the [[GEOS-Chem_v11-02_benchmark_history#v11-02d-Run2|1-year benchmark v11-02d-Run2]] and has been spun up for many years. The updated restart file also allowed us to [[UCX_chemistry_mechanism#v11-02|remove the initial stratospheric 2D mixing ratio option from input.geos]] and use the values provided by the restart file to initialize species concentrations for the relevant UCX species.
 +
 
 +
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 18:26, 16 March 2018 (UTC)
 +
 
 +
=== Restore use of prod/loss rates from GMI in the mesosphere for UCX-based simulations ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018).'''''</span>
 +
 
 +
Following the implementation of the [[Stratospheric_chemistry#Prod.2Floss_rates_from_UCX|prod/loss rates from UCX]] in [[GEOS-Chem v11-02#v11-02e|v11-02e]], it was discovered that we need to restore the use of prod/loss rates from GMI in the mesosphere for UCX-based simulations. The prod/loss rates saved out from a UCX simulation via KPP are only saved out on 59 levels (i.e. the maximum number of levels in the chemistry grid). Values in levels 59-72 (~1 hPa to ~ 0.01 hPa according to our vertical grids wiki page) are set to zero.
 +
 
 +
'''''Chris Holmes wrote:'''''
 +
 
 +
:I think it’s important to do some kind of mesospheric P/L in the UCX simulations. If we’re going to have a stratospheric chemistry capability in GEOS-Chem, we need to have a reasonable upper boundary condition (and zero P/L in the mesosphere is not reasonable). It seems to me that the easiest thing is to retain the GMI P/L in the mesosphere in UCX simulations.
 +
 
 +
'''''Sebastian Eastham wrote:'''''
 +
 
 +
:I think that's the most sensible approach for now. One additional point that's worth raising is that, in the course of the evaluation of the UCX by the GMI group, they actually ended up extending the UCX up into the mesosphere. This removed the need for any mesospheric paramaterizations, and did not involve adding any new chemistry or code. They found that this gave the same results as the GMI model - no big surprise since the GMI and UCX models are very similar, and the GMI model is already run up to the model top. As a result we may want to consider (in v11-03) adding an "extended UCX" option which runs UCX directly up to the model top which will resolve this problem when generating P/L data. It will also allow interested users to avoid reading in all the P/L data (a not-insignificant startup cost) or using mesospheric paramaterizations.
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 15:56, 17 May 2018 (UTC)
 +
 
 +
=== Fix bug in units of UCX P/L rates===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018).'''''</span>
 +
 
 +
'''''Christoph Keller wrote:'''''
 +
 
 +
:I think I found a bug in the implementation of the UCX P/L rates: in subroutine Set_PLVEC (in strat_chem_mod.F90), there is a unit conversion of the UCX P/L fields (from molec/cm3/s to v/v/s and 1/s, respectively), which will be erased after the first simulation month because the arrays in the PLVEC(:) structure are only *pointers* to the 3D fields managed by HEMCO (PLVEC(N)%PROD and PLVEC(N)%LOSS). After the first simulation month, HEMCO automatically updates the UCX P/L arrays (because they change monthly), thus overwriting the unit conversion performed in Set_PLVEC. So I'd expect the first month to be ok but then after the second month, GEOS-Chem will actually apply UCX P/L loss rates that are in molec/cm3/s instead of v/v/s and 1/s.
 +
 
 +
:In short, the arrays in PLVEC (PLVEC(N)%PROD and PLVEC(N)%LOSS) should not be modified because HEMCO changes them automatically.
 +
 
 +
The solution for this issue is to apply the unit conversion to the prod/loss rates at the time that we copy them into the loss frequency (k) and production term (P) variables in routine DO_STRAT_CHEM.
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 15:56, 17 May 2018 (UTC)
 +
 
 +
=== Use 1x1 GEIA NH3 emissions ===
 +
 
 +
<span style="color:red">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018), but was later removed in [[#GEOS-Chem v11-02 release candidate|v11-02-release-candidate]]. See below for details.'''''</span>
 +
 
 +
'''''Eleanor Morris wrote:'''''
 +
 
 +
:We have been looking at the ammonia emissions in the model in our comparisons to observational data from Africa and have a couple of questions.
 +
 
 +
:Although the GEIA anthropogenic emissions are overwritten by other inventories, for Africa at least, the model still uses the GEIA Natural and Biofuel emission of NH3:
 +
 
 +
    GEIA_NH3_BIOF $ROOT/NH3/v2014-07/NH3_biofuel.geos.4x5.nc NH3_BIOFUEL 1998/1-12/1/ C xy  kg/m2/s NH3 - 2 1
 +
    GEIA_NH3_NATU $ROOT/NH3/v2014-07/NH3_geos.4x5.nc        NH3_NATU    1990/1-12/1/ C xy  kg/m2/s NH3 - 3 1
 +
 
 +
:The standard HEMCO_Config.rc file uses the 4x5 degree resolution emission files (NH3_biofuel.geos.4x5.nc and NH3_geos.4x5.nc) although there are files at 1x1 degree resolution in the same directory (NH3_biofuel.geos.1x1.nc and NH3_geos.1x1.nc). We were wondering why the default setting is to use the 4x5 degree resolution rather than the 1x1 degree resolution? It looks like the original GEIA NH3 emissions are at 1x1 degree resolution.
 +
 
 +
:There are large differences between the files at different resolutions. Attached are some plots of the NH3 natural emissions for the three resolutions of files in the $ROOT/NH3/v2014-07 folder. When we calculate the total global burden, the 1x1 degree emissions give a total of ~26 Tg NH3 yr-1 which is relatively similar to the Park et al. Transboundary Pollution paper which quotes 20 Tg N yr-1 global natural ammonia emissions. The 4x5 degree emissions are a factor of 20 lower than this and we wondered whether perhaps there was an error in the re-gridding from 1x1 degree to 4x5 degree?
 +
 
 +
:When we run our simulation using the default 4x5 degree emissions, the modelled ammonia concentrations are too low compared to observations, whereas when we run with the 1x1 degree resolution emissions we get much better agreement.
 +
 
 +
:We just wanted to know whether there was a reason that 4x5 is used as the standard emissions instead of 1x1 degree, and if not whether the HEMCO_Config.rc file might need updating to use the 1x1 degree emissions instead?
 +
 
 +
'''''Bob Yantosca wrote:'''''
 +
 
 +
:I am not 100% sure why we still use the GEIA NH3 emissions at 4x5 in the HEMCO.  It may have been that while we were developing HEMCO, we needed to keep those NH3 emissions at 4x5 in order to better match the emissions in the pre-HEMCO code, so that we could get “apples-to-apples” comparisions.  It could have been as simple as that.  I would recommend bringing this to the attention of the Emissions and Deposition Working Group. 
  
 
'''''Eloise Marais wrote:'''''
 
'''''Eloise Marais wrote:'''''
  
:The [[GEOS-Chem_v11-02#Add_bug_fixes_for_GLYX.2C_MGLY.2C_IEPOX.2C_and_IMAE_heterogeneous_chemistry_added_in_v11-02c|bug fixes that Seb found in the KPP implementation of the heterogeneous chemistry]] mean that the [[Secondary_organic_aerosols#SOA_formation_from_aqueous_isoprene_uptake|isoprene SOA]] is now no longer consistent with the version that I used to implement this chemistry. This is concerning, as I will no doubt receive emails asking why isoprene SOA in the public release isn’t consistent with what’s reported in my paper. This now requires that I provide updated effective Henry’s Law constants and possibly also yields of SOA precursors.
+
:GCST, in summary, can we use the GEIA natural NH3 emissions file at 1x1 resolution rather than 4x5? Currently the default is 4x5 and Eleanor from U York found large differences between these two files, where the 1x1 emissions are more consistent with observations. It looks like GEIA natural NH3 will eventually be retired and replaced with emissions that Amos Tai is working on, so this would be an interim fix.
  
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 17:12, 5 January 2018 (UTC)
+
<span style="color:red">'''''NOTE: The 1-year benchmark v11-02f-Run0 showed that the 1&deg; x 1&deg; NH3 emissions file was problematic, producing 20X too much emissions as expected.  This is probably due to a regridding issue.  We have since removed this update from the [[#GEOS-Chem v11-02 release candidate|v11-02-release-candidate]] version and this issue remains open. (Bob Y., 07 Jun 2018).'''''</span>
 +
 
 +
==== GEIA anthro and biofuel NH3 are skipped if either CEDS or EDGARv43 is used ====
 +
 
 +
The [[GCST]] has confirmed that HEMCO's default behavior is to skip the following GEIA anthropogenic (<tt>GEIA_NH3_ANTH</tt>) and biofuel (<tt>GEIA_NH3_BIOF</tt>) if either the CEDS or EDGARv43 is selected as the global inventory.  This is because these GEIA emissions are listed in the <tt>HEMCO_Config.rc</tt> files with a lower priority than CEDS or EDGARv43.
 +
 
 +
The GEIA NH3 natural emissions (<tt>GEIA_NH3_NATU</tt>) and emissions from arctic seabirds (<tt>SEABIRD_NH3</tt>) will always be read, regardless of the choice of global base inventory.
 +
 
 +
For more information, please see [[GEOS-Chem_12#Avoid_double_counting_CEDS_and_GEIA_NH3_biofuel_emissions|this post on our ''GEOS-Chem 12'' wiki page]].
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 16:09, 3 July 2018 (UTC)
 +
 
 +
=== Fix for organic nitrate chemistry ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018).'''''</span>
 +
 
 +
'''''Xuan Wang wrote:'''''
 +
 
 +
:For MONIT and HONIT = INDOL + HNO3, I see from wiki we should use uptake rate of 2.78e-04 s-1 (lifetime = 1hour). However, we are setting gamma = 2.78e-04 s-1, not the uptake rate. I change it to manually setting the rate to 2.78e-04 s-1.
 +
 
 +
'''''Jenny Fisher wrote:'''''
 +
 
 +
:In my version there were 3 gas-phase monoterpene nitrate species (MONITU, MONITS, HONIT) and one monoterpene nitrate aerosol species (let’s call it pONITm since it had a confusing name in my code), plus INDIOL.
 +
:In my version, for each gas phase species we had e.g. MONITU -> pONITm with gamma = 0.01
 +
:Then pONITm -> HNO3 + INDIOL with a lifetime of 1 hour (set as a simple first order rate constant, not as an uptake rate).
 +
 
 +
'''''Xuan Wang wrote:'''''
 +
 
 +
:In v11-02, your pONITm seems to be IONITA and MONITA. For each gas it is like MONITU = MONITA with gamma = 0.01. For MONITA/IONITA =HNO3 + INDIOL, it is like this:
 +
 
 +
        ! Define gamma for MONITA
 +
        ! Imposed lifetime = 1 hour (Fisher et al., 2016)
 +
        XSTKCF = 2.78e-4_fp
 +
 
 +
:And later:
 +
        ! Reaction rate for surface of aerosol
 +
        ADJUSTEDRATE=ARSL1K(XAREA(N),XRADI(N),XDENA,XSTKCF,XTEMP, &
 +
                              (A**0.5_FP))
 +
 
 +
'''''Jenny Fisher wrote:'''''
 +
 
 +
:Ok, this does seem like an error I think, although I don’t really know how to fix it in the framework of v11-02. Should we even be using XSTKCF at all? Is there a way to do this in the v11-02 chemistry framework to just set the rate constant to 2.78e-4 s^-1?
 +
 
 +
'''''Melissa Sulprizio wrote:'''''
 +
 
 +
:If you would like to just set the rate constant for this reaction to 2.78e-4 then we can change it from a heterogeneous reaction in the .eqn files (i.e. utilizing the HET function in gckpp_HetRates.F90) to a standard reaction utilizing the Arrhenius equation (i.e. function GCARR). For example:
 +
 
 +
    IONITA = INDIOL + HNO3 :  HET(ind_IONITA,1);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}
 +
    MONITA = INDIOL + HNO3 :  HET(ind_MONITA,1);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}
 +
 
 +
:becomes
 +
 
 +
    IONITA = INDIOL + HNO3 :  GCARR(2.78E-04, 0.0E+00, 0.0);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}
 +
    MONITA = INDIOL + HNO3 :  GCARR(2.78E-04, 0.0E+00, 0.0);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 15:57, 17 May 2018 (UTC)
 +
 
 +
=== Fix compatibility issue with ifort17 ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018).'''''</span>
 +
 
 +
GEOS-Chem v11-02e with the netCDF diagnostics option enabled is not compatible with Intel Fortran compiler 17.0.4 and likely other ifort17 versions as well. The source code compiles but a string concatenation error is encountered in state_diag_mod.F90 during run-time. It appears that ifort17 is stricter than earlier versions of intel compilers with regards to concatenating indexed character arrays. Setting the indexed character array to a temporary variable prior to concatenation prevents the issue.
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 16:09, 17 May 2018 (UTC)
 +
 
 +
=== KPP errors with halogen chemistry and tropchem-based mechanisms ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018).'''''</span>
 +
 
 +
'''''[[User:Melissa_Payer|Melissa Sulprizio]] wrote:'''''
 +
 
 +
:Katie Travis discovered an issue starting with [[GEOS-Chem v11-02#v11-02d|v11-02d]] when the halogen chemistry updates were introduced. When running 2x2.5 tropchem simulations for an extended period (usually a month will do), GEOS-Chem crashes with a KPP integrator error. This happens for both GEOS-FP and MERRA-2 met fields. We don’t seem to encounter these errors for 2x2.5 standard (i.e. UCX-based) simulations. I’m wondering if we need to change the tolerance values in v11-02? Right now, in flexchem_mod.F we have:
 +
 
 +
    !%%%%% CONVERGENCE CRITERIA %%%%%
 +
 +
    ! Absolute tolerance
 +
    ATOL      = 1e-2_dp 
 +
 +
    ! Relative tolerance
 +
    IF ( Input_Opt%LUCX  ) THEN
 +
      ! UCX-based mechanisms
 +
      !RTOL      = 2e-2_dp
 +
      !RTOL      = 1e-2_dp
 +
      RTOL      = 0.5e-2_dp
 +
    ELSE
 +
      ! Non-UCX mechanisms
 +
      RTOL      = 2e-1_dp
 +
    ENDIF
 +
 
 +
:I ran two simulations to see if we can resolve the issue in the 2x2.5 tropchem simulations. In the first simulation I changed the convergence criteria (RTOL) to 1e-1_dp following Mat's suggestion and in the second I changed the value to 1e-2_dp. In both simulations, we still get the integrator errors (<tt>Forced exit from Rosenbrock due to the following error: Step size too small...</tt>) but the 1-month runs finish successfully because the error is not encountered twice for a given grid box in a given timestep. Also, comparing the two simulations, we see the integrator errors at different times and for different grid boxes.
 +
 
 +
'''''Barron Henderson wrote:'''''
 +
 
 +
:I have a couple thoughts that you may have already tried. They basically boil down to finding the few species at fault and tightening tolerances for everything else.
 +
 
 +
:First, add a print statement for SPC_NAMES and VAR and in flexchem_mod.F90 when the IERR is -7. That will help us to quickly see if any concentrations at cells are unexpected.
 +
 
 +
:Second, add a print statement for Yerr after 535 in gckpp_Integrator.F90. That should further isolate the source of error and will only print when the error is triggered.
 +
 
 +
:Hopefully, a few species will jump out. Then, the tolerances can targeted. Right now, it looks like tolerances are set as a vector, but all values are identical. The values are pretty loose for many species, but I suspect they were chosen to allow integration for a few radical species. One common approach for this situation is to identify a few key species and set the ATOL/RTOL separately for those species.
 +
 
 +
    # set default
 +
    ATOL = 1e-3
 +
    RTOL = 1e-5
 +
    # set species specific
 +
    ATOL(ind_O1D) = 1e-2
 +
    RTOL(ind_O1D) = 1e-2
 +
    ATOL(ind_XXX) = 1e-2
 +
    RTOL(ind_XXX) = 1e-2
 +
 
 +
:Tightening the tolerances will make it more stable, and oddly often faster.
 +
 
 +
'''''Mat Evans wrote:'''''
 +
 
 +
:So even if we set the tolerance to be 1e-2_dp we still get a convergence error.  Its hard to know what to advise here. It sounds like we need a bit of work to be done looking at this and looking at the situation that the convergence fails and to get an understanding of when we get the problem. However given the pressure at the moment to get v11.2 done I think my advise would be for us to move to using 1e-2 in the interim and we come back to this.
 +
 
 +
:The major problem I can see from this is in the run time. I don’t know how much longer the simulation is to run with this as the convergence criteria.
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 16:09, 17 May 2018 (UTC)
 +
 
 +
=== Disable Xiao propane emissions to avoid double-counting with CEDS ===
 +
 
 +
<span style="color:green">'''''This fix was incldued in [[#GEOS-Chem v11-02 release candidate|GEOS-Chem v11-02-rc]] (issued 22 Jun 2018).'''''</span>
 +
 
 +
'''''[[User:Lizzie Lundgren|Lizzie Lundgren]] wrote:'''''
 +
 
 +
<blockquote>I am turning off the XIAO inventory in v11-02 but just noticed Daniel said "for propane" while Christoph said to turn off the entire inventory. Christoph pointed out this would avoid double-counting C2H6 since we also get it from [the C2H6_2010 inventory of] Tzompa-Sosa et al. (2017)</blockquote>
 +
 
 +
'''''[[User:Christoph Keller|Christoph Keller]] wrote:'''''
 +
 
 +
<blockquote>I should have been more specific: in the current configuration, XIAO emissions are only used for propane. For ethane, the Tzompa-Sosa emissions always have priority over XIAO (hierarchy of 1000 vs. 5) and XIAO ethane emissions are thus ignored. The problem with propane is that CEDS and XIAO have the same hierarchy (=5) and the emissions are therefore added. Just disabling XIAO will fix this - and won't have an impact on ethane.</blockquote>
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 18:54, 29 June 2018 (UTC)
 +
 
 +
=== Bug fix: Specifying NO_REDUCED=no now compiles GEOS-Chem for reduced grids ===
 +
 
 +
<span style="color:green">'''''This fix was included in [[GEOS-Chem v11-02#v11-02a|v11-02a]] and approved on 12 May 2017.'''''</span>
 +
 
 +
'''''[[User:Zhuangjw|Jiawei Zhang]] wrote:'''''
 +
 
 +
<blockquote>I found a small but annoying bug in <tt>Makefile_header.mk</tt>.  The problem is that if I use <tt>NO_REDUCED=n</tt> (i.e. use 47 levels) as a compile option, then the makefile will treat it as <tt>NO_REDUCED=yes</tt> (i.e. use 72 levels). The <tt>LLPAR</tt> parameter will then be set to 72 instead of 47, leading to unexpected behavior.
 +
 
 +
That's because <tt>Makefile_header.mk</tt> contains something like:</blockquote>
 +
 
 +
      # %%%%% REDUCED VERTICAL GRID (default, unless specified otherwise) %%%%
 +
      ifndef NO_REDUCED
 +
        USER_DEFS        += -DGRIDREDUCED
 +
      else
 +
        REGEXP          :=(^[Yy]|^[Yy][Ee][Ss])
 +
        ifeq ($(shell [[ <nowiki>"$(NO_REDUCED)" =~ $(REGEXP)</nowiki> ]] && echo true),true)
 +
        endif
 +
      endif
 +
 
 +
<blockquote>The <tt>GRIDREDUCED</tt> C-preprocessor switch will only be activated if <tt>NO_REDUCED</tt> is not provided. The following "yes" or "no" just doesn't matter because the nested <tt>ifeq</tt> statement doesn't do anything. With this problem, the only way to compile a 47L model is NOT providing the <tt>NO_REDUCED</tt> option. It also affects compiling GEOS-Chem from the <tt>geosfp_4x5_tropchem</tt> rundir, because <tt>NO_REDUCED=n</tt> is one of the default options.</blockquote>
 +
 
 +
We have now fixed this issue by setting the <tt>NO_REDUCED</tt> Makefile variable to "no" by default, and then evaluating its value in a separate IF statement.  This code has now been added to the <tt>Makefile_header.mk</tt>.
 +
 
 +
  # %%%%% REDUCED VERTICAL GRID (default, unless specified otherwise) %%%%
 +
  ifndef NO_REDUCED
 +
    NO_REDUCED      :=no
 +
  endif
 +
  REGEXP              :=(^[Nn]|^[Nn][Oo])
 +
  ifeq ($(shell [[ <nowiki>"$(NO_REDUCED)" =~ $(REGEXP)</nowiki> ]] && echo true),true)
 +
    USER_DEFS        += -DGRIDREDUCED
 +
  endif
 +
 
 +
This fix now will compile GEOS-Chem for 47 levels when <tt>NO_REDUCED=no</tt> is specified, or if <tt>NO_REDUCED</tt> is omitted.
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 16:09, 1 March 2017 (UTC)
  
 
== Outstanding issues not yet resolved in GEOS-Chem v11-02 ==
 
== Outstanding issues not yet resolved in GEOS-Chem v11-02 ==
  
TBD
+
=== GCHP issues ===
 +
 
 +
For a complete list of outstanding issues not yet resolved in GEOS-Chem v11-02 with High Performance option, see the [[GCHP_v11-02#Outstanding_HP_issues_not_yet_resolved_in_GEOS-Chem_v11-02|outstanding issues section]] of the GCHP v11-02 wiki page.
 +
 
 +
=== Discrepancies in single-stage simulations vs multi-stage simulations ===
 +
 
 +
'''''Jack Kodros wrote'''''
 +
 
 +
<blockquote>I am having an issue re-starting simulations. I was trying to restart a simulation midway through the year using the restart file the sim had created, yet I noticed I started getting odd results.
 +
 
 +
To test things out I ran a tropchem sim starting in February and ran straight through until April. I then  ran a second, separate, but identical sim with the March restart file and HEMCO restart file from the first sim. I then compared the two sim outputs for the month of March. I had expected to get identical results, but instead I get small differences. I am attaching a percent change plot of OC, BC, NH3, and SO4 from this simulation showing the changes.
 +
 
 +
This is a concerning issue for me because I will now get different results if my simulation runs straight through or if I have to restart the simulation part way through.</blockquote>
 +
 
 +
'''''[[User:Christoph Keller|Christoph Keller]] wrote:'''''
 +
 
 +
<blockquote>At least for standard GEOS-Chem the updates should be pretty simple. It should be as easy as adding a few more variables to the GEOS-Chem and HEMCO restart files:
 +
 
 +
#Move PSC_RESTART variable from HEMCO restart file to GEOS-Chem restart file.
 +
#Archive fields DRY_TOTN, WET_TOTN, into HEMCO restart file.
 +
#Archive H2O2s, SO2s, and HSAVE_FOR_KPP as GEOS-Chem restart file variables.  (NOTE: HSAVE_FOR_KPP might be obsolete in v11-02.)
 +
</blockquote>
  
For a complete list of outstanding issues not yet resolved in GEOS-Chem v11-02 with High Performance option, see the [[GEOS-Chem_HP_v11-02#Outstanding_HP_issues_not_yet_resolved_in_GEOS-Chem_v11-02|outstanding issues section]] of the GCHP v11-02 wiki page.
+
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 18:19, 24 May 2018 (UTC)

Latest revision as of 21:07, 20 September 2022



GEOS-Chem v11-02-final will also carry the designation GEOS-Chem 12.0.0. We are migrating to a purely numeric versioning system in order to adhere more closely to software development best practices. For a complete description of the new versioning system, please see our GEOS-Chem version numbering system wiki page.




Contents

Version history

The table below shows the previous, current, and successive versions of GEOS-Chem:

Previous version This version Next version
GEOS-Chem v11-01-public GEOS-Chem 12.0.0
(aka v11-02-final)

GEOS-Chem 12.1.0
Released 01 Feb 2017 Released 10 Aug 2018 Release TBD
GC v11-01 online manual GC 12 online manual GC 12 benchmark history

The sections below describe each GEOS-Chem version in the v11-02 series. Versions v11-02a, v11-02b, etc. were validated with 1-month benchmark simulations (and 1-year benchmark simulations if necessary).

GEOS-Chem v11-02-final (aka 12.0.0)

This version was released on 10 Aug 2018

For a list of the features that were added to GEOS-Chem during the v11-02 "release candidate" testing period, please see visit this section of our GEOS-Chem 12 wiki page.

--Bob Yantosca (talk) 14:33, 21 June 2018 (UTC)

GEOS-Chem v11-02 release candidate

NOTE: These updates were not received in time for the v11-02f benchmarks, but were included in the v11-02-release-candidate version (aka v11-02-rc).

Feature Submitted by Type Status
Features affecting the full-chemistry simulation:
Only add ISOAAQ species to the PM2.5 diagnostics when using the Complex SOA option Aerosols WG Bug fix Released 22 Jun 2018
Only add ISOAAQ species to OCPISOA when using the complex SOA option Aerosols WG Bug fix " "
Disable Xiao propane emissions to avoid double-counting with CEDS Christoph Keller (NASA/GMAO)
GCST
Bug fix " "
Grid offset bug fix in high-resolution MEGAN files at southern mid-latitudes Jenny Fisher (Wollongong)
Dylan Millet (UMN)
Bug fix " "
Features NOT affecting the full-chemistry simulation:
Minor fixes (including updated parameters) for the Mercury simulation. Colin Thackray (Harvard) Bug Fix
Science
" "

--Bob Yantosca (talk) 19:16, 21 June 2018 (UTC)

v11-02f

Please see the following link for complete information about the validation of GEOS-Chem v11-02f:

  1. Approval form for 1-month benchmark simulation v11-02f
  2. Results for 1-year benchmark simulation v11-02f-Run0
Feature Submitted by Type Status
Features affecting the full-chemistry simulation:
Update HEMCO from v2.1.001 to v2.1.005 Christoph Keller (NASA GMAO) Structural 1-month benchmark:
Approved on 10 May 2018

1-year benchmark:
Approved on 17 May 2018

Use global 0.25x0.25 files for MEGAN input and MODIS LAI Jenny Fisher (U. Wollongong)
Eloise Marais (U. Birmingham)
Dylan Millet (UMN)
Bug fix/
Science
" "
Update lightning BETA values and OTD/LIS files for consistency with latest met fields Lee Murray (Rochester) Science " "
CEDS global anthropogenic emissions (1950-2014)
  • This inventory will be the default global anthropogenic emissions inventory
  • Includes a fix for NOx unit conversion
Lu Shen (Harvard)
Melissa Sulprizio (Harvard)
Science " "
EDGAR v4.3.1 emissions
  • This inventory will be included as an option
Chi Li (Dalhousie) Science " "
Use 1x1 GEIA NH3 emissions Eleanor Morris (York) Bug fix " "
Historical Canadian emissions Chi Li (Dalhousie) Science " "
Global ethane emissions Zitely Tzompa-Sosa (CSU) Science " "
Inefficient combustion emissions in Africa (DICE-Africa inventory) Eloise Marais (U. Birmingham)
Christine Wiedinmyer (NCAR)
Science " "
Trash burning emissions Eloise Marais (U. Birmingham)
Christine Wiedinmyer (NCAR)
Science " "
Update GFED4 emission factors Katie Travis (MIT) Bug fix " "
Retire GFED3 biomass burning emissions Katie Travis (MIT) Structural " "
Ammonia emissions from Arctic seabirds Betty Croft (Dalhousie) Science " "
Fixes for prod/loss rates used in linearized chemistry routines: Chris Holmes (FSU)
Christoph Keller (NASA GMAO)
Bug fix " "
Add fix for UCX mesospheric NOy routine Sebastian Eastham (MIT) Bug fix " "
Fix Henry's Law parameters for HCl Xuan Wang (Harvard) Bug fix " "
Fix for N2O5 uptake Xuan Wang (Harvard) Bug fix " "
Add fix for organic nitrate chemistry Xuan Wang (Harvard) Bug fix " "
Features not affecting the full-chemistry simulation:
Updated CH4 simulation Bram Maasakkers (Harvard)
Melissa Sulprizio (Harvard)
Science " "
Updated Tagged CO simulation Jenny Fisher (Wollongong) Science " "
Updates to Hg chemistry from Horowitz et al. (2017) Hannah Horowitz (U. Washington)
Colin Thackray (Harvard)
Rebecca Stern (Harvard)
Science " "
Change default timestep for MERRA-2 nested grid simulations Bram Maasakkers (Harvard) Science " "
Fix bug in biomass burning emissions when GFED is turned off Corey Trujillo (U. Colorado)
Therese Carter (MIT)
Bug fix " "
Fix missing HCOOH factors for RCP in HEMCO_Config.rc Seb Eastham (MIT) Bug fix " "
Fix bug in calculation of L(Ox) for tagged Ox simulations Lu Shen (Harvard) Bug fix " "
Change tolerance value as a temporary fix for convergence errors in 2x2.5 tropchem simulations Katie Travis (MIT)
Chemistry WG
Bug fix " "
Fix compatibility issue with ifort17 Bob Yantosca (GCST) Bug fix " "

v11-02e

Please see the following link for complete information about the validation of GEOS-Chem v11-02e:

  1. Approval form for 1-month benchmark simulation v11-02e
  2. Approval form for 1-month benchmark simulation v11-02e using a spun-up restart file
  3. Results for 1-year benchmark simulation v11-02e-Run0
  4. Results for 1-year Rn-Pb-Be benchmark simulation with fixed GEOS-FP met fields
  5. Results for 1-year Rn-Pb-Be benchmark simulation using 72 levels

We also completed a second set of benchmark simulations for GEOS-Chem v11-02e using 2016 GEOS-FP met fields. Please see the following links:

  1. Approval form for 1-month benchmark simulation v11-02e using 2016 GEOS-FP met fields
  2. Results for 1-year benchmark simulation v11-02e-Run1
  3. Results for 1-year Rn-Pb-Be benchmark simulation using 2016 GEOS-FP met fields
Feature Submitted by Type Status
Features affecting the full-chemistry simulation:
Monthly mean surface methane distributions Lee Murray (Rochester) Science 1-month benchmark:
Approved on 15 Mar 2018

1-year benchmark:
Approved 24 Mar 2018

UCX updates, including:
Seb Eastham (Harvard)
Chris Holmes (Florida State)
Melissa Sulprizio (GCST)
Science " "
Use monthly mean stratospheric prod/loss rates from UCX simulation for all species Melissa Sulprizio (GCST) Science " "
Update density of BC to 1.8 and add absorption enhancement factor in input.geos Xuan Wang (MIT) Science " "
Metal catalyzed oxidation of SO2 as a switch (on by default) Viral Shah (UW)
Becky Alexander (UW)
Science " "
Add SOA updates and bug fixes, including: Sam Silva (MIT)
Chemistry Working Group
Aerosols Working Group
GCST
Eloise Marais (U. Birmingham)
Science & Bug fix " "
Fix dust rainout settings to match BCPO Colette Heald (MIT) Bug fix " "
Fixes in gckpp_HetRates.F90 for halogen chemistry Lei Zhu (Harvard)
Seb Eastham (MIT)
Bug fix " "
Add bug fix in sulfate_mod.F involving reorder IF statements for Fahey and Seinfeld decision algorithm Chris Holmes (FSU)
Becky Alexander (UW)
Bug fix " "
Sync GEOS-FP files on Harvard ftp with files at Dalhousie GCST Bug fix " "
Features not affecting the full-chemistry simulation:
Spatially varying OM/OC as a switch (off by default) Sajeev Philip (Dalhousie/NASA Ames) Science " "
Update Rn-Pb-Be-Pasv simulations to use 72 levels by default Melissa Sulprizio (GCST) Science " "
Bug fixes for ISORROPIA II stable mode Jingyuan Shao (UW)
Shaojie Song (Harvard)
Bug fix " "
Change input.geos timestep units from minutes to seconds Lizzie Lundgren (GCST) Structural " "
Fixes for ND21 binary diagnostics Lizzie Lundgren (GCST) Bug fix " "
Fix bug in HEMCO reference time code Barron Henderson (US EPA) Bug fix " "

v11-02d

Please see the following link for complete information about the validation of GEOS-Chem v11-02d:

  1. Approval form for 1-month benchmark simulation v11-02d
  2. Results for 1-year benchmark simulation v11-02d-Run2
Feature Submitted by Type Status
Features affecting the full-chemistry simulation:
Halogen chemistry updates Tomás Sherwen (York)
Johan Schmidt (Harvard)
Seb Eastham (MIT)
Lei Zhu (Harvard)
Oxidants and Chemistry WG
Science Approved 12 Feb 2018
Sulfur oxidation by reactive halogens Qianjie Chen (UW) Science
Add bug fixes for GLYX, MGLY, IEPOX, and IMAE heterogeneous chemistry Sebastian Eastham (MIT) Bug fix
Fix bug impacting dust radius obtained in RDUST_ONLINE for UCX simulations Melissa Sulprizio (GCST) Bug fix
Add fix for IPMN reactions contributing to high P(Ox) and L(Ox) Melissa Sulprizio (GCST)
Eloïse Marais (U. Birmingham)
Bug fix
Features not affecting the full-chemistry simulation:
Update CH4 latitude bands for 2014-2016 Katie Travis (MIT) Science
Remove support for GCAP, GEOS-4, MERRA, GEOS-5 Melissa Sulprizio (GCST) Structural
Timestamp HEMCO diagnostics filename with start time rather than end time Lizzie Lundgren (GCST) Structural
Structural updates for netCDF diagnostics GCST Structural

v11-02c

Please see the following link for complete information about the validation of GEOS-Chem v11-02c:

  1. Approval form for 1-month benchmark simulation v11-02c
  2. Results for 1-year benchmark simulation v11-02c-Run0
Feature Submitted by Type Status
Features affecting the full-chemistry simulation:
Enhance default GEOS-Chem simple SOA Sal Farina (Colorado State)
Aerosols Working Group Chairs
Science 1-month benchmark:
Approved on 07 Sep 2017

1-year benchmark:
Approved 21 Sep 2017

Updates to isoprene and monoterpene chemistry Katie Travis (MIT)
Jenny Fisher (U. Wollongong)
Eloïse Marais (U. Birmingham)
Christopher Chan Miller (Harvard)
Kelvin Bates (Caltech)
Rebecca Schwantes(Caltech)
Science " "
Add aqueous isoprene uptake to SOA scheme Eloise Marais (U. Birmingham) Science " "
Carbon balance (fix C creation) Sarah Safieddine (MIT) Bug fix " "
Fix bugs for EOH and MGLY following implementation of PAN updates in v11-02a Melissa Sulprizio (GCST) Bug fix " "
Update HEMCO from v2.0.004 to v2.1.001 Christoph Keller (NASA GMAO) Structural " "
Features not affecting the full-chemistry simulation:
HEMCO updates:
Christoph Keller (NASA GMAO)
Andy Jacobson (NOAA)
Paulo Tuccella (L'Aquila)
Structural " "
Fixes for several minor issues:
Chris Holmes (Florida State)
Chris Holmes (Florida State)
Jiawei Zhuang (Harvard)
Amanda Giang (MIT)
Katie Travis (MIT)
Katie Travis(MIT)
Bug fix " "
Fixes for the TOMAS simulation:
Sal Farina (Colorado State)
Jack Kodros (Colorado State)
Bob Yantosca (GCST)
Bug fix " "
Fix STE flux diagnostics and add to benchmark procedure Melissa Sulprizio (GCST) Bug fix &
Benchmarking
" "
Initial structural modifications for netCDF diagnostics:
  • Introduce Headers/State_Diag as a stub module (for now)
  • Add a registry object into State_Met, State_Diag, and State_Chm in order to obtain a pointer to any module variable (or slice) by looking up its name
  • Add new module Headers/registry_mod.F90 which contains derived types and routines for registering module variables.
Bob Yantosca (GCST) Structural " "
Combine timestep settings in input.geos in a Timesteps menu Melissa Sulprizio (GCST) Structural " "
Update CO data used in 1-year benchmark plots Jenny Fisher (U. Wollongong) Benchmarking " "

v11-02b

Please see the following link for complete information about the validation of GEOS-Chem v11-02b:

  1. Approval form for 1-month benchmark simulation v11-02b
  2. Approval form for 1-month benchmark simulation v11-02b with high performance option (v11-02b-HP)
  3. Results for 1-year standard benchmark simulation of v11-02b-HP Run0
  4. Results for 1-year standard benchmark simulation of v11-02b-HP Run1 (uses offline archived lightning and dust emissions)
  5. Results for 1-year RnPbBe benchmark simulation of v11-02b-HP
Feature Submitted by Type Status
Features not affecting the full-chemistry simulation:
Source code updates for high performance option (GCHP) Seb Eastham (Harvard)
Lizzie Lundgren (Harvard)
Mike Long (Harvard)
Jiawei Zhuang (Harvard)
Bob Yantosca (Harvard)
Structural Approved 16 Jun 2017
Bug fixes for diagnostics:
Ilya Stanevic (Toronto)
Lee Murray (Rochester)
Bug fix " "
Write initial and final Ox mass to file when using the tagged O3 simulation Bob Yantosca (GCST) Structural " "
Convert CO2 emissions units using dry pressure used in advection Meemong Lee (JPL) Science " "
Add QFED emissions for 2014-2016 Christoph Keller (NASA GMAO) Science " "

v11-02a

Please see the following link for complete information about the validation of GEOS-Chem v11-02a:

  1. Approval form for 1-month benchmark simulation v11-02a
  2. Results for 1-year benchmark simulation v11-02a-Run0
  3. Results for 1-year benchmark simulation v11-02a-Run1 (final version)
Feature Submitted by Type Status
Features affecting the full-chemistry simulation:
Update chemistry rate constants based on JPL Publication 15-10 Barron Henderson (US EPA),
Mat Evans (U. York), &
Oxidants and Chemistry WG
Science Approved 12 May 2017
Fixes to correct ALK4 lumping issue Barron Henderson (US EPA) Science " "
PAN updates, including:
  • Updates to emissions, deposition, species, chemistry
  • Ocean+biogenic ALD2 emissions from Millet al. [2010]
Emily Fischer (CSU) Science " "
Monthly mean NEI2011 emissions GCST &
Katie Travis (Harvard)
Science " "
Bug fixes in the GEOS-Chem sulfate module:
Prasad Kasibhatla (Duke)
Prasad Kasibhatla (Duke)
Qianjie Chen (UW)
Viral Shah (UW)
Bug fix " "
Fix bug in dry deposition aerodynamic resistance Brian Boys (Dalhousie) Bug fix " "
Fix acetone parameterization in hcox_seaflux_mod.F90 GCST Bug fix " "
Bug fix in cos(SZA) for start of timestep Lizzie Lundgren (GCST) Bug fix " "
Include TOA pressure when calculating dry pressure edges Seb Eastham (Harvard) Science " "
Features not affecting the full-chemistry simulation:
Implement ISORROPIA v2.0 as a Fortran module Seb Eastham (Harvard) &
GCST
Structural " "
Updates to the HEMCO emissions component:
Christoph Keller (NASA GMAO)
Seb Eastham (Harvard)
Jessica Morena (Dalhousie)
Brian Boys (Dalhousie) & GCST
Structural
& Bug fix
" "
Updates to gain computational speedup:
Mike Long (Harvard)
Bob Yantosca (GCST)
Structural
& Bug fix
" "
Fixed typo in INIT_WINDOW Bob Yantosca (GCST) Bug fix " "
netCDF file I/O updates:
Chris Holmes (Florida State)
Andy Jacobson (NOAA)
GCST
Bug fix &
Structural
" "
Makefile and build sequence updates:
Jiawei Zhang (Harvard)
GCST
Bug fix &
Structural
" "
Bug fixes for running UCX in ESMF environment Christoph Keller (NASA GMAO) Bug fix " "
Bug fixes for diagnostics:
Aaron van Donkelaar (Dalhousie)
Jenny Fisher (U. Wollongong)
GCST
Jenny Fisher (U. Wollongong)
Chris Holmes (Florida State)
Bug fix " "
Removal of obsolete variables: GCST Structural " "

New data directories

The following HEMCO data directories have been added or updated for GEOS-Chem v11-02. You will have to download the directories relevant to your simulation.

Inventory What was added? When added? Directory
MEGAN Global 0.25x0.3125 CLM4 PFT and emission factor input files for MEGAN v11-02f HEMCO/MEGAN/v2017-07/
MODIS LAI Updated 0.25 x 0.25 MODIS LAI files (2005-2011) v11-02f CHEM_INPUTS/MODIS_LAI_201707/
Lightning NOx GEOS-FP OTD-LIS redistribution factors for Apr 2012 - Jul 2017 and MERRA-2 OTD-LIS redistribution factors for any date v11-02f HEMCO/LIGHTNOX/v2017-09/
Volcanic SO2 emissions OMI-based volcanic emissions for 2005-2012 v11-02f HEMCO/VOLCANO/v2018-03/
APEI Historical Canadian emissions v11-02f HEMCO/APEI/v2016-11/
EDGAR v4.3 Annual global anthropogenic emissions for 1970-2010 at 0.1x0.1 resolution v11-02f HEMCO/EDGARv43/v2016-11/
C2H6_2010 Global fossil fuel and biofuel emissions of C2H6 for 2010 v11-02f HEMCO/C2H6_2010/v2017-05/
DICE-Africa Biofuel and diffuse anthropogenic emissions for Africa v11-02f HEMCO/DICE_Africa/v2016-10/
Trash burning emissions Global trash burning emissions from Wiedinmyer et al. (2014) v11-02f HEMCO/TrashEmis/v2015-03/
SEABIRD NH3 NH3 emissions from arctic sea birds v11-02f HEMCO/NH3/v2018-04
OMOC Spatially varying OM/OC ratio v11-02e HEMCO/OMOC/v2018-01
NOAA GMD CH4 Monthly mean surface methane distributions v11-02e HEMCO/NOAA_GMD/v2018-01
UCX Input files for applying linearized stratospheric chemistry production and loss rates saved out from a UCX simulation using GEOS-Chem v11-02d v11-02e HEMCO/UCX/v2018-02/
IODINE Monthly emission fluxes of organic iodine compounds (CH3I, CH2I2, CH2ICl, CH2IBr) from Carlos Ordonez. v11-02d HEMCO/IODINE/v2017-03/
GMI New symbolic links:
  • gmi.clim.RIPA.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25
  • gmi.clim.RIPB.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25
  • gmi.clim.RIPD.geos5.2x25.nc -> gmi.clim.RIP.geos5.2x25
  • gmi.clim.IPMN.geos5.2x25.nc -> gmi.clim.PMN.geos5.2x25
  • gmi.clim.NPMN.geos5.2x25.nc -> gmi.clim.PMN.geos5.2x25
v11-02c HEMCO/GMI/v2015-02/
QFED Data for 2014-2016 v11-02b HEMCO/QFED/v2014-09/
NEI2011 North American emissions Monthly mean data files (created from the existing hourly data in HEMCO/NEI2011/v2015-03) v11-02a HEMCO/NEI2011/v2017-02-MM/
ALD2 Fields for ALD2 emissions, including:
  • Seawater concentration of acetaldehyde
  • Heterotrophic respiration rates, used to compute biogenic emissions of ALD2 and EOH

(added with the PAN updates)

v11-02a HEMCO/ALD2/v2017-03/
POET Anthropogenic emissions of ethanol from the POET inventory
(added with the PAN updates)
v11-02a HEMCO/POET/v2017-03/
Soil NOx New file: DepReservoirDefault.nc
Read default DEP_RESERVOIR fields from file when not found in HEMCO restart file
v11-02a HEMCO/SOILNOX/v2014-07/

Validation

In this section we provide information about the benchmarks and tests that we have done to validate GEOS-Chem v11-02.

1-month and 1-year benchmarks

For complete information about the benchmark simulations used to validate GEOS-Chem v11-02, please see our GEOS-Chem v11-02 benchmark history wiki page.

Unit tests for GEOS-Chem v11-02

We perform a full set of unit tests on GEOS-Chem on an almost daily basis. These tests ensure that each new feature being added does not break any existing functionality.

Here are the unit tests that we ran before submitting 1-month benchmark simulation for v11-02. Click on the links below to see the results from each set of unit tests:

Unit test Date
v11-02a 17 Apr 2017
v11-02b 18 May 2017
v11-02c 5 Sep 2017
v11-02d 5 Jan 2018
v11-02e 07 Mar 2018
v11-02f 02 May 2018

Previous issues now resolved in GEOS-Chem v11-02

The following bugs and/or technical issues have been resolved in v11-02.

Removed the NNPAR parameter from CMN_SIZE_mod.F

This update was included in v11-02a and approved on 12 May 2017.

The NNPAR parameter is a holdover from when we used COMMON blocks to store GEOS-Chem arrays. With the advent of the GEOS-Chem species database, we now use State_Chm%nSpecies to represent the total number of species in the chemistry mechanism, and State_Chm%nAdvect to represent the total number of advected species.

We have therefore removed the following lines (in RED) from CMN_SIZE_mod.F:

      !=================================================================
      ! TRACER & EMISSION SPECIES PARAMETERS
      !
      ! NNPAR   = max number of tracers
      ! NEMPARA = max number of anthropogenic emission species
      ! NEMPARB = max number of biogenic      emission species
      !=================================================================
      ! increase NNPAR and NEMPARA an extra amount (hotp 7/31/09)
#if   defined( TOMAS )
# if    defined( TOMAS40 )
      INTEGER,    PARAMETER :: NNPAR   = 430   ! For TOMAS40 (sfarina 6/11/13)
# elif  defined( TOMAS15 )
      INTEGER,    PARAMETER :: NNPAR   = 205   ! For TOMAS15 (sfarina 6/11/13)
# elif  defined( TOMAS12 )
      INTEGER,    PARAMETER :: NNPAR   = 178   ! For TOMAS12 (sfarina 6/11/13)
# else
      INTEGER,    PARAMETER :: NNPAR   = 340   ! For TOMAS (win, bmy, sfarina 6/11/13)f
# endif
#else
      INTEGER,    PARAMETER :: NNPAR   = 150   ! For non-TOMAS simulations
#endif

--Bob Yantosca (talk) 20:27, 17 March 2017 (UTC)

Print the version number to the log file

This update was included in v11-02a and approved on 12 May 2017.

In routine DISPLAY_GRID_AND_MODEL of GeosCore/main.F, we now print out the version number, e.g.

Using GEOS-Chem version: v11-02

to the log file (i.e. stdout output stream). This will facilitate debugging.

--Bob Yantosca (talk) 18:21, 16 March 2017 (UTC)

Include TOA pressure when calculating dry pressure edges

This fix was included in v11-02a and approved on 12 May 2017.

For physicality and consistency with GCHP, the top-of-atmosphere (TOA) pressure is now included when calculating dry pressure edges. The following lines in GeosCore/dao_mod.F were changed from:

     ! Reset dry surface pressure to zero
     PS_DRY = 0.0e+0_fp

to:

     ! Reset dry surface pressure to TOA value
     PS_DRY = GET_AP(LLPAR+1)

--Melissa Sulprizio (talk) 20:47, 14 April 2017 (UTC)

ND26 vertical flux diagnostic does not yield mass balance in advection

This update was included in v11-02b and approved on 16 Jun 2017.

To achieve mass balance in the advection diagnostic output (ND24, ND25, and ND26) Ilya Stanevic (U. Toronto) suggested that the FZ array be used to compute the ND26 vertical mass flux diagnostic for all levels. This is similar to how FX is used for ND24 and FY for ND25. Note that this fix is only for global advection since the FZ array is not defined for the nested model. With this update, the global ND26 vertical flux diagnostic is defined as the flux into the grid box above (e.g. ND26(:,:,1) is flux from level 1 to level 2 and ND26(:,:,LLPAR) is all zeros).

--Lizzie Lundgren (talk) 21:41, 9 May 2017 (UTC)

ND40 planeflight diagnostic not always written out when turned on

This update was included in v11-02b and approved on 16 Jun 2017.

We have moved the discussion of this issue to our Planeflight diagnostics wiki page.

CO2 emissions unit conversion not consistent with new advection pressure

This update was included in v11-02b and approved on 16 Jun 2017.

We have moved the discussion of this issue to our CO2 simulation wiki page.

--Lizzie Lundgren (talk) 19:44, 10 May 2017 (UTC)

Fixed an incorrect format statement in input_mod.F

This update was included in v11-02c and approved on 21 Sep 2017.

Chris Holmes wrote:

An incorrect format specification caused GEOS-Chem to crash rather than write a helpful error message when the transport timestep is too big. The code in RED was removed and the lines in GREEN was added at approximately line 3324 of input_mod.F.
        IF ( Input_Opt%TS_DYN > MAX_DYN .and. LTRAN ) THEN
           write(MSG,'(I8,a,L)') '<>', Input_Opt%TS_DYN, Input_Opt%HPC
           WRITE( MSG, 300 ) 'Transport timestep exceeds max:',
    &                         Input_Opt%TS_DYN, MAX_DYN
300        FORMAT( a, i8, ' >', i8 ) 
           CALL ERROR_STOP( MSG, LOCATION )
        ENDIF

--Bob Yantosca (talk) 17:55, 18 April 2017 (UTC)

Bug fixes for the ND50 timeseries diagnostic

This update was included in v11-02c and approved on 21 Sep 2017.

Chris Holmes wrote:

I have found what I suspect are two bugs in ND50.

1. Function ITS_TIME_FOR_WRITE_DIAG50() checks if the integer part of the hour is zero to decide whether to write to disk. With the 20 minute time steps recommended in v11-01, Diag50 gets written and overwritten multiple times in the first hour of every day! The fix is simply to have this function return TRUE if the GMT time is 00:00.

2. Within WRITE_DIAG50 there is the following block of code
        ! Pick the proper divisor, depending on whether or not the
        ! species in question is archived only each chem timestep
        SELECT CASE ( ND50_TRACERS(W) )
           CASE (91, 92, 76, 77 )
              DIVISOR = COUNT_CHEM3D
           CASE DEFAULT
              DIVISOR = COUNT
        END SELECT

I see nothing special about tracers 76, 77, 91, and 92 that should make them divided by a different number. The special tracers numbers are all >= 150 for ND50. I suspect this is left over from a previous code revision. I think it should be DIVISOR=COUNT for all ND50_Tracers.

--Bob Yantosca (talk) 14:42, 14 April 2017 (UTC)

Add MERRA2 to #if and #elif statements where it had been omitted

This update was included in v11-02c and approved on 21 Sep 2017.

Jiawei Zhuang wrote:

In v11-02, (commit 8d27b9a), I am trying to unify defined( GEOS_FP ) and defined( MERRA ) by a adding a new tag defined( FlexGrid ). However, I find some model codes that treat them differently without explanation:
(1) In GeosCore/convection_mod.F: At the #if statement at line 1443, MERRA2 will execute the #else block. But it should execute the #if block:
      #if defined( GEOS_5 ) || defined( GEOS_FP ) || defined( MERRA )
          ! Do nothing, these are all hybrid grids
      #else
          ... etc ...
(2) In HEMCO/Core/hco_interp_mod.F90, MERRA2 will skip the #elif block at line 531:
      #elif defined( GEOS_5 ) || defined( MERRA ) || defined( GEOS_FP )
          ! Full grid
          ... etc ...
(3) At line 1193 of GeosCore/hcoi_gc_main_mod.F90, I noticed what appears to be a typo:
      #if defined ( GEOS_FP ) || ( MERRA2 )
(4) At line 779 of GeosCore/land_mercury_mod.F, I noticed this #if block, which could be written better:
            IF ( FIRSTTIME ) THEN
      #if !defined( MERRA2 ) || !defined( GEOS_FP )
               WRITE( 6, * ) 'WARNING: Snow emission parameterization '//
         &       'designed for use with incident '
               WRITE( 6, * ) 'shortwave radiation. Only net shortwave '//
         &       'radiation available.'
      #endif
               FIRSTTIME = .FALSE.
           ENDIF
(5) In KPP/Standard/gckpp_HetRates.F90, there are 3 places that seem to miss defined (MERRA2):
      ! continental or marine clouds only...
#if defined( GEOS_5 ) || defined( MERRA ) || defined( GEOS_FP )
      IF ( (FRLAND (I,J) > 0) .or. (FROCEAN(I,J) > 0) ) THEN
#else
      ! Above line is to skip over land ice (Greenland and Antartica). This
      ! should do the same (and also work for GEOS-5, but leave above for now).
      IF ( IS_LAND .and. .not. IS_ICE  ) THEN
#endif

#if defined( GEOS_5 ) || defined( MERRA ) || defined( GEOS_FP )
      IF ( FRLAND(I,J) > FROCEAN(I,J) ) THEN
#else
      IF ( IS_LAND ) THEN
#endif

#if defined( GEOS_5 ) || defined( MERRA ) || defined( GEOS_FP )
      ! QL is [g/g]
      Vc = CLDF(I,J,L) * QL * AD(I,J,L) / dens_h2o
#else
      ! QL is [cm3/cm3]
      Vc = CLDF(I,J,L) * QL * XAIRM3
#endif

These appear to have been oversights. To fix the issues, simply add the code in GREEN to the #if blocks mentioned above:

 At line 1443 of convection_mod.F: 
 #if defined( GEOS_5 ) || defined( GEOS_FP ) || defined( MERRA ) || defined( MERRA2 )

 At line 531 of hco_interp_mod.F90:
 #elif defined( GEOS_5 ) || defined( MERRA ) || defined( GEOS_FP ) || defined( MERRA2 )
 
 At line 1193 of hcoi_gc_main_mod.F90
 #if defined( GEOS_FP ) || defined( MERRA2 )

 At line 779 of land_mercury_mod.F:
 #if defined( GEOS_FP ) || defined( MERRA2 )
    ! do nothing
 #else
     WRITE( 6, * ) ... etc ...

 At lines 1953, 2029, and 2047 of KPP/Standard/gckpp_HetRates.F90
 #if defined( GEOS_5 ) || defined( MERRA ) || defined( GEOS_FP ) || defined( MERRA2 )

--Bob Yantosca (talk) 20:28, 22 May 2017 (UTC)

Code cleanup suggestions

Suggestion 1 was included in v11-02c and approved on 21 Sep 2017.

Suggestion 2 (Git ID: d06ae04a) was included in GEOS-Chem 12.1.0, which was released on 26 Nov 2018.

Chris Holmes wrote:

I have a few housekeeping suggestions that I think would be helpful to GC users.
  1. Combine all of the timestep settings in input.geos in a single “Timesteps” menu. With the default timesteps being much shorter now, there seems to be much more interest in changing these timesteps (e.g. do model development with large time steps, then publication runs with small time steps), myself included. It would be nice to change all of the time steps in one place. Besides making it easier to double or halve them all, a single menu would also help clarify whether the time steps for each process have the recommended 1:2 ratio for operator splitting.
  2. Move all of the photolysis and aerosol optics data files out of the run directory and into a folder in ExtData. That would allow them to be documented, tracked and updated like emissions. As far as I have seen, these files are rarely modified by users, so they don’t need to be in the run directory. That would include the following files...
      brc.dat
      dust.dat
      FJX_j2j.dat
      FJX_spec.dat
      h2so4.dat
      jv_spec_mie.dat
      org.dat
      so4.dat
      soot.dat
      ssa.dat
      ssc.dat

The .dat files required for FAST-JX can now be found in CHEM_INPUTS/FAST_JX/v2018-09/.

--Melissa Sulprizio (talk) 16:37, 19 July 2017 (UTC)

Update CH4 latitude bands for 2014-2016

This update was included in v11-02d (approved 12 Feb 2018).

The last year of latitudinal CH4 concentration data (based on CMDL flask observations) was 2013. (Data is provided in 4 latitude bands: 90°S-30°S, 30°S-0°N:, 0°N-30°N, 30°N-90°N.) Katherine Travis has updated the subroutine get_global_ch4.F with additional years of latitudinal CH4 data from 2014 thru 2016.

Lee Murray will provide additional latitudinal CH4 concentration data through present day, once it is available. See this post on our NOx-Ox-HC-aerosol wiki page for more information.

--Melissa Sulprizio (talk) 15:16, 5 September 2017 (UTC)
--Bob Yantosca (talk) 21:59, 20 February 2018 (UTC)

Add bug fixes for GLYX, MGLY, IEPOX, and IMAE heterogeneous chemistry added in v11-02c

This update was included in v11-02d (approved 12 Feb 2018).

Sebastian Eastham wrote:

There’s a bug in the v11-02c het rate code. In KPP/Standard/gckpp_HetRates.F90, many of the dedicated routines (e.g. HetIMAE) loop over NAERO. The loop looks like
   Do N=1,NAERO
      ! [Comment]
      If (N==8) Then
         XSTKCF = [number] 
      End If
      If (N == 13) Then
         ADJUSTEDRATE = XAREA(N)*XSTKCF
      Else
         ADJUSTEDRATE = ARSL1K( XAREA(N), XRADI(N), XDENA, XSTKCF, XTEMP, (A**0.5_FP) )
      End If
      Het_IMAE = Het_IMAE + ADJUSTEDRATE
   End Do
The problem here is that every aerosol after aerosol 8 is inheriting the sticking coefficient from aerosol 8. Crucially, this includes SOAS – so the fact that the simple SOA tracer is missing could affect this. GCC is artificially increasing the het rate on all aerosols because of this bug, whereas GCHP is not. Fixing this would be as simple as inserting either XSTKCF = 0.0 or XSTKCF = B at the start of the loop, where I’ve written ! [Comment]. The question of which is correct is really a science question (probably XSTKCF = B, although I’d ask the original authors). However, what’s definitely true is that one of these two options needs to be put in or the result of the calculation will be wrong.
Fixing this bug is likely to change SOA burdens in both GCC and GCHP. In gckpp_HetRates, the following routines are definitely affected:
  • HetGLYX
  • HetMGLY
  • HetIEPOX (forms SOAIE)
  • HetIMAE (forms SOAME)
The other het routines look like they have the correct form (either by a "header" of XSTKCF = B, or because they have an If...ElseIf...Else...End clause that catches all possible values of N and sets XSTKCF accordingly).

The solution is to set XSTKCF = TINY(1e0) where ! [Comment] is in the code snippet above, following what is done in Eloise Marais's version of calcrate.F.

--Melissa Sulprizio (talk) 17:57, 20 September 2017 (UTC)

Add fix for IPMN reactions contributing to high P(Ox) and L(Ox)

This update was included in v11-02d (approved 12 Feb 2018).

During the 1-year benchmark for v11-02d, Lu Hu pointed out an issue in the O3 budget where P(Ox) and L(Ox) were too high (~40000 Tg/yr compared to the expected 4000-5000 Tg/yr). This issue was tracked to the following reactions in the chemical mechanism:

IMAO3 + NO2 {+M} = IPMN  # Impacts L(Ox)
IPMN = IMAO3 + NO2       # Impacts P(Ox)

Eloise Marais wrote:

I'd suggest just using the same reaction rate values as for NPMN: non-isoprene PMN, as the reaction rates I used come from a rather speculative chemical mechanism anyway proposed by Lin et al., 2013 (http://www.pnas.org/content/suppl/2013/04/03/1221150110.DCSupplemental/pnas.201221150SI.pdf) and the isoprene SOA that these reactions form is relatively minor.

This fix involves changing the reactions rates for the following IPMN reactions in the .eqn files from:

IMAO3 + NO2 {+M} = IPMN : GCJPLPR(1.21E-11, 1.1E+00, 0.0, 7.7E-12, 0.2, 0.0, 0.6, 0.0, 0.0);
NMAO3 + NO2 {+M} = NPMN : GCJPLPR(9.00E-28, 8.9E+00, 0.0, 7.7E-12, 0.2, 0.0, 0.6, 0.0, 0.0);
IPMN = IMAO3 + NO2 : GCJPLEQ(1.60E+16, 0.0E+00, -13486.0, 9.00E-28, 8.9E0, 0.0, 7.7E-12, 0.2, 0., 0.6, 0., 0.);
NPMN = NMAO3 + NO2 : GCJPLEQ(9.00E-29, 0.0E+00, 14000.0, 9.00E-28, 8.9E0, 0.0, 7.7E-12, 0.2, 0., 0.6, 0., 0.);

to:

IMAO3 + NO2 {+M} = IPMN : GCJPLPR(9.00E-28, 8.9E+00, 0.0, 7.7E-12, 0.2, 0.0, 0.6, 0.0, 0.0);
NMAO3 + NO2 {+M} = NPMN : GCJPLPR(9.00E-28, 8.9E+00, 0.0, 7.7E-12, 0.2, 0.0, 0.6, 0.0, 0.0);
IPMN = IMAO3 + NO2 : GCJPLEQ(9.00E-29, 0.0E+00, 14000.0, 9.00E-28, 8.9E0, 0.0, 7.7E-12, 0.2, 0., 0.6, 0., 0.);
NPMN = NMAO3 + NO2 : GCJPLEQ(9.00E-29, 0.0E+00, 14000.0, 9.00E-28, 8.9E0, 0.0, 7.7E-12, 0.2, 0., 0.6, 0., 0.);

In addition IPMN should be added to the Ox family in gckpp.kpp:

POx : O3 + NO2 + 2NO3 + PAN + NPMN + IPMN + PPN + HNO4 + 3N2O5 + HNO3 + BrO + HOBr + BrNO2 + 2BrNO3 + MPN + ETHLN + ISN1 + ISOPNB + ISOPND + MACRN + MVKN + PROPNN + R4N2 + INPN + ISNP + INO2 + ISNOOA + ISNOOB + ISNOHOO + MAN2 + PRN1 + PRPN + R4N1 + PMNN + MACRNO2 + ClO + HOCl + ClNO2 + 2ClNO3 + 2Cl2O2 + 2OClO + O + O1D + IO + HOI + IONO + 2IONO2 + 2OIO + 2I2O2 + 3I2O3 + 4I2O4;
LOx : O3 + NO2 + 2NO3 + PAN + NPMN + IPMN + PPN + HNO4 + 3N2O5 + HNO3 + BrO + HOBr + BrNO2 + 2BrNO3 + MPN + ETHLN + ISN1 + ISOPNB + ISOPND + MACRN + MVKN + PROPNN + R4N2 + INPN + ISNP + INO2 + ISNOOA + ISNOOB + ISNOHOO + MAN2 + PRN1 + PRPN + R4N1 + PMNN + MACRNO2 + ClO + HOCl + ClNO2 + 2ClNO3 + 2Cl2O2 + 2OClO + O + O1D + IO + HOI + IONO + 2IONO2 + 2OIO + 2I2O2 + 3I2O3 + 4I2O4;

After implementing those fixes, the chemical mechanisms will need to be rebuilt with KPP.

--Melissa Sulprizio (talk) 21:42, 25 January 2018 (UTC)

Remove support for GCAP, GEOS-4, MERRA, GEOS-5

This update was included in v11-02d (approved 12 Feb 2018).

Daniel Jacob wrote to the GEOS-Chem Steering Committee:

I met today with the GCST and they brought up again the issue of discontinuing support for old met fields (GEOS-4, GEOS-5, MERRA) in GEOS-Chem. They told me that it would benefit the code structurally to stop supporting these old met fields, also that GCHP will not be able to use them since the data are in binary format. We would also reclaim some disk storage, and we would satisfy GMAO who doesn’t like that we’re using these old met fields.
We will discontinue support for the old met fields (MERRA, GEOS-4, GEOS-5) in future model versions BUT we will keep the data archives [at Dalhousie] at least for now.

Loretta Mickley consulted with Lee Murray, Amos Tai, and Jeff Geddes of the Chemistry-Ecosystems-Climate Working Group and wrote:

It looks like we can safely retire the GCAP met fields and #ifdef GCAP sections of the code now. As Lee describes, there are other capabilities in use by the chemistry-climate community that we will eventually implement into standard code.
I believe we currently have GCAP2 met fields stored at Harvard, and these should be maintained. GCAP2 consists of met fields from the GISS ModelE2.1 present-day and future. Lee Murray and Eric Leibensperger created these files, after a lot of work!

--Melissa Sulprizio (talk) 16:11, 2 November 2017 (UTC)
--Bob Yantosca (talk) 21:59, 20 February 2018 (UTC)

Bug fixes for aerosols in v11-01

This fix was included in v11-02c.

Jingyuan Shao (UW) wrote:

When I use GEOS-Chem v11-01, I found two small bugs in the codes.
The first is described in more detail at this post on our ISORROPIA II wiki page.
Second, In aerosol_mod.F(line 2077) from
     ! Mass of hydrophilic aerosol from Mian Chin
     ALLOCATE( WAERSL(IIPAR,JJPAR,LLPAR,NAER), STAT=AS)
     IF ( AS /= 0) CALL ALLOC_ERR( 'WAERSL' )
     DAERSL = 0d0
to
     ! Mass of hydrophilic aerosol from Mian Chin
     ALLOCATE( WAERSL(IIPAR,JJPAR,LLPAR,NAER), STAT=AS)
     IF ( AS /= 0) CALL ALLOC_ERR( 'WAERSL' )
     WAERSL = 0d0

The fix in aerosol_mod.F was included in v11-02c with the implementation of aqueous isoprene uptake updates from Eloise Marais.

--Melissa Sulprizio (talk) 16:28, 15 November 2017 (UTC)

Remove MONX from chemical mechanisms

This update was included in v11-02e (approved 24 Mar 2018).

Sam Silva (MIT) wrote:

I was wondering why there were three monoterpene species in the chemical mechanism: MTPA (lumped alpha pinene), LIMO (lumped limonene) and then MONX (all monoterpenes).
MONX doesn’t play with any of the monoterpene nitrate chemistry, it just is converted into 3 formaldehyde:
   MONX = 3.000CH2O :                          GCARR(1.07E-05, 0.0E+00, 0.0);
It seems like having the MONX group (all monoterpenes) and the speciated lumped monoterpenes is an inconsistency. Ultimately, I’m interested in adding Glyoxal formation from monoterpenes, and it’s not clear to me how to properly do that, given the three different species.
Do you know of a reason why MONX is still in the mechanism?

Melissa Sulprizio responded:

I had a look through the Git history of the globchem.dat file. The MONX -> 3CH2O reaction was originally included in the SOA mechanism. Just prior to the v10-01 public release, the GEOS-Chem Steering Committee requested that we include SOA chemistry in the benchmark simulations. At that time, MONX was made an active species and the MONX -> 3CH2O reaction was added to the globchem.dat file for the benchmark simulations. No other monoterpene chemistry was included in that chemical mechanism at the time.
In v11-02a, we introduced monoterpene chemistry during Emily’s PAN updates. At that time MTPA was introduced. The monoterpene chemistry was further modified and expanded upon in v11-02c following Jenny’s updates originally made for SEAC4RS.
Looking back at Emily’s original PAN updates that she sent to us (based on v9-02h), it looks like MONX is a dead species in globchem.dat and the MONX -> 3CH2O reaction is not included. That may have been an oversight during the PAN updates implementation. Since Emily’s globchem.dat did not include the SOA chemistry by default and our standard/benchmark mechanism did, it wasn’t clear that we should remove the MONX species and reaction.

The MONX species will be removed from GEOS-Chem mechanisms in v11-02e, upon the recommendation given to us by Emily Fischer. This update also requires the addition of the MTPO species to represents all the monoterpenes not covered by species MTPA. See this post on our Secondary organic aerosols wiki page for more details.

--Melissa Sulprizio (talk) 16:28, 15 November 2017 (UTC)

Fix dust rainout settings to match BCPO

This update was included in v11-02e (approved 24 Mar 2018).

Colette Heald (MIT) wrote:

BCPO/DST1-4: are all "insoluble" but can act as IN, so I think the DST RainEff should be modified to [1,1,0] to match BCPO

--Lizzie Lundgren (talk) 15:44, 30 November 2017 (UTC)

Change timestep units from minutes to seconds

This update was included in v11-02e (approved 24 Mar 2018).

Andrea Molod (GMAO) wrote:

As part of the side by side testing of GMI and GEOSChem chemistry in our coupled framework, we have run across a small but annoying issue. The "annoyance" is that the geoschem code does not accomodate a call at a time step that is not a multiple of a minute. our current model "heartbeat" is 450 seconds (for many reasons) at what is now a wide range of resolutions, and we end up having to run with geoschem at 300 seconds (for other reasons internal to GEOS).

To address this issue, the timestep duration specifications in the configuration file input.geos will be changed from minutes to seconds. The data type will stay as integer. This change has implications for the GEOS-Chem source code. Subroutines of form GET_TS_{component} now return seconds instead of minutes and usages of time units since a reference time must incorporate seconds as well as minutes in the calculation.

--Lizzie Lundgren (talk) 23:19, 29 January 2018 (UTC)

Fixes for two new ND21 binary diagnostics introduced in v11-02

This update was included in v11-02e (approved 24 Mar 2018).

Two new ND21 diagnostics introduced in v11-02c are assigned to the wrong tracer numbers in the AD21 array written to binary output. The optical depth of SOA from aqueous isoprene is written to tracer 34 but should be written to tracer 58. The aqueous aerosol volume is written to tracer 35 but should be written to tracer 59.

These tracer number errors resulted in an inconsistency with metadata supplied by GAMAP and prevented the diagnostics from being written to file during a full chemistry simulation. The new diagnostics also overwrite the dust AOD wavelength 2 diagnostics for bins 2 and 3 if run with RRTMG. A parallelization issue in setting one of the diagnostics is also present.

To get around these issues in v11-02c and v11-02d, replace AD21 index 34 with 58 and 35 with 59 where the diagnostics are set in aerosol_mod.F and declare local variable TAERVOL as private in the OMP loop in subroutine RDAER.

--Lizzie Lundgren (talk) 23:19, 29 January 2018 (UTC)

Use spun-up restart file to initialize benchmark simulations

This update was included in v11-02e (approved 24 Mar 2018).

Starting in GEOS-Chem v11-02e, we have replaced the default restart file (initial_GEOSChem_rst.4x5_standard.nc) used to initialize 1-month benchmark simulations and UCX-based simulations in general. The previous default restart file had not been spun up in many years, so several species concentrations did not reflect updates to the model. The new restart file was obtained from the 1-year benchmark v11-02d-Run2 and has been spun up for many years. The updated restart file also allowed us to remove the initial stratospheric 2D mixing ratio option from input.geos and use the values provided by the restart file to initialize species concentrations for the relevant UCX species.

--Melissa Sulprizio (talk) 18:26, 16 March 2018 (UTC)

Restore use of prod/loss rates from GMI in the mesosphere for UCX-based simulations

This update was included in v11-02f (approved 17 May 2018).

Following the implementation of the prod/loss rates from UCX in v11-02e, it was discovered that we need to restore the use of prod/loss rates from GMI in the mesosphere for UCX-based simulations. The prod/loss rates saved out from a UCX simulation via KPP are only saved out on 59 levels (i.e. the maximum number of levels in the chemistry grid). Values in levels 59-72 (~1 hPa to ~ 0.01 hPa according to our vertical grids wiki page) are set to zero.

Chris Holmes wrote:

I think it’s important to do some kind of mesospheric P/L in the UCX simulations. If we’re going to have a stratospheric chemistry capability in GEOS-Chem, we need to have a reasonable upper boundary condition (and zero P/L in the mesosphere is not reasonable). It seems to me that the easiest thing is to retain the GMI P/L in the mesosphere in UCX simulations.

Sebastian Eastham wrote:

I think that's the most sensible approach for now. One additional point that's worth raising is that, in the course of the evaluation of the UCX by the GMI group, they actually ended up extending the UCX up into the mesosphere. This removed the need for any mesospheric paramaterizations, and did not involve adding any new chemistry or code. They found that this gave the same results as the GMI model - no big surprise since the GMI and UCX models are very similar, and the GMI model is already run up to the model top. As a result we may want to consider (in v11-03) adding an "extended UCX" option which runs UCX directly up to the model top which will resolve this problem when generating P/L data. It will also allow interested users to avoid reading in all the P/L data (a not-insignificant startup cost) or using mesospheric paramaterizations.

--Bob Yantosca (talk) 15:56, 17 May 2018 (UTC)

Fix bug in units of UCX P/L rates

This update was included in v11-02f (approved 17 May 2018).

Christoph Keller wrote:

I think I found a bug in the implementation of the UCX P/L rates: in subroutine Set_PLVEC (in strat_chem_mod.F90), there is a unit conversion of the UCX P/L fields (from molec/cm3/s to v/v/s and 1/s, respectively), which will be erased after the first simulation month because the arrays in the PLVEC(:) structure are only *pointers* to the 3D fields managed by HEMCO (PLVEC(N)%PROD and PLVEC(N)%LOSS). After the first simulation month, HEMCO automatically updates the UCX P/L arrays (because they change monthly), thus overwriting the unit conversion performed in Set_PLVEC. So I'd expect the first month to be ok but then after the second month, GEOS-Chem will actually apply UCX P/L loss rates that are in molec/cm3/s instead of v/v/s and 1/s.
In short, the arrays in PLVEC (PLVEC(N)%PROD and PLVEC(N)%LOSS) should not be modified because HEMCO changes them automatically.

The solution for this issue is to apply the unit conversion to the prod/loss rates at the time that we copy them into the loss frequency (k) and production term (P) variables in routine DO_STRAT_CHEM.

--Bob Yantosca (talk) 15:56, 17 May 2018 (UTC)

Use 1x1 GEIA NH3 emissions

This update was included in v11-02f (approved 17 May 2018), but was later removed in v11-02-release-candidate. See below for details.

Eleanor Morris wrote:

We have been looking at the ammonia emissions in the model in our comparisons to observational data from Africa and have a couple of questions.
Although the GEIA anthropogenic emissions are overwritten by other inventories, for Africa at least, the model still uses the GEIA Natural and Biofuel emission of NH3:
    GEIA_NH3_BIOF $ROOT/NH3/v2014-07/NH3_biofuel.geos.4x5.nc NH3_BIOFUEL 1998/1-12/1/ C xy  kg/m2/s NH3 - 2 1
    GEIA_NH3_NATU $ROOT/NH3/v2014-07/NH3_geos.4x5.nc         NH3_NATU    1990/1-12/1/ C xy  kg/m2/s NH3 - 3 1
The standard HEMCO_Config.rc file uses the 4x5 degree resolution emission files (NH3_biofuel.geos.4x5.nc and NH3_geos.4x5.nc) although there are files at 1x1 degree resolution in the same directory (NH3_biofuel.geos.1x1.nc and NH3_geos.1x1.nc). We were wondering why the default setting is to use the 4x5 degree resolution rather than the 1x1 degree resolution? It looks like the original GEIA NH3 emissions are at 1x1 degree resolution.
There are large differences between the files at different resolutions. Attached are some plots of the NH3 natural emissions for the three resolutions of files in the $ROOT/NH3/v2014-07 folder. When we calculate the total global burden, the 1x1 degree emissions give a total of ~26 Tg NH3 yr-1 which is relatively similar to the Park et al. Transboundary Pollution paper which quotes 20 Tg N yr-1 global natural ammonia emissions. The 4x5 degree emissions are a factor of 20 lower than this and we wondered whether perhaps there was an error in the re-gridding from 1x1 degree to 4x5 degree?
When we run our simulation using the default 4x5 degree emissions, the modelled ammonia concentrations are too low compared to observations, whereas when we run with the 1x1 degree resolution emissions we get much better agreement.
We just wanted to know whether there was a reason that 4x5 is used as the standard emissions instead of 1x1 degree, and if not whether the HEMCO_Config.rc file might need updating to use the 1x1 degree emissions instead?

Bob Yantosca wrote:

I am not 100% sure why we still use the GEIA NH3 emissions at 4x5 in the HEMCO. It may have been that while we were developing HEMCO, we needed to keep those NH3 emissions at 4x5 in order to better match the emissions in the pre-HEMCO code, so that we could get “apples-to-apples” comparisions. It could have been as simple as that. I would recommend bringing this to the attention of the Emissions and Deposition Working Group.

Eloise Marais wrote:

GCST, in summary, can we use the GEIA natural NH3 emissions file at 1x1 resolution rather than 4x5? Currently the default is 4x5 and Eleanor from U York found large differences between these two files, where the 1x1 emissions are more consistent with observations. It looks like GEIA natural NH3 will eventually be retired and replaced with emissions that Amos Tai is working on, so this would be an interim fix.

NOTE: The 1-year benchmark v11-02f-Run0 showed that the 1° x 1° NH3 emissions file was problematic, producing 20X too much emissions as expected. This is probably due to a regridding issue. We have since removed this update from the v11-02-release-candidate version and this issue remains open. (Bob Y., 07 Jun 2018).

GEIA anthro and biofuel NH3 are skipped if either CEDS or EDGARv43 is used

The GCST has confirmed that HEMCO's default behavior is to skip the following GEIA anthropogenic (GEIA_NH3_ANTH) and biofuel (GEIA_NH3_BIOF) if either the CEDS or EDGARv43 is selected as the global inventory. This is because these GEIA emissions are listed in the HEMCO_Config.rc files with a lower priority than CEDS or EDGARv43.

The GEIA NH3 natural emissions (GEIA_NH3_NATU) and emissions from arctic seabirds (SEABIRD_NH3) will always be read, regardless of the choice of global base inventory.

For more information, please see this post on our GEOS-Chem 12 wiki page.

--Bob Yantosca (talk) 16:09, 3 July 2018 (UTC)

Fix for organic nitrate chemistry

This update was included in v11-02f (approved 17 May 2018).

Xuan Wang wrote:

For MONIT and HONIT = INDOL + HNO3, I see from wiki we should use uptake rate of 2.78e-04 s-1 (lifetime = 1hour). However, we are setting gamma = 2.78e-04 s-1, not the uptake rate. I change it to manually setting the rate to 2.78e-04 s-1.

Jenny Fisher wrote:

In my version there were 3 gas-phase monoterpene nitrate species (MONITU, MONITS, HONIT) and one monoterpene nitrate aerosol species (let’s call it pONITm since it had a confusing name in my code), plus INDIOL.
In my version, for each gas phase species we had e.g. MONITU -> pONITm with gamma = 0.01
Then pONITm -> HNO3 + INDIOL with a lifetime of 1 hour (set as a simple first order rate constant, not as an uptake rate).

Xuan Wang wrote:

In v11-02, your pONITm seems to be IONITA and MONITA. For each gas it is like MONITU = MONITA with gamma = 0.01. For MONITA/IONITA =HNO3 + INDIOL, it is like this:
        ! Define gamma for MONITA 
        ! Imposed lifetime = 1 hour (Fisher et al., 2016)
        XSTKCF = 2.78e-4_fp
And later:
        ! Reaction rate for surface of aerosol
        ADJUSTEDRATE=ARSL1K(XAREA(N),XRADI(N),XDENA,XSTKCF,XTEMP, &
                              (A**0.5_FP))

Jenny Fisher wrote:

Ok, this does seem like an error I think, although I don’t really know how to fix it in the framework of v11-02. Should we even be using XSTKCF at all? Is there a way to do this in the v11-02 chemistry framework to just set the rate constant to 2.78e-4 s^-1?

Melissa Sulprizio wrote:

If you would like to just set the rate constant for this reaction to 2.78e-4 then we can change it from a heterogeneous reaction in the .eqn files (i.e. utilizing the HET function in gckpp_HetRates.F90) to a standard reaction utilizing the Arrhenius equation (i.e. function GCARR). For example:
    IONITA = INDIOL + HNO3 :  HET(ind_IONITA,1);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}
    MONITA = INDIOL + HNO3 :  HET(ind_MONITA,1);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}
becomes
    IONITA = INDIOL + HNO3 :  GCARR(2.78E-04, 0.0E+00, 0.0);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}
    MONITA = INDIOL + HNO3 :  GCARR(2.78E-04, 0.0E+00, 0.0);  {2017/07/14; Fisher2016; KRT,JAF,CCM,EAM,KHB,RHS}

--Bob Yantosca (talk) 15:57, 17 May 2018 (UTC)

Fix compatibility issue with ifort17

This update was included in v11-02f (approved 17 May 2018).

GEOS-Chem v11-02e with the netCDF diagnostics option enabled is not compatible with Intel Fortran compiler 17.0.4 and likely other ifort17 versions as well. The source code compiles but a string concatenation error is encountered in state_diag_mod.F90 during run-time. It appears that ifort17 is stricter than earlier versions of intel compilers with regards to concatenating indexed character arrays. Setting the indexed character array to a temporary variable prior to concatenation prevents the issue.

--Bob Yantosca (talk) 16:09, 17 May 2018 (UTC)

KPP errors with halogen chemistry and tropchem-based mechanisms

This update was included in v11-02f (approved 17 May 2018).

Melissa Sulprizio wrote:

Katie Travis discovered an issue starting with v11-02d when the halogen chemistry updates were introduced. When running 2x2.5 tropchem simulations for an extended period (usually a month will do), GEOS-Chem crashes with a KPP integrator error. This happens for both GEOS-FP and MERRA-2 met fields. We don’t seem to encounter these errors for 2x2.5 standard (i.e. UCX-based) simulations. I’m wondering if we need to change the tolerance values in v11-02? Right now, in flexchem_mod.F we have:
   !%%%%% CONVERGENCE CRITERIA %%%%%

   ! Absolute tolerance
   ATOL      = 1e-2_dp  

   ! Relative tolerance
   IF ( Input_Opt%LUCX  ) THEN
      ! UCX-based mechanisms
      !RTOL      = 2e-2_dp
      !RTOL      = 1e-2_dp
      RTOL      = 0.5e-2_dp
   ELSE
      ! Non-UCX mechanisms
      RTOL      = 2e-1_dp
   ENDIF
I ran two simulations to see if we can resolve the issue in the 2x2.5 tropchem simulations. In the first simulation I changed the convergence criteria (RTOL) to 1e-1_dp following Mat's suggestion and in the second I changed the value to 1e-2_dp. In both simulations, we still get the integrator errors (Forced exit from Rosenbrock due to the following error: Step size too small...) but the 1-month runs finish successfully because the error is not encountered twice for a given grid box in a given timestep. Also, comparing the two simulations, we see the integrator errors at different times and for different grid boxes.

Barron Henderson wrote:

I have a couple thoughts that you may have already tried. They basically boil down to finding the few species at fault and tightening tolerances for everything else.
First, add a print statement for SPC_NAMES and VAR and in flexchem_mod.F90 when the IERR is -7. That will help us to quickly see if any concentrations at cells are unexpected.
Second, add a print statement for Yerr after 535 in gckpp_Integrator.F90. That should further isolate the source of error and will only print when the error is triggered.
Hopefully, a few species will jump out. Then, the tolerances can targeted. Right now, it looks like tolerances are set as a vector, but all values are identical. The values are pretty loose for many species, but I suspect they were chosen to allow integration for a few radical species. One common approach for this situation is to identify a few key species and set the ATOL/RTOL separately for those species.
    # set default
    ATOL = 1e-3
    RTOL = 1e-5
    # set species specific
    ATOL(ind_O1D) = 1e-2
    RTOL(ind_O1D) = 1e-2
    ATOL(ind_XXX) = 1e-2
    RTOL(ind_XXX) = 1e-2
Tightening the tolerances will make it more stable, and oddly often faster.

Mat Evans wrote:

So even if we set the tolerance to be 1e-2_dp we still get a convergence error. Its hard to know what to advise here. It sounds like we need a bit of work to be done looking at this and looking at the situation that the convergence fails and to get an understanding of when we get the problem. However given the pressure at the moment to get v11.2 done I think my advise would be for us to move to using 1e-2 in the interim and we come back to this.
The major problem I can see from this is in the run time. I don’t know how much longer the simulation is to run with this as the convergence criteria.

--Bob Yantosca (talk) 16:09, 17 May 2018 (UTC)

Disable Xiao propane emissions to avoid double-counting with CEDS

This fix was incldued in GEOS-Chem v11-02-rc (issued 22 Jun 2018).

Lizzie Lundgren wrote:

I am turning off the XIAO inventory in v11-02 but just noticed Daniel said "for propane" while Christoph said to turn off the entire inventory. Christoph pointed out this would avoid double-counting C2H6 since we also get it from [the C2H6_2010 inventory of] Tzompa-Sosa et al. (2017)

Christoph Keller wrote:

I should have been more specific: in the current configuration, XIAO emissions are only used for propane. For ethane, the Tzompa-Sosa emissions always have priority over XIAO (hierarchy of 1000 vs. 5) and XIAO ethane emissions are thus ignored. The problem with propane is that CEDS and XIAO have the same hierarchy (=5) and the emissions are therefore added. Just disabling XIAO will fix this - and won't have an impact on ethane.

--Bob Yantosca (talk) 18:54, 29 June 2018 (UTC)

Bug fix: Specifying NO_REDUCED=no now compiles GEOS-Chem for reduced grids

This fix was included in v11-02a and approved on 12 May 2017.

Jiawei Zhang wrote:

I found a small but annoying bug in Makefile_header.mk. The problem is that if I use NO_REDUCED=n (i.e. use 47 levels) as a compile option, then the makefile will treat it as NO_REDUCED=yes (i.e. use 72 levels). The LLPAR parameter will then be set to 72 instead of 47, leading to unexpected behavior. That's because Makefile_header.mk contains something like:
      # %%%%% REDUCED VERTICAL GRID (default, unless specified otherwise) %%%%
      ifndef NO_REDUCED
        USER_DEFS        += -DGRIDREDUCED
      else
        REGEXP           :=(^[Yy]|^[Yy][Ee][Ss])
        ifeq ($(shell [[ "$(NO_REDUCED)" =~ $(REGEXP) ]] && echo true),true)
        endif
      endif
The GRIDREDUCED C-preprocessor switch will only be activated if NO_REDUCED is not provided. The following "yes" or "no" just doesn't matter because the nested ifeq statement doesn't do anything. With this problem, the only way to compile a 47L model is NOT providing the NO_REDUCED option. It also affects compiling GEOS-Chem from the geosfp_4x5_tropchem rundir, because NO_REDUCED=n is one of the default options.

We have now fixed this issue by setting the NO_REDUCED Makefile variable to "no" by default, and then evaluating its value in a separate IF statement. This code has now been added to the Makefile_header.mk.

  # %%%%% REDUCED VERTICAL GRID (default, unless specified otherwise) %%%%
  ifndef NO_REDUCED
    NO_REDUCED       :=no
  endif
  REGEXP              :=(^[Nn]|^[Nn][Oo])
  ifeq ($(shell [[ "$(NO_REDUCED)" =~ $(REGEXP) ]] && echo true),true)
    USER_DEFS        += -DGRIDREDUCED
  endif

This fix now will compile GEOS-Chem for 47 levels when NO_REDUCED=no is specified, or if NO_REDUCED is omitted.

--Bob Yantosca (talk) 16:09, 1 March 2017 (UTC)

Outstanding issues not yet resolved in GEOS-Chem v11-02

GCHP issues

For a complete list of outstanding issues not yet resolved in GEOS-Chem v11-02 with High Performance option, see the outstanding issues section of the GCHP v11-02 wiki page.

Discrepancies in single-stage simulations vs multi-stage simulations

Jack Kodros wrote

I am having an issue re-starting simulations. I was trying to restart a simulation midway through the year using the restart file the sim had created, yet I noticed I started getting odd results.

To test things out I ran a tropchem sim starting in February and ran straight through until April. I then ran a second, separate, but identical sim with the March restart file and HEMCO restart file from the first sim. I then compared the two sim outputs for the month of March. I had expected to get identical results, but instead I get small differences. I am attaching a percent change plot of OC, BC, NH3, and SO4 from this simulation showing the changes.

This is a concerning issue for me because I will now get different results if my simulation runs straight through or if I have to restart the simulation part way through.

Christoph Keller wrote:

At least for standard GEOS-Chem the updates should be pretty simple. It should be as easy as adding a few more variables to the GEOS-Chem and HEMCO restart files:
  1. Move PSC_RESTART variable from HEMCO restart file to GEOS-Chem restart file.
  2. Archive fields DRY_TOTN, WET_TOTN, into HEMCO restart file.
  3. Archive H2O2s, SO2s, and HSAVE_FOR_KPP as GEOS-Chem restart file variables. (NOTE: HSAVE_FOR_KPP might be obsolete in v11-02.)

--Bob Yantosca (talk) 18:19, 24 May 2018 (UTC)