aqmbc

Air Quality Model Boundary Condition (aqmbc) Tools

Air Quality Model Boundary Conditions is a tool to create time and space boundary conditions for an Air Quality Model from one or many data sources.

aqmbc is a python tool to make boundary conditions for CMAQ.

The key value of aqmbc is to provide a single interface to multiple models. This makes maintaining a system easy to acheive and hopefully easy to use.

Currently, aqmbc supports:

• GEOS-Chem v12 and v14
• CMAQ (Hemispheric or Regional)
• NASA GMAO's GEOS-CF
• RAQMS

To do:

• aqmbc would like to support WACCM, AM4, and CAMS
• A future version will allow for explicit output of CAMx readable boundary conditions. For now, use the cmaq2camx tools to convert these files.

Documentation

https://barronh.github.io/aqmbc

Installation

pip install git+https://github.com/barronh/aqmbc.git

or

pip install https://github.com/barronh/aqmbc/archive/refs/heads/main.zip

Example

# Create a working directory with standard configurations
python -m aqmbc --template examples
cd examples
# Edit a standard configuration file to point to your inputs/outputs
# - gcncv12.cfg
# - gcnd49v11.cfg
# - raqms.cfg
# Or make your own cfg
python -m aqmbc gcncv12.cfg
# Right now, the raqms reader/interpolator must be defined at runtime
# If using raqms, run python raqms.py

This process makes files that mirror the global air quality model structure. However, CMAQ requires hourly lateral boundary condition files with instantaneous values for every hour edges (24/day + 1).

It is often convenient to concatenate all days in a month using ncrcat or other tools. Below is an example using PseudoNetCDF.

import PseudoNetCDF as pnc
import glob

paths = (
    sorted(glob.glob('BCON/*2016-01-*.12US1.BCON.nc'))
    + sorted(glob.glob('BCON/*2016-02-01*.12US1.BCON.nc'))
)
mf = pnc.pncmfopen(paths, format='ioapi', stackdim='TSTEP')
mf.save('BCON_2016-01.12US1.BCON.nc')

If your data has fewer than hourly values, you will also need to add interpolation. Currently, these steps are not supported because the possibilities are too varied. I recommend xarray for this type of processing.

How to Contribute

• If you add a new set of definitions (new version of GEOS-Chem or a new model), consider opening an issue on github and posting it there. It would be great to include these in future versions.
• If you find a bug, open an issue.

Prerequisites

• PseudoNetCDF (github.com/barronh/PseudoNetCDF).
• pyproj