NCAR / hrldas

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

correct specific humidity magnitude for ERA5-Land forcing #204

Closed CharlesZheZhang closed 3 months ago

CharlesZheZhang commented 3 months ago

This bug was reported in this issue: https://github.com/NCAR/hrldas/issues/203

This pull request is to correct the specific humidity magnitude in preprocessing ERA5-Land data in bash code https://github.com/NCAR/hrldas/blob/master/hrldas/HRLDAS_forcing/run/examples/ERA5/Prepare_ERA5-Land_Forcing.sh https://github.com/NCAR/hrldas/blob/087f1fd78601b37eaeead93da9f30253441f7cc2/hrldas/HRLDAS_forcing/run/examples/ERA5/Prepare_ERA5-Land_Forcing.sh#L165

According to Bolton 1980: https://archive.eol.ucar.edu/projects/ceop/dm/documents/refdata_report/eqns.html The output specific humidity unit is already in [kg/kg], thus no need to divide it again by 1000.

cenlinhe commented 3 months ago

According to https://github.com/NCAR/hrldas/issues/203, the range of specific humidity is about 0~0.025 [kg/kg] after this fix, which is typical range for global values: https://geography.name/how-does-specific-humidity-vary-globally-and-seasonally/