Difference between revisions of "GEOS-Chem v8-01-04"
(→Run2) |
|||
Line 181: | Line 181: | ||
|-align="center" | |-align="center" | ||
| Red | | Red | ||
− | | [[#Run0 | + | | [[#Run0|v8-01-04 Run0]] |
| GEOS-5 4x5 version 5.0.1<br>with [[GEOS-5_issues#GEOS-5_update_and_clarification|"quick fix" for optical depth]]<br><br>spinup: 2005<br>run: 2005 | | GEOS-5 4x5 version 5.0.1<br>with [[GEOS-5_issues#GEOS-5_update_and_clarification|"quick fix" for optical depth]]<br><br>spinup: 2005<br>run: 2005 | ||
| GEIA/Piccot (scale year 1998)<br>EPA/NEI99 w/o ICARTT fix<br>EDGAR ship emissions | | GEIA/Piccot (scale year 1998)<br>EPA/NEI99 w/o ICARTT fix<br>EDGAR ship emissions | ||
Line 188: | Line 188: | ||
|-align="center" | |-align="center" | ||
| Green | | Green | ||
− | | [[#Run1 | + | | [[#Run1|v8-01-04 Run1]] |
| " " | | " " | ||
| GEIA/Piccot<br>EDGAR emissions<br>EMEP European emissions<br>BRAVO Mexican emissions<br>David Streets 2006 emissions<br>CAC Canadian emissions<br>EPA/NEI99 with ICARTT fix<br>EDGAR ship emissions<br>ARCTAS ship SO2 emissions<br>Anthro scale year 2005 | | GEIA/Piccot<br>EDGAR emissions<br>EMEP European emissions<br>BRAVO Mexican emissions<br>David Streets 2006 emissions<br>CAC Canadian emissions<br>EPA/NEI99 with ICARTT fix<br>EDGAR ship emissions<br>ARCTAS ship SO2 emissions<br>Anthro scale year 2005 | ||
Line 195: | Line 195: | ||
|-align="center" | |-align="center" | ||
| Green | | Green | ||
− | | [[# | + | | [[#Run2|v8-01-04 Run2]] |
| GEOS-5 4x5 version 5.1.0,<br>"reprocessed" met fields<br><br>spinup: 2004<br>run: 2005 | | GEOS-5 4x5 version 5.1.0,<br>"reprocessed" met fields<br><br>spinup: 2004<br>run: 2005 | ||
| " " | | " " |
Revision as of 16:11, 6 May 2009
Contents
Overview
EMISSIONS RELEASE date March 10, 2009
Contains everything in v8-01-03 plus the following updates:
- Annual scale factors for NOx-CO-SOx from, applied to ALL inventory if needed (Aaron)
- Daily NOx variability from EDGAR applied to all others anthro NOx inventories
- 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
- new EPA mask (Overlap Bravo/Epa and CAC/Epa)
- ship NOx emitted as HNO3 + 10*O3
- pre-Arctas SO2 ship emissions from Streets
- 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:
- v8-01-01 Run0 vs. v8-01-01 Run1 is a clean comparison for the switch from ISORROPIA to RPMARES.
- 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.
- 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:
- 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.
- 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.
- 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:
- 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.
- 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.
- All runs used the same:
- A.T.E. Algorithm = RPMARES
- Advection routine: "new" tpcore
- Biogenic emissions: "new" MEGAN AEF's for isoprene
- GFED2 biomass emissions
- Aircraft NOx emissions
- Photolysis: approximate random cloud overlap
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)