MERRA implementation details
NOTE: Information subject to change!
On this page we list the necessary modifications to GEOS-Chem for MERRA compatibility.
- 1 Overview
- 2 Required coding changes in GEOS-Chem
- 3 Resolved issues
- 4 Outstanding issues
The GMAO MERRA data product is a 30-year reanalysis done with the GEOS-5.2.0 assimilation system. Therefore, within GEOS-Chem we can (for the most part) treat MERRA in the same way as we do for the GEOS-5 meteorology.
However, in order to reduce the volume of data produced (and for other technical concerns), GMAO chose to reduce the vertical and/or horizontal resolution of some important MERRA data fields. Also, some MERRA data fields have different units than in the GEOS-5 product. Because of this, some additional special handling may be required within GEOS-Chem.
Key similarities to GEOS-5
The following aspects of MERRA are identical to the GEOS-5 operational data product.
- The MERRA high-resolution grid (0.5° x 0.666°) is identical to that of GEOS-5.
- The MERRAvertical grid (72 hybrid levels) is identical to that of GEOS-5.
- By and large, the same data fields that were in the GEOS-5 product are also contained in MERRA.
Key differences from GEOS-5
The following aspects of MERRA differ from the GEOS-5 operational data product:
- The MERRA temporal resolution differs from GEOS-5 as follows:
- Surface data: hourly averages (MERRA) vs. 3-hour averages (GEOS-5)
- 3-D data: 3-hour averages (MERRA) vs. 6-hour averages (GEOS-5)
- Differences in individual MERRA met fields:
- Most of the releveant MERRA 3-D data for GEOS-Chem are only available at 1° x 1.25° horizontal resolution.
- OMEGA, PV, RH are 6-hr instantaneous data in MERRA but are 6-hr averaged data in GEOS-5
- SNOMAS (snow mass) has units of kg/m2 in MERRA but mm in GEOS-5
- Many of the precipitation fields (e.g. DQIDTMST, DQVDTMST, DQLDTMST, etc.) are saved on a 1.25° x 1.25° x 42 pressure levels instead of 72 hybrid levels. P -> hybrid interpolation is necessary.
- The MERRA fields FRLAKE, FRLAND, FRLANDIC, FROCEAN, PHIS are time-invariant.
Required coding changes in GEOS-Chem
As noted in the section above, the similarities between GEOS-5 and MERRA (especially of the vertical grids) means that in most cases we can just treat MERRA in the same way as for GEOS-5. In many areas of the code, the coding changes are straightforward and simply involve an extension of C-preprocessor statements such as:
#if defined( GEOS_5 )
#if defined( GEOS_5 ) || defined( MERRA )
In other areas of GEOS-Chem the changes are more involved. Below we provide a description of the areas of GEOS-Chem that were modified for MERRA:
The following modifications for MERRA were made in the Headers/ directory:
The following modifications for MERRA were made in the GeosUtil/ directory:
The following modifications for MERRA were made in the GeosCore/ directory:
The following modifications for MERRA were made in the GeosTomas/ directory:
--Bob Y. 11:18, 25 August 2010 (EDT)
For the time being, GEOS-Chem/MERRA will read in the same the OTD/LIS redistribution factors as for GEOS-5. MERRA is based on GEOS-5.2.0, so this should be OK for a first-cut.
However, for proper treatment, the OTD/LIS redistribution factors for MERRA will be rescaled for the entire OTD/LIS observation period. This cannot be done until we have processed and regridded several years of MERRA met data.
Wet deposition and scavenging
We are currently looking into how some of the MERRA precipitation fields can be used to improve the existing GEOS-Chem wet scavenging algorithm.
Acetone ocean source
As a first cut, the ocean source of acetone for MERRA (computed in routine OCEAN_SOURCE_ACET in acetone_mod.f) is computed in the same way as for GEOS-5. The following #if block was based on the GEOS-5 code:
#elif defined( MERRA ) !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% !%%% NOTE: Treat MERRA in the same way as GEOS-5, although %%% !%%% we may want to rescale this in the future (bmy, 8/13/10) %%% !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% #if defined( GRID4x5 ) REAL*8, PARAMETER :: SCALE_FACTOR = 0.9551d0 ! MERRA 4x5 #elif defined( GRID2x25 ) REAL*8, PARAMETER :: SCALE_FACTOR = 0.2212d0 ! MERRA 2x25 #endif
At some point we may need to rescale this. However this cannot be done until we have sufficient MERRA met fields processed and regridded.
--Bob Y. 15:03, 16 August 2010 (EDT)