tests/testthat/test-sensible.R

In jonasViehweger/energyBudget: Calculates Fluxes for Weather Stations

context("Sensible")

hum1 <- 89     # Relative Humidity Hoehe 1
hum2 <- 88     # Relative Humidity Hoehe 2
p0 <- 1013.25  # Referenzdruck Meeresniveau in hPa
p <- 1000      # Gemessener Druck in hPa
t1 <- 22        # Lufttemperatur auf 2 m
t2<- 21        # Lufttemperatur auf 10 m
z1 <- 2        # Erste Messhoehe Meter
z2 <- 10       # 2te MEsshoehe
v1 <- 1        # m/Sec Windgeschwindigkeit Hoehe 1
v2 <- 2.3      # m/Sec Windgeschwindigkeit Hoehe 2
ludi <- 1.225  # Luftdiche Meeresniveau kg/m?
rad_bal <- 400      # STrahlungsbilanz W/m?
soil_flux <- 40       # Bodenw?rmestrom W/m?

test_that("sensible_bowen", {
  expect_equal(sensible_bowen(t1 = t1, t2 = t2,
                              hum1 = hum1, hum2 = hum2,
                              p1 = p, p2 = p,
                              z1 = z1, z2 = z2,
                              rad_bal = rad_bal, soil_flux = soil_flux),
               -128.3694, tolerance = 1e-3)
})

grad_rich_no <- -0.1573365
monin <- -63.55804
ustar <- 0.5867643

test_that("sensible_monin", {
  expect_equal(sensible_monin(t1, t2, p1 = p, p2 = p, z1, z2,
                              monin, ustar, grad_rich_no),
               237.7302, tolerance = 1e-3)
})


# Priestley taylor test data

t <- 22  # Temperatur
rad_bal <- 400
soil_flux <- -40

test_that("sensible_priestley_taylor", {
  expect_equal(sensible_priestley_taylor(t, rad_bal, soil_flux), -35.11673, tolerance = 1e-3)
})