-
Notifications
You must be signed in to change notification settings - Fork 53
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Update ERA5-Land pre-processing code #121
Update ERA5-Land pre-processing code #121
Conversation
@CharlesZheZhang maybe you can also put the discussion # and link in your description above so that people can cross-check them. |
hrldas/docs/README.ERA5
Outdated
@@ -238,6 +232,14 @@ See the full documentation for a complete description of items, but these are th | |||
|
|||
ZLVL = forcing height | |||
|
|||
For ERA5-Land data, will also need to change the four layer depth, as they are different in ERA5-Land data: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think I fixed this bug in the latest HRLDAS (ba09bf0). Users can still use the original Noah-MP soil thickness and the ERA-5 thickness will be interpolated into Noah-MP values. Could you please do a quick test of this?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks Cenlin, I did test the model can run with the original soil thickness (0.1, 0.3, 0.6 and 1.0). Will modify this part in the documentation.
@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. |
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. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Some typos found (NLDAS2) in the tutorial/Note5_Regional_modeling_ERA5-Land_forcing.ipynb
- Generate geo_em domain file
This step is the same as previous tutorial with NLDAS2 forcing
3.3.2 Run the create_forcing.exe
After modifying the namelist for create_forcing.exe, we can now start to generate NLDAS2 forcing data for our 2D domain:
- Running CONUS domain with NLDAS2 forcing
! YEND = 203 | ||
|
||
NSOIL=4 | ||
soil_thick_input(1) = 0.07 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
you may also want to change this back to Noah-MP soil thickness
Looks good to me. I approve it. @tslin2 What do you think? |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me. I approve it.
Thank you all for your time to review this update. I will merge it now. |
We sincerely thank our users and developers for contributing code to pre-process ERA5-Land forcing and in the 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 #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 #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. (#111)