tests/testthat/test_model_runs.R

In stineb/rpmodel: P-Model

context("test model runs")

test_that("default model run",{
  skip_on_cran()
  
  # no atmospheric pressure given
  expect_warning(rpmodel( 
    tc             = 20,
    vpd            = 1000,
    co2            = 400,
    fapar          = 1,
    ppfd           = 300,
    elv            = 0,
    kphio          = 0.049977,
    beta           = 146,
    c4             = FALSE,
    method_jmaxlim = "none",
    do_ftemp_kphio = FALSE,
    do_soilmstress = FALSE,
    verbose        = TRUE
  ))
  
  out_pmodel <- rpmodel( 
    tc             = 20,
    vpd            = 1000,
    co2            = 400,
    fapar          = 1,
    ppfd           = 300,
    elv            = 0,
    kphio          = 0.049977,
    beta           = 146,
    patm           = 1024,
    c4             = FALSE,
    method_jmaxlim = "none",
    do_ftemp_kphio = FALSE,
    do_soilmstress = FALSE,
    verbose        = TRUE
  )
  
  # output must be a list
  expect_type(out_pmodel, "list")
  
  # kphio temp
  out_pmodel_kphio <- rpmodel( 
    tc             = 20,
    vpd            = 1000,
    co2            = 400,
    fapar          = 1,
    ppfd           = 300,
    elv            = 0,
    kphio          = 0.049977,
    beta           = 146,
    patm           = 1024,
    c4             = FALSE,
    method_jmaxlim = "none",
    do_ftemp_kphio = TRUE,
    do_soilmstress = FALSE,
    verbose        = TRUE
  )
  
  # output must be a list
  expect_type(out_pmodel_kphio, "list")
  
  # for c4 run, assimilation and GPP are not identical
  # results in error (tests routine, not sure if
  # settings are meaningful)
  expect_warning(
    rpmodel( 
      tc             = 20,
      vpd            = 1000,
      co2            = 400,
      fapar          = 1,
      ppfd           = 300,
      elv            = 0,
      kphio          = 0.049977,
      beta           = 146,
      patm           = 1024,
      c4             = TRUE,
      method_jmaxlim = "none",
      do_ftemp_kphio = FALSE,
      do_soilmstress = FALSE,
      verbose        = TRUE
    )
  )
  
  # enables soil moisture stress routine
  out_rpmodel_soilm <- rpmodel( 
      tc             = 20,
      vpd            = 1000,
      co2            = 400,
      fapar          = 1,
      ppfd           = 300,
      elv            = 0,
      kphio          = 0.049977,
      beta           = 146,
      patm           = 1024,
      c4             = FALSE,
      soilm          = 1,
      method_jmaxlim = "wang17",
      do_ftemp_kphio = FALSE,
      do_soilmstress = TRUE,
      verbose        = TRUE
    )
  
  # check model output
  expect_type(out_rpmodel_soilm, "list")
})

test_that("jmax wang17",{
  skip_on_cran()
  
  
  out_pmodel <- rpmodel( 
    tc             = 20,
    vpd            = 1000,
    co2            = 400,
    fapar          = 1,
    ppfd           = 300,
    elv            = 0,
    kphio          = 0.049977,
    beta           = 146,
    patm           = 1024,
    c4             = FALSE,
    method_jmaxlim = "wang17",
    do_ftemp_kphio = FALSE,
    do_soilmstress = FALSE,
    verbose        = TRUE
  )
  
  # output must be a list (no atmosphere warning)
  expect_type(out_pmodel, "list")
})

test_that("jmax smith19",{
  skip_on_cran()
  
  
  out_pmodel <- rpmodel( 
    tc             = 20,
    vpd            = 1000,
    co2            = 400,
    fapar          = 1,
    ppfd           = 300,
    elv            = 0,
    kphio          = 0.049977,
    beta           = 146,
    patm           = 1024,
    c4             = FALSE,
    method_jmaxlim = "smith19",
    do_ftemp_kphio = FALSE,
    do_soilmstress = FALSE,
    verbose        = TRUE
  )
  
  # output must be a list (no atmosphere warning)
  expect_type(out_pmodel, "list")
})