GEOS-Chem v9-01-03
Contents
- 1 Overview
- 2 Previous issues now resolved in v9-01-03
- 3 Outstanding issues not yet resolved in v9-01-03
Overview
BETA RELEASE -- TBD
- Previous version: GEOS-Chem v9-01-02
- Next version: ???
What's new in this version
GEOS-Chem v9-01-03 shall contain the following major updates and improvements. (NOTE: The order of inclusion is subject to change):
| Feature | Type | Submitted by | Benchmark (1 Month) |
Status |
|---|---|---|---|---|
| Fix for GLCO3/GLPAN bug in standard chemistry mechanism | Bug Fix | Fabien Paulot (Harvard) | v9-01-03a | Approved 08 Dec 2011 |
| Updated acetone chemistry (ocean source & sink) | Science | Emily Fischer (Harvard) | v9-01-03b | Approved 14 Dec 2011 |
| Retirement of obsolete GEIA biogenic emissions | Structural | Emily Fischer (Harvard) Dylan Millet (U. Minnesota) |
" " | " " |
| OVOC dry deposition updates | Science | Jingqiu Mao (NOAA/GFDL) | v9-01-03c | Approved 19 Dec 2011 |
| Improved sea salt emission & deposition | Science | Lyatt Jaeglé (U. Washington) | v9-01-03d | Approved 12 Jan 2012 |
| Update dry deposition to use local surface pressure | Science | Lyatt Jaeglé (U. Washington); GEOS-Chem Support Team | " " | " " |
| Dust submicron size distribution for optics | Science | David Ridley (Colo. State) | " " | " " |
| Improved snow scavenging and washout parameterization | Science | Qiaoqiao Wang (Harvard) | v9-01-03e + Rn-Pb-Be benchmark |
Pending approval |
| Structural changes to clarify science behind washout | Structural | Helen Amos (Harvard) | " " | " " |
| Updates for Grid-independent GEOS-Chem | Structural | GEOS-Chem Support Team | v9-01-03f | GI code has reached a stage where we should sync with the standard code |
| PARANOX ship plume emissions | Science | Geert Vinken (Eindhoven) | v9-01-03g | TBD |
| SOA simulation with semi-volatile POA | Science | Havala Pye (EPA) GEOS-Chem Support Team |
v9-01-03h | TBD |
| Stratospheric P and k (monthly climatological) | Science | Lee Murray (Harvard) | v9-01-03i | TBD |
| EPA/NEI05 North American Hg emissions | Science | Yanxu Zhang (U. Washington) Lyatt Jaeglé (U. Washington) |
v9-01-03j | TBD |
| Streets future Hg emissions | Science | Bess Corbitt (Harvard) | " " | " " |
| Satellite-based NOx emission trends | Science | Lok Lamsal (NASA GSFC) | " " | " " |
| Hg(II) gas-aerosol partitioning | Science | Helen Amos (Harvard) | " " | " " |
| Nested-grid Hg simulation | Science | Yanxu Zhang (U. Washington) Lyatt Jaeglé (U. Washington) |
" " | " " |
| Nested-grid CO2 | Science | Yuxuan Wang (Tsinghua) | " " | " " |
| Historical emission inventories of SO2, NOx, BC, and POA | Science | Eric Leibensperger (MIT) | " " | " " |
| Add scaling of lightning NOx for 2x2.5 MERRA simulations | Science | Lee Murray (Harvard) Matt Cooper (Dalhousie) |
" " | " " |
| Bug fix in emep_mod.F for offline simulations | Bug fix | GEOS-Chem Support Team | " " | " " |
| Bug fix to prevent crash in wetdep when encountering anomalously low PRECTOT values | Bug fix | GEOS-Chem Support Team | " " | " " |
| Minor fix in Rn-Pb-Be simulation to prevent missing drydep diagnostics | Bug fix | GEOS-Chem Support Team | " " | " " |
| Correct out-of-bounds error in offline aerosol simulation | Bug fix | GEOS-Chem Support Team | " " | " " |
| Bug fixes in diag3.F | Bug fix | Mat Evans (Leeds) GEOS-Chem Support Team |
" " | " " |
| Bug fixes for nested-grid simulation | Bug fix | Dylan Millet (U. Minnesota) | " " | " " |
| Updated tagged Ox simulation | Structural | Hongyu Liu (Natl' Institute of Aerospace) | TBD | TBD |
| Updating photolysis from FAST-J to FAST-JX | Science | Jingqiu Mao (Princeton) | TBD | TBD |
| Daily and hourly GFED3 biomass emissions | Science | Prasad Kasibhatla (Duke) | TBD | We will include this as soon as it is completed |
--Melissa Payer 16:04, 23 January 2012 (EST)
--Bob Y. 10:46, 25 January 2012 (EST)
Previous issues now resolved in v9-01-03
Bug fix for GLCO3 and GLPAN in standard mechanism
This update was tested in the 1-month benchmark simulation v9-01-03a and approved on 08 Dec 2011.
Prior to v9-01-02, some species in the globchem.dat and ratj.d files were left as inactive, but should have been declared as dead. For a full description of this issue, please see this wiki post on the Chemistry Issues page.
--Bob Y. 17:07, 14 December 2011 (EST)
Correct out-of-bounds error in offline aerosol simulation
NOTE: This bug fix is slated for inclusion into v9-01-03j.
When using the offline aerosol simulation, an out-of-bounds error was found to occur in array WTAERA within subroutine RDAER (aerosol_mod.F). To fix the error, add the following lines to INIT_COMODE (comode_mod.F) below the offline aerosol simulation IF statement:
ALLOCATE( WTAREA( ITLOOP, NDUST+NAER ), STAT=AS )
IF ( AS /= 0 ) CALL ALLOC_ERR( 'WTAREA' )
WTAREA = 0d0
ALLOCATE( WERADIUS( ITLOOP, NDUST+NAER ), STAT=AS )
IF ( AS /= 0 ) CALL ALLOC_ERR( 'WERADIUS' )
WERADIUS = 0d0
--Melissa Payer 17:47, 11 January 2012 (EST)
Bug fixes in diag3.F
Bug fix for mass transport flux diagnostics
NOTE: This bug fix is slated for inclusion into v9-01-03i.
Mat Evans wrote:
- In diag3.f line 1740
ARRAY(:,:,1:LD24) = MASSFLEW(:,:,LD24:1:-1,N) / SCALEDYN
- Should be
ARRAY(:,:,1:LD24) = MASSFLEW(:,:,LLPAR:LLPAR-LD24+1:-1,N) / SCALEDYN
- And the same for LD25 Line 1773
- And for same for LD26 line 1806
- Everything is fine if LD24 = LLPAR but it seems to go wrong if you try to output a subset of the model levels..
--Melissa Payer 17:28, 12 January 2012 (EST)
Bug fix in ND67 sea level pressure diagnostic
When plotting SLP in GAMAP for GEOS-5 at both 4°x5° and 2°x2.5° resolutions, the values appeared to be half of what they should be (SLP ~500 hPa globally).
In diag3.F, the lines:
CASE ( 18 )
SCALEX = SCALE_ND67
UNIT = 'hPa'
should be changed to:
CASE ( 18 )
#if defined( MERRA )
SCALEX = SCALE_A1 ! MERRA SLP is hourly
#else
SCALEX = SCALE_I6 ! Otherwise SLP is 6-h inst.
#endif
UNIT = 'hPa'
As it is now, SCALE_ND67 = SCALE_A3 but SLP is an I6 field for GCAP, GEOS-3, GEOS-4, and GEOS-5.
--Melissa Payer 17:42, 12 January 2012 (EST)
Bug fixes for nested-grid simulation
Dylan Millet wrote:
- My student Lu is running N. American nested grid with the v9-3-1b that Bob provided. We found a few minor glitches, mostly specific to the nested grid, that we thought we'd pass along. Perhaps they've already been fixed on your end.
- in land_mercury_mod.f, the variable "SOIL_EMS_FAC" is undefined for the 05x0666 case.
- the Makefiles in the GeosApm and GeosTomas directories each have an intermediate "else" that needs to be changed to an "elif"
- in tpcore_bc_mod.F, when reading in the boundary conditions if ITS_A_NEW_DAY, the call to ITS_A_NEW_DAY() needs to be ITS_A_NEW_DAY(.TRUE.) so that it will read in the boundary conditions actually at midnight. Otherwise it won't read them in at the start of the run.
--Melissa Payer 11:11, 18 January 2012 (EST)
Structural changes to clarify science behind washout
Helen Amos wrote:
- Qiaoqiao's aerosol scavenging updates should *not* be applied to HNO3. HNO3 should still be scavenged according to Eq (14) from Jacob (2000). This issue came up a while back and as a result I wrote code to split washout into three cases: aerosols, HNO3, and soluble gases. I restructured the washout code specifically to handle Qiaoqiao's updates and make sure HNO3 didn't change, so if no one objects I suggest these structure changes go into the std code with Qiaoqiao's updates.
--Melissa Payer 16:01, 23 January 2012 (EST)
Outstanding issues not yet resolved in v9-01-03
TBD
