Dew Point Temperature from Specific Humidity

Usage

td_from_hus(hus, pa, nthreads = 1)

Arguments

hus: specific humidity [lon, lat, lev, time] (kg/kg).
pa: air pressure [lon, lat, lev, time](hPa).
nthreads: an integer specifying the number of threads to use for computation (OpenMP). Default is 1.

Value

the computed dew point temperature [lon, lat, lev, time] (K).

Note

air pressure variable should match the specific humidity correctly. If hus is a three dimensional array on pressure levels [lon, lat, lev, time], pa values should the same on each [lon, lat,,time] slice. In case you are using model levels, estimate pressure array beforehand.

Examples


data("ERA5_pl")
hus <- ERA5_pl$hus
hus_d <- dim(hus)
pa <- ERA5_pl$plev / 100 # convert to hPa
pa <- convert_plev_to_array(pa, hus_d)
td <- td_from_hus(hus, pa, nthreads = 1)