README
M. Sulprizio (Aug 20, 2026)

Overview
========

This folder contains scripts needed to crop global GEOS-FP meteorology fields
to regional domains supported by GEOS-Chem (geos-chem.org) and the IMI
(carboninversion.com).

GEOS-Chem and the IMI support regions within the following domains:

    Africa (AF)        : 37°S-40°N, 20°W-53°E

    Asia (AS)          : 11°S-55°N, 60°E-150°E

    Europe (EU)        : 33°N-61°N, 30°W-70°E

    Middle East (ME)   : 12°N-44°N, 20°W-70°E

    North America (NA) : 10°N-70°N, 140°W-40°W

    South America (SA) : 59°S-16°N, 88°W-31°W

    Oceania (OC)       : 50°S-5°N, 110°E-180°E 

    Russia (RU)        : 41°N-83°N, 19°E-180°E

These are the regional windows for which pre-cut meteorological files are
available at 0.25° x 0.3125° resolution. You may utilize other regions, but
this requires either using global meteorological fields which can be
computationally expensive (not recommended) or by cropping global meteorological
fields via a pre-processing step.

To facilite cropping global meteorological fields at at 0.25° x 0.3125°
resolution, a sample script (crop_met.sh) is included here. This script
utilizes the Climate Data Operators (CDO). It also includes an option to first
download global meteorological fields at 0.25° x 0.3125° resolution. The global
files are large (approx. 300G per month), so when using that option it is
recommend that you process short periods at a time and delete the global files
before processing additional periods.

Instructions
============

In a text editor, modify the user settings section in crop_met.sh.

   - RegionID          : Two-letter string to identify region (e.g. SA for South
                         America). See regions listed above.
   - LonMin            : Minimum longitude edge of the region of interest.
   - LonMax            : Maximum longitude edge of the region of interest.
   - LatMin            : Minimum latitude edge of the region of interest.
   - LatMax            : Maximum latitude edge of the region of interest.
   - DownloadGlobalMet : Boolean for downloading global 0.25° x 0.3125°
                         meteorology fields for cropping. Default is true.
                         Set to false if you already have these files on your
                         system.
   - RemoveGlobalMet   : Boolean for deleting global meteorology files after
                         cropping. Default is false. Set to true once you are
		       	 sure the cropped meteorology files are working properly
                         in the model and you do not plan to generate cropped
	                 files for additional regions.
   - InDir             : Directory containing the global high-resolution
                         meteorology fields.
   - OutDir            : Directory where the cropped meteorology files will be
                         placed.

The cropped meteorology files can be generated by then executing ./crop_met.sh
at the command line or submitting the script to your cluster’s scheduler if
available. Headers for the SLURM scheduler are included at the top of the
script, but you can modify or remove those as needed.