Skip to content

Compute aerosol properties maps from ECMWF IFS model data

License

Notifications You must be signed in to change notification settings

AugustinMortier/aer-ifs

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

88 Commits
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

aer-IFS

Compute aerosol optical properties maps based on ECMWF IFS forecasts.

  • Lidar Ratio (LR)
  • Mass to Extinction Coefficient (MEC)

LR computation methodology:

  • We get the IFS forecasts data of aerosol optical depth forecasts for the different chemical components.
  • We get the IFS forecasts data of relative humidity.
  • For each aerosol species, we calculate a LR weighted by the relative optical depth contribution for different RH values.
  • The total LR is determined as the sum of the weighted LR of the different aerosols species.

Individual aerosol LR are taken from the IFS 49R1 input configuration file. For each species, a specific column is selected, which might correspond to a specific range bins or according to the given literature reference.

example

aer-ifs --date 2024-10-15

LR

IFS LR - 2024-10-15 - 1064 nm

MEC

IFS MEC - 2024-10-15 - 1064 nm

get started

1. clone repo

git clone https://github.com/AugustinMortier/aer-ifs.git

2. install

  • via poetry
poetry install
  • via pip/pipx
pip install .

Note

You can also install lr-ifs with this one liner:

pip install "git+ssh://[email protected]/AugustinMortier/aer-ifs.git"

how to use

Compute LR and MEC for the 2024-09-26

aer-ifs --date 2024-09-26

This will create into output_path (default: ./data/):

  • {yyyy}/{mm}/lr_ifs-{yyyymmdd}.nc: netcdf file which contains the computed LR and MEC.
  • {yyyy}/{mm}/lr_ifs-{yyyymmdd}.json (if aprofiles option enabled (default)): json file which contains, for each E-PROFILE station available and for the selected day, the corresponding IFS-LR and IFS-MEC.

limitations

  • Only one value per day is computed at the moment and corresponds to 00:00:00Z.
  • Only one value per column is computed at the moment. Only the surface RH is considered for the calculation.