FINN2CMAQ Processor

Overview

This processor make Fire Emissions ready for CMAQ.

FINN makes global predictions of fire emissions that are quickly available as a part of the WACMM system. This processor takes the WACMM FINN input and converts it to a 2D daily ioapi file. It then converts the 2D daily file into a 3D hourly file for use in CMAQ. It should be noted that this 3D product is a very simplistic version of fires.

The processor can work with retrospective tar.gz files from bai.acom.ucar.edu or daily near-real-time tz files from www.acom.ucar.edu.

Download a daily text file (./scripts/get.py)
Translate to CMAQ input (./scripts/txt2daily.py)
Add LAY and hourly (./daily2hourly3d.py)

Example:

See examples/FINN2CMAQ.ipynb

Configuration Notes

daily2hourly3d.py has three notable configuration options:

--tpropath specifies how all fires will be allocated in time.
--layerpath specifies how all fires will be allocated to model layers.
--exprpath specifies how FINN file columns will be translated to CMAQ-ready species.

For more details, see run ./scripts/daily2hourly3d.py --help

Notable Updates

2022-09-16
- Added support for PMC from v2 FINN files.
- Removed run_retro.sh because the download part should be manual because the path structures change too much.
- Changed default hourly temporal allocation from 210 to 1008. Select either by moving to the first position in aux/tpro.txt
- Modified aux/layerfrac.csv to extend to 44 layers. The EPA 35 and 44 layer structures are very similar in the first 19 levels, so you can use the aux/layerfrac.csv for either by removing the number of layers you are not using.
2022-08-04
- Versions have been added to txt2daily.py and daily2hourly3d.py scripts. (0.1.0)
- Histories including versions have been added to outputs.
- Additional documentation has been added to functions.
2021-06-24
- An example Jupyter Notebook has been added for testing purposes.
- The PM speciation has been updated to use g/s definitions.
- The scripts/daily2hourly3d.py has been updated to output a CMAQ-ready unit and to have the option to process individual days.

To Do

Add speciation translations for MOZART or SAPRC99. This would provide access to different explicit species (e.g., HONO).

Prerequisites

Command-Line (i.e, Linux, Mac)
Python3
- pandas
- numpy
- pyproj
- PseudoNetCDF

Annotated Directory Structure

.
|-- README.md
|   # This file
|-- run_nrt.sh
|   # A wrapper to run 1 day for the near-real-time FINN
|-- examples
|   `-- FINN2CMAQ.ipynb
|-- scripts
|   |-- get_retro.py
|   |-- get_nrt.py
|   |-- txt2daily.py
|   `-- daily2hourly3d.py
|-- aux
|   |-- README.md
|   |   # More details about aux files
|   |-- gc2_to_cb6r3_ae7.txt
|   |   # speciation definition for GEOS-Chem in v1.5
|   |-- gc2_to_cb6r3_ae7_v2.txt
|   |   # speciation definition for GEOS-Chem in v2
|   |-- GRIDDESC
|   |   # example grid definitions
|   |-- layerfrac.csv
|   |   # example layer fractions
|   `-- tpro.txt
|       # example hour fractions
|-- daily
|   `-- %Y
|       |-- FINNv1.5_%Y.GEOSCHEM.nc
|       |   # file and folder only exist if run_retro.sh has been run
|       `-- FINNv1.5_%Y-%m-%d.GEOSCHEM.NRT.108NHEMI2.nc
|           # file and folder only exist if run_nrt.sh has been run
|-- hourly
|   `-- %Y
|       `-- %m
|           |-- FINNv1.5_%Y.CB6r3.3D.%Y-%m-%d.nc
|           |   # file and folders only exist if run_retro.sh have been run
|           `-- FINNv1.5_%Y.CB6r3.3D.%Y-%m-%d.nc
|               # file and folders only exist if run_nrt.sh have been run
`-- www.acom.ucar.edu
    |-- Data
    |   `-- fire
    |       `-- data
    |           `-- finn1
    |               `-- FINNv1.5_2016.GEOSCHEM.tar.gz
    |                   # file and folders only exist if run_retro.sh has been run
    `-- acresp
        `-- MODELING
            `-- finn_emis_txt
                `-- GLOB_GEOSchem_%Y%j.txt.gz
                    # file and folders only exist if run_nrt.sh has been run

barronh / finn2cmaq