Ensure that allo_find() matches equations by site (not only species) about fgeo.biomass HOT 4 CLOSED

library(tidyverse)
library(fgeo.biomass)

census <- fgeo.biomass::scbi_tree1
species <- fgeo.biomass::scbi_species



census_species <- census %>%
  add_species(species, "scbi")

result <- fgeo.biomass::scbi_tree1 %>%
  add_species(species, site = "scbi") %>%
  allo_find()

# The result has equations of multiple types
result
#> # A tibble: 5 x 2
#>   eqn_type       data                 
#>   <chr>          <list>               
#> 1 species        <tibble [8,930 x 8]> 
#> 2 genus          <tibble [5,642 x 8]> 
#> 3 mixed_hardwood <tibble [5,516 x 8]> 
#> 4 family         <tibble [10,141 x 8]>
#> 5 woody_species  <tibble [0 x 8]>

# But they all come from SCBI
result %>%
  unnest() %>%
  pull(site) %>%
  unique()
#> [1] "scbi"