Changes upper bounds of logkappa2

using VAST_3.9.0 the upper bound of logkappa2 is reached
using the default settings b/c of the small number of mesh points.
This commit changes the upper bound to 10 using
run_model = FALSE to get the correct map.
This trick can be found in
James-Thorson-NOAA/VAST#300
which was started by @kellijohnson-NOAA b/c of similar issues
when running the model for West Coast species.

sandbox/pcod_delta_gamma.R

@@ -67,7 +67,21 @@ fit_vast <- fit_model(
   b_i = pcod[, "density"],
   a_i = pcod[, "effort"],
   input_grid = input_grid,
-  working_dir = paste0(here("doc", "appendix-VAST"), "/")
+  working_dir = paste0(here("doc", "appendix-VAST"), "/"),
+  run_model = FALSE
+)
+upperbounds <- fit_vast$Upper
+upperbounds["logkappa2"] <- 10
+fit_vast <- fit_model(
+  settings = settings,
+  Lat_i = pcod[, "lat"],
+  Lon_i = pcod[, "lon"],
+  t_i = pcod[, "year"],
+  b_i = pcod[, "density"],
+  a_i = pcod[, "effort"],
+  input_grid = input_grid,
+  working_dir = paste0(here("doc", "appendix-VAST"), "/"),
+  upper = upperbounds
 )
 #   saveRDS(fit_vast, file = f)
 # } else {