Difference between revisions of "GEOS-Chem v9-01-03"
(→What's new in this version) |
(→Bug fix in streets_anthro_mod.F for offline simulations) |
||
Line 360: | Line 360: | ||
=== Bug fix in streets_anthro_mod.F for offline simulations === | === Bug fix in streets_anthro_mod.F for offline simulations === | ||
+ | |||
+ | '''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-03_benchmark_history#v9-01-03h|v9-01-03h]] and approved on 09 Mar 2012.''''' | ||
'''''[mailto:fabienpaulot@gmail.com Fabien Paulot] wrote:''''' | '''''[mailto:fabienpaulot@gmail.com Fabien Paulot] wrote:''''' |
Revision as of 14:15, 12 March 2012
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 |
Approved 02 Feb 2012 |
Structural changes to clarify science behind washout | Structural | Helen Amos (Harvard) | " " | " " |
PARANOX ship plume emissions | Science | Geert Vinken (Eindhoven) | v9-01-03f | Approved 17 Feb 2012 |
Bug fixes for David Streets emissions for years 2005 and earlier | Bug fix | Jenny Fisher (Harvard) | " " | " " |
Bug fix in emep_mod.F for offline simulations | Bug fix | GEOS-Chem Support Team | v9-01-03g | Approved 27 Feb 2012 |
Bug fix in nei2005_anthro_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) | " " | " " |
Bug fix for reading Hg emissions | Bug fix | Jenny Fisher (Harvard) | " " | " " |
Add scaling of lightning NOx for 2x2.5 MERRA simulations | Science | Lee Murray (Harvard) Matt Cooper (Dalhousie) |
" " | " " |
Initial source code modifications for GEOS-5.7.2 | Structural | GEOS-Chem Support Team | " " | " " |
Hg(II) gas-aerosol partitioning | Science | Helen Amos (Harvard) | v9-01-03h | Approved 09 March 2012 |
Nested-grid CO2 | Science | Yuxuan Wang (Tsinghua) | " " | " " |
Historical emission inventories of SO2, NOx, BC, and POA | Science | Eric Leibensperger (MIT) | " " | " " |
CH4 simulation updates | Science | Kevin Wecht (Harvard) | " " | " " |
Remove support for GEOS-3 meteorology | Structural | GEOS-Chem Support Team | " " | " " |
Bug fix in streets anthro mod.F for offline simulations | Bug fix | Fabien Paulot (Harvard) | " " | " " |
Updates from the grid-independent GEOS-Chem model code | Structural | GEOS-Chem Support Team | v9-01-03i | GI code was tested in a CTM in Feb 2012. The validated code will be merged back into the standard model repository at some point thereafter (date as yet unknown) |
Daily and 3-hourly GFED3 biomass emissions | Science | Prasad Kasibhatla (Duke) | v9-01-03j | TBD |
SOA simulation with semi-volatile POA | Science | Havala Pye (EPA) GEOS-Chem Support Team |
v9-01-03k | TBD |
Stratospheric P and k (monthly climatological) | Science | Lee Murray (Harvard) | TBD | TBD |
EPA/NEI05 North American Hg emissions | Science | Yanxu Zhang (U. Washington) Lyatt Jaeglé (U. Washington) |
TBD | TBD |
Nested-grid Hg simulation | Science | Yanxu Zhang (U. Washington) Lyatt Jaeglé (U. Washington) |
" " | " " |
Streets future Hg emissions | Science | Bess Corbitt (Harvard) | " " | " " |
Satellite-based NOx emission trends | Science | Lok Lamsal (NASA GSFC) | " " | " " |
Fix to jv_spec_aod.dat | Bug fix | Gabriele Curci (University of L'Aquila) | TBD | TBD |
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 |
--Bob Y. 10:46, 25 January 2012 (EST)
--Melissa Payer 14:36, 8 March 2012 (EST)
New data directories
The following new data directories have been added for GEOS-Chem v9-01-03. You will have to download the directories relevant to your simulation.
You can download these directories with anonymous FTP or the Wget utility. For instructions, please see Chapter 2.4, Downloading the GEOS-Chem shared data directories in the GEOS-Chem Online User's Guide.
GEOS_1x1/PARANOX_201202 # For PARANOX ship plume model
--Melissa Payer 16:53, 22 February 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)
Structural changes to clarify science behind washout
This update was tested in the 1-month benchmark simulation v9-01-03e and approved on 02 Feb 2012.
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)
Correct out-of-bounds error in offline aerosol simulation
This update was tested in the 1-month benchmark simulation v9-01-03g and approved on 27 Feb 2012.
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
This update was tested in the 1-month benchmark simulation v9-01-03g and approved on 27 Feb 2012.
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
This update was tested in the 1-month benchmark simulation v9-01-03g and approved on 27 Feb 2012.
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
This update was tested in the 1-month benchmark simulation v9-01-03g and approved on 27 Feb 2012.
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)
Bug fix in streets_anthro_mod.F for offline simulations
This update was tested in the 1-month benchmark simulation v9-01-03h and approved on 09 Mar 2012.
Fabien Paulot wrote:
- I just found a little bug in the streets emission routine
!================================================================= ! Apply seasonal variation to NH3 based on seasonality from ! Lex Bouwman. Follow methodology in emep_mod.f (jaf, 3/2/11) !================================================================= ! Get TAU value for 1995, since the data is timestamped w/ this TAU1995 = GET_TAU0( GET_MONTH(), 1, 1995 ) ! For NH3 only ... IF ( TRACERN == IDTNH3 ) THEN ! File name containing scaling factors FILENAME_S = TRIM( DATA_DIR_1x1 ) // & 'Streets_200607/NH3-Streets-SeasonalScalar.generic.1x1' ! Tracer number for scale factor data RATIOID = 74 ! Echo info WRITE( 6, 101 ) TRIM( FILENAME_S ) 101 FORMAT( ' - READ_STREETS: Reading ', a ) ! Read scaling factors CALL READ_BPCH2( FILENAME_S, 'RATIO-2D', RATIOID, & TAU1995, I1x1, J1x1-1, & 1, SCALAR_1x1, QUIET=.TRUE. ) ! Apply seasonal scalar to NH3 emissions ARRAY(:,:,1) = ARRAY(:,:,1) * SCALAR_1x1(:,:,1) ENDIF
- The instruction TRACERN == IDTNH3 should be replaced by TRACER == 30 to avoid issues with the offline chemistry where idtnh3 is not equal to 30.
- This routine is being called by
CALL READ_STREETS( FILENAME, 'ANTHSRCE', 30, TAU2000, NH3 )
- for NH3
- So it seems to me, you are not going to scale NH3 if IDTNH3 is not equal to 30.
--Melissa Payer 11:51, 2 March 2012 (EST)
Outstanding issues not yet resolved in v9-01-03
TBD