Feature/aqm canopy

[drnimbusrain]

@BrianCurtis-NOAA @JianpingHuang-NOAA @KaiWang-NOAA

This PR adds the inline canopy photolysis effects to AQM/Online-CMAQ, and references issue #16 , which provides some more information on the canopy photolysis effect.

Note that this is only an option and does not change results from the base configuration. When turned on, this option also only affects grid cells that are satisfy contiguous canopy conditions in the code.

Canopy inputs are added to the aqm.rc like any other external input sources. Required variables are the forest canopy height (FCH), forest fraction (FRT), clumping index (CLU), population density (POPU), leaf area index from ECCC (LAIE), and four in-canopy level specific cumulative LAI fractions (C1R, C2R, C3R, and C4R). I have created these canopy inputs based on collaboration with ECCC. The differences are both AQM and CMAQ specific.

I have created the canopy input files for Online-CMAQ at the RRFS_CONUScompact_13km and NA13km domains here on Hera: /scratch2/NCEPDEV/naqfc/RRFS_CMAQ/canopy/

Here are the example canopy additions to the aqm.rc template that should be added and then populated based on updated sample configuration files in UFS-SRW-App/online-cmaq for EMCs testing.

canopy_yn: true

canopy_type: canopy

canopy_format: netcdf

canopy_file: {{ aqm_rc_canopy_file_fp }}

canopy_frequency: static

canopy_species::
FCH 1.00000 FCH m
FRT 1.00000 FRT 1
CLU 1.00000 CLU 1
POPU 1.00000 POPU 10000_people/10km2
LAIE 1.00000 LAIE 1
C1R 1.00000 C1R 1
C2R 1.00000 C2R 1
C3R 1.00000 C3R 1
C4R 1.00000 C4R 1
::

These effects have been tested in both the offline NAQFC ("NACC-CMAQ") and Online-CMAQ for the RRFS_CONUScompact_13km domain during August 2019 on Hera. Results show improved model performance and reduced ozone overpredictions in regions of contiguous canopies. There are negligible impacts on summer PM2.5. Please see my DSRA/EMC presentation for results.

Also, I have tested the canopy photolysis codes in the latest workflow with latest EMC/AQM codes, and ran the model for a 6-hr engineering test over the target larger NA13km domain on WCOSS2. Engineering tests of the photolysis effect over the NA domain seem OK:

@chan-hoo If/when this PR is accepted, can we rework the UFS-SRW-App to eventually remove the --canopy option and the canopy specific workflow changes/options as the integrated/default? Thank you.

[BrianCurtis-NOAA]

@drnimbusrain If there's still work ongoing, it might be helpful to change this to a draft PR until the work is ready.