GEOS-Chem v9-01-01

Overview

PUBLIC RELEASE -- TBD

Will contain everything in GEOS-Chem v8-03-02, plus:

1. Capability to run with GMAO MERRA met fields (R. Yantosca)
2. Improvements to wet scavenging algorithms for MERRA met data (H. Amos, B. Corbitt)

--Bob Y. 10:14, 22 September 2010 (EDT)

Previous issues now resolved in v8-03-02

Patches from v8-03-02

The following bug fix patches are now standardized in v9-01-01:

1. Transport fix for GEOS-5 nested grid simulations at 0.5° x 0.666°
2. Typo in globchem.dat affecting yield of isoprene nitrates
3. Mike Barkley noted that the KPP solver files in the KPP/standard/ and the KPP/SOA directories had to be rebuilt to incorporate the changes in the globchem.dat file. This was done on 17 Sep 10 and added to the GEOS-Chem Git repository as a patch in GEOS-Chem v8-03-02 (name "New KPP scheme for SOA and standard chemistry", tag "v8-03-02-Patch-KPP").

--Bob Y. 10:14, 22 September 2010 (EDT)

Duplicate diagnostics

Chris Holmes wrote:

I turned on diagnostics ND14-15 (convection and BL mixing) and ND24-26 (U/V/W transport fluxes), plus the other diagnostics that are on by default. I ran the simulation for 1 day. The ctm.bpch output contains *two* copies of ND14-15 and ND24-26 outputs. The other diagnostics appear normal.

In fact, the diagnostics are not written twice in the ctm.bpch file. The problem is a double definition of the categories of these diagnostics in gamap_mod.f. These diagnostics appear then twice in diaginfo.dat with two different offsets (3000 and 4000). GAMAP will then print these diagnostics twice. The fix is to comment one definition of these diagnostics in gamap_mod.f. Please comment the lines with the 4000 offset (see the SPACING * 4) in gamap_mod.f:

!      N           = N + 1
!      CATEGORY(N) = 'CV-FLX-$'
!      DESCRIPT(N) = 'Upward flux from wet conv'
!      OFFSET(N)   = SPACING * 4
!
!      N           = N + 1
!      CATEGORY(N) = 'TURBMC-$'
!      DESCRIPT(N) = 'Upward flux from PBL mixing'
!      OFFSET(N)   = SPACING * 4
!
!      N           = N + 1
!      CATEGORY(N) = 'EW-FLX-$'
!      DESCRIPT(N) = 'E/W transport flux'
!      OFFSET(N)   = SPACING * 4
!
!      N           = N + 1
!      CATEGORY(N) = 'NS-FLX-$'
!      DESCRIPT(N) = 'N/S transport flux'
!      OFFSET(N)   = SPACING * 4
!
!      N           = N + 1
!      CATEGORY(N) = 'UP-FLX-$'
!      DESCRIPT(N) = 'Up/down transport flux'
!      OFFSET(N)   = SPACING * 4

--Ccarouge 14:58, 21 September 2010 (EDT)