Difference between revisions of "GEOS-Chem v7-04-13"

From Geos-chem
Jump to: navigation, search
Line 53: Line 53:
 
# 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.
 
# 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.
  
--[[User:Bmy|Bob Y.]] 11:18, 6 May 2009 (EDT)
+
== Previous issues now resolved in v7-04-13 ==
 +
 
 +
=== regrid_1x1_mod.f ===
 +
 
 +
This bug has now been fixed in [[wiki:geos-chem:dev_forum|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
 +
 
 +
--[[User:Bmy|Bob Y.]] 16:59, 7 May 2009 (EDT)

Revision as of 20:59, 7 May 2009

Overview

BETA RELEASE November 2007

Contains everything in GEOS-Chem v7-04-12, with the following additions:

  1. Offline H2/HD simulation (cf. Lyatt Jaegle et al)
  2. 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.
  3. Bug fix -- replace GEMISNOX w/ module arrays to avoid common block errors in lightning & aircraft NOx emissions
  4. Bug fix -- cap variable tropopause at 200hPa for polar regions (90S-60S and 60N-90N)
  5. Lee Murray's local lightning redistribution of lightning flashes
  6. 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:

  1. Results from v7-04-13-Run can be downloaded from FTP directory
    • pub/geos/1yr_benchmarks/v7-04-13/geos4/2005/Run/output/pdf
  2. Output plots from v7-04-13-Run can be viewed on the web HERE.
  3. 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)