GEOS-Chem v11-02 benchmark history

From Geos-chem
Jump to: navigation, search


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.




On this page we have posted complete information about all benchmark simulations (both 1-month and 1-year) for GEOS-Chem v11-02.

1-year full-chemistry benchmarks

v11-02-final (aka 12.0.0)

Please see visit this page for details on 1-year benchmark simulation 12.0.0-Run0.

v11-02f-Run1

This 1-year benchmark simulation was approved by the GEOS-Chem Steering Committee on 11 Jun 2018.

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02f-Run1 GEOS-FP,
72L, 4x5
2016 Updates introduced in this benchmark: 11.532
Green v11-02e-Run1 GEOS-FP,
72L, 4x5
2016 Updates introduced in this benchmark:
  • Use GEOS-FP 4x5 met fields for 2016
11.541
Red v11-02e-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02e: 11.686

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02f/Run1/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02f/Run1/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

--Bob Yantosca (talk) 21:00, 7 June 2018 (UTC)

v11-02f-Run0

1-year benchmark simulation v11-02f-Run0 contained all of the same updates described in v11-02f-Run1 above, PLUS

Subsequent analysis revealed that the switch from the 4° x 5° GEIA NH3 emissions data to the 1° x 1° data produced 20X too much NH3 emissions as expected. We believe that this is likely attributed to an regridding error during the preparation of the original file.

We have since reverted to using the 1° x 1° GEIA natural NH3 emissions data, which is benchmarked with v11-02f-Run1. We have also de-listed the v11-02f-Run0 benchmark from the wiki, given that it contained clearly incorrect emissions.

--Bob Yantosca (talk) 21:11, 7 June 2018 (UTC)

v11-02e-Run1

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02e-Run1 GEOS-FP,
72L, 4x5
2016 Updates introduced in this benchmark:
  • Use GEOS-FP 4x5 met fields for 2016
11.541
Green v11-02e-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02e: 11.686
Red v11-02d-Run2 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02d: 11.857

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/Run1/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/Run1/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

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

v11-02e-Run0

This 1-year benchmark simulation was approved by the GEOS-Chem Steering Committee on 24 Mar 2018.

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02e-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02e: 11.686
Green v11-02d-Run2 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02d: 11.857
Red v11-02c-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02c: 11.596

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/Run0/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/Run0/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

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

v11-02d-Run2

This 1-year benchmark simulation was approved by the GEOS-Chem Steering Committee on 12 Feb 2018.

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02d-Run2 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02d: 11.857
Green v11-02c-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02c: 11.596
Red v11-02a-Run1 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02a: 11.829

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02d/Run2/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02d/Run2/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

--Melissa Sulprizio (talk) 21:18, 8 February 2018 (UTC)

v11-02d-Run1

This 1-year benchmark simulation was not approved by the GEOS-Chem Steering Committee due to concerns about high P(Ox) and L(Ox) values in the O3 budget table. Please see v11-02d-Run2 for the approved 1-year benchmark for this version.

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02d-Run1 GEOS-FP,
72L, 4x5
2013

Updates introduced in v11-02d:

11.979
Green v11-02c-Run0 GEOS-FP,
72L, 4x5
2013

Updates introduced in v11-02c:

11.596
Red v11-02a-Run1 GEOS-FP,
72L, 4x5
2013

Updates introduced in v11-02a:

11.829

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02d/Run1/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02d/Run1/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

--Melissa Sulprizio (talk) 19:56, 8 January 2018 (UTC)

v11-02d-Run0

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02d-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02d: 11.662
Green v11-02c-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02c: 11.596
Red v11-02a-Run1 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02a: 11.829

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02d/Run0/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02d/Run0/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

--Melissa Sulprizio (talk) 16:34, 13 December 2017 (UTC)

v11-02c-Run0

This 1-year benchmark simulation was approved by the GEOS-Chem Steering Committee on 21 Sep 2017.

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02c-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02c: 11.596
Green v11-02a-Run1 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02a: 11.829
Red v11-01-public-release-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-01 public release: 12.001

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02c/Run0/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02c/Run0/output/

New plots were added this version for evaluating the two SOA schemes, including:

Please also view the following pages comparing this version to past 1-year benchmarks:

--Melissa Sulprizio (talk) 16:52, 15 September 2017 (UTC)

Comments about the 1-year benchmark v11-02c-Run0

Computing PM2.5 and AOD

Jeff Pierce wrote:

Can you tell us which SOA is being included into the PM2.5 and AOD calculations for the three different simulations? I assume the SOA scheme changes between these different simulations?

Melissa Sulprizio responded:

Following the recommendation of the Aerosol WG leaders, simple SOA is the default in all full-chemistry simulations for v11-02c and later. However, for the benchmark simulations we turn on both complex SOA and simple SOA so that we can continually evaluate the two schemes. That means, when a user downloads v11-02 and creates a run directory for the full-chemistry simulations (tropchem, standard w/ UCX, etc.), input.geos will automatically be set up to use simple SOA. The complex SOA option will only be turned on if a user creates a run directory for the “complexSOA”, “complexSOA_SVPOA”, or “benchmark” simulation. As I mentioned below, by default in aerosol_mod.F, when complex SOA is turned on the complex SOA species will be used in the calculation of PM2.5 and AOD over the simple SOA species to avoid double counting. This is only an issue in the benchmark simulations when both options are turned on – in all other simulations, it will be clear whether complex SOA or simple SOA species are used in the calculations based on the settings in input.geos.
Ocean ALD2

Dylan Millet wrote:

I missed it in the prior benchmark but it seems as part of the PAN update the ocean acetaldehyde source went into the model. That seems fine to me but it is a quite large global reactivity source (can someone tell me if the 66TgC for ALD2osr in the emission summary file is gross or net? I hope gross) so there needs to be traceability for that in the model narrative (e.g. if it's based on Millet et al. 2010 or on something else).

Melissa Sulprizio responded:

The ALD2 ocean source (ALD2osr) in the emissions summary file is gross. HEMCO handles the ALD2 ocean sink and we can’t currently output that as a diagnostic. This is similar to how we currently handle the ACET ocean source/sink. Here are comments from Christoph added during the implementation from HEMCO:
    !  (1 ) in HEMCO, Acetone sink from OCEAN is handled by drydep and the
    !       corresponding diagnostics. If needed, we can write a wrapper that
    !       explicitly calculates the Acetone ocean sink. (ckeller, 08/04/14)
Halogens

Dylan Millet wrote:

There are some large relative changes for many of the halogen species, which I guess is linked to OH/NOx feedbacks from the isoprene chemistry. It would be great if someone close to the halogen stuff can confirm that they're comfortable with the changes.

Mat Evans wrote:

So as Dylan says the halogens change significantly with these updated. I think most of these changes relate to the increases in the oceanic ALD2 source. ALD2 reacts pretty quickly with Br resulting in HBr. So BrOx concentrations tend to decrease and HBr increases.
This is useful as we currently don’t include the sea-salt debromination source. It seems likely that the additional oceanic oVOCs leads to a shorter BrOx lifetime and so might allow for an increase in the Br source consistent with the observational BrO constraint we’ll have to see in future releases.

--Melissa Sulprizio (talk) 14:53, 25 September 2017 (UTC)

v11-02a-Run1

This 1-year benchmark simulation was an unofficial benchmark to validate removal partitioning of NOx biomass burning emissions directly to PAN and HNO3 from v11-02a.

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02a-Run1 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02a: 11.829
Green v11-01-public-release-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-01 public release: 12.001
Red v11-01k-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-01h:

Updates introduced in v11-01i:

Updates introduced in v11-01j:

Updates introduced in v11-01k:

12.004

The output plots for this run may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02a/Run1/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02a/Run1/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

--Melissa Sulprizio (talk) 16:37, 1 June 2017 (UTC)

v11-02a-Run0

This 1-year benchmark simulation was approved by the GEOS-Chem Steering Committee on 12 May 2017, with the recommendation to remove partitioning of NOx biomass burning emissions directly to PAN and HNO3.

Three GEOS-Chem model versions were compared to each other:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Black Observations        
Blue v11-02a-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-02a: 11.750
Green v11-01-public-release-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-01 public release: 12.001
Red v11-01k-Run0 GEOS-FP,
72L, 4x5
2013 Updates introduced in v11-01h:

Updates introduced in v11-01i:

Updates introduced in v11-01j:

Updates introduced in v11-01k:

12.004

The output plots for Run0 may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02a/Run0/output
mget *

You may also view the PDF files online by pointing your browser to

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02a/Run0/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

--Melissa Sulprizio (talk) 18:14, 10 May 2017 (UTC)

1-month benchmarks

v11-02-final (aka 12.0.0)

Please see visit this page for the assessment form for the 1-month benchmark simulation of 12.0.0.

v11-02f

This section includes the assessment form for the 1-month benchmark simulation of v11-02f.

Description
New features added into GEOS-Chem
Feature Submitted by
Features affecting the full-chemistry simulation:
Update HEMCO from v2.1.001 to v2.1.005 Christoph Keller (NASA GMAO)
Use global 0.25x0.25 files for MEGAN input and MODIS LAI Jenny Fisher (U. Wollongong)
Eloise Marais (U. Birmingham)
Dylan Millet (UMN)
Update lightning BETA values and OTD/LIS files for consistency with latest met fields Lee Murray (Rochester)
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)
EDGAR v4.3.1 emissions
  • This inventory will be included as an option
Chi Li (Dalhousie)
Use 1x1 GEIA NH3 emissions Eleanor Morris (York)
Historical Canadian emissions Chi Li (Dalhousie)
Global ethane emissions Zitely Tzompa-Sosa (CSU)
Inefficient combustion emissions in Africa (DICE-Africa inventory) Eloise Marais (U. Birmingham)
Christine Wiedinmyer (NCAR)
Trash burning emissions Eloise Marais (U. Birmingham)
Christine Wiedinmyer (NCAR)
Update GFED4 emission factors Katie Travis (MIT)
Retire GFED3 biomass burning emissions Katie Travis (MIT)
Ammonia emissions from Arctic seabirds Betty Croft (Dalhousie)
Fixes for prod/loss rates used in linearized chemistry routines: Chris Holmes (FSU)
Christoph Keller (NASA GMAO)
Add fix for UCX mesospheric NOy routine Sebastian Eastham (MIT)
Fix Henry's Law parameters for HCl Xuan Wang (Harvard)
Fix for N2O5 uptake Xuan Wang (Harvard)
Add fix for organic nitrate chemistry Xuan Wang (Harvard)
Features not affecting the full-chemistry simulation:
Updated CH4 simulation Bram Maasakkers (Harvard)
Melissa Sulprizio (Harvard)
Updated Tagged CO simulation Jenny Fisher (Wollongong)
Updates to Hg chemistry from Horowitz et al. (2017) Hannah Horowitz (U. Washington)
Colin Thackray (Harvard)
Rebecca Stern (Harvard)
Change default timestep for MERRA-2 nested grid simulations Bram Maasakkers (Harvard)
Add ND49 timeseries diagnostic for correcting ozone from the height of the lowest model level to 10m Katie Travis (MIT)
Fix bug in biomass burning emissions when GFED is turned off Corey Trujillo (U. Colorado)
Therese Carter (MIT)
Fix missing HCOOH factors for RCP in HEMCO_Config.rc Seb Eastham (MIT)
Fix bug in calculation of L(Ox) for tagged Ox simulations Lu Shen (Harvard)
Change tolerance value as a temporary fix for convergence errors in 2x2.5 tropchem simulations Katie Travis (MIT)
Chemistry WG
Fix compatibility issue with ifort17 Bob Yantosca (GCST)
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2016
1-month benchmark finished on: 2018-05-08 06:17:34
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 "Haswell" @ 2.50 GHz)
  • Wall time: 06:25:56
  • CPU time / wall time: 22.2501
  • % of ideal performance: 92.71%
  • Memory: 6.737 GB
Compared to previous benchmark: v11-02e using 2016 GEOS-FP met fields
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify): Several new emissions inventories added
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02f/output/v11-02f.results.html
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02f/output/
Metrics
Global mean OH from log file (% change): 12.3128742319025 x 105 molec/cm3 (-1.18%)
Methyl chloroform lifetime (% change): 5.0338 years (+1.28%)
Methane lifetime (% change): 8.4770 years (+1.23%)
O3 STE flux (% change): 496.2694 Tg O3/year (-0.522%)
At the SURFACE, list all species that changed by 10% or more: The updates to emissions and chemistry impact several species:
  • NO, O3, PAN, CO, ALK4, ISOP (and species that react with it), HNO3, ACET, H2O2, MEK, ALD2, RCHO, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, HNO4, MP, DMS, SO2, SO4, SO4s, NH3, NH4, NIT, NITs, BCPI, BCPO, OCPI, OCPO, BrO, HOBr, HBr, BrNO2, BrNO3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, ISOA*, ASOA*, TSOA*, SOAP, SOAS, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOX, IEPOX*, ISN1, RIP*, , IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1O*, MONITS, MONITU, HONIT, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2O2, HOI, I2, IBr, ICl, I, IO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI, OH

Species that are largely unaffected:

  • Dust and sea-salt species (these are not affected by the emissions updates)
  • The remaining Br, Cl, I species not listed above
  • CFC species
  • N2O, OCS, CH4
  • IONITA
Comments on SURFACE differences:
  • Most of these differences can be attributed to the new emissions inventories that were added (C2H6, DICE-Africa, Trash burning, etc.)
  • Large differences in NH3 can be attributed to using the GEIA natural NH3 emissions file at 1° x 1° resolution, plus some new emission inventories (such as Arctic Seabirds).
  • N2O5 uptake was modified in this version leading to an increase in concentration over South America but a decrease in mostly N mid-latitudes.
  • Some of the changes in the NOy species could be attributed to the fix that was made in UCX for mesospheric NOy.
  • Differences in HCl can be attributed to the fix of an incorrect Henry's law parameter.
  • Species for which there are biogenic emissions are affected by the switch to the 0.25° x 0.25° MODIS LAI data. This might also affect dry deposition.
  • NOTE: Some large ratios (e.g. of isoprene or SOA species over oceans) might be due to small-number statistics, where concentrations are already very low.
At 500 hPa, list all species that changed by 10% or more: Mostly the same species as listed in SURFACE DIFFERENCES
Comments on 500 hPa differences:
  • See comments for surface differences
  • CO is more or less unchanged at 500hPa
In the ZONAL MEAN differences, list all species that changed by 10% or more: Mostly the same species as listed above in SURFACE differences
Comments on ZONAL MEAN differences: See comments on SURFACE differences
In the EMISSION RATIO maps, list all species that changed by 10% or more:
  • All emitted species except CH2Br2
Comments on EMISSION RATIO differences:
  • Differences are largely caused by the introduction of new inventories as listed in the table above.
  • Species for which there are biogenic emissions are affected by the switch to the 0.25° x 0.25° MODIS LAI data.
Additional or summary comments:
  • While the plots show large regional differences, the overall metrics w/r/t the v11-02e_2016 run did not change very much. Mean OH decreased by 1.18% and MCF lifetime increased by about 1.2%.
Approval
Requires further investigation: No
Approved by: Daniel Jacob
Date of approval: 10 May 2018

--Bob Yantosca (talk) 21:32, 8 May 2018 (UTC)

v11-02e using 2016 GEOS-FP met fields

This section includes the assessment form for the 1-month benchmark simulation of v11-02e.

Description
New features added into GEOS-Chem
Feature Submitted by
Features affecting the full-chemistry simulation:
Use GEOS-FP 4x5 met fields for July 2016 GCST
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2016
1-month benchmark finished on: TBD
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 "Haswell" @ 2.50 GHz)
  • Wall time: 06:08:11
  • CPU time / wall time: 22.1148
  • % of ideal performance: 92.15%
  • Memory: 6.9433 GB
Compared to previous benchmark: v11-02e using spun up restart file
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify):
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02e_2016.tar.gz
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02e_2016.tar.gz
Metrics
Global mean OH from log file (% change): 12.4599571735462 x 105 molec/cm3 (-0.19%)
Methyl chloroform lifetime (% change): 4.9701 years (-0.73%)
Methane lifetime (% change): 8.3737 years (-0.81%)
O3 STE flux (% change): 498.8760 Tg O3/year (+8.12%)
At the SURFACE, list all species that changed by 10% or more: All species
Comments on SURFACE differences: The difference in met field year (2016 vs. 2013) affects all operations (transport, convection, PBL mixing, chemistry, drydep, wetdep).
At 500 hPa, list all species that changed by 10% or more: All species
Comments on 500 hPa differences: See comments on surface differences
In the ZONAL MEAN differences, list all species that changed by 10% or more: All species
Comments on ZONAL MEAN differences: See comments on surface differences
In the EMISSION RATIO maps, list all species that changed by 10% or more:
  • All species with biomass burning emissions
  • All species with biogenic emissions
    • EXCEPT for the Bromine species (Br2, CHBr2), since these have a different biogenic source
  • All species with oceanic emissions
  • ACET direct emissions
  • ACET emissions from MBO
  • NO lightning emissions
  • NO fertilizer emissions
  • NO ship emissions
  • NO soil emissions
  • All dust emissions
  • All sea salt (accum & coarse mode) emissions
Comments on EMISSION RATIO differences: Switching the met field year from 2013 to 2016 had the greatest impact on emissions that rely on physical parameters such as temperature, pressure, PAR, 10-m wind speed, etc. Anthropogenic and biofuel emissions inventories (such as EDGAR, NEI, EMEP, etc.) that are read directly from disk are not affected.
Additional or summary comments:
  • There are many differences caused by the switch of the met field year, almost too many to summarize. We invite you to look at the plots for your species of interest.
  • GEOS-FP met data for 2013 was produced with version 5.11.0 of the GEOS-DAS system.
  • GEOS-FP met data for 2016 was produced with version 5.16.5 of the GEOS-DAS system. There were many updates made to the GEOS-DAS system since v5.11.0. See GMAO's GEOS Near-Real Time Data Products page for more information.
Approval
Requires further investigation: No
Approved by: Daniel Jacob
Date of approval: 02 Apr 2018

--Bob Yantosca (talk) 16:18, 9 April 2018 (UTC)

v11-02e using spun-up restart file

This section includes the assessment form for the 1-month benchmark simulation of v11-02e.

Description
New features added into GEOS-Chem
Feature Submitted by
Features affecting the full-chemistry simulation:
Initialize simulation with a new, spun-up restart file saved out from 1-year benchmark v11-02d-Run2 GCST
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2013
1-month benchmark finished on: Wed Mar 7 19:41:24 EST 2018
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50 GHz)
  • Wall time: 05:59:35
  • CPU time / wall time: 21.4413
  • % of ideal performance: 89.34%
  • Memory: 6.9177 GB
Compared to previous benchmark: v11-02e
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify): Initial conditions
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02e_newRst.tar.gz
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02e_newRst.tar.gz
Metrics
Global mean OH from log file (% change): 12.4840921614489 x 105 molec/cm3 (6.01%)
Methyl chloroform lifetime (% change): 5.0067 years (-4.77%)
Methane lifetime (% change): 8.4425 years (-4.76%)
O3 STE flux (% change): 461.3745 Tg O3/year (-21.24%)
At the SURFACE, list all species that changed by 10% or more: NO, O3, PAN, CO, ALK4, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, MSA, NH3, NH4, NIT, NITs, BCPI, OCPI, BCPO, OCPO, DST1, DST2, DST3, DST4, SALA, SALC, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CHBr3, CH2Br2, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG0-3, TSOA0-3, ISOG1-3, ISOA1-3, ASOG1-3, ASOAN, ASOA1-3, SOAP, SOAS, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, CFC115, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, BrSALA, BrSALC, CH3I, CH2I2, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI, OH, HO2
Comments on SURFACE differences: We expect using initial species concentrations from a properly spun-up restart file would have a large impact on most of the species concentrations.
At 500 hPa, list all species that changed by 10% or more: NO, O3, PAN, CO, ALK4, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, MSA, NH3, NH4, NIT, NITs, BCPI, OCPI, BCPO, OCPO, DST1, DST2, DST3, DST4, SALA, SALC, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CHBr3, CH2Br2, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG0-3, TSOA0-3, ISOG1-3, ISOA1-3, ASOG1-3, ASOAN, ASOA1-3, SOAP, SOAS, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, CCl4, CH3CCl3, CFC113, CFC14, CFC115, HCFC141b, HCFC142b, H1211, H2402, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, BrSALA, BrSALC, CH3I, CH2ICl, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI, OH, HO2
Comments on 500 hPa differences: See comments on surface differences
In the ZONAL MEAN differences, list all species that changed by 10% or more: NO, O3, PAN, CO, ALK4, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, MSA, NH3, NH4, NIT, NITs, BCPI, OCPI, BCPO, OCPO, DST1, DST2, DST3, DST4, SALA, SALC, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CHBr3, CH2Br2, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG0-3, TSOA0-3, ISOG1-3, ISOA1-3, ASOG1-3, ASOAN, ASOA1-3, SOAP, SOAS, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, N2O, OCS, CH4, BrCl, HCl, CCl4, CH3Cl, CH3CCl3, CFC113, CFC114, CFC115, HCFC141b, HCFC142b, CFC11, CFC12, HCFC22, H1211, H1301, H2402, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, BrSALA, BrSALC, CHCl3, CH2Cl2, CH3I, CH2I2, CH2ICl, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI, OH, HO2
Comments on ZONAL MEAN differences: See comments on surface differences
In the EMISSION RATIO maps, list all species that changed by 10% or more:
  • NO anthropogenic+biofuel emissions
  • NO fertilizer emissions
  • NO soil emissions
Comments on EMISSION RATIO differences:
  • NO anthropogenic, soil, and fertilizer emissions change by more than 10% over the oceans. This is due to changing NOx concentrations, which are used as input in the PARANOX and soil NOx extensions.
Additional or summary comments:
Approval
Requires further investigation: No
Approved by: GCST
Date of approval: 15 Mar 2018

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

v11-02e

This section includes the assessment form for the 1-month benchmark simulation of v11-02e.

Description
New features added into GEOS-Chem
Feature Submitted by
Features affecting the full-chemistry simulation:
Monthly mean surface methane distributions Lee Murray (NASA GISS/LDEO)
UCX updates, including:
Seb Eastham (Harvard)
Chris Holmes (Florida State)
Melissa Sulprizio (GCST)
Use monthly mean stratospheric prod/loss rates from UCX simulation for all species Melissa Sulprizio (GCST)
Update density of BC to 1.8 and add absorption enhancement factor in input.geos Xuan Wang (MIT)
Metal catalyzed oxidation of SO2 as a switch (on by default) Viral Shah (UW)
Becky Alexander (UW)
Add SOA updates and bug fixes, including: Sam Silva (MIT)
Chemistry Working Group
Aerosols Working Group
GCST
Eloise Marais (U. Birmingham)
Fix dust rainout settings to match BCPO Colette Heald (MIT)
Fixes in gckpp_HetRates.F90 for halogen chemistry Lei Zhu (Harvard)
Seb Eastham (MIT)
Add bug fix in sulfate_mod.F involving reorder IF statements for Fahey and Seinfeld decision algorithm Chris Holmes (FSU)
Becky Alexander (UW)
Sync GEOS-FP files on Harvard ftp with files at Dalhousie GCST
Features not affecting the full-chemistry simulation:
Spatially varying OM/OC as a switch (off by default) Sajeev Philip (Dalhousie/NASA Ames)
Update Rn-Pb-Be-Pasv simulations to use 72 levels by default Melissa Sulprizio (GCST)
Bug fixes for ISORROPIA II stable mode Jingyuan Shao (UW)
Shaojie Song (Harvard)
Change input.geos timestep units from minutes to seconds Lizzie Lundgren (GCST)
Fixes for ND21 binary diagnostics Lizzie Lundgren (GCST)
Fix bug in HEMCO reference time code Barron Henderson (US EPA)
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2013
1-month benchmark finished on: Wed Mar 7 00:38:16 EST 2018
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50 GHz)
  • Wall time: 06:15:26
  • CPU time / wall time: 21.880
  • % of ideal performance: 91.17
  • Memory: 6.9255 GB
Compared to previous benchmark: v11-02d
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify):
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02e.tar.gz
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02e.tar.gz
Metrics
Global mean OH from log file (% change): 11.7765580248825 x 105 molec/cm3 (-6.43%)

Updates that had the largest contribution to changing global mean OH:

Methyl chloroform lifetime (% change): 5.2574 years (5.44%)
Methane lifetime (% change): 8.8642 years (5.16%)
O3 STE flux (% change): 585.7695 Tg O3/year (0.86%)
At the SURFACE, list all species that changed by 10% or more: NO, PAN, ALK4, ISOP, HNO3, H2O2, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, CH2O, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, MSA, NH3, NH4, NIT, NITs, DST1, DST2, DST3, DST4, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG1, TSOG2, TSOG3, TSOG0, TSOA1, TSOA2, TSOA3, TSOA0, ASOG1, ASOG2, ASOG3, ASOAN, ASOA1, ASOA2, ASOA3, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, BrSALA, BrSALC, CH2I2, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI, OH, HO2
Comments on SURFACE differences: Below we summarize the notable changes caused by specific updates.

Monthly mean surface methane distributions

  • This update contributes to slight >10% decreases in several species concentrations at the surface over oceans where concentrations are low.

Remove initial stratospheric 2D mixing ratio option

  • This update contributes to >10% increases or decreases in many species concentrations over the oceans, poles, and aloft where species concentrations are low.

Use monthly mean stratospheric prod/loss rates from UCX simulation for all species

  • This update causes large % changes in the stratosphere for many species, as observed in the zonal mean plots.

Metal catalyzed oxidation of SO2

  • This update contributes to slight (~0.005 ppbv) decreases in SO2 concentrations at the surface, primarily over Africa, the Middle East, India, and China.

Add MTPO as a species in full-chemistry simulations

  • The addition of MTPO in the chemical mechanism caused large % decreases in MTPO concentrations globally. This also contributes to slight decreases in the monoterpene products present in the complex SOA scheme (TSOG*, TSOA*).

Fixes for isoprene SOA for consistency with Marais et al. (2016)

  • Updating the yield of LVOC in several reactions results in increased LVOC and LVOCOA as expected.
  • Changing daytime and nighttime gamma in function HETGLYX increases SOAGX globally as expected and causes slight decreases in GLYX.
  • Changing the Henry’s law constant in routine HSTAR_EPOX increases SOAIE and SOAME globally as expected and causes large % decreases globally in IEPOXA, IEPOXB, IEPOXD.
  • Removing INDIOL from AOD and aerosol mass calculations causes mainly small (<10%) changes in species concentrations.

Fix dust rainout settings to match BCPO

  • This fix results in large global increases in the dust species concentrations at the surface and aloft.

Fixes in gckpp_HetRates.F90 for halogen chemistry

  • This update causes large % changes globally in the halogen species, including Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, BrSALA, BrSALC, I2, IBr, ICl, I, IO, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4.
  • This update also contributes to increase in concentrations over southern oceans and decreases in concentrations at northern high latitudes for many species (NO, PAN, HNO3, N2O5, HNO4, MPN, NO2, NO3, HNO2).

Sync GEOS-FP files on Harvard ftp with files at Dalhousie

  • This fix, which addresses issues A3cld met fields, primarily causes changes in species concentrations over the oceans. Several species show increases (e.g. NO, H2O2, PPN, R4N2, N2O5, HNO4, SO2, NIT, MPN, NO2, NO3, OH, HO2), while some show decreases (e.g. ALK4, ISOP, MEK, ALD2, RCHO, MVK, MACR, NPMN, PRPE, DMS, SO4, NH4, many Br species, ISN1, RIPA, RIPB, RIPD, IPMN).

The following updates impacted species concentrations slightly, but the changes did not amount to >10%.

At 500 hPa, list all species that changed by 10% or more: NO, PAN, ALK4, ISOP, HNO3, H2O2, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, MSA, NH3, NIT, NITs, DST1, DST2, DST3, DST4, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CHBr3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG0, TSOA1, TSOA2, TSOA3, TSOA0, ISOA1, ISOA3, ASOG1, ASOG2, ASOG3, ASOAN, ASOA1, ASOA2, ASOA3, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, CCl4, CH3CCl3, CFC113, CFC114, CFC115, HCFC141b, HCFC142b, CFC11, H1211, H1301, H2402, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, BrSALA, BrSALC, CH3I, CH2I2, CH2ICl, CH2IBr, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI, OH, HO2
Comments on 500 hPa differences: See comments on surface differences
In the ZONAL MEAN differences, list all species that changed by 10% or more: NO, O3, PAN, CO, ALK4, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, NH3, NH4, NIT, NITs, DST1, DST2, DST3, DST4, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CG2Br2, CH3Br, MPN, ISOPND, ISOPNB, MOBA, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG0, TSOA0, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, LVOC, LVOCOA, ISN1OG, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, N2O, OCS, CH4, BrCl, HCl, CCl4, CH3Cl, CH3CCl3, CFC113, CFC114, CFC115, HCFC123, HCFC141b, HCFC142b, CFC11, CFC12, HCFC22, H1211, H1301, H2402, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, BrSALA, BrSALC, CH3I, CH2I2, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, OH, HO2
Comments on ZONAL MEAN differences: See comments on surface differences
In the EMISSION RATIO maps, list all species that changed by 10% or more:
  • NO anthropogenic+biofuel emissions
Comments on EMISSION RATIO differences:
  • NO anthropogenic emissions change by more than 10% over the oceans. This is due to changing NOx concentrations, which are used by the PARANOX ship plume model to compute ship emissions.
Additional or summary comments:

Special tropchem benchmark:

http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02e_Tropchem.tar.gz
Approval
Requires further investigation: No
Approved by: Daniel Jacob
Date of approval: 15 Mar 2018

--Melissa Sulprizio (talk) 15:39, 7 March 2018 (UTC)

v11-02d

This section includes the assessment form for the 1-month benchmark simulation of v11-02d.

Description
New features added into GEOS-Chem
Feature Submitted by
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
Sulfur oxidation by reactive halogens Qianjie Chen (UW)
Add bug fixes for GLYX, MGLY, IEPOX, and IMAE heterogeneous chemistry Sebastian Eastham (MIT)
Fix bug impacting dust radius obtained in RDUST_ONLINE for UCX simulations Melissa Sulprizio (GCST)
Add fix for IPMN reactions contributing to high P(Ox) and L(Ox) Melissa Sulprizio (GCST)
Eloïse Marais (U. Birmingham)
Features not affecting the full-chemistry simulation:
Update CH4 latitude bands for 2014-2016 Katie Travis (MIT)
Remove support for GCAP, GEOS-4, MERRA, GEOS-5 Melissa Sulprizio (GCST)
Timestamp HEMCO diagnostics filename with start time rather than end time Lizzie Lundgren (GCST)
Structural updates for netCDF diagnostics GCST
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2013
1-month benchmark finished on: Wed Jan 24 20:03:24 EST 2018
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50 GHz)
  • Wall time: 05:51:49
  • CPU time / wall time: 22.2282
  • % of ideal performance: 92.62
  • Memory: 6.5056 GB
Compared to previous benchmark: v11-02c
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify):
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02d.tar.gz
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02d.tar.gz
Metrics
Global mean OH from log file (% change): 12.5858002257716 x 105 molec/cm3 (2.60 %)
Methyl chloroform lifetime (% change): 4.9862 years (-3.15 %)
Methane lifetime (% change): 8.4289 years (New in this version)
O3 STE flux (% change): 580.7734 Tg O3/year (-0.71%)
At the SURFACE, list all species that changed by 10% or more: NO, O3, PAN, ALK4, ISOP, HNO3, H2O2, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, NH3, NH4, NIT, NITs, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, ASOG2, ASOAN, ASOA2, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, OH, HO2

New species added in this version (these show up uniform >100% changes): BrSALA, BrSALC, CHCl3, CH2Cl2, CHClI, CH2I2, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI

Comments on SURFACE differences: Below we summarize the notable changes caused by specific updates.

Halogen chemistry updates

  • This update caused large changes in many of the bromine and chlorine species. This is the direct result of updating the bromine and chlorine chemistry following Sherwen et al (2016). This update also introduced several new iodine species and reactions.

Sulfur oxidation by reactive halogens

  • This update causes large percent changes globally in SO2, SO4, SO4s, NH3, NH4, NIT, NITs, and may of the bromine species. The introduction of HOBr+S(IV) chemistry reduces HOBr (and several bromine species) over mid- and high latitudes in the Northern and Southern Hemisphere. These changes also impact NOx species, with large percent differences showing up in bands at ~60S and >60N where concentrations of these species are low.

Add bug fixes for GLYX, MGLY, IEPOX, and IMAE heterogeneous chemistry

  • As the description of this fix implies, this update causes large percent changes in GLYX, MGLY, IEPOX, and IMAE and related species. A bug was artificially increasing the het rate on all aerosols, so with this fix the concentrations of these species increases.

Fix bug impacting dust radius obtained in RDUST_ONLINE for UCX simulations

  • This fix impacts many species concentrations over northern Africa and the Middle East, where dust concentrations are high.
At 500 hPa, list all species that changed by 10% or more: NO, ALK4, ISOP, H2O2, MEK, ALD2, RCHO, MVK, MACR, NPMN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, NH3, NIT, NITs, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CH3Br3, MPN, ISOPND, ISOPNB, MOBA, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, ASOAN, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, OH, HO2

New species added in this version (these show up uniform >100% changes): BrSALA, BrSALC, CHCl3, CH2Cl2, CHClI, CH2I2, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI

Comments on 500 hPa differences: See comments on surface differences
In the ZONAL MEAN differences, list all species that changed by 10% or more: NO, O3, PAN, ALK4, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, NH3, NIT, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CHBr3, MPN, ISOPND, ISOPNB, MOBA, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, TOLU, XYLE, MTPA, LIMO, MTPO, EOH, MGLY, GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEOPXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, ISN1OG, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187, BrCl, HCl, H1211, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, OH, HO2

New species added in this version (these show up uniform >100% changes): BrSALA, BrSALC, CHCl3, CH2Cl2, CHClI, CH2I2, CH2ICl, CH2IBr, HOI, I2, IBr, ICl, I, IO, HI, OIO, INO, IONO, IONO2, I2O2, I2O3, I2O4, ISALA, ISALC, AERI

Comments on ZONAL MEAN differences: See comments on surface differences
In the EMISSION RATIO maps, list all species that changed by 10% or more:
  • NO anthropogenic+biofuel emissions
Comments on EMISSION RATIO differences:
  • NO anthropogenic emissions change by more than 10% over the oceans. This is due to changing NOx concentrations, which are used by the PARANOX ship plume model to compute ship emissions.
Additional or summary comments:

Special tropchem benchmark:

  • A 1-month benchmark for the tropchem mechanism with halogen chemistry was also completed. Those benchmark plots and metrics can be viewed at:
http://ftp.as.harvard.edu/pub/index/geos-chem/1mo_benchmarks/v11-02/v11-02d_Tropchem.tar.gz
Large differences in J-values between the standard benchmark and the tropchem simulation can be attributed to using online O3 for Fast-JX in the benchmark simulation and O3 columns from met in the tropchem simulation.
Approval
Requires further investigation: Mat Evans wrote:
We’d like a little bit more confirmation that there are no obvious bugs lurking in the code before we sign off on things from our perspective. The HOBR+S(IV) and the Br+ALD2 both have potentially large impacts and I think the problem with IPMN is indicative of what can be lurking if we don’t take some care to check on things.

Daniel Jacob wrote:

There are still some fixes to be done to halogen and SOA chemistry, but v11-2d gives reasonable results and greatly enhances our halogen simulation capability. I am approving v11-2d, and we will implement the residual halogen and SOA fixes in in v11-2e along with some other developments.
Approved by: Daniel Jacob
Date of approval: 12 Feb 2018

--Melissa Sulprizio (talk) 18:52, 7 December 2017 (UTC)

v11-02c

This section includes the assessment form for the 1-month benchmark simulation of v11-02c.

Description
New features added into GEOS-Chem
Feature Submitted by
Features affecting the full-chemistry simulation:
Enhance default GEOS-Chem simple SOA Sal Farina (Colorado State)
Aerosols Working Group Chairs
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)
Add aqueous isoprene uptake to SOA scheme Eloïse Marais (U. Birmingham)
Carbon balance (fix C creation) Sarah Safieddine (MIT)
Fix bugs for EOH and MGLY following implementation of PAN updates in v11-02a Melissa Sulprizio (GCST)
Update HEMCO from v2.0.004 to v2.1.001 Christoph Keller (NASA GMAO)
Features not affecting the full-chemistry simulation:
HEMCO updates:
Christoph Keller (NASA GMAO)
Andy Jacobson (NOAA)
Paulo Tuccella (L'Aquila)
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)
Fixes for the TOMAS simulation:
Sal Farina (Colorado State)
Jack Kodros (Colorado State)
Bob Yantosca (GCST)
Fix STE flux diagnostic and add to benchmark procedure Melissa Sulprizio (GCST)
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)
Combine timestep settings in input.geos in a Timesteps menu Melissa Sulprizio (GCST)
Update CO data used in 1-year benchmark plots Jenny Fisher (U. Wollongong)
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2013
1-month benchmark finished on: Wed Sep 6 00:02:51 EDT 2017
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50 GHz)
  • Wall time: 05:23:05
  • CPU time / wall time: 22.3511
  • % of ideal performance: 93.13%
  • Memory: 5.8683 GB
Compared to previous benchmark: v11-02b
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify):
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02c.tar.gz
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02c.tar.gz
Metrics
Global mean OH from log file (% change): 12.3036776026986 x 105 molec/cm3 (-0.63 %)
Methyl chloroform lifetime (% change): 5.1288 years (1.03 %)
O3 STE flux (% change): 585.2005 Tg O3/year (-42.77%)
At the SURFACE, list all species that changed by 10% or more: NO, O3, PAN, ALK4, ISOP, HNO3, H2O2, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, C3H8, CH2O, N2O5, HNO4, MP, DMS, SO2, SO4, NH3, NH4, NIT, NITs, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG0-3, TSOA0-4, ISOG1-3, ISOA1-3, ASOG1-3, ASOAN, ASOA1-3, EOH, MGLY, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, OH, HO2
Comments on SURFACE differences: Below we summarize the notable changes caused by specific updates.

Enhance default GEOS-Chem simple SOA

  • The following species are new in this version and show up as large positive ratios: SOAP, SOAS. See this table for species descriptions.

Updates to isoprene and monoterpene chemistry, which includes aqueous isoprene uptake

  • The majority of the large differences in species concentrations are caused by these chemistry updates. These updates included adjusting reaction rates, modifying products, and adding heterogeneous chemistry and greatly impacted species like ISOPND, ISOPNB, PROPNN, HAC, GLYC, MVKN, MACRN, R4N2.
  • PMN is now split into NPMN and IPMN (PMN from non-isoprene and isoprene sources). This is shown as a large decrease in NPMN in the ratio plots.
  • Changes in MTPA and MTPO concentrations are caused by updates to reaction rates in the monoterpene chemistry originally introduced in v11-02a.
  • The large decrease in LIMO is caused by the addition of gas-phase reactions of limonene with OH, O3, and NO3. This decrease also contributes to large decreases in SOA products (TSOA/G*) produced by oxidation of LIMO.
  • The decrease in ISOG2 and ISOA2 may be explained by the removal of the ISOP lost to reaction with NO3. This pathway was removed from the SOA code in the aqueous isoprene uptake updates.
  • The following species are new in this version and show up as large positive ratios: GLYX, ACTA, HPALD, DHDN, ETHLN, HCOOH, IEPOXA, IEPOXB, IEPOXD, ISN1, RIPA, RIPB, RIPD, IMAE, SOAIE, SOAME, SOAGX, SOAMG, LVOC, LVOCOA, ISN1OG, ISN1OA, MONITS, MONITU, HONIT, IONITA, MONITA, INDIOL, IPMN, HC187. See this table for species descriptions.

Carbon balance (fix C creation)

  • The large decrease in R4N2 globally is caused by (1) replacing R4N2 by PROPNN in the ISNOOB + NO3 and ISNOOB + NO reactions and (2) removing R4N2 as product from the ATO2 + NO reaction.

Fix bugs for EOH and MGLY following implementation of PAN updates in v11-02a

  • A bug in the species definitions of EOH and MGLY meant the concentrations for these species were not getting properly updated after gas-phase chemistry. This is now fixed and results in large % differences in the EOH and MGLY species concentrations

Update HEMCO from v2.0.004 to v2.1.001

  • This update contributes to >10% changes in some species north of 70°N.
  • Species such as DST1-4, SALA, and SALC have small (<10%) but structured changes at the surface. These changes are caused by the bug fixes, I/O updates, and regridding updates added in HEMCO v2.1.
At 500 hPa, list all species that changed by 10% or more: NO, PAN, ALK4, ISOP, HNO3, H2O2, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, CH2O, N2O5, HNO4, MP, DMS, SO2, NH3, NIT, NITs, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, BENZ, TOLU, XYLE, MTPA, LIMO, MTPO, TSOG1, TSOG3, TSOA1, TSOA3, TSOA0, ISOG1-3, ISOA1-3, EOH, MGLY, BrCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, OH, HO2
Comments on 500 hPa differences: See comments on surface differences
In the ZONAL MEAN differences, list all species that changed by 10% or more: NO, PAN, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, NPMN, PPN, R4N2, PRPE, CH2O, N2O5, HNO4, MP, DMS, SO2, NH3, NIT, NITs, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, MAP, NO2, NO3, HNO2, TOLU, XYLE, MTPA, LIMP, MTPO, TSOG1, TSOA1, TSOA3, ISOG2, ISOA2, EOH, MGLY, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClNO2, ClOO, OClO, Cl2, Cl2O2, OH, HO2
Comments on ZONAL MEAN differences: See comments on surface differences
In the EMISSION RATIO maps, list all species that changed by 10% or more:
  • Anthropogenic emissions: CO, NH3, NO, SO2, SO4
  • Ship emissions: CO, SO2
Comments on EMISSION RATIO differences: The band of positive differences ~70°N over Europe is caused by the HEMCO update from v2.0.004 to v2.1.001. Christoph Keller wrote:
I added an update to the regridding routines (regrid_a2a_mod.F90) to support non-global grids as well as regridding of missing values. Grid boxes with missing values are now ignored. The old code did not support missing values and would assign a value of 0.0 to all grid boxes with no input value. For regional emission inventories such as EMEP, this meant that all values outside of the EMEP domain were treated as zero's and the regridding did essentially dilute the emissions at the grid boundaries. The new code now ignores the grid boxes outside of the boundaries, which I find a better approach (and which explains the increase at 70N - right where the EMEP domain ends).
Additional or summary comments: Following an initial 1-month benchmark, the GCSC advised we restore the Henry's law constant for H2O2 wet deposition. This value was changed with the updates to isoprene and monoterpene chemistry. Daniel Jacob wrote:
For wetdep of H2O2 we should restore the old Henry’s law constant of 8.3E4exp[7400(1/T – 1/298)] because as Dylan points out that’s the physical value. For drydep of H2O2 we can keep the value of 5E7 as parameterized by Nguyen to fit his drydep data.
Approval
Requires further investigation:
  • Diagnostic issues will be resolved in the 1-year benchmark. This involves adding simple SOA to the PM2.5 diagnostic and fixing total SOA concentrations that are too high (likely from calling OASAVE too often). Both of these issues are manifested in the ND42 diagnostic output.
  • Open item is tracking of SOA in complex scheme and avoiding double-counting. This will be discussed in the Aerosol WG and will be resolved at a later date.
Approved by: Jeff Pierce, Eloise Marais, Jenny Fisher, Katie Travis, Daniel Jacob
Date of approval: 07 Sep 2017

--Melissa Sulprizio (talk) 17:43, 6 September 2017 (UTC)

v11-02b

This section includes the assessment form for the 1-month benchmark simulation of v11-02b.
New in this version: See also the assessment form for the 1-month benchmark of v11-02b with high performance option (GCHP).

Description
New features added into GEOS-Chem
Feature Submitted by
Features not affecting the full-chemistry simulation:
Source code updates for compatibility with high performance option (GCHP) Seb Eastham (Harvard)
Lizzie Lundgren (GCST)
Mike Long (GCST)
Jiawei Zhuang (Harvard)
Bob Yantosca (GCST)
Bug fixes for diagnostics:
Ilya Stanevic (Toronto)
Lee Murray (Rochester)
Write initial and final Ox mass to file when using the tagged O3 simulation Bob Yantosca (GCST)
Convert CO2 emissions units using dry pressure used in advection Meemong Lee (JPL)
Add QFED emissions for 2014-2016 Christoph Keller (NASA GMAO)
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2013
1-month benchmark finished on: Thurs June 07 23:08:02 EDT 2017
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50 GHz)
  • Wall time: 04:25:40
  • CPU time / wall time: 22.5723
  • % of ideal performance: 94.05%
  • Memory: 5.0852 GB
Compared to previous benchmark: v11-02a
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify):
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02b.tar.gz
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02b.tar.gz
Comparison with high performance option enabled:
(new in this version)

We now include a plots subdirectory called gchp_comparison containing 1-month standard simulation output comparisons using the standard "classic" capabilities (GCC) and the high performance option (GCHP). GCHP features the same science as GCC but operates on a cubed-sphere grid and is parallelized using a message-passing interface (MPI) implementation. GCHP improves upon GCC by (1) enabling more accurate transport, and (2) providing efficient scaling making finer resolution global simulations possible.

There is a separate benchmark approval form for the comparison of GCHP and GCC. See the Assessment form for the 1-month benchmark of v11-02b with high performance option.

Metrics
Global mean OH (from log file): 12.3818817518568 x 105 molec/cm3
Methyl chloroform lifetime: 5.0764 years
Did either of these change by more than 5%? No, both values are identical to values of the previous version.
At the SURFACE, list all species that changed by 10% or more: none
Comments on SURFACE differences: none
At 500 hPa, list all species that changed by 10% or more: none
Comments on 500 hPa differences: No differences.
In the ZONAL MEAN differences, list all species that changed by 10% or more: none
Comments on ZONAL MEAN differences: No differences.
In the EMISSION RATIO maps, list all species that changed by 10% or more: none
Comments on EMISSION RATIO differences: none
Additional or summary comments: As expected, merging in the modifications for GCHP has zero impact on the GEOS-Chem "classic" standard simulation.
Approval
Requires further investigation: No
Approved by: Lizzie Lundgren (GCST)
Date of approval: 13 June 2017

--Lizzie Lundgren (talk) 18:54, 13 June 2017 (UTC)

v11-02a

Here is the assessment form for 1-month benchmark simulation v11-02a.

Description
New features added into GEOS-Chem
Feature Submitted by
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
Fixes to correct ALK4 lumping issue Barron Henderson (US EPA)
PAN updates, including: Emily Fischer (CSU)
Monthly mean NEI2011 emissions GCST &
Katie Travis (Harvard)
Bug fixes in the GEOS-Chem sulfate module:
Prasad Kasibhatla (Duke)
Prasad Kasibhatla (Duke)
Qianjie Chen (UW)
Viral Shah (UW)
Fix bug in dry deposition aerodynamic resistance Brian Boys (Dalhousie)
Fix acetone parameterization in hcox_seaflux_mod.F90 GCST
Bug fix in cos(SZA) for start of timestep Lizzie Lundgren (GCST)
Include TOA pressure when calculating dry pressure edges Seb Eastham (Harvard)
Features not affecting the full-chemistry simulation:
Implement ISORROPIA v2.0 as a Fortran module Seb Eastham (Harvard) &
GCST
Updates to the HEMCO emissions component:
Christoph Keller (NASA GMAO)
Seb Eastham (Harvard)
Jessica Morena (Dalhousie)
Brian Boys (Dalhousie) & GCST
Updates to gain computational speedup:
Mike Long (Harvard)
Bob Yantosca (GCST)
Fixed typo in INIT_WINDOW Bob Yantosca (GCST)
netCDF file I/O updates:
Chris Holmes (Florida State)
Andy Jacobson (NOAA)
GCST
Makefile and build sequence updates:
Jiawei Zhang (Harvard)
GCST
Bug fixes for running UCX in ESMF environment Christoph Keller (NASA GMAO)
Bug fixes for diagnostics:
Aaron van Donkelaar (Dalhousie)
Jenny Fisher (U. Wollongong)
GCST
GCST
Jenny Fisher (U. Wollongong)
Chris Holmes (Florida State)
Removal of obsolete variables: GCST
Version, resolution, met fields used: v11-02, GEOS-FP (72L), 4x5, July 2013
1-month benchmark finished on: Sun Apr 16 01:37:36 EDT 2017
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50 GHz)
  • Wall time: 04:22:54
  • CPU time / wall time: 22.5896
  • % of ideal performance: 94.12%
  • Memory: 5.2038 GB
Compared to previous benchmark: v11-01 public release
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Anthro Emiss, Biogenic Emiss, Biomass Emiss, Photolysis, Chemistry, Other (please specify):
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02a.tar.gz
Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02a.tar.gz
Metrics
Global mean OH (from log file): 12.3128108012973 x 105 molec/cm3
Methyl chloroform lifetime: 5.1159 years
Did either of these change by more than 5%? No. Mean OH changed by -2.57% and MCF lifetime changed by 2.94%. These changes were primarily caused by two updates:
  • The JPL 15-10 update changed mean OH by -1.03% and MCF lifetime by 0.58%.
  • The PAN updates changed mean OH by -1.62% and MCF lifetime by 2.41%
At the SURFACE, list all species that changed by 10% or more: NO, O3, PAN, CO, ALK4, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, PMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, MSA, NH3, NH4, NIT, NITs, BCPI, OCPI, BCPO, OCPO, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CHBr3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, RIP, IEPOX, MAP, NO2, NO3, HNO2, BrCl, HCl, Cl, ClO, HOCl, ClNO3, ClOO, OClO, Cl2, Cl2O2, MTPA, LIMO, MTPO, TSOG1, TSOG2, TSOG2, TSOG0, TSOA1, TSOA2, TSOA3, TSOA0, ISOG1, ISOG2, ISOG3, ISOA1, ISOA2, ISOA3, BENZ, TOLU, XYLE, ASOG1, ASOG2, ASO3G3, ASOAN, ASOA1, ASOA2, ASOA3, OH, HO2
Comments on SURFACE differences:

Below we summarize the notable changes caused by specific updates.

Update chemistry rate constants based on JPL Publication 15-10

  • Notable changes caused by this update:
    • >10% increase in MVK, PRPE, N2O5, Cl2O2
    • >10% decrease in MACR, PMN, CHBr3
  • Matt Evans wrote
    The only major changes [from JPL 15-10] are to do with the MVK / MACR + OH changes we expected, changes in the PRPE +OH and changes in the N2O5 / HO2NO2 formation / loss rates. Unless I’ve missed something I think we are happy with all of this.

PAN updates

Monthly mean NEI2011 emissions

  • This update caused numerical noise differences over the US in NO, ALK4, RCHO, N2O5, DMS, SO2, NH3, NIT, BCPO, OCPO, NO2, NO3, HNO2, BENZ, TOLU, XYLE.

Bug fixes in the GEOS-Chem sulfate module:

  • The bug fix in SEASALT_CHEM caused a decrease in SO4s and NITs globally.
  • The bug fix in CHEM_NIT did not cause any changes that amount to >10%.
  • The fix for sulfate production in HET_DROP_CHEM caused an increase in SO2, NH4, and NITs over the oceans and a decrease in SO4, SO4s, NH3, and NITs over the oceans
  • The bug fixes in sulfate chemistry routines caused changes in SO2, SO4, SO4s, NH3, NH4, NIT, and NITs.
    • Viral Shah wrote:
      The effects of my proposed fixes are a bit more pronounced during the NH winter. The plots show significant increases in sulfate concentrations over large areas of the NH. The SO4 increases are accompanied by decreases in NIT in some areas. The increase in SO4 in the cold areas of NH is mostly because of the correction to HNO3 effective Henry's law coefficient that leads to higher cloud water pH (between 4 and 5) and higher production by in-cloud oxidation by O3. The lower SO4 concentrations in the SH are probably due to the fix to ∆H/R parameter. I expect the effects [of the LWC units fix] to be relatively small.

Fix bug in dry deposition aerodynamic resistance

  • This update caused slight decreases in many species at the surface. Species changed by >10% from this update include HNO3, H2O2, N2O5, and NIT.

The following updates impacted species concentrations slightly, but the changes did not amount to >10%.

At 500 hPa, list all species that changed by 10% or more: NO, PAN, CO, ALK4, ISOP, HNO3, H2O2, ACET, MEK, ALD2, RCHO, MVK, MACR, PMN, PPN, R4N2, PRPE, C3H8, CH2O, C2H6, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, MSA, NH3, NH4, NIT, NITs, BCPI, OCPI, BCPO, OCPO, Br2, Br, BrO, HOBr, HBr, BrNO2, BrNO3, CHBr3, MPN, ISOPND, ISOPNB, MOBA, PROPNN, HAC, GLYC, MVKN, MACRN, RIP, IEPOX, MAP, NO2, NO3, HNO2, BrCl, Cl, ClO, HOCl, ClNO3, ClOO, OClO, Cl2, Cl2O2, MTPA, LIMO, MTPO, TSOG1, TSOG2, TSOG2, TSOG0, TSOA1, TSOA2, TSOA3, TSOA0, ISOG1, ISOG2, ISOG3, ISOA1, ISOA2, ISOA3, BENZ, TOLU, XYLE, ASOG1, ASOG2, ASO3G3, ASOAN, ASOA1, ASOA2, ASOA3, OH, HO2
Comments on 500 hPa differences: See comments on surface differences.
In the ZONAL MEAN differences, list all species that changed by 10% or more: NO, PAN, ALK4, ISOP, HNO3, H2O2, MEK, ALD2, RCHO, MVK, MACR, PMN, PPN, PRPE, C3H8, CH2O, N2O5, HNO4, MP, DMS, SO2, SO4, SO4s, NH3, NH4, NIT, NITs, BCPI, OCPI, BCPO, OCPO, Br2, Br, BrO, HOBr, HBr, BrNO2, CHBr3, MPN, ISOPND, ISOPNB, MOBA, HAC, GLYC, MVKN, MACRN, RIP, IEPOX, MAP, NO2, NO3, HNO2, BrCl, Cl, ClO, HOCl, ClNO3, ClOO, OClO, Cl2, Cl2O2, MTPA, LIMO, MTPO, TSOG1, TSOG2, TSOG3, TSOG0, TSOA1, TSOA2, TSOA3, TSOA0, ISOG1, ISOG2, ISOG3, ISOA1, ISOA2, ISOA3, BENZ, TOLU, XYLE, ASOG1, ASOG2, ASOG3, ASOAN, ASOA1, ASOA2, ASOA3, OH, HO2
Comments on ZONAL MEAN differences: See comments on surface differences.
In the EMISSION RATIO maps, list all species that changed by 10% or more:
  • NO: Anthropogenic+biofuel, biomass burning, and fertilizer
Comments on EMISSION RATIO differences:
  • The differences in NO emissions can be explained by the following updates:
    • NO biomass burning emissions decrease because NO is partitioned directly to 40% PAN and 20% HNO3.
    • Changes in NO anthropogenic emissions are isolated to the oceans, indicating a change in ship emissions. Changes in ship emissions may result from slight changes in NO or NO2 concentrations.
    • Differences in fertilizer emissions are caused by changes in NOx concentrations.
  • Anthropogenic emissions over US change slightly (<10%) because of the change from hourly NEI2011 emissions to monthly mean NEI2011 emissions. This impacts the following species: ALD2, ALK4, BC, C2H6, C3H8, CH2O, CO, MEK, NH3, NO, OC, PRPE, SO2, SO4.
  • Slight (<10%) differences can be seen in the biomass burning emission difference plots. This is numerical noise caused by injecting biomass burning emissions at different levels (65% to the boundary layer, 35% to the free troposphere).
  • ACET ocean source decreases slightly (<10%) because of the fix for acetone parameterization in hcox_seaflux_mod.F90.
Additional or summary comments: Removal of biomass burning updates

After the 1-month benchmark for v11-02a, the GCSC decided we should remove the code updates for vertically distributed biomass burning emissions (see this wiki post for more details). In addition, after the 1-year benchmark v11-02a-Run0, the GCSC asked that we comment out the partitioning of NOx biomass burning emissions directly to PAN and HNO3 (see this wiki post for more details). An unofficial 1-month benchmark was run for v11-02a to evaluate the model changes with these two biomass burning updates removed. The plots for that run may be viewed at:

http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02a_final/output/

Results from our v11-02a benchmark using the GNU Fortran Compiler

In order to evaluate the performance of GEOS-Chem using the the free and open source GNU Fortran compiler, we performed an additional 1-month benchmark for v11-02a. For this additional benchamark, we used GNU Fortran v6.2.0 to compile GEOS-Chem instead of our usual Intel Fortran Compiler version 11.1.069.

For more details, see this post on our GNU Fortran compiler wiki page.

Summary: The benchmark using GNU Fortran yielded essentially identical results to the benchmark using Intel Fortran. This is very encouraging, as it will allow GEOS-Chem development to take place on computational platforms that do not have proprietary compilers (such as Intel Fortran or PGI Fortran), which can be prohibitively expensive to purchase.

Approval
Requires further investigation: Yes, a 1-year benchmark simulation
Approved by: Emily Fischer, Prasad Kasibhatla, Brian Boys, Daniel Jacob
Date of approval: 24 Apr 2017

--Melissa Sulprizio (talk) 19:27, 20 April 2017 (UTC)

1-year Rn-Pb-Be benchmarks

v11-02e-RnPbBePasv-Run2

Two 1-year Rn-Pb-Be simulations were performed using GEOS-Chem v11-02e. The simulations utilized 4° x 5° GEOS-FP met fields for the year 2016, with a 4-year spinup.

For comparison of the Pb-210 and Be-7 budgets to previous versions, please see the following posts on the Rn-Pb-Be simulation wiki page:

  1. Budget of Pb210 from 1-year benchmark simulations
  2. Budget of Be7 from 1-year benchmark simulations

Please note that we include a passive species to test for mass conservation in GEOS-Chem transport. The passive tracer continues to show that mass is not conserved when the non-local PBL mixing scheme (VDIFF) is turned on. We first discovered this issue in v11-01 and continue to complete the 1-year Rn-Pb-Be simulations with and without non-local PBL mixing to track this issue.

You may view the benchmark plots for the simulation by pointing your browser to the following links.

Using the non-local PBL mixing scheme (VDIFF):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/RnPbBePasv-Run2/NLPBL/output/

Using the full PBL mixing scheme (TURBDAY):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/RnPbBePasv-Run2/TURBDAY/output/

--Melissa Sulprizio (talk) 17:59, 16 March 2018 (UTC)

v11-02e-RnPbBePasv-Run1

Two 1-year Rn-Pb-Be simulations were performed using GEOS-Chem v11-02e. The simulations utilized 4° x 5° GEOS-FP met fields for the year 2013, with a 4-year spinup. This simulation marks the first Rn-Pb-Be benchmark to utilize 72 levels.

For comparison of the Pb-210 and Be-7 budgets to previous versions, please see the following posts on the Rn-Pb-Be simulation wiki page:

  1. Budget of Pb210 from 1-year benchmark simulations
  2. Budget of Be7 from 1-year benchmark simulations

Please note that we include a passive species to test for mass conservation in GEOS-Chem transport. The passive tracer continues to show that mass is not conserved when the non-local PBL mixing scheme (VDIFF) is turned on. We first discovered this issue in v11-01 and continue to complete the 1-year Rn-Pb-Be simulations with and without non-local PBL mixing to track this issue.

You may view the benchmark plots for the simulation by pointing your browser to the following links.

Using the non-local PBL mixing scheme (VDIFF):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/RnPbBePasv-Run1/NLPBL/output/

Using the full PBL mixing scheme (TURBDAY):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/RnPbBePasv-Run1/TURBDAY/output/

--Melissa Sulprizio (talk) 17:59, 16 March 2018 (UTC)

v11-02e-RnPbBePasv-Run0

Two 1-year Rn-Pb-Be simulations were performed using GEOS-Chem v11-02e. The simulations utilized 4° x 5° GEOS-FP met fields for the year 2013, with a 4-year spinup.

For comparison of the Pb-210 and Be-7 budgets to previous versions, please see the following posts on the Rn-Pb-Be simulation wiki page:

  1. Budget of Pb210 from 1-year benchmark simulations
  2. Budget of Be7 from 1-year benchmark simulations

Please note that we include a passive species to test for mass conservation in GEOS-Chem transport. The passive tracer continues to show that mass is not conserved when the non-local PBL mixing scheme (VDIFF) is turned on. We first discovered this issue in v11-01 and continue to complete the 1-year Rn-Pb-Be simulations with and without non-local PBL mixing to track this issue.

You may view the benchmark plots for the simulation by pointing your browser to the following links.

Using the non-local PBL mixing scheme (VDIFF):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/RnPbBePasv-Run0/NLPBL/output/

Using the full PBL mixing scheme (TURBDAY):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02e/RnPbBePasv-Run0/TURBDAY/output/

--Melissa Sulprizio (talk) 17:59, 16 March 2018 (UTC)

v11-02b-RnPbBePasv

Two 1-year Rn-Pb-Be simulations were performed using GEOS-Chem v11-02b. The simulations utilized 4° x 5° GEOS-FP met fields for the year 2013, with a 4-year spinup.

For comparison of the Pb-210 and Be-7 budgets to previous versions, please see the following posts on the Rn-Pb-Be simulation wiki page:

  1. Budget of Pb210 from 1-year benchmark simulations
  2. Budget of Be7 from 1-year benchmark simulations

Please note that we include a passive species to test for mass conservation in GEOS-Chem transport. The passive tracer continues to show that mass is not conserved when the non-local PBL mixing scheme (VDIFF) is turned on. We first discovered this issue in v11-01 and continue to complete the 1-year Rn-Pb-Be simulations with and without non-local PBL mixing to track this issue.

You may view the benchmark plots for the simulation by pointing your browser to the following links.

Using the non-local PBL mixing scheme (VDIFF):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b/RnPbBePasv/RnPbBePasv_VDIFF/output/

Using the full PBL mixing scheme (TURBDAY):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b/RnPbBePasv/RnPbBePasv_TURBDAY/output/

--Lizzie Lundgren (talk) 16:49, 24 July 2017 (UTC)

High Performance Option benchmarks

v11-02b HP 1-month

This section includes the assessment form for the 1-month benchmark simulation GEOS-Chem v11-02b with high performance v11-02b. GCHP features the same science as GCC but operates on a cubed-sphere grid and is parallelized using a message-passing interface (MPI) implementation. GCHP improves upon GCC by (1) enabling more accurate transport, and (2) providing efficient scaling making finer resolution global simulations possible.
See also the assessment form for the 1-month v11-02b benchmark simulation using the standard "classic" capability.

Description
New features added:

As this is the first formal GCHP benchmark, new features include the high performance model capability in its entirety. This includes all recent feature updates included in GCHP v11-02b as well as updates in v1.0.0 that were introduced over many years and are too numerous to list.

GEOS-Chem version: GCHP v11-02b
Meteorology fields: GEOS-FP (72L), 2x2.5, July 2013

NOTE: Input 2x2.5 met-fields are regridded to cubed-sphere grid resolution c24 (approximately equivalent to 4x5) within GCHP. In contrast, the GCC simulations in this comparison use input 4x5 met-fields.

Cubed-sphere resolution: c24
1-month benchmark finished on: Thurs June 12 22:01:08 EDT 2017
Performance statistics:
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50 GHz)
  • Wall time: 03:33:49
  • CPU time / wall time: 23.0338
  • % of ideal performance: 95.97%
  • Memory: 5.0852 GB
Compared to benchmark: GEOS-Chem Classic v11-02b
Additional comments:

We compared GEOS-Chem with and without the high performance option (GCHP and GCC respectively) for two different scenarios:

  1. Out-of-the-box
  2. Custom

The out-of-the-box comparison uses as-is code and run directory settings for both GCHP and GCC. Each simulation uses restart files regridded from higher resolution 2x2.5. GCHP uses input 2x2.5 met-fields which are regridded to cubed-sphere grid resolution c24 (approximately equivalent to 4x5) within the model. In contrast, the GCC simulations in this comparison use input 4x5 met-fields. GCC is compiled with netcdf diagnostics turned on and binary diagnostics turned off. GCHP output data are regridded from c24 (cubed-sphere grid) to 1x1.25 and GCC output data are regridded from 4x5 to 1x1.25 for comparison plotting.

The custom comparison includes additional modifications to the GCC code base and run directory to remove certain remaining known differences between GCHP and GCC. These include the following changes to GCC:

  1. Turn off OTD-LIS factors and set OTD-LIS scaling to 0.5568042545
  2. Turn off all 2D initial mixing ratios except for stratospheric H2O
  3. Turn off usage of the HEMCO restart file
  4. Force cos(SZA) at timestep mid-point to equal cos(SZA) calculated at timestep start

NOTE: Non-concentration diagnostics are not yet implemented in GCHP. The v11-02b HP benchmark therefore includes only comparisons of species concentrations. A complete set of diagnostics will be available pending implementation of netCDF diagnostics in GEOS-Chem (in progress).

Plots may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/v11-02/v11-02b/gchp_comparisons/
Metrics
Global mean OH: Diagnostic not yet implemented. See note in additional comments section above.
Methyl chloroform lifetime: Diagnostic not yet implemented. See note in additional comments section above.
Did either of these change by more than 5%? N/A
At the SURFACE, list all species that changed by 10% or more:
  • Out-of-the-box: All species changed by 10% or more except CH2Br2, CH3Br, N2O, OCS, CH4, CCl4, CH3Cl, CH3CCl3, CFC113, CFC114, CFC115, HCFC141b, HCFC142b, CFC11, CFC12, HCFC22, H1211, H1301, and H2404.
  • Custom: All species except those listed above and HCFC123.
Comments on SURFACE differences:

Sources of general differences include the following:

  1. Many species are impacted by different met-fields: GCC uses 4x5 lat/lon while GCHP uses 2x2.5 regridded to cubed-sphere c24 (approximately 4x5). Differences due to regridding also extends to other inputs such as emissions.
  2. Differences are greatest along coastlines for many species. Gradients are generally greater along coastlines and steep gradients impart the strongest differences in regridding. Differences in land cover will also impact the differences along coastlines.
  3. Large fractional differences occur where concentrations are very low and can be attributed to numerical noise.
  4. All differences are further amplified by regridding GCHP output from cubed-sphere resolution to lat/lon for comparison with GCC. Characteristics of output regridding differences include banded patterns of alternating positive and negative differences.
  5. Several large differences between out-of-the-box GCC and GCHP are not present or reduced in the custom run comparison due to turning off initial mixing ratios in the GCC custom run. This includes HCFC123 and differences in chlorinated species away from the poles (e.g. BrCl, ClNO2, and Cl2). Removal of the initial mixing ratios is slated for inclusion in GCC v11-02e.
  • Differences specific to particular species are as follows:
  1. Numerical noise differences in NH3 and NIT are due to ISORROPIA.
  2. Differences in dust species (DST1, DST2, DST3, DST4) are especially impacted by the effects of differences in met fields. Dust tuning factors and emissions patterns are both dependent on met fields. Furthermore, edge cases where dust is subject to small local emissions based on met conditions are very sensitive to regridding, whereas consistent emissions locations (e.g. over the Middle East) will produce dust in similar quantities on any grid.
  3. Differences in isoprene (ISOP) and related biogenic species (e.g. MEK, MVK, MACR, PMN, and PRPE) are also especially impacted by differences in met fields. The MEGAN biogenic emissions module is sensitive to grid resolution because emissions are estimated based on several met fields.
  4. Differences in NO are influenced by differences in PARANOX ship emissions, soil NOx emissions, and fertilizer NOx emissions, all of which use NOx concentrations as input. This will also be affected by resolution-specific lightning factors.
  5. There are notable differences in chlorinated species (e.g. BrCl, ClNO2, Cl2) at the poles which are due to the differences in vertical transport between lat/lon and cubed-sphere grids, specifically since GCC estimates faster vertical mixing at the poles.
At 500 hPa, list all species that changed by 10% or more: All species changed by 10% or more except O3, CO, ACET, and all species listed above for the surface.
Comments on 500 hPa differences: See comments on surface differences above.
In the ZONAL MEAN differences, list all species that changed by 10% or more: All species changed by 10% or more.
Comments on ZONAL MEAN differences:
  • See comments on surface differences above.
  • H2O is consistently lower in GCHP in the stratosphere. This could be the result of differences in stratosphere/troposphere exchange, differences in convection, or GMAO mass scaling used for mass conservation. These would also all affect stratospheric CFC115.
Additional or summary comments:

Overall, the largest fractional differences between GCC and GCHP occur where concentrations are very low. Differences between GCHP and GCC are primarily due to the effects of the following:

  1. Different met fields (2x2.5 lat/lon regridded to cubed-sphere in GCHP versus 4x5 lat/lon in GCC)
  2. Higher-order vertical transport in GCHP due to a more recent version of tpcore
  3. Transport on a cubed-sphere grid versus a cartesian grid
  4. Regridding GCHP inputs and outputs to and from the cubed-sphere grid

The concentration differences we see are inevitable given the fundamental differences between GEOS-Chem with and without high performance capability.

Approval
Requires further investigation: Yes, a 1-year benchmark simulation
Approved by: Daniel Jacob
Date of approval: 16 June 2017

--Lizzie Lundgren (talk) 18:05, 13 June 2017 (UTC)

v11-02b HP 1-year Standard

v11-02b-HP-Run1

This 1-year benchmark simulation is an unofficial benchmark to help validate GEOS-Chem with the high performance option enabled. It uses offline lightning and dust emissions developed to minimize differences due to lightning and dust and provide a better understanding of differences observed in v11-02b-HP-Run0.

We compared two GEOS-Chem model versions for this 1-year benchmark:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Green v11-02b-HP-Run0 (GCHP) GEOS-FP,
72L, 2° x 2.5° regridded online to c24
2013 Updates introduced in GCHP v11-02b:
  • The high performance model capability in its entirety.
  • Differences from v11-02b-HP-Run0:
    • v11-02b-HP-Run1 differs from v11-02b-HP-Run0 by using offline dust and lightning emissions developed by Seb Eastham (Harvard) to reduce error associated with lightning and dust. This serves to better understand the sources of differences observed in v11-02b-HP-Run0.
    • Lightning and dust emissions are pre-calculated by storing the hourly average flux from a previous 4x5 simulation using GEOS-FP meteorology. Lightning NO mass flux is stored as an hourly 3D field at 4x5 and dust emissions are stored as four hourly 2D fields, one for each of the four dust tracers. Since GCHP conservatively regrids the emissions rates to the C24 grid, using these 4x5 emissions ensures that the temporal and spatial distributions are well matched albeit with some diffusion due to the regridding process.
    • See also the resolution dependency of emissions open issue description on the GEOS-Chem HP v11-02 wiki page.
N/A
Red v11-02a-Run1 (GEOS-Chem Classic) GEOS-FP,
72L, 4° x 5°
2013 Updates introduced in GEOS-Chem Classic v11-02a: 11.829
Black Observations        

The output plots for Run1 may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b-HP/Run1/output
mget *

You may also view the PDF files online by pointing your browser to:

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b-HP/Run1/output/

--Lizzie Lundgren (talk) 16:40, 27 July 2017 (UTC)

v11-02b-HP-Run0

This 1-year benchmark simulation is the first official benchmark to validate GEOS-Chem with the high performance option enabled.

We compared two GEOS-Chem model versions for this 1-year benchmark:

Color Model Version Met Type Year Updates affecting the benchmark simulation Annual Mean OH
[105 molec/cm3]
Green v11-02b-HP-Run0 (GCHP) GEOS-FP,
72L, 2° x 2.5° regridded online to c24
2013 Updates introduced in GCHP v11-02b:
  • The high performance model capability in its entirety. As discussed in the v11-02b HP 1-month benchmark assessment form, notable features of GEOS-Chem HP that result in expected differences with the GCC benchmark results include:
    • Use of different met fields, specifically 2° x 2.5° lat/lon regridded to cubed-sphere in GCHP versus 4° x 5° lat/lon in GCC
    • Higher-order vertical transport in GCHP due to a more recent version of tpcore
    • Transport on a cubed-sphere grid versus a cartesian grid
    • Regridding GCHP inputs and outputs to and from the cubed-sphere grid
    • Lack of availability of OTD-LIS redistribution grids for cubed-sphere resulting in lightning NOX distribution differences. This is part of a larger issue of the resolution dependency of emissions.
  • There are no updates introduced in v11-02b that impact the full-chemistry simulation.
N/A
Red v11-02a-Run1 (GEOS-Chem Classic) GEOS-FP,
72L, 4° x 5°
2013 Updates introduced in GEOS-Chem Classic v11-02a: 11.829
Black Observations        

The output plots for Run0 may be downloaded from:

ftp ftp.as.harvard.edu
cd gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b-HP/Run0/output
mget *

You may also view the PDF files online by pointing your browser to:

http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b-HP/Run0/output/

Please also view the following pages comparing this version to past 1-year benchmarks:

--Lizzie Lundgren (talk) 15:30, 12 July 2017 (UTC)

v11-02b HP 1-year Rn-Pb-Be

This 1-year benchmark simulation is the first official benchmark to validate GEOS-Chem with the high performance option enabled.

We compared 1-year Rn-Pb-Be simulation output of GEOS-Chem Classic v11-02b and GCHP v11-02b for this 1-year benchmark:

Versions Compared Internal Resolution Met Source Met Resolution Year Notes
GEOS-Chem Classic v11-02b 72L, 4° x 5° GEOS-FP 72L, 4° x 5° 2013 See also the Rn-Pb-Be 1-year benchmark for GEOS-Chem v11-02b which compares v11-02b Rn-Pb-Be results with previous version v11-01i
GCHP v11-02b 72L, c24 (cubed-sphere) GEOS-FP 72L, 2° x 2.5° regridded online to c24 2013 Notable features of v11-02b HP that result in expected differences with the v11-02b Rn-Pb-Be simulation include:
  • Use of different met fields, specifically 2° x 2.5° lat/lon regridded to cubed-sphere in GCHP versus 4° x 5° lat/lon in GCC
  • Higher-order vertical transport in GCHP due to a more recent version of tpcore
  • Transport on a cubed-sphere grid versus a cartesian grid
  • Regridding GCHP inputs and outputs to and from the cubed-sphere grid

You may view the benchmark comparison plots for the simulation by pointing your browser to the following links.

Using the non-local PBL mixing scheme (VDIFF):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b-HP/RnPbBePasv/VDIFF/output/

Using the full PBL mixing scheme (TURBDAY):
http://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-02/v11-02b-HP/RnPbBePasv/TURBDAY/output/

Please note that this benchmark does not include comparisons of Pb-210 and Be-7 budgets or passive tracer mass conservation because the required diagnostics are not yet output by GCHP.

--Lizzie Lundgren (talk) 17:57, 25 July 2017 (UTC)