Adding Fengsha dust emission scheme

[bbakernoaa]

This PR add the fengsha windblown dust emission model. Note that this is only an option and does not change results from the base configuration. Inputs are added to the aqm.rc like any other emission sources. Required variables are the clay and sand fraction, threshold friction velocity, and drag partition. I have created the inputs based on the latest inputs from UFS-Aerosols which is running the same flavor of FENGSHA. The only differences are CMAQ specific.

I have created inputs for Online-CMAQ at the CONUS 25 and 13 km grids as well as the new large 13 km domain here: /scratch1/RDARCH/rda-arl-gpu/Barry.Baker/emissions/nexus/FENGSHA/

Example aqm.rc additions:

#
# FENGSHA
#
fengsha_yn: true

fengsha_type: fengsha

fengsha_format: netcdf

fengsha_file: /scratch1/RDARCH/rda-arl-gpu/Barry.Baker/emissions/nexus/FENGSHA/test_fengsha.nc

fengsha_frequency: monthly

fengsha_species::
    CLAYF       1.00000   clayfrac     1
    SANDF       1.00000   sandfrac     1
    DRAG        1.00000   albedo_drag  1
    UTHR        1.00000   uthres       1
::

[BrianCurtis-NOAA]

A brief read through the code changes, and I'm curious if you know that anyone uses the makefiles any more. If not, with CMake we shouldn't need to keep those Makefiles in the repository any more.

[BrianCurtis-NOAA]

Talked with Raffaele and agreed removing Makefiles not a good idea at this time.

[BrianCurtis-NOAA]

@bbakernoaa a couple questions:
Will Fengsha be used in the operational implementation?
Have you been testing with CCPP_SUITES: FV3_GFS_V16 or FV3_GFS_v15p2 or both?

Potential concern:
It looks like fensgha is enabled by default in aqm_config_mod?

AQM/src/shr/aqm_config_mod.F90

Line 38 in 058973c

logical :: fengsha_yn = .true.

I believe this will break the current RT in UFS. Depending on your answers to the questions we can figure out how to get this into UFS properly.

[BrianCurtis-NOAA]

@bbakernoaa I've created a UFSWM branch here: https://github.com/BrianCurtis-NOAA/ufs-weather-model/tree/feature/fengsha2 I will need to get an official FENGSHA_FILES.nc from you for the date of the RT:

export SYEAR=2019
export SMONTH=08
export SDAY=1
export SHOUR=12

I looked in your Hera dir you mentioned earlier and I did not see FENGSHA_FILES.nc to use.

If you are using FV3_GFS_v16 then I will need to get a run_dir from you to get the updated input files for that version. Then we can create a new RT for that configuration.

[bbakernoaa]

Brian these are agnostic from the physics anything that works with AQM will work with with FENGSHA. It will not change with different physics routines as they all provide the necessary inputs. What grid is the RT on?

On Fri, Aug 19, 2022 at 3:55 PM Brian Curtis ***@***.***> wrote: @bbakernoaa <https://github.com/bbakernoaa> I've created a UFSWM branch here: https://github.com/BrianCurtis-NOAA/ufs-weather-model/tree/feature/fengsha2 I will need to get an official FENGSHA_FILES.nc from you for the date of the RT: export SYEAR=2019 export SMONTH=08 export SDAY=1 export SHOUR=12 I looked in your Hera dir you mentioned earlier and I did not see FENGSHA_FILES.nc to use. If you are using FV3_GFS_v16 then I will need to get a run_dir from you to get the updated input files for that version. Then we can create a new RT for that configuration. — Reply to this email directly, view it on GitHub <#6 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AFIUVN5ZYJZMKZJFQ5W45WTVZ7RANANCNFSM56WAPFFQ> . You are receiving this because you were mentioned.Message ID: ***@***.***>

-- Barry Baker National Oceanic and Atmospheric Administration Air Resources Laboratory Physical Research Scientist Chemical Modeling and Emissions Group Leader NCWCP, R/ARL, Rm. 4204 5830 University Research Court College Park, Maryland 20740 Phone: ‪(301) 683-1395‬

[BrianCurtis-NOAA]

Brian these are agnostic from the physics anything that works with AQM will work with with FENGSHA. It will not change with different physics routines as they all provide the necessary inputs. What grid is the RT on?

export OUTPUT_GRID='lambert_conformal'

[bbakernoaa]

@BrianCurtis-NOAA This isn't quite the question I asked. The input fields need to be created on the model grid (with halo=0).

Also note that the fengsha input fields do not exist for a specific date. Rather they are a climatology. There isn't a file for a specific date for this emission input as the emissions are weather driven unlike the anthropogenic emissions.

[BrianCurtis-NOAA]

@bbakernoaa Can you update this with the latest AQM develop branch code?

[bbakernoaa]

@BrianCurtis-NOAA What is the status of this?

[BrianCurtis-NOAA]

Its close to merge status. I'm out but I will try late tonight if a merge is ready.