GEOS-Chem v9-01-01

From Geos-chem
Revision as of 13:40, 23 September 2010 by Bmy (Talk | contribs) (Duplicate diagnostics)

Jump to: navigation, search

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)

Starting runs at times other than 00 GMT

The MEGAN biogenic emissions reads in several days of A3 data before the first timestep of a G-C simulation. In routine DO_OPEN_A3 (which determines if it is the right time to open a new file), we have this IF statement:

     LOGICAL, SAVE :: FIRST    = .TRUE.
        ...   
     ! Open A-3 file if it's 01:30 GMT,  or on the first call
     IF ( NHMS == 013000 .or. FIRST ) THEN
        DO_OPEN = .TRUE.
        ...
   

Therefore, when MEGAN was turned on, FIRST was set to .FALSE. because we were reading more than one day of A3 data. This caused the IF statement above to fail if we wanted to start a G-C simulation at any other time than 00 GMT.

The quick fix was to introduce an optional logical flag RESET, such that:

     ! Reset the FIRST flag if
     IF ( PRESENT( RESET ) ) THEN
        IF ( RESET ) FIRST = .TRUE.
     ENDIF
   

This makes sure that the FIRST flag is not affected by the prior A3 reads by the MEGAN emissions.

We now call OPEN_A3_FIELDS with RESET=.TRUE. from main.f for the very first time we read in the A3 data (before the time-stepping loop).

The analogous fix was also made for the MERRA reading code in merra_a1_mod.f.

--Bob Y. 09:40, 23 September 2010 (EDT)

Outstanding issues not yet resolved in v8-03-02