test_that("combining cfp_dat works", {
  mod1 <- filter(base_dat, site == "site_a")
  mod2 <- filter(base_dat, site == "site_b")

  y <- combine_models(list(mod1, mod2))

  expect_equal(n_profiles(base_dat), n_profiles(y))
  expect_equal(cfp_id_cols(y),
               c(cfp_id_cols(base_dat), "cmb_id"))


})