Difference between revisions of "GEOS-Chem v9-01-02"

Revision as of 17:24, 31 March 2011

Overview

BETA RELEASE -- Spring 2011

• Previous version: GEOS-Chem v9-01-01
• Next version: GEOS-Chem v9-01-03

What's new in this version

GEOS-Chem v9-01-02 contains the following major updates and improvements:

1. APM aerosol microphysics (F. Yu & G. Luo, SUNY Albany)
2. Updated dry deposition velocities for aerosols over ice and snow (J. Fisher)
3. Various emissions updates
   1. Imposed seasonality on NH3 emissions from David Streets 2000 inventory (J. Fisher)
   2. GFED3 biomass emissions
   3. RETRO anthropogenic VOC emissions
   4. FAA aircraft emissions
   5. Updated soil NOx emissions
4. Fixes for several unresolved issues from v9-01-01

New data directories

The following new data directories have been added for GEOS-Chem v9-01-02. 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/APM_201011                                             # For APM microphysics only

GEOS_2x2.5/CO_201003/fossilfuel_andres/annual_v2010             # CO2 annual emissions from CDIAC
GEOS_2x2.5/CO_201003/fossilfuel_andres/monthly_v2010            # CO2 monthly emissions from CDIAC

GEOS_4x5/CO_201003/fossilfuel_andres/annual_v2010               # CO2 annual emissions from CDIAC
GEOS_4x5/CO_201003/fossilfuel_andres/monthly_v2010              # CO2 monthly emissions from CDIAC

GEOS_1x1/Streets_200607/NH3-Streets-SeasonalScalar.generic.1x1  # For NH3 seasonality in Streets 2000

--Bob Y. 11:25, 15 March 2011 (EDT)

Previous issues now resolved in v9-01-02

Patches from v9-01-01

The following patches that were issued after the release of GEOS-Chem v9-01-01 are now standardized in v9-01-02:

1. Double counting of biofuel emissions over Asia
2. Updated annual and monthly CO2 emissions from CDIAC
3. Typos in vdiff_mod.f90

--Bob Y. 11:36, 25 March 2011 (EDT)

Minor bug fixes in ISORROPIA II

Please see this wiki post for a description of two minor bug fixes that were implemented in source code file ISOROPIA/isorropiaIIcode.f. This bug does not affect the current GEOS-Chem ISORROPIA configuration, but may affect future development.

--Bob Y. 13:13, 1 March 2011 (EST)

Diagnostic updates

The following GEOS-Chem diagnostics have been updated:

Aerosol optical depth

Patrick Kim discovered some inconsistencies in the way aerosol optical depth diagnostics were being computed in GEOS-Chem v9-01-01 and prior versions. We have now corrected these inconsistencies.

Sungshik Patrick Kim wrote:

Here is an example of a daily (first page) and monthly (second page) average AOD for January 2008 from the public release of v9-01-01 and a newer version that addresses the changes we discussed earlier vs. the corresponding average AERONET data, a difference plot, and a histogram of the differences by grid box. [T]he largest differences appear to be due to the change in met field rather than changing the timing of the AOD calculation and it may be largest in the Arctic due to uncertainties in the previous calculation of RH.

Planeflight diagnostic

The MERRA SEAICExx fields are now saved to the plane flight diagnostic, per request of Jenny Fisher.

--Bob Y. 10:17, 7 March 2011 (EST)

Bug fixes for compatibility with the PGI compiler

Gan Luo from SUNY/Albany pointed out several minor typos that cause the PGI compiler to choke when building the GEOS-Chem executable. These will all be fixed in v9-01-02.

(1) Line 303 of GeosCore/depo_mercury_mod.f:

Change this:
 #include 'define.h'

to this:
 #include "define.h"

(2) Line 169 of GeosCore/ocean_mercury_mod.f:

Change this:
 REAL*8, PARAMETER     :: Kd_part   = 10**(5.5)

to this:
 REAL*8, PARAMETER     :: Kd_part   = 10d0**(5.5d0)

(3) Line 1281 of GeosCore/ocean_mercury_mod.f:

Change this:
 HgPaq(I,J)   = MAX ( HgPaq(I,J) , 0.0 )

to this:
 HgPaq(I,J)   = MAX ( HgPaq(I,J) , 0d0 )

(4) In subroutine DO_MERRA_CONVECTION of convection_mod.f

Change this:
 REAL*8            :: ALPHA,       ALPHA2,   CLDBASE, CMFMC_BELOW

to this:
 INTEGER           :: CLDBASE
 REAL*8            :: CMFMC_BELOW, ALPHA,    ALPHA2

(5) In module vdiff_mod.f90: Several numerical constants had to be rewritten so that explicitly use "D" exponents. In other words:

Change numbers like these:
 1.
 1.e-12

to these:
 1.d0
 1.d-12

In short, the PGI compiler expects constants to agree strictly with the type of the variable that they are being assigned to (i.e. if A is a REAL*8, then you must use A = 1.0d0 instead of A = 1.0). IFORT will correct for this situation during compilation, but PGI will throw an error.

--Bob Y. 13:27, 2 March 2011 (EST)

Out-of-bounds error in readchem.f

It was discovered that this code near the end of subroutine readchem.f:

                 ELSEIF ( J == MAPPL(ILBRO2H,NCS) .or.
    &                     J == MAPPL(ILBRO2N,NCS) .or.
    &                     J == MAPPL(ILTRO2H,NCS) .or.
    &                     J == MAPPL(ILTRO2N,NCS) .or.
    &                     J == MAPPL(ILXRO2H,NCS) .or.
    &                     J == MAPPL(ILXRO2N,NCS) ) THEN
                    ITS_NOT_A_ND65_FAMILY(J) = .FALSE.
                    EXIT

will cause an out-of-bounds error if any of the IL* indices are zero. This error manifested itself in a standard 43-tracer simulation.

We have corrected this error by first making sure that each IL* index is nonzero before doing the corresponding IF ( J = ... ) test:

                 ! Avoid out-of-bounds-errors by making sure that each
                 ! IL* index is nonzero before doing the IF test for J. 
                 ! (bmy, 3/1/11)
                 IF ( ILBRO2H > 0 ) THEN 
                    IF ( J == MAPPL(ILBRO2H,NCS) ) THEN
                       ITS_NOT_A_ND65_FAMILY(J) = .FALSE.
                       EXIT
                    ENDIF
                 ENDIF

                 ... etc ...

--Bob Y. 14:03, 2 March 2011 (EST)

NaN output in mercury diagnostics

Please see this wiki post for a description of a bug in the Hg simulation that was causing NaN to appear in some diagnostic output. This is now fixed in v9-01-02.

--Bob Y. 10:45, 9 March 2011 (EST)

Code.v9-01-01

Bug fix for leap years with GCAP

Recall that the GCAP met fields are climatological GCM output, and do not have any leap years (i.e. there are no Feb 29ths).

In GEOS-Chem v9-01-01 and prior versions, routine SET_CURRENT_TIME in time_mod.f) uses the Astronomical Julian Date algorithms (located in julday_mod.f), which assume that there are leap-years. If you ran GEOS-Chem with the GCAP met fields for more than one year, the current code caused the time & date to be computed incorrectly at the start of the second year. GEOS-Chem would die because the wrong met fields would be read in from disk.

We have added an improved algorithm fix into SET_CURRENT_TIME for GCAP met fields only. Whenever a leap-year day (Feb 29th) is encountered, we shall just skip ahead a day to March 1st and compute the corresponding time & date values. However, we must also keep track of how many leap years we have encountered since the start of the GEOS-Chem simulation. This count is needed so that we can increment the current Astronomical Julian Date so as to "skip over" from Feb 29th to Mar 1st for each new leap year that we encounter.

This fix now skips from Feb 28 to Mar 1st for all leap years. We have tested this for multiple leap years and it works fine.

--Bob Y. 13:58, 14 March 2011 (EDT)

Outstanding issues not yet resolved in v9-01-02

TBA