Error generating extrapolation grid with combined IBTS Regions

[AndyCampbell]

Combining IBTS survey regions raises the following error

With region as c("EVHOE","IE-IGFS","ATL-IBTS-Q1") in the call to make_settings (or any subset of these) the extrapolation grid is successfully created during fit_model().
However, if I attempt to add either or both of "SP-PORC" and "SP-NORTH" to the vector of regions the following error results

Reading shapefile with projargs: +proj=longlat +ellps=WGS84 +no_defs
Error in combine_extrapolation_info(Return, Extrapolation_List, create_strata_per_region = create_strata_per_region) :
Must use same projargs for projection of all extrapolation grids

SessionInfo

R version 4.3.2 (2023-10-31 ucrt)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 17763)

Matrix products: default

locale:
[1] LC_COLLATE=English_Ireland.1252 LC_CTYPE=English_Ireland.1252 LC_MONETARY=English_Ireland.1252
[4] LC_NUMERIC=C LC_TIME=English_Ireland.1252

time zone: Europe/Dublin
tzcode source: internal

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] Cairo_1.6-1 lubridate_1.9.2 forcats_1.0.0 stringr_1.5.1
[5] dplyr_1.1.4 purrr_1.0.2 readr_2.1.4 tidyr_1.3.0
[9] tibble_3.2.1 ggplot2_3.4.4 tidyverse_2.0.0 VAST_3.10.1
[13] FishStatsUtils_2.12.1 marginaleffects_0.16.0 units_0.8-5 TMB_1.9.9

settings = FishStatsUtils::make_settings(n_x = knots,
Region = c("EVHOE","IE-IGFS","ATL-IBTS-Q1"),
FieldConfig = c("Omega1" = "IID", "Epsilon1" = "IID", "Omega2" = "IID", "Epsilon2" = "IID"),
RhoConfig = c("Beta1"=0, "Beta2"=0, "Epsilon1"=0, "Epsilon2"=0),
OverdispersionConfig = c("eta1"=0,"eta2"=0),
ObsModel = c(PosDist = 1, Link = 0),
purpose = "index2",
bias.correct = FALSE,
max_cells = Max_cells)