search

NCAR / hrldas

HRLDAS (High Resolution Land Data Assimilation System)
Other
50 stars 53 forks source link

Update ERA5-Land pre-processing code #121

Closed CharlesZheZhang closed 1 year ago

CharlesZheZhang commented 1 year ago

We sincerely thank our users and developers for contributing code to pre-process ERA5-Land forcing and in the discussions: (https://github.com/NCAR/hrldas/discussions/111)

Alessandro Anav (alessandro.anav@enea.it) from Italian National Agency for New Technologies, Energy and Sustainable Economic Development, Ehsan Jalilvand (ehsanj@msu.edu) from Michigan State University, Tzu-Shun Lin (tslin2@ucar.edu) and Zhe Zhang (zhezhang@ucar.edu) from NCAR.

Please see /hrldas/tutorialNote5_2D_ERA5-Land_forcing.ipynb for a jupyter notebook on how to pre-process and run hrldas with ERA5-Land data step-by-step, also see /hrldas/docs/README.ERA5 for documentations. (not ERA5 atmoshpere https://github.com/NCAR/hrldas/discussions/88)

Prepare_ERA5-Land_Forcing.sh: this extracting program was originally written with python, while users might have encountered issues with discrepancies between python version 2 and 3 https://github.com/NCAR/hrldas/discussions/101. Thus, we decide to adapt to a shell code Prepare_ERA5-Land_Forcing.sh provided by Alessandro Anav (alessandro.anav@enea.it), that is also faster and simpler in extracting ERA5-land forcing data in one step. We really appreciate the user's contribution!

extract_era5_init.py: extract the initial-only fields (SKT, SD, STL1-4, SWV1-4) from the ERA5-Land GRIB output files is still used in python. Note: in the ERA5-Land data documentation (https://codes.ecmwf.int/grib/param-db/?id=141), the variable snow depth water equivalent (SWE) is with a short name sd. It is easy to get confused with another variable snow depth that has a short name sde. (https://github.com/NCAR/hrldas/discussions/111)

cenlinhe commented 1 year ago

@CharlesZheZhang maybe you can also put the discussion # and link in your description above so that people can cross-check them.

cenlinhe commented 1 year ago

@CharlesZheZhang maybe also rename the "Note5_2D_ERA5-Land_forcing.ipynb" to "Note5_Regional_modeling_ERA5-Land_forcing.ipynb" since you also included model simulation instructions in this file. It is not limited to forcing pre-processing.

CharlesZheZhang commented 1 year ago

Update documentation on soil thickness, now the latest HRLDAS code can cope with different soil layer thickness, so users don't need to change them in the namelist.

cenlinhe commented 1 year ago

Looks good to me. I approve it. @tslin2 What do you think?

CharlesZheZhang commented 1 year ago

Thank you all for your time to review this update. I will merge it now.