tests/testthat/test-humidity.R

In FabianMitze/FieldClimRevised: Calculates Fluxes for Weather Stations

context("Humidity")

# test data
t1 <- 22
t2 <- 21
t_pot <- 23.1132
hum1 <- 89
hum2 <- 88
elev <- 200
z1 <- 2
z2 <- 10
p_vapor <- 23.83956
p1 <- 1000
p2 <- 1000


test_that("hum_sat_vapor_pres", {
  expect_equal(hum_sat_vapor_pres(t = t1), 26.78602, tolerance = 1e-3)
})

test_that("hum_vapor_pres", {
  expect_equal(hum_vapor_pres(hum = hum1, t = t1), 23.83956, tolerance = 1e-3)
})

test_that("hum_specific", {
  expect_equal(hum_specific(p_vapor = p_vapor, p = 1000), 0.0148282, tolerance = 1e-4)
})

test_that("hum_absolute", {
  expect_equal(hum_absolute(p_vapor = p_vapor, t_pot = t_pot), 0.01743569, tolerance = 1e-4)
})

test_that("hum_evap_heat", {
  expect_equal(hum_evap_heat(t = t1), 2448643, tolerance = 1e-4)
})

test_that("hum_precipitable_water", {
  expect_equal(hum_precipitable_water(p = p1, t = t1, elev = elev),
               3.322423, tolerance = 1e-4)
})

test_that("hum_moisture_gradient", {
  expect_equal(hum_moisture_gradient(hum1, hum2, t1, t2, p1, p2),
               -0.00013045, tolerance = 1e-6)
})