GEOS-Chem v7-04-13
From Geos-chem
Contents
Overview
BETA RELEASE November 2007
Contains everything in GEOS-Chem v7-04-12, with the following additions:
- Offline H2/HD simulation (cf. Lyatt Jaegle et al)
- Bug fix in "smvgear.f" -- due to a design flaw, you could have a situation where you could have negative or zero tracer in a box even if the global & local error criteria were satisfied. If this condition is encountered, we now revert to the concentrations at the end of the previous internal timestep, and then reduce the internal timestep and try again until we get convergence.
- Bug fix -- replace GEMISNOX w/ module arrays to avoid common block errors in lightning & aircraft NOx emissions
- Bug fix -- cap variable tropopause at 200hPa for polar regions (90S-60S and 60N-90N)
- Lee Murray's local lightning redistribution of lightning flashes
- Minor bug fixes for GCAP and diagnostics
Also note that the Last year of GFED2 biomass emissions is still 2005 in this version. If you want to use the 2006 emissions you will have to edit the source code module "gfed2_biomass_mod.f".
1-Year Benchmarks
Run
The 3 models compared (and their colors on the plots) were:
Color | Quantity Plotted | Met Type | Lightning algorithm |
---|---|---|---|
Red | v7-04-12 Run1 | GEOS-4 4x5 | ancient |
Green | v7-04-12 Run2 | GEOS-4 4x5 | near-land w/ OTD-LIS regional scaling |
Blue | v7-04-13 | GEOS-4 4x5 | updated non-near-land w/ OTD-LIS local scaling |
Black | Observations |
NOTES:
- Results from v7-04-13-Run can be downloaded from FTP directory
- pub/geos/1yr_benchmarks/v7-04-13/geos4/2005/Run/output/pdf
- Output plots from v7-04-13-Run can be viewed on the web HERE.
- v7-04-13-Run uses the same input conditions and emissions as v7-04-12-Run2, except for the updated lightning emissions as described above.
Previous issues now resolved in v7-04-13
regrid_1x1_mod.f
This bug has now been fixed in GEOS-Chem v7-04-13.
Mike Barkley (mbarkley@staffmail.ed.ac.uk) wrote:
- I think I've found an error in the regrid_1x1_mod.f subroutine (attached in the text file):
SUBROUTINE REGRID_MASS_TO_2x25( I1, J1, L1, IN, I2, J2, OUT )
- There is a do loop over longitude with the upper limit defined as the input latitude (J1) instead what should (?) be the output longitude (I2) - I've indicated where this in the program, Which is correct? We didn't notice this until we were running multi-processor 2x2.5 simulations on different servers.
- The bug was:
!----------------------- ! Non-polar latitudes !----------------------- DO J = 2, J2-1 ... DO I = 1, J1
- which needs to be replaced by:
!----------------------- ! Non-polar latitudes !----------------------- DO J = 2, J2-1 ... DO I = 1, I1
--Bob Y. 16:59, 7 May 2009 (EDT)