GEOS-Chem v8-01-04

From Geos-chem
Revision as of 14:59, 6 May 2009 by Bmy (Talk | contribs)

Jump to: navigation, search

Overview

EMISSIONS RELEASE date March 10, 2009

Contains everything in v8-01-03 plus the following updates:

  1. Annual scale factors for NOx-CO-SOx from, applied to ALL inventory if needed (Aaron)
  2. Daily NOx variability from EDGAR applied to all others anthro NOx inventories
  3. Several emissions updates
    • Streets 2006, with monthly variability
    • EMEP 1990-2005 w/ ship emissions separated
    • Canadian CAC,
    • EPA fixes = CA transport + Rynda scaling,
    • VISTAS/ARP for NOx in the USA
    • MEGAN2
  4. new EPA mask (Overlap Bravo/Epa and CAC/Epa)
  5. ship NOx emitted as HNO3 + 10*O3
  6. pre-Arctas SO2 ship emissions from Streets
  7. GFED2 biomass options:
    • Monthly mean (as was in the code previously)
    • 8-day mean
    • 3-hourly mean (selected months only)
    • 3-hourly synoptic (selected months only)

Philippe Le Sager has prepared a document (PDF format) which describes the anthropogenic emissions sources in some detail.

1-year benchmarks

The following benchmark simulations have been done with a few different versions of GEOS-5 meteorology:

  • GEOS-5.0.1
  • GEOS-5.1.0
  • GEOS-5.2.0

At present, all of the GEOS-5.0.1 met data files for GEOS-Chem are being replaced with the newer GEOS-5.1.0 data product. GEOS-5.1.0 will have data available from late 2003 thru September 2008. Data after September 2008 will be the newest GEOS-5.2.0 data product. For more information about these GEOS-5 data versions, please see the GEOS-5 documentation from GMAO.

Run0

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

Color Quantity Plotted Met Field Type A.T.E. algorithm Anthro Emissions Biogenic Emissions Advection Annual Mean OH
[105 molec/cm3]
Red v8-01-01 Run0 GEOS-5 4x5
version 5.0.1

spinup: 2005
run: 2005
ISORROPIA GEIA/Piccot (scale year 1998)
EPA/NEI99 w/o ICARTT fix
EDGAR ship emissions
"old" MEGAN AEF's for isoprene "old" tpcore 12.006
Green v8-01-01 Run1 " " RPMARES " " " " " " 12.040
Blue v8-01-04-Run0 GEOS-5 4x5
version 5.0.1
with "quick fix"
for optical depth


spinup: 2005
run: 2005
" " " " "new" MEGAN AEF's for isoprene "new" tpcore 10.689
Black Observations            

SUMMARY:

  1. v8-01-01 Run0 vs. v8-01-01 Run1 is a clean comparison for the switch from ISORROPIA to RPMARES.
  2. v8-01-04 Run1 used the same emissions as v8-01-01 Run1, but also implemented the the "quick fix" for the GEOS-5 optical depth. Dylan Millet reported that the the "quick fix" was responsible for lowering the mean OH concentration from ~12 x 105 to ~11 x 105 molecules/cm3 in simulations that he had performed. Therefore the "quick fix" is probably the major factor driving the drop in mean OH in the v8-01-04 runs.
  3. All runs listed here used the same:
    • biomass emissions (GFED2 monthly)
    • biofuel emissions
    • lightning NOx emissions
    • aircraft emissions

The output plots for Run0 may be downloaded from:

ftp ftp.as.harvard.edu
cd pub/geos-chem/1yr_benchmarks/v8-01-04/geos5/2005/Run0/output

Run1

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

Color Quantity Plotted Met Field Type Photolysis Anthro Emissions Biogenic Emissions Advection Annual Mean OH
[105 molec/cm3]
Red v8-01-01 Run1 GEOS-5 4x5
version 5.0.1

spinup: 2005
run: 2005
linear
cloud overlap
GEIA/Piccot (scale year 1998)
EPA/NEI99 w/o ICARTT fix
EDGAR ship emissions
"old" MEGAN AEF's for isoprene "old" tpcore 12.040
Green v8-01-04 Run0 GEOS-5 4x5
version 5.0.1
with "quick fix"
for optical depth


spinup: 2005
run: 2005
" " " " "new" MEGAN AEF's for isoprene "new" tpcore 10.689
Blue v8-01-04 Run1 " " approximate
random
cloud overlap
GEIA/Piccot
EDGAR emissions
EMEP European emissions
BRAVO Mexican emissions
David Streets 2006 emissions
CAC Canadian emissions
EPA/NEI99 with ICARTT fix
EDGAR ship emissions
ARCTAS ship SO2 emissions
Anthro scale year 2005
"new" MEGAN AEF's for isoprene " " 10.294
Black Observations            

SUMMARY:

  1. v8-01-04 Run1 used the same emissions as v8-01-01 Run1, but also implemented the the "quick fix" for the GEOS-5 optical depth. Dylan Millet reported that the the "quick fix" was responsible for lowering the mean OH concentration from ~12 x 105 to ~11 x 105 molecules/cm3 in simulations that he had performed. Therefore the "quick fix" is probably the major factor driving the drop in mean OH in the v8-01-04 runs.
  2. v8-01-04 Run1 vs. v8-01-04 Run0 is a clean comparison for the new anthropogenic emissions options that were implemented into v8-01-04.
  3. All runs used the same:
    • A.T.E. Algorithm = RPMARES
    • GFED2 biomass emissions
    • Aircraft NOx emissions
    • Lightning NOx emissions

The output plots for Run1 may be downloaded from:

ftp ftp.as.harvard.edu
cd pub/geos-chem/1yr_benchmarks/v8-01-04/geos5/2005/Run1/output

Run2

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

Color Quantity Plotted Met Field Type Anthro Emissions Lightning NOx
[Tg N/yr]
Annual Mean OH
[105 molec/cm3]
Red v8-01-04 Run0 GEOS-5 4x5 version 5.0.1
with "quick fix" for optical depth

spinup: 2005
run: 2005
GEIA/Piccot (scale year 1998)
EPA/NEI99 w/o ICARTT fix
EDGAR ship emissions
5.979 10.689
Green v8-01-04 Run1 " " GEIA/Piccot
EDGAR emissions
EMEP European emissions
BRAVO Mexican emissions
David Streets 2006 emissions
CAC Canadian emissions
EPA/NEI99 with ICARTT fix
EDGAR ship emissions
ARCTAS ship SO2 emissions
Anthro scale year 2005
5.979 10.294
Green v8-01-04 Run2 GEOS-5 4x5 version 5.1.0,
"reprocessed" met fields

spinup: 2004
run: 2005
" " 5.936 11.099
Black Observations        

SUMMARY:

  1. v8-01-04 Run1 vs. v8-01-04 Run0 is a clean comparison for the new anthropogenic emissions options that were implemented into v8-01-04.
  2. v8-01-04 Run2 vs. v8-01-04 Run1 is a clean comparison between the GEOS-5.0.1 and GEOS-5.1.0 met products.
  3. All runs used the same:

The output plots for Run1 may be downloaded from:

ftp ftp.as.harvard.edu
cd pub/geos-chem/1yr_benchmarks/v8-01-04/geos5/2005/Run2/output

Previous issues now resolved in v8-01-04

Post-Release Patches

A couple of issues were reported after the release of v8-01-04:

Fabien Paulot (fabienpaulot@gmail.com) wrote:

The tracer id is hardcoded in the sulfate module (sulfate_mod.f which causes error when calling the ship emission module (GET_EMEP_ANTHRO) if the user wants to modify the tracer list. I just replaced 26 by IDTSO2 and 30 by IDTNH3 to fix the problem. See line 6594.

This issue was due to a typo in sulfate_mod.f. The error only occurred during aerosol-only offline simulations.


Prasad Kasibhatla (psk9@duke.edu) wrote:

I also noticed in the standard distribution, that the CPU time for the one month full chem simulation at Harvard went from 56 minutes for v8-01-02 to 75 minutes for v8-01-03 to 121 minutes for v8-01-04, an increase of a factor of 2.2. I was wondering if you have looked at this slowdown?

This issue was traced to the way the !$OMP PARALLEL DO loops were implemented in the TPCORE advection module tpcore_fvdas_mod.f. Claire Carouge has fixed this by changing the order of some of the parallel loops in the code. The updated TPCORE now runs much faster than before!


Hongyu Liu (hyl@nianet.org) reported:

The following references need to be deleted from the OBJS section in the Makefile:
   arctas_ship_emiss_mod.o       \
   cac_anthro_mod.o              \
   scale_anthro_mod.o            \
   tpcore_fvdas_mod.o            \
   vistas_anthro_mod.o           \
since these are already listed in the MODS section. Having these listed twice will cause a link-time error.


Colette Heald (heald@atmos.colostate.edu) wrote:

I wanted to let you know that to get v8-01-04 to compile [with the SunStudio compiler], I needed to modify the code in Makefile.sparc on L354 (i.e. under the line for gamap_mod.o) from:
   f90 -fpp -stackvar -O0 -xfilebyteorder=big16:%all
to the following:
   f90 -fpp -O0 -stackvar -xfilebyteorder=big16:%all $(FARCH) -c $*.f
otherwise it fails in input_mod.f compilation with messages that it did not recognize gamap_mod.f.
I also wanted to point out that there is a syntax error on L1040 of rpmares_mod.f (a comma after the write(6,*) but before the statement), which can cause compile to fail depending on the compiler options.


Fabien Paulot (fabienpaulot@gmail.com) also reported...

..an overestimate in the North-South mass flux computed in the new TPCORE. A geometric scale factor was applied twice when computing the flux. The updated TPCORE now correctly computes the NS mass flux.


Claire Carouge (ccarouge@seas.harvard.edu) wrote:

A GEOS-Chem user reported a slowdown also in the mass conservation calculation in transport_mod.f. I've changed the loops to parallelize over the tracers (it previously parallelized over levels) and it runs faster now.

All users of GEOS-Chem v8-01-04 should replace the existing versions of sulfate_mod.f, tpcore_fvdas_mod.f, and transport_mod.f with the fixed files. If you are using the PGI compiler, then also get the new Makefile.pgi file. If you are using the SunStudio compiler, then get the new Makefile.sparc file.

The updated files are available for download from:

ftp ftp.as.harvard.edu
cd pub/geos-chem/patches/v8-01-04
get sulfate_mod.f
get tpcore_fvdas_mod.f90
get transport_mod.f
get Makefile.sparc
get Makefile.pgi

Note

All users of GEOS-Chem v8-01-04 should turn on the mass flux diagnostics (ND24/ND25/ND26) in input.geos to be safe. See this bug report.

--Bob Y. 10:58, 6 May 2009 (EDT)