Rewrote rh2qair to account for air pressure. Now true inverse of qair2rh

[ankurdesai]

Still need to change units to be same for both functions and met2CF Ameriflux/CSV/ALMA and temporal downscaling need to pass pressure; two of these functions don’t read in pressure currently. qair2rh is called in met2model for BIOCRO/SIPNET and temporal downscaling, so if units change, these will need to change too. refs #518

[dlebauer]

@ankurdesai Per #518, shouldn't the following return the input value of qair (0.01)?

rh2qair(rh = qair2rh(qair = 0.01, temp = 20, press = 1013.25), T = 20 + 272.15, press = 1013.25)

[dlebauer]

(note that I've started to draft a test for this in modules/data.atmosphere/tests/testthat/test.metutils.R@qair

[ankurdesai]

It does if you do this:
rh2qair(rh = qair2rh(qair = 0.01, temp = 20, press = 1013.25), T = 20 + 273.15, press=101325)

(replace 1013.25 with 101325, replace 272.15 with 273.15)

0.01

rh2qair expects pressure in Pa, temperature in Kelvin. qair2rh expects pressure in hPa, temperature in C. Once the functions calls in various places are fixed, we can change these to be consistent.

---

Ankur R Desai, Associate Professor
University of Wisconsin - Madison, Atmospheric and Oceanic Sciences
http://flux.aos.wisc.edu http://flux.aos.wisc.edu/ desai@aos.wisc.edu mailto:desai@aos.wisc.edu
O: +1-608-520-0305 / M: +1-608-218-4208

On Jun 30, 2015, at 3:40 PM, David LeBauer notifications@github.com wrote:

@ankurdesai https://github.com/ankurdesai Per #518 #518, shouldn't the following return the input value of qair (0.01)?

rh2qair(rh = qair2rh(qair = 0.01, temp = 20, press = 1013.25), T = 20 + 272.15, press = 1013.25)
—
Reply to this email directly or view it on GitHub #528 (comment).