diff --git a/tests/testthat/test-evap.R b/tests/testthat/test-evap.R
new file mode 100644
index 0000000000000000000000000000000000000000..0bf33e2b16c125d8fb18cece1c1050aaf98d3c2a
--- /dev/null
+++ b/tests/testthat/test-evap.R
@@ -0,0 +1,15 @@
+context("Test evaporation")
+
+test_that("PEdaily_Oudin works", {
+  skip_on_cran()
+  rm(list = ls())
+  
+  data(L0123001)
+  PotEvap <- PEdaily_Oudin(JD = as.POSIXlt(BasinObs$DatesR)$yday + 1,
+                           Temp = BasinObs$T,
+                           Lat = 0.8, LatUnit = "rad")
+  PotEvapFor <- PEdaily_Oudin(JD = as.POSIXlt(BasinObs$DatesR)$yday + 1,
+                              Temp = BasinObs$T,
+                              Lat = 0.8, LatUnit = "rad", run_fortran = T)
+  expect_true(all(range(PotEvap - PotEvapFor) < 0.001))
+})