Difference between revisions of "GEOS-Chem 12"
(→TOMAS40 simulations exit with a floating-point exception) |
(→12.0.0) |
||
Line 50: | Line 50: | ||
=== 12.0.0 === | === 12.0.0 === | ||
− | [[GEOS-Chem | + | The following updates and fixes are slated to be added to GEOS-Chem 12.0.0 (aka v11-02-final). Many of these updates fix issues that were identified during the [[GEOS-Chem v11-02#v11-02 GEOS-Chem v11-02 release candidate|GEOS-Chem v11-02 release candidate]] testing period. |
+ | |||
+ | Please see the following link for complete information about the validation of GEOS-Chem 12.0.0: | ||
+ | #[[GEOS-Chem_12_benchmark_history#12.0.0|''Approval form for 1-month benchmark simulation 12.0.0'']] | ||
+ | #[[GEOS-Chem_12_benchmark_history#12.0.0-Run0|''Results for 1-year benchmark simulation 12.0.0-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 (benchmark) simulation: | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[Trash_burning_emissions#Avoid_double_counting_trash_emissions|Turn off TRASH emissions by default to avoid double counting]] | ||
+ | |[[GCST]]<br>Eloise Marais (U. Birmingham) | ||
+ | |Bug fix | ||
+ | |Implemented 27 Jul 2018 | ||
+ | *This update introduces relatively small changes: | ||
+ | **Mean OH decreased by 0.12% | ||
+ | **MCF lifetime increased by 0.23% | ||
+ | **CH4 lifetime increased by 0.24% | ||
+ | **O3 STE flux increased by 0.042% | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem 12#Add_entries_to_HEMCO_Config.rc_for_SOAP_and_pFe_in_CEDS_inventory|Add entries to HEMCO_Config.rc for SOAP and pFe in CEDS inventory]] | ||
+ | |[[GCST]] | ||
+ | |Bug fix | ||
+ | |Implemented 27 Jul 2018 | ||
+ | *This update introduces relatively small changes: | ||
+ | **Mean OH decreased by 0.15% | ||
+ | **MCF lifetime increased by 0.28% | ||
+ | **CH4 lifetime increased by 0.29% | ||
+ | **O3 STE flux increased by 0.043% | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[#Bug fix for liquid water content in gckpp_HetRates.F90|Bug fix for liquid water content in gckpp_HetRates.F90]] | ||
+ | |Jiayue Huang (U. Washington)<br>Quanjie Chen (U. Washington) | ||
+ | |Bug fix | ||
+ | |Implemented 19 Jul 2018 | ||
+ | *This update introduces relatively small changes: | ||
+ | **Mean OH increased by 1.91% | ||
+ | **MCF lifetime decreased by 2.41% | ||
+ | **CH4 lifetime decreased by 2.41% | ||
+ | **O3 STE flux decreased by 0.18% | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem 12#Bug_fix_for_CH4_boundary_conditions|Bug fix for CH4 boundary conditions]] | ||
+ | |Lee Murray (Rochester) | ||
+ | |Bug fix | ||
+ | |Implemented 19 Jul 2018 | ||
+ | *Results approved by Lee Murray | ||
+ | *This update introduces small changes: | ||
+ | **Mean OH decreased by 0.48% | ||
+ | **MCF lifetime increased by 0.82% | ||
+ | **O3 STE decreased by: 0.022% | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem 12#Fixed typo in call to GAMMA_HOBr_AER|Fixed typo in call to GAMMA_HOBr_AER]] | ||
+ | |Xuan Wang<br>[[GCST]] | ||
+ | |Bug fix | ||
+ | |Implemented 18 Jul 2018 | ||
+ | *NOTE: This update produced only very small differences w/r/t internal run v11-02-rc-DICEnewrc | ||
+ | **Mean OH decreased by -0.058% | ||
+ | **MCF lifetime increased by 0.082% | ||
+ | **CH4 lifetime increased by 0.084% | ||
+ | **O3 STE flux decreased by 0.005% | ||
+ | |||
+ | |-valign="top" | ||
+ | |Correct issues when specifying [[DICE-Africa anthropogenic emissions inventory|DICE-Africa emissions]]: | ||
+ | *[[DICE-Africa_anthropogenic_emissions_inventory#Avoid_double-counting_between_DICE-Africa_and_other_inventories|Avoid double-counting between DICE-Africa and other inventories]] | ||
+ | *[[DICE-Africa_anthropogenic_emissions_inventory#Use_unique_names_for_DICE-Africa_entries_in_HEMCO_Config.rc|Use unique names for DICE-Africa entries in the <tt>HEMCO_Config.rc</tt> file]] | ||
+ | |Eloise Marais (U. Birmingham)<br>[[GCST]] | ||
+ | |Bug Fix | ||
+ | |Implemented 16 Jul 2018. | ||
+ | *Results approved by Eloise Marais. | ||
+ | *NOTE: This update produced only very small differences w/r/t internal run v11-02-rc-arctasship: | ||
+ | **Mean OH increased by 0.07% | ||
+ | **MCF lifetime decreased by 0.12% | ||
+ | **CH4 lifetime decreased by 0.13% | ||
+ | **O3 STE flux decreased by -0.0023% | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[Organics_Working_Group#Add_missing_Criegee_intermediate_reaction|Add missing Criegee intermediate reaction]] | ||
+ | |Xin Chen (UMN) | ||
+ | |Bug fix | ||
+ | |Implemented 11 Jul 2018, pending validation | ||
+ | *NOTE: This update produced only very small differences w/r/t internal run v11-02-rc-arctasship: | ||
+ | **Mean OH decreased by 0.00002% | ||
+ | **MCF lifetime was unchanged | ||
+ | **CH4 lifetime increased by 0.001% | ||
+ | **O3 STE flux increased by 0.0001% | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[Ship_emissions#Fix_bug_in_ARCTAS_SHIP_SO2_scale_factors|Bug fix for ARCTAS_SHIP in <tt>HEMCO_Config.rc</tt> files]] | ||
+ | |Barron Henderson (EPA) | ||
+ | |Bug Fix | ||
+ | |Implemented 03 Jul 2018 | ||
+ | *NOTE: This update only produced very small differences w/r/t internal run v11-02-rc-gfed2016: | ||
+ | **Mean OH decreased by 0.0004% | ||
+ | **O3 STE flux decreased by 0.0007% | ||
+ | **MCF and CH4 lifetimes were unchanged | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[FlexChem#FlexChem bug fix: do not zero ACTA, EOH, HCOOH|Bug fix in FlexChem: do not zero out ACTA, EOH, HCOOH concentrations]] | ||
+ | |Katie Travis (MIT) | ||
+ | |Bug Fix | ||
+ | |Implemented 27 Jun 2018 | ||
+ | *NOTE: This update only produced very small differences w/r/t v11-02-rc: | ||
+ | **Mean OH decreased by 0.083% | ||
+ | **MCF lifefime increased by 0.112% | ||
+ | **CH4 lifetime increased by 0.112% | ||
+ | **O3 STE flux decreased by 0.011% | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GFED4_biomass_burning_emissions#Years_2015_and_2016_have_been_added|Add GFED4.1 data for 2015 and 2016; Use GFED4.1 2016 in benchmarks]] | ||
+ | |Kate O'Dell (CSU)<br>[[GCST]] | ||
+ | |Science | ||
+ | |Data files downloaded 27 Jun 2018 | ||
+ | *NOTE: This update only produced small differences: | ||
+ | **Mean OH increased by 0.60% | ||
+ | **MCF lifefime decreased by 0.65% | ||
+ | **CH4 lifetime decreased by 0.64% | ||
+ | **O3 STE flux decreased by 0.0093% | ||
+ | |||
+ | |- | ||
+ | !colspan="4" bgcolor="#CCFFFF"|Features NOT affecting the full-chemistry simulation: | ||
+ | |||
+ | |- | ||
+ | |[[Mercury#Added_chemistry_and_ocean_from_Horowitz_et_al._2017|Minor fixes for mercury chemistry updates implemented in v11-02f]] | ||
+ | |Colin Thackray (Harvard) | ||
+ | |Bug fix | ||
+ | |Implemented 02 Aug 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[MERRA-2_implementation_details#Default_timestep_for_MERRA-2_nested_grid_simulations|Restore default timestep for MERRA-2 nested grid simulations to 10/20]] | ||
+ | |Ke Li (Harvard) | ||
+ | |Bug fix | ||
+ | |Implemented 30 Jul 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem 12#Remove_duplicate_C3HI_spectrum_in_FJX_spec.dat_and_FJX_j2j.dat|Remove duplicate C3HI spectrum in FJX_spec.dat and FJX_j2j.dat]] | ||
+ | |Tomas Sherwen (York) | ||
+ | |Bug fix | ||
+ | |Implemented 27 Jul 2018 | ||
+ | *A difference test w/r/t the previous commit yielded identical results. | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem 12#Fixed typo that prevented PM2.5 from being archived to netCDF output|Fixed typo that prevented PM2.5 from being archived to netCDF output]] | ||
+ | |Jonathan Moch (Harvard)<br>[[GCST]] | ||
+ | |Bug fix | ||
+ | |Implemented 24 Jul 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[Two-way_coupling_between_global_and_nested_GEOS-Chem_models|Bring two-way nested-grid code option up-to-date]] | ||
+ | |Yingying Yan (PKU)<br>Jintai Lin (PKU) | ||
+ | |Structural | ||
+ | |Implemented 24 Jul 2018 | ||
+ | *This option is invoked by compiling with <tt>EXCHANGE=y</tt> | ||
+ | *A difference test w/r/t the previous commit yielded identical results. | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[Planeflight_diagnostic#Do_not_force_L.3D1_by_default_in_planeflight_mod.F|Do not force L=1 by default in planeflight_mod.F]] | ||
+ | |Katie Travis (MIT)<br>[[GCST]] | ||
+ | |Bug fix | ||
+ | |Implemented 23 Jul 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[#Add PH2O2 to the default list of prod/loss families in the KPP mechanisms|Add PH2O2 to the default list of prod/loss families in the KPP mechanisms]] | ||
+ | |[[GCST]] | ||
+ | |Structural | ||
+ | |Implemented 19 Jul 2018 | ||
+ | *This update facilitates saving out oxidant fields from a full-chemistry simulation that can be used to drive an aerosol-only simulation. | ||
+ | |||
+ | |-valign="top" | ||
+ | |Update to HEMCO v2.1.007, which includes minor updates: | ||
+ | #[[Implementation of HEMCO in GEOS-Chem#Stop with error if multiple containers have the same name|Stop with error if multiple containers have the same name]] | ||
+ | #[[Implementation of HEMCO in GEOS-Chem#Bug fix for distributing emissions in the vertical dimension|Bug fix for distributing emissions in the vertical dimension]] | ||
+ | #[[Implementation of HEMCO in GEOS-Chem#Add extra error checks in HEMCO standalone module|Add extra error checks in HEMCO standalone module]] | ||
+ | #[[Implementation of HEMCO in GEOS-Chem#Remove null string character from netCDF unit string|Remove null string character from netCDF unit string]] | ||
+ | #[[Implementation_of_HEMCO_in_GEOS-Chem#Bug fix for HEMCO soil NOx error with ifort 17|Bug fix for HEMCO soil NOx error with ifort 17]] | ||
+ | |[[GCST]]<br>Christoph Keller<br>Jenny Fisher | ||
+ | |Structural &<br>Bug fix | ||
+ | |Update history: | ||
+ | #Implemented 20 Jul 2018 | ||
+ | #Implemented 19 Jul 2018 | ||
+ | #Implemented 18 Jul 2018 | ||
+ | #Implemented 17 Jul 2018 | ||
+ | #Implemented 26 Jun 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem 12#Bug fix for reference time in netCDF diagnostic files|Bug fix for reference time in netCDF diagnostic files]] | ||
+ | |[[GCST]] | ||
+ | |Bug fix | ||
+ | |Implemented 12 Jul 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |Updates for [[QFED biomass burning emissions]]: | ||
+ | #[[QFED biomass burning emissions#Update to QFED v2.5r1|Update to QFED v2.5r1]] | ||
+ | #[[QFED_biomass_burning_emissions#Fix_QFED_files_for_2016.2F10.2F13|Fix QFED v2.4r6 files for 2016/10/13]] | ||
+ | |Christoph Keller (NASA GMAO) | ||
+ | |Bug fix | ||
+ | |Update history: | ||
+ | #Implemented 19 Jul 2018 | ||
+ | #Implemented 11 Jul 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem_12#Fix_issues_in_complexSOA_and_complexSOA_SVPOA_restart_files|Fix issues in complexSOA and complexSOA_SVPOA restart files]] | ||
+ | |Eloise Marais (U. Birmingham) | ||
+ | |Bug fix | ||
+ | |Implemented 05 Jul 2018 | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem 12#Avoid double counting CEDS and GEIA NH3 biofuel emissions|Avoid double counting CEDS and GEIA NH3 biofuel emissions]] | ||
+ | |Christoph Keller (NASA/GMAO)<br>[[GCST]] | ||
+ | |Bug fix | ||
+ | |Implemented 29 Jun 2018 | ||
+ | *NOTE: The default HEMCO behavior is to prioritize CEDS over GEIA NH3 anthro and biofuel emissions. But to avoid confusion, we will explicitly request this behavior in the <tt>HEMCO_Config.rc</tt> files that ship with the [[GEOS-Chem Unit Tester]]. | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[Coupling_GEOS-Chem_with_RRTMG#Bug_fixes_for_RRTMG_in_v11-02|Bug fixes for v11-02 simulations using the RRTMG radiative transfer model]] | ||
+ | |Chris Holmes (U. Florida) | ||
+ | |Bug Fix | ||
+ | |Implemented 26 Jun 2018 | ||
+ | *NOTE: This only affects the RRTMG simulations. Full-chemistry difference-tests yielded identical results w/r/t v11-02-rc. | ||
+ | |||
+ | |-valign="top" | ||
+ | |Now always declare the netCDF time variable with an unlimited dimension | ||
+ | |Chris Holmes (U. Florida) | ||
+ | |Structural | ||
+ | |Implemented 26 Jun 2018 | ||
+ | *NOTE: Full-chemistry difference-tests yielded identical results w/r/t v11-02-rc. | ||
+ | |||
+ | |} | ||
== Current stable version == | == Current stable version == |
Revision as of 15:39, 7 August 2018
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
- 1 Version history summary
- 2 Current version in development
- 3 Current stable version
- 4 Previously released stable versions
- 5 New data directories
- 6 Validation
- 7 Previous issues now resolved in GEOS-Chem 12
- 7.1 Bug fix for HEMCO soil NOx error with ifort 17
- 7.2 FlexChem bug fix: do not zero ACTA, EOH, HCOOH
- 7.3 Avoid double counting CEDS and GEIA NH3 biofuel emissions
- 7.4 Fix bug in ARCTAS_SHIP_SO2 scale factors
- 7.5 Fix issues in complexSOA and complexSOA_SVPOA restart files
- 7.6 Bug fix for reference time in netCDF diagnostic files
- 7.7 Bug fix for CH4 boundary conditions
- 7.8 Fixed typo in call to GAMMA_HOBr_AER
- 7.9 Bug fix for liquid water content in gckpp_HetRates.F90
- 7.10 Add PH2O2 to the default list of prod/loss families in the KPP mechanisms
- 7.11 Fixed typo that prevented PM2.5 from being archived to netCDF output
- 7.12 Add entries to HEMCO_Config.rc for SOAP and pFe in CEDS inventory
- 7.13 Remove duplicate C3HI spectrum in FJX_spec.dat and FJX_j2j.dat
- 8 Outstanding issues not yet resolved in GEOS-Chem 12
Version history summary
The table below shows the previous, current, and successive major GEOS-Chem versions:
Previous major release | Current release candidate | Pending major release | Next stable version | Successive major release |
---|---|---|---|---|
GEOS-Chem v11-01-public | v11-02-release-candidate (aka v11-02-rc) |
GEOS-Chem 12.0.0 aka v11-02-final in the old version number system |
GEOS-Chem 12.1.0 | GEOS-Chem 13.0.0 |
RELEASED 01 Feb 2017 | RELEASED 22 Jun 2018 | RELEASE TBD | RELEASE TBD | RELEASE TBD 2019 |
GC v11-01 online manual | GC v11-02 benchmark history | GC 12 online manual | GC 12 benchmark history |
The sections below describe each GEOS-Chem feature version and GEOS-Chem patch version belonging to the GEOS-Chem 12 series. Each feature version carries the tag 12.Y.0, which indicates that the version will be (or has been) validated with a 1-month benchmark simulation (and one or more 1-year benchmark simulations if necessary). Each patch version carries the tag 12.Y.Z, which denotes that the version preserves backwards-compatibility with the prior feature version 12.Y.0. For more information, please see our GEOS-Chem version numbering system wiki page.
Several GEOS-Chem updates have been proposed by the relevant Working Groups, but have not yet been assigned to a specific GEOS-Chem version. For a complete list of these updates, please visit our GEOS-Chem model development priorities wiki page.
Current version in development
12.0.0
The following updates and fixes are slated to be added to GEOS-Chem 12.0.0 (aka v11-02-final). Many of these updates fix issues that were identified during the GEOS-Chem v11-02 release candidate testing period.
Please see the following link for complete information about the validation of GEOS-Chem 12.0.0:
- Approval form for 1-month benchmark simulation 12.0.0
- Results for 1-year benchmark simulation 12.0.0-Run0
Feature | Submitted by | Type | Status |
---|---|---|---|
Features affecting the full-chemistry (benchmark) simulation: | |||
Turn off TRASH emissions by default to avoid double counting | GCST Eloise Marais (U. Birmingham) |
Bug fix | Implemented 27 Jul 2018
|
Add entries to HEMCO_Config.rc for SOAP and pFe in CEDS inventory | GCST | Bug fix | Implemented 27 Jul 2018
|
Bug fix for liquid water content in gckpp_HetRates.F90 | Jiayue Huang (U. Washington) Quanjie Chen (U. Washington) |
Bug fix | Implemented 19 Jul 2018
|
Bug fix for CH4 boundary conditions | Lee Murray (Rochester) | Bug fix | Implemented 19 Jul 2018
|
Fixed typo in call to GAMMA_HOBr_AER | Xuan Wang GCST |
Bug fix | Implemented 18 Jul 2018
|
Correct issues when specifying DICE-Africa emissions: | Eloise Marais (U. Birmingham) GCST |
Bug Fix | Implemented 16 Jul 2018.
|
Add missing Criegee intermediate reaction | Xin Chen (UMN) | Bug fix | Implemented 11 Jul 2018, pending validation
|
Bug fix for ARCTAS_SHIP in HEMCO_Config.rc files | Barron Henderson (EPA) | Bug Fix | Implemented 03 Jul 2018
|
Bug fix in FlexChem: do not zero out ACTA, EOH, HCOOH concentrations | Katie Travis (MIT) | Bug Fix | Implemented 27 Jun 2018
|
Add GFED4.1 data for 2015 and 2016; Use GFED4.1 2016 in benchmarks | Kate O'Dell (CSU) GCST |
Science | Data files downloaded 27 Jun 2018
|
Features NOT affecting the full-chemistry simulation: | |||
Minor fixes for mercury chemistry updates implemented in v11-02f | Colin Thackray (Harvard) | Bug fix | Implemented 02 Aug 2018 |
Restore default timestep for MERRA-2 nested grid simulations to 10/20 | Ke Li (Harvard) | Bug fix | Implemented 30 Jul 2018 |
Remove duplicate C3HI spectrum in FJX_spec.dat and FJX_j2j.dat | Tomas Sherwen (York) | Bug fix | Implemented 27 Jul 2018
|
Fixed typo that prevented PM2.5 from being archived to netCDF output | Jonathan Moch (Harvard) GCST |
Bug fix | Implemented 24 Jul 2018 |
Bring two-way nested-grid code option up-to-date | Yingying Yan (PKU) Jintai Lin (PKU) |
Structural | Implemented 24 Jul 2018
|
Do not force L=1 by default in planeflight_mod.F | Katie Travis (MIT) GCST |
Bug fix | Implemented 23 Jul 2018 |
Add PH2O2 to the default list of prod/loss families in the KPP mechanisms | GCST | Structural | Implemented 19 Jul 2018
|
Update to HEMCO v2.1.007, which includes minor updates: | GCST Christoph Keller Jenny Fisher |
Structural & Bug fix |
Update history:
|
Bug fix for reference time in netCDF diagnostic files | GCST | Bug fix | Implemented 12 Jul 2018 |
Updates for QFED biomass burning emissions: | Christoph Keller (NASA GMAO) | Bug fix | Update history:
|
Fix issues in complexSOA and complexSOA_SVPOA restart files | Eloise Marais (U. Birmingham) | Bug fix | Implemented 05 Jul 2018 |
Avoid double counting CEDS and GEIA NH3 biofuel emissions | Christoph Keller (NASA/GMAO) GCST |
Bug fix | Implemented 29 Jun 2018
|
Bug fixes for v11-02 simulations using the RRTMG radiative transfer model | Chris Holmes (U. Florida) | Bug Fix | Implemented 26 Jun 2018
|
Now always declare the netCDF time variable with an unlimited dimension | Chris Holmes (U. Florida) | Structural | Implemented 26 Jun 2018
|
Current stable version
The GEOS-Chem version listed below is considered stable and can be used to start a new research project.
Version TBD
You can use the following Git commands to download the GEOS-Chem source code and GEOS-Chem Unit Tester package.
# Download the GEOS-Chem source code git clone https://github.com/geoschem/geos-chem Code.X.Y.Z cd Code.X.Y.Z git checkout master cd .. # Download the GEOS-Chem Unit Tester, from which run directories can be generated git clone https://github.com/geoschem/geos-chem-unittest UT cd UT git checkout master
For more information about getting set up with GEOS-Chem, please see the GEOS-Chem Online User's Guide.
Previously released stable versions
TBD
New data directories
The following HEMCO data directories have been added or updated for GEOS-Chem 12. You will have to download the directories relevant to your simulation.
Inventory | Type | What was added? | When added? | Directory |
---|---|---|---|---|
QFED biomass burning | OPTIONAL | New data files for QFED v2.5r1 | 12.0.0 | HEMCO/QFED/v2018-07 |
GFED4 biomass inventory | DEFAULT | GFED4.1 data files for 2015 and 2016 have been added | 12.0.0 | HEMCO/GFED4/v2015-10/2015 HEMCO/GFED4/v2015-10/2016 |
NOTES:
- DEFAULT denotes that the inventory is part of the default HEMCO emissions configuration that you get when you run GEOS-Chem "out-of-the-box".
- OPTIONAL denotes optional inventories that you may use (or not) according to your research needs.
--Bob Yantosca (talk) 20:02, 29 June 2018 (UTC)
Validation
In this section we provide information about the benchmarks and tests that we have done to validate GEOS-Chem v11-03.
1-month and 1-year benchmarks
For complete information about the benchmark simulations used to validate GEOS-Chem v11-03, please see our GEOS-Chem 12 benchmark history wiki page.
Unit tests for GEOS-Chem 12
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-03. Click on the links below to see the results from each set of unit tests:
Unit test | Date |
---|---|
TBD |
Previous issues now resolved in GEOS-Chem 12
The following bugs and/or technical issues have been resolved in GEOS-Chem 12.
Bug fix for HEMCO soil NOx error with ifort 17
This fix will be included in GEOS-Chem 12.0.0.
For complete information about this issue and how it was resolved, please see this post on our Implementation of HEMCO in GEOS-Chem wiki page..
--Bob Yantosca (talk) 13:54, 27 June 2018 (UTC)
FlexChem bug fix: do not zero ACTA, EOH, HCOOH
This fix will be included in GEOS-Chem 12.0.0.
For complete information about this bug and how it was fixed, please see this post on our FlexChem wiki page.
--Bob Yantosca (talk) 14:42, 27 June 2018 (UTC)
Avoid double counting CEDS and GEIA NH3 biofuel emissions
This fix will be included GEOS-Chem 12.0.0.
Eloise Marais wrote:
If CEDS includes NH3 emissions, then this will be double counting. I recommend using CEDS, as GEIA is very out of date (1998 according to the HEMCO configuration file I have).
Daniel Jacob replied:
We should use the GEIA NH3 natural emission inventory (1° x 1° file from Eleanor Morris) but not the GEIA NH3 biofuel inventory.
UPDATE 29 Jun 2018:
We discovered that HEMCO's default behavior is to ignore the GEIA NH3 anthro and biofuel emissions if CEDS is turned on. This is shown in the HEMCO.log file when VERBOSE = 3:
Register_Base: Checking GEIA_NH3_ANTH HEMCO: Entering Get_targetID (hco_config_mod.F90) ( 4) Skip container GEIA_NH3_ANTH because of CEDS_NH3_WST . . . Register_Base: Checking GEIA_NH3_BIOF HEMCO: Entering Get_targetID (hco_config_mod.F90) ( 4) Skip container GEIA_NH3_BIOF because of CEDS_NH3_WST_Cat00002
But for the sake of clarity, we will modify the HEMCO_Config.rc files to explicitly request this behavior:
(((.not.CEDS.or.EDGARv43 # GEIA NH3 anthro (skip if CEDS or EDGARv43 is on) 0 GEIA_NH3_ANTH $ROOT/NH3/v2014-07/NH3_geos.4x5.nc NH3_ANTH 1990/1-12/1/0 C xy kg/m2/s NH3 - 1 1 # GEIA NH3 biofuel (skip if CEDS or EDGARv43 is on) (((BIOFUEL 0 GEIA_NH3_BIOF $ROOT/NH3/v2014-07/NH3_biofuel.geos.1x1.nc NH3_BIOFUEL 1998/1-12/1/0 C xy kg/m2/s NH3 - 2 1 )))BIOFUEL ))).not.CEDS.or.EDGARv43
--Bob Yantosca (talk) 21:02, 29 June 2018 (UTC)
Fix bug in ARCTAS_SHIP_SO2 scale factors
This fix will be included GEOS-Chem 12.0.0.
For a complete description of this issue and how it was resolved, please see this post on our Ship Emissions wiki page.
--Bob Yantosca (talk) 15:54, 3 July 2018 (UTC)
Fix issues in complexSOA and complexSOA_SVPOA restart files
This fix will be included GEOS-Chem 12.0.0.
Eloise Marais wrote:
I just wanted to draw your attention to a few issues with the restart files provided on the ftp site (http://ftp.as.harvard.edu/gcgrid/data/ExtData/SPC_RESTARTS/). The complexSOA and complexSOA_SVPOA restart files (4° x 5° and 2° x 2.5°) both have a species entry for IEPOX and also IEPOXD, IEPOXB, and IEPOXA, but should only have the latter 3. This is the same for RIP. The number of species in the 4° x 5° and 2° x 2.5° files for complexSOA is also different.
To fix these issues, we have generated new restart files for the complexSOA and complexSOA_SVPOA simulations. You can download these to your SPC_RESTARTS folder with the following commands:
wget "ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/SPC_RESTARTS/initial_GEOSChem_rst.2x25_complexSOA.nc" wget "ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/SPC_RESTARTS/initial_GEOSChem_rst.2x25_complexSOA_SVPOA.nc" wget "ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/SPC_RESTARTS/initial_GEOSChem_rst.4x5_complexSOA.nc" wget "ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/SPC_RESTARTS/initial_GEOSChem_rst.4x5_complexSOA_SVPOA.nc"
--Bob Yantosca (talk) 18:01, 5 July 2018 (UTC)
Bug fix for reference time in netCDF diagnostic files
This fix will be included GEOS-Chem 12.0.0.
Maria Zatko wrote:
For the sake of run time, we would like to run GEOS-Chem with the old (15/30 minute timesteps) instead of the recommended timesteps (10/20 minutes). I noticed that the output netcdf files have different timestamps depending on the timestep (e.g., 10/20 minutes, 15/30 minutes) used. For example, when 10/20 minute timesteps are used, the time units in the following netcdf output files range from 0 to 23 UTC:
GEOSChem3hrInst.SpeciesConcWithMet.20150701_0000z.nc4 (file size: 4.6 GB) GEOSChem3hrInst.SpeciesConcWithMet.20150702_0000z.nc4 (file size: 5.3 GB) time:units = "minutes since 2015-07-01 00:00:00 UTC"
When 15/30 minute timesteps are used, the time units in the following output files include hour 24:
GEOSChem3hrInst.SpeciesConcWithMet.20150701_0000z.nc4 (file size: 4.6 GB) GEOSChem3hrInst.SpeciesConcWithMet.20150701_2400z.nc4 (file size: 5.3 GB) GEOSChem3hrInst.SpeciesConcWithMet.20150702_2400z.nc4 (file size: 5.3 GB) time:units = "minutes since 2015-07-01 24:00:00 UTC"
I wanted to check in to see if others have run into this issue as well.
Bob Yantosca replied:
This issue is caused by floating-point roundoff error in converting from the Julian date to the calendar date. I have created a more robust fix for this issue, which now computes the correct netCDF file reference time independent of the model "heartbeat" timestep.To apply the fix:
The fix will be standardized in the GEOS-Chem 12.0.0 release.
- Copy the patch file ftp://ftp.as.harvard.edu/pub/exchange/bmy/nc_ref_time_fix.txt to your computational server.
- Change to your GEOS-Chem v11-02 release candidate code directory and check out the "master" branch.
- Type: git am < nc_ref_time_fix.txt
--Bob Yantosca (talk) 15:10, 12 July 2018 (UTC)
Bug fix for CH4 boundary conditions
This fix will be included GEOS-Chem 12.0.0.
Lee Murray wrote:
I’ve been testing the release candidate, and believe that we need to change the default input.geos options for the Standard and Benchmark simulations in the unit tester to :
%%% EMISSIONS MENU %%% :
Turn on emissions? : T
HEMCO Input file : HEMCO_Config.rc
=> 1ppt MBL BRO Sim.? : F
Switches for UCX :---
=> Use CH4 emissions? : F
=> Turn on surface BCs :---
=> CH4? : F
=> OCS? : T
=> CFCs? : T
=> Cl species? : T
=> Br species? : F
=> N2O? : T
=> Set initial glob MRs:---
=> strat. H2O? : T
=> CFC emission year : 0
By default in the unit tester, => CH4 is still enabled. This appears to set up conflict between the older UCX surface layer and the newer HEMCO methane PBL boundary conditions, causing air to be a mix of preindustrial and present-day values (see below for Jul 1, 2013 surface methane in a benchmark run). Setting it to F and/or removing LCH4SBC causes the methane to behave as expected, i.e., being set by HEMCO in all simulations, with the expected spatial pattern. Also, we can slate Online CH4 chemistry? : T in the Chemistry Menu and Input_Opt%LCH4CHEM for future removal, as it no longer does anything.
--Melissa Sulprizio (talk) 18:16, 16 July 2018 (UTC)
Fixed typo in call to GAMMA_HOBr_AER
This fix will be included GEOS-Chem 12.0.0.
Xuan Wang wrote:
When looking back at the standard GC codes of v11-02d, I found a bug in the function HETHOBr_SS_JS (in module KPP/Standard/gckpp_HetRates.F90). When calling the gamma calculation, in v11-02d, around line 3514, we have:
CALL Gamma_HOBr_AER(rAer, denAir, 2, TK, clConc, brConc, &
hConc, GAM_HOBr, r_gp)
This will cause both HOBr+Cl- (X=1) and HOBr+Br- (X=2) on sea salt use the gamma of HOBr + Br- (X=2). It should be:
CALL Gamma_HOBr_AER(rAer, denAir, X, TK, clConc, brConc, &
hConc, GAM_HOBr, r_gp)
It looks like a typo which could lead to quite different result for halogen chemistry.
--Bob Yantosca (talk) 20:26, 17 July 2018 (UTC)
Bug fix for liquid water content in gckpp_HetRates.F90
This fix will be included GEOS-Chem 12.0.0.
Jiayue Huang wrote:
Qianjie Chen and I think that there is potential a bug in KPP/Standard/gckpp_HetRates.F90 for halogen chemistry in the model v11-02d/e/f.In the original code, we were using the grid-box average liquid water content, which likely underestimates the in-cloud Br- and Cl- concentration. So we think that we should change it to in-cloud liquid water content by dividing the cloud fraction. This results in increases the S(IV)+Br reaction in cloud, and thus reduces the BrO abundance, mainly in the 30-90N region.
This bug fix is to use cloud fraction when calculating liquid water content [cm3 H2O/cm3 air] in routine GET_HALIDE_CLDConc, which is used to calculate the in-cloud concentrations of Br- and Cl-.
--Bob Yantosca (talk) 15:48, 20 July 2018 (UTC)
Add PH2O2 to the default list of prod/loss families in the KPP mechanisms
This fix will be included GEOS-Chem 12.0.0.
The P(H2O2) field is required as input in the aerosol-only simulations. To simplify saving out new oxidant fields, we now include P(H2O2) by default in the prod/loss families so that it is saved out when users turn on the prod/loss diagnostic in either bpch or netCDF format.
--Bob Yantosca (talk) 16:16, 20 July 2018 (UTC)
Fixed typo that prevented PM2.5 from being archived to netCDF output
This fix will be included GEOS-Chem 12.0.0.
In routine Set_AerMass_Diagnostic of module GeosCore/diagnostics_mod.F90, there was a typo in the IF block where the PM2.5 diagnostic values get copied into State_Diag%PM25. The code in RED below:
IF ( Archive_BetaNO ) THEN
State_Diag%PM25(I,J,L) = ...etc...
needs to be replaced by the code in GREEN:
IF ( Archive_PM25 ) THEN
State_Diag%PM25(I,J,L) = ...etc...
--Bob Yantosca (talk) 20:33, 24 July 2018 (UTC)
Add entries to HEMCO_Config.rc for SOAP and pFe in CEDS inventory
This fix will be included GEOS-Chem 12.0.0.
The CEDS emissions entries in HEMCO_Config.rc were missing for species SOAP and pFe. To fix this, add the following lines in GREEN:
0 CEDS_CO_AGR $ROOT/CEDS/v2018-04/CO-em-anthro_CMIP_CEDS_195001-201412.nc CO_agr 1950-2014/1-12/1/0 C xy kg/m2/s CO 26 1/2/12 5 0 CEDS_SOAP_AGR - - - - - - SOAP 26/280 1/2/12 5 0 CEDS_CO_ENE $ROOT/CEDS/v2018-04/CO-em-anthro_CMIP_CEDS_195001-201412.nc CO_ene 1950-2014/1-12/1/0 C xy kg/m2/s CO 26 1/2/12 5 0 CEDS_SOAP_ENE - - - - - - SOAP 26/280 1/2/12 5 0 CEDS_CO_IND $ROOT/CEDS/v2018-04/CO-em-anthro_CMIP_CEDS_195001-201412.nc CO_ind 1950-2014/1-12/1/0 C xy kg/m2/s CO 26 1/2/12 5 0 CEDS_SOAP_IND - - - - - - SOAP 26/280 1/2/12 5 0 CEDS_CO_TRA $ROOT/CEDS/v2018-04/CO-em-anthro_CMIP_CEDS_195001-201412.nc CO_tra 1950-2014/1-12/1/0 C xy kg/m2/s CO 26 1/2/12 5 0 CEDS_SOAP_TRA - - - - - - SOAP 26/280 1/2/12 5 0 CEDS_CO_SLV $ROOT/CEDS/v2018-04/CO-em-anthro_CMIP_CEDS_195001-201412.nc CO_slv 1950-2014/1-12/1/0 C xy kg/m2/s CO 26 1/2/12 5 0 CEDS_SOAP_SLV - - - - - - SOAP 26/280 1/2/12 5 0 CEDS_CO_WST $ROOT/CEDS/v2018-04/CO-em-anthro_CMIP_CEDS_195001-201412.nc CO_wst 1950-2014/1-12/1/0 C xy kg/m2/s CO 26 1/2/12 5 0 CEDS_SOAP_WST - - - - - - SOAP 26/280 1/2/12 5 0 CEDS_SO2_AGR $ROOT/CEDS/v2018-04/SO2-em-anthro_CMIP_CEDS_195001-201412.nc SO2_agr 1950-2014/1-12/1/0 C xy kg/m2/s SO2 - 1/2/12 5 0 CEDS_SO4_AGR - - - - - - SO4 63 1/2/12 5 0 CEDS_pFe_AGR - - - - - - pFe 66 1/2/12 5 0 CEDS_SO2_ENE $ROOT/CEDS/v2018-04/SO2-em-anthro_CMIP_CEDS_195001-201412.nc SO2_ene 1950-2014/1-12/1/0 C xy kg/m2/s SO2 - 1/2/12 5 0 CEDS_SO4_ENE - - - - - - SO4 63 1/2/12 5 0 CEDS_pFe_ENE - - - - - - pFe 66 1/2/12 5 0 CEDS_SO2_IND $ROOT/CEDS/v2018-04/SO2-em-anthro_CMIP_CEDS_195001-201412.nc SO2_ind 1950-2014/1-12/1/0 C xy kg/m2/s SO2 - 1/2/12 5 0 CEDS_SO4_IND - - - - - - SO4 63 1/2/12 5 0 CEDS_pFe_IND - - - - - - pFe 66 1/2/12 5 0 CEDS_SO2_TRA $ROOT/CEDS/v2018-04/SO2-em-anthro_CMIP_CEDS_195001-201412.nc SO2_tra 1950-2014/1-12/1/0 C xy kg/m2/s SO2 - 1/2/12 5 0 CEDS_SO4_TRA - - - - - - SO4 63 1/2/12 5 0 CEDS_pFe_TRA - - - - - - pFe 66 1/2/12 5 0 CEDS_SO2_SLV $ROOT/CEDS/v2018-04/SO2-em-anthro_CMIP_CEDS_195001-201412.nc SO2_slv 1950-2014/1-12/1/0 C xy kg/m2/s SO2 - 1/2/12 5 0 CEDS_SO4_SLV - - - - - - SO4 63 1/2/12 5 0 CEDS_pFe_SLV - - - - - - pFe 66 1/2/12 5 0 CEDS_SO2_WST $ROOT/CEDS/v2018-04/SO2-em-anthro_CMIP_CEDS_195001-201412.nc SO2_wst 1950-2014/1-12/1/0 C xy kg/m2/s SO2 - 1/2/12 5 0 CEDS_SO4_WST - - - - - - SO4 63 1/2/12 5 0 CEDS_pFe_WST - - - - - - pFe 66 1/2/12 5
--Melissa Sulprizio (talk) 16:07, 27 July 2018 (UTC)
Remove duplicate C3HI spectrum in FJX_spec.dat and FJX_j2j.dat
This fix will be included GEOS-Chem 12.0.0.
Tomas Sherwen wrote:
- I've just found that both FJX_spec.dat the redundant CH3I spectrum and the new spectrum are both present in FJX_spec.dat and FJX_j2j.dat (in v12.0.0). This will have been due to the spectrum in my branch being under the species alias of "CH3IT" to preserve (and ignore) the use of "CH3I" in the namespace of the code for the existing methyl iodide simulation when the iodine simulation was developed. The "CH3IT" variables were replaced "CH3I" in the v11-1-Patches branches with halogen chemistry provided to GCST, but the existing methyl iodide simulation variable was not removed. There are no double ups in the KPP mechanism itself.
- Please remove lines 276-281 in FJX_spec.dat (copied below).
" CH3I 243 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 J10 3.958E-20 4.769E-20 8.257E-19 2.154E-19 7.349E-20 3.351E-20 1.034E-20 4.435E-21 2.272E-21 4.393E-22 4.052E-24 0.000E+00 CH3I 300 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 J10 3.958E-20 4.769E-20 8.376E-19 2.438E-19 9.355E-20 4.680E-20 1.529E-20 6.532E-21 3.328E-21 6.783E-22 1.038E-23 0.000E+00 "
- Please also remove the line 124 in FJX_j2j.dat (copied below).
" 123 CH3I PHOTON I PRODUCTS 1.000 /CH3I / "
--Melissa Sulprizio (talk) 17:44, 27 July 2018 (UTC)
Outstanding issues not yet resolved in GEOS-Chem 12
RRTMG does not compile with GNU Fortran
For a complete description of this issue, please see this post on our Coupling GEOS-Chem with RRTMG wiki page.
--Bob Yantosca (talk) 18:55, 3 August 2018 (UTC)
TOMAS40 simulations exit with a floating-point exception
We were unit testing the TOMAS40 code in the GEOS-Chem v11-02-release-candidate (soon to be finalized as 12.0.0) with the debugging flags:
make –j8 DEBUG=y BOUNDS=y FPEX=y NO_ISO=y ut &
which revealed a a NaN error. This was thrown at line 903 of GeosCore/tomas_mod.F:
ERR_VAR = 'Gcout'
ERR_MSG = 'After COND_NUC'
! check for NaN and Inf (win, 10/4/08)
do jc = 1, icomp-1
ERR_IND(1) = I
ERR_IND(2) = J
ERR_IND(3) = L
ERR_IND(4) = 0
call check_value( Gcout(jc), ERR_IND, ERR_VAR, ERR_MSG )
! if( IT_IS_FINITE(Gcout(jc))) then
! print *,'xxxxxxxxx Found Inf in Gcout xxxxxxxxxxxxxx'
! print *,'Location ',I,J,L, 'comp',jc
! call debugprint( Nkout, Mkout, i,j,l,'After COND_NUC')
! stop
! endif
enddo
This error might be caused by the several chemistry and emissions updates that have been added to GEOS-Chem. The GEOS-Chem Support Team has referred the error to the TOMAS development team, who will be investigating further.
--Bob Yantosca (talk) 16:56, 6 August 2018 (UTC)