Skip to content

LBC Tutorials

Jump to bottom
barronh edited this page Dec 18, 2017 · 4 revisions

Boundary Conditions from GEOS-Chem ND49 files

Full Process described on GEOS-Chem wiki

Steps:

  1. Run GEOS-Chem with ND49 enabled with all chemical tracers, pressure, and temperature
  2. Prepare CMAQ MCIP inputs
  3. Download Mapping files from GEOS-Chem to CMAQ mechanisms
  4. Run pncglobal2cmaq.py

Boundary Conditions from GEOS-Chem Benchmark

Making boundary conditions from GEOS-Chem Benchmark simulations is easy and accessible, though not ideal.

In Issue #7, I was asked to create functionality to produce lateral boundary conditions from GEOS-Chem Benchmarks. These benchmarks are run for a single year (w/ 1-year spinup) and so are only year-specific for the benchmark year, and only output monthly averages. Ideally, boundary conditions would come from day-specific ND49 results. That being said, representative month-specific inputs would be better than a constant and easy to obtain.

Steps:

  1. Download Mapping files from GEOS-Chem to CMAQ mechanisms
  2. Download GEOS-Chem Benchmark
  3. Run pncglobal2cmaq.py

Download mappings.json and update

Download mappings.json or request one -- finding the right home for downloading these now.

This requires one change to the typical mapping JSON file. Replace PSURF in AIRMASSDEN and AIRDEN definitions with (hyam[:].reshape(1, -1).T + hybm[:].reshape(1, -1).T * PSURF[:][:, [0]].T).T as in AIRMOLDEN.

Download GEOS-Chem benchmark

wget ftp://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-01/v11-01-public-release/Run0/bpch.tar.gz
wget ftp://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-01/v11-01-public-release/Run0/tracerinfo.dat
wget ftp://ftp.as.harvard.edu/gcgrid/geos-chem/1yr_benchmarks/v11-01/v11-01-public-release/Run0/diaginfo.dat

tar xzf bpch.tar.gz

Run pncglobal2cmaq.py

The command below requires METCRO3D, METBDY3D and ctm.bpch.v11-01-public-Run0.2013*00 will work. METCRO3D and METBDY3D are meteorological inputs to CMAQ and must have the VGLVLS property. No meteorological variables are required. ctm.bpch* files are from the benchmark.

pncglobal2cmaq.py --sigma --stack=time -f benchmark ctm.bpch.v11-01-public-Run0.2013*00 --mapping=mapping.json --METCRO3D /path/to/cmaq/data/mcip/METCRO3D_Benchmark --METBDY3D /path/to/cmaq/data/mcip/METBDY3D_Benchmark
Clone this wiki locally