README -- Contents of ann_mean_trop/
Bob Yantosca
bmy@io.harvard.edu
20 Feb 2004


Files:
------
ann_mean_trop.geos1.2x25 : Annual mean tropopause heights and flux flags
ann_mean_trop.geos3.2x25 :  for GEOS-1, GEOS-3, GEOS-STRAT, and GEOS-4 grids
ann_mean_trop.geoss.2x25 :  at 2 x 2.5 horizontal resolution.
ann_mean_trop.geos4.2x25 :

Files are in binary punch format.


Source:
-------

For GEOS-1, GEOS-STRAT, and GEOS-3:
===================================
To create the "ann_mean_trop.geos*.2x25" files, you need to first run the
GEOS-CHEM code for a full year (or for as long as you have met fields on
disk, if less than one year), and save out the ND55 diagnostic to the binary
punch file. 

Then you can use the IDL code

	~bmy/IDL/regrid/create_amt.pro

which will create the "ann_mean_trop.geos*.2x25" files.  Tracer 1 in this file
is the annual mean tropopause height in levels, and tracer 2 is the flux
diagnostic flags (which determine if the Strat/Trop exchange is coming from
the top, left, right, or a combination of these).


For GEOS_4:
===========
The GEOS-4 annual mean tropopause were computed as follows:
  
(a) Bob Yantosca downloaded the *mis_x* HDF-EOS files from GMAO's machine 
    mintz.gsfc.nasa.gov for each day of 2003.  These files contain the 
    TROPP (tropopause pressure in hPa0 met field. 

(b) Routine ~bmy/IDL/regrid/create_geos4_tropp.pro was used to read
    the HDF-EOS files and save them to BPCH format.  

(c) Routine ~bmy/IDL/regrid/create_geos4_amt.pro was used to create
    an annual mean tropopause based on the surface pressure and
    tropopause pressure for the year 2003.  This should be OK for
    other years as well...it's an average.


Notes:
------
(1) The "ann_mean_trop.geos*.2x25" files are read into GEOS_CHEM by routine
    "read_tropopause.f".