tests/testthat/test-turbulence.R

In jonasViehweger/energyBudget: Calculates Fluxes for Weather Stations

context("Turbulence")

# Testdata
ah <- 10             # Measurement height
h <- 1.2   # Obstacle height in m
vh <- 1.2  # Vegetation height in m
v <- 3.5   # wind velocity

test_that("turb_roughness_length", {
  expect_equal(turb_roughness_length(obs_height = h), 0.12)
  expect_equal(turb_roughness_length(surface_type = "Meadow"), 0.02)
})

test_that("turb_displacement", {
  expect_equal(turb_displacement(vh), 0.8)
})

test_that("turb_ustar", {
  expect_equal(turb_ustar(v, ah, 0.12), 0.3165381, tolerance = 1e-5)
})