test_that("whether ind_gw_recharge works", {
  expect_equal(
    ind_gw_recharge(
      B_LU_BRP = c(265,265,265,233,233,233),
      D_PSP = c(100,200,300,400,500,600),
      D_WRI_K = c(50,90,50,100,50,90),
      B_DRAIN = c(TRUE,FALSE,TRUE,FALSE,TRUE,FALSE),
      B_GWL_CLASS = c('III','III','III','IV','III','IV'),
      B_SC_WENR = rep(c('Zeer groot', 'Zeer beperkt'), 3),
      D_SE = c(1, 1, 1, 6.049445, 6.049445, 6.049445)
    ),
    expected = c(0.1514698, 0.3802220, 0.5656138, 0.9614681, 0.6772030, 0.9557414),
    tolerance = 0.01
  )
 
})

test_that("B_DRAIN affects scores output for IIIb and IV", {
  expect_all_false(
    ind_gw_recharge(
      B_LU_BRP = c(265,265,265,233,233,233),
      D_PSP = c(100,200,300,400,500,600),
      D_WRI_K = c(50,90,50,100,50,90),
      B_DRAIN = rep(TRUE, 6),
      B_GWL_CLASS = c('IIIb','IIIb','IIIb','IV','IIIb','IV'),
      B_SC_WENR = rep(c('Zeer groot', 'Zeer beperkt'), 3),
      D_SE = c(1, 1, 1, 6.049445, 6.049445, 6.049445)
    ) ==
      ind_gw_recharge(
        B_LU_BRP = c(265,265,265,233,233,233),
        D_PSP = c(100,200,300,400,500,600),
        D_WRI_K = c(50,90,50,100,50,90),
        B_DRAIN = rep(FALSE, 6),
        B_GWL_CLASS = c('IIIb','IIIb','IIIb','IV','IIIb','IV'),
        B_SC_WENR = rep(c('Zeer groot', 'Zeer beperkt'), 3),
        D_SE = c(1, 1, 1, 6.049445, 6.049445, 6.049445)
        )
      )
})

test_that("Groundwater recharge is correctly modified with a target by ind_gw_target", {
  dt <- data.table(
    D_RISK_GWR = c(0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7),
    B_SOILTYPE_AGR = c('veen',
                       'zeeklei', 'zeeklei', 'zeeklei', 'zeeklei',
                       'loess', 'loess'),
    B_GWL_CLASS = c('IIa',
                    'IIa', 'Va', 'Vb', 'IVu',
                    'IIa', 'VIIIo'),
    B_AREA_DROUGHT = c(FALSE,
                       FALSE, FALSE, FALSE, FALSE,
                       TRUE, TRUE)
  )
  
  dt[,S_H_GWR := ind_gw_target(D_RISK_GWR, B_SOILTYPE_AGR, B_GWL_CLASS, B_AREA_DROUGHT = B_AREA_DROUGHT)]
  
  expect_gt(dt[B_SOILTYPE_AGR == 'veen' & B_GWL_CLASS == 'IIa', S_H_GWR], 
            dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'IIa', S_H_GWR])
  expect_gt(dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'IIa', S_H_GWR], 
            dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'Va', S_H_GWR])
  expect_gt(dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'Va', S_H_GWR], 
            dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'Vb', S_H_GWR])
  expect_gt(dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'Vb', S_H_GWR], 
            dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'IVu', S_H_GWR])
  expect_gt(dt[B_SOILTYPE_AGR == 'zeeklei' & B_GWL_CLASS == 'IIa' & B_AREA_DROUGHT == FALSE, S_H_GWR],
            dt[B_SOILTYPE_AGR == 'loess' & B_GWL_CLASS == 'IIa' & B_AREA_DROUGHT == TRUE, S_H_GWR])
  expect_gt(dt[B_SOILTYPE_AGR == 'loess' & B_GWL_CLASS == 'IIa', S_H_GWR], 
            dt[B_SOILTYPE_AGR == 'loess' & B_GWL_CLASS == 'VIIIo', S_H_GWR])

})

test_that('A warning is thrown when deprecated argument I_P_SE is used in ind_gw_recharge',{
  expect_warning(
    ind_gw_recharge(
      B_LU_BRP = 265,
      D_PSP = 400,
      D_WRI_K = 50,
      I_P_SE = 0.8,
      B_DRAIN = TRUE, 
      B_GWL_CLASS = 'IIIb',
      B_SC_WENR = 'Matig'
    )
  )
})

test_that('A warning is thrown when deprecated argument I_P_CO is used in ind_gw_recharge',{
  expect_warning(
    ind_gw_recharge(
      B_LU_BRP = 265,
      D_PSP = 400,
      D_WRI_K = 50,
      I_P_CO = 0.1,
      B_DRAIN = TRUE, 
      B_GWL_CLASS = 'IIIb',
      D_SE = 25
    )
  )
})

test_that('A warning is thrown when both D_SE and I_P_SE are supplied to ind_gw_recharge',{
  expect_warning(
    ind_gw_recharge(
      B_LU_BRP = 265,
      D_PSP = 400,
      D_WRI_K = 50,
      B_DRAIN = TRUE,
      I_P_SE = 0.8,
      B_GWL_CLASS = 'IIIb',
      B_SC_WENR = 'Matig',
      D_SE = 25
    )
  )
})

test_that('A warning is thrown when both I_P_CO and B_SC_WENR or D_P_CO are suppolied to ind_gw_recharge',{
  expect_warning(
    ind_gw_recharge(
      B_LU_BRP = 265,
      D_PSP = 400,
      D_WRI_K = 50,
      B_DRAIN = TRUE,
      I_P_CO = 0.1,
      B_GWL_CLASS = 'IIIb',
      B_SC_WENR = 'Matig',
      D_SE = 25
    )
  )
  
  expect_warning(
    ind_gw_recharge(
      B_LU_BRP = 265,
      D_PSP = 400,
      D_WRI_K = 50,
      B_DRAIN = TRUE,
      I_P_CO = 0.1,
      B_GWL_CLASS = 'IIIb',
      D_SE = 25,
      D_P_CO = 0.5
    )
  )
})


test_that('D_P_CO is favoured over B_SC_WENR',{
  expect_gt(
    object = 
      ind_gw_recharge(
        B_LU_BRP = 2014,
        D_PSP = 400,
        D_WRI_K = 50,
        B_DRAIN = TRUE,
        B_GWL_CLASS = 'VI',
        B_SC_WENR = 'Groot',
        D_SE = 25,
        D_P_CO = 1 # this should be used instead of B_SC_WENR
      ),
    expected = 
      ind_gw_recharge(
        B_LU_BRP = 2014,
        D_PSP = 400,
        D_WRI_K = 50,
        B_DRAIN = TRUE,
        B_GWL_CLASS = 'VI',
        B_SC_WENR = 'Groot',
        D_SE = 25
      )
  )
})

test_that('D_SE is used instead of I_P_SE when both are provided and a warning is given',{
  expect_warning(
    ind_gw_recharge(
      B_LU_BRP = 2014,
      D_PSP = 400,
      D_WRI_K = 50,
      B_DRAIN = TRUE,
      B_GWL_CLASS = 'VI',
      B_SC_WENR = 'Groot',
      D_SE = 5,
      I_P_SE = 0.95
    )
  )
  
  suppressWarnings(
  expect_gt(
    object = 
      ind_gw_recharge(
        B_LU_BRP = 2014,
        D_PSP = 400,
        D_WRI_K = 50,
        B_DRAIN = TRUE,
        B_GWL_CLASS = 'VI',
        B_SC_WENR = 'Groot',
        D_SE = 25,
        I_P_SE = 0.15
      ),
    expected = 
      ind_gw_recharge(
        B_LU_BRP = 2014,
        D_PSP = 400,
        D_WRI_K = 50,
        B_DRAIN = TRUE,
        B_GWL_CLASS = 'VI',
        B_SC_WENR = 'Groot',
        D_SE = 2,
        I_P_SE = 0.95
      )
  )
  )
})