Wet-bulb temperature

Usage

theta_w(hus, pa, ta, nthreads = 1)

Arguments

hus: specific humidity (kg/kg).
pa: air pressure (hPa).
ta: air temperate (K).
nthreads: an integer specifying the number of threads to use for computation (OpenMP). Default is 1.

Value

the computed wet-bulb temperature (K).

Note

air pressure variable should match the specific humidity correctly. If hus is a three dimensional array on pressure levels [lon, lat, plev], pa values should the same on each [lon, lat] slice. See convert_plev_to_array.

In case you are using model levels, estimate pressure array beforehand.

References

Stull, R., 2011: Wet-Bulb Temperature from Relative Humidity and Air Temperature. J. Appl. Meteor. Climatol., 50, 2267–2269, doi:10.1175/JAMC-D-11-0143.1 .

Examples

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