README -- describes contents of HEMCO/OH/v2022-11/v5-07-08 directory
01 Dec 2022
GEOS-Chem Support Team
geos-chem-support@g.harvard.edu


Files:
===============================================================================

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%  NOTE: The HEMCO/OH/v2014-09/v5-07-08/OH_3Dglobal.geos5*nc files       %%%
%%%  have time:units = "hours since 1985-1-1 00:00:0.0".  The 0.0 causes   %%%
%%%  GCHP simulations to die with error.  We have now edited time:units    %%%
%%%  for compatibility with GCHP using the following commands:             %%%
%%%                                                                        %%%
%%%  for f in ./*.nc; do                                                   %%% 
%%%    ncatted -a units,time,o,c,"hours since 1985-01-01 00:00:00" $f      %%%
%%%    nc_chunk.pl $f 1                                                    %%%
%%%  done                                                                  %%%
%%%                                                                        %%%
%%%   -- Bob Yantosca (01 Dec 2022)                                        %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

(1) OH_3Dglobal.geos5.47L.2x25.nc  : 3D OH fields for 2x2.5 GEOS-5  grid (47L)
(2) OH_3Dglobal.geos5.47L.4x5.nc   : 3D OH fields for 4x5   GEOS-5  grid (47L)
(3) OH_3Dglobal.geos5.72L.2x25.nc  : 3D OH fields for 2x2.5 GEOS-5  grid (72L)
(4) OH_3Dglobal.geos5.72L.4x5.nc   : 3D OH fields for 4x5   GEOS-5  grid (72L)
(5) OH_3Dglobal.geosfp.47L.2x25.nc : Links to OH_3Dglobal.geos5.47L.2x25.nc
(6) OH_3Dglobal.geosfp.47L.4x5.nc  : Links to OH_3Dglobal.geos5.47L.4x5.nc
(7) OH_3Dglobal.geosfp.72L.2x25.nc : Links to OH_3Dglobal.geos5.72L.2x25.nc
(8) OH_3Dglobal.geosfp.72L.4x5.nc  : Links to OH_3Dglobal.geos5.72L.4x5.nc

 -- For all files:
    Units       : kg/m3
    Timestamp   : 1985/01 thru 1985/12
    Compression : Level 1 (e.g. nccopy -d1)

###############################################################################
#####  NOTE: For most simulations, you can specify that you want to read  #####
#####  OH data from the OH_3Dglobal.geos5.72L.2x25.nc file, and HEMCO     #####
#####  will perform the vertical and horizontal regridding.               #####
#####                                                                     #####
#####  We will keep the OH files defined on other other vertical and      #####
#####  horizontal grids for backwards compatibility.                      #####
#####                                                                     #####
#####    -- Bob Yantosca (12 Mar 2015)                                    #####
###############################################################################


Source:
===============================================================================

The original bpch files used to create these files were from the
GEOS_MEAN/OHmerge/v5-07-08/ data directory. The bpch files were converted to
netCDF using GAMAP routine bpch2coards.pro and the data was converted from
molec/cm3 to kg/m3 using NCL4GC/hemco/hco_adj_global_oh.ncl.


From the original README file qin GEOS_MEAN/OHmerge/v5-07-08:

  The monthly mean OH data was taken from the GEOS-CHEM v5-07-08 1-year 
  benchmark run using GEOS-3 winds at 4 x 5 resolution for model year 2001.  
  This run was submitted by Mat Evans and was completed in September 2003.  
  
  The punch files from the v5-07-08 benchmark run are stored at Harvard at:
  
  (1) /data/ctm/GEOS_MEAN/2001/v5-07-08.mje/ctm.bpch.v5-07-08.Jan-Mar-2001
  (2) /data/ctm/GEOS_MEAN/2001/v5-07-08.mje/ctm.bpch.v5-07-08.Apr-Jun-2001
  (3) /data/ctm/GEOS_MEAN/2001/v5-07-08.mje/ctm.bpch.v5-07-08.Jul-Sep-2001
  (4) /data/ctm/GEOS_MEAN/2001/v5-07-08.mje/ctm.bpch.v5-07-08.Oct-Dec-2001
  
  Tropospheric OH data (computed by the SMVGEAR II solver) is stored in these 
  punch files under GAMAP diagnostic category "CHEM-L=$" and tracer=1.  The OH 
  data is saved from levels 1-20, since the SMVGEAR II solver cannot operate
  in the stratosphere. 


How the netCDF files were created:
===============================================================================

The files in this directory were created from the original data files in
binary punch format.  NetCDF files were generated from the bpch files with
GAMAP routine BPCH2COARDS.

Bob Yantosca wrote the following script to adjust the netCDF attributes for 
better COARDS compliance.

   #!/bin/tcsh
   
   # Fixes attributes
   
   ncatted -a axis,time,o,c,"T" $1
   
   ncatted -a long_name,lev,o,c,"GEOS-Chem levels" $1
   ncatted -a axis,lev,o,c,"Z" $1
   ncatted -a positive,lev,o,c,"up" $1
   
   ncatted -a axis,lat,o,c,"Y" $1
   
   ncatted -a axis,lon,o,c,"X" $1
   
   ncatted -a References,global,o,c,"www.geos-chem.org; wiki.geos-chem.org" $1
   ncatted -a Format,global,o,c,"NetCDF-3" $1
   
   exit(0)

NOTE: Using "GEOS-Chem levels" for the level dimension will allow HEMCO to
perform the vertical regridding.