tests/testthat/test-radiation.R

In jonasViehweger/fieldClim: Calculates Fluxes for Weather Stations

context("Radiation")

# Test data
t <- 20.8
st <- 15 #surface temp
p <- 1000
elev <- 200
datetime <- as.POSIXlt("2018-9-29 11:12:00", tz = "CET")
lat <- 50.84050277777778
lon <- 8.683303333333333
sky_view <- 0.8
em <- 0.95
tges <- 0.6058589
albedo <- 0.15
slope <- 20
sol_azimuth <- 179.1576
sol_elevation <- 36.98267
hex <- 180

test_that("rad_emissivity_air", {
  expect_equal(rad_emissivity_air(t, elev, p), 0.2090643, tolerance = 1e-5)
})

test_that("rad_lw_out", {
  expect_equal(rad_lw_out(st, em), 371.3726, tolerance = 1e-5)
})

test_that("rad_lw_in", {
  expect_equal(rad_lw_in(0.2090643, t), 88.50861, tolerance = 1e-4)
})

test_that("rad_sw_toa", {
  expect_equal(rad_sw_toa(datetime, lat, lon), 819.5255, tolerance = 1e-3)
})

test_that("rad_sw_in", {
  expect_equal(rad_sw_in(819.5255, tges), 496.5168, tolerance = 1e-3)
})

test_that("rad_sw_out", {
  expect_equal(rad_sw_out(496.5168, albedo), 74.47752, tolerance = 1e-3)
})

test_that("rad_sw_radiation_balance", {
  expect_equal(rad_sw_radiation_balance(496.5168, 74.47752), 422.0393, tolerance = 1e-3)
})

test_that("rad_bal_total", {
  expect_equal(rad_bal_total(rad_sw_radiation_balance = 422.0393,
                             rad_lw_out = 371.3726,
                             rad_lw_in = 88.50861), 139.1754, tolerance = 1e-3)
})