data_ss3.R

# data.R - condition OM(s)
# WKMFOA_toolset/data.R

# Copyright (c) WUR, 2023.
# Author: Iago MOSQUEIRA (WMR) <iago.mosqueira@wur.nl>
#
# Distributed under the terms of the EUPL-1.2


library(icesTAF)
mkdir("data")

library(mse)
library(FLSRTMB)

# CHOOSE number of cores for doFuture
cores <- 2

source("utilities.R")

# LOAD SS3 SA results, 2024 ICES WGNSSK sol.27.4
load('boot/data/sol274.rda', verbose=TRUE)

# DATA year
dy <- dims(run)$maxyear

# INTERMEDIATE year
iy <- dy + 1

# HINDCAST years
hy <- iy - dims(run)$age
hys <- ac(seq(hy, dy))

# FINAL year
fy <- 2050

# NUMBER of iterations
it <- 500

# RNG seed
set.seed(987)

# - CONSTRUCT OM

# GENERATE lognormal rec deviances, sigma and rho from SS3
srdevs <- rlnormar1(n=500, sdlog=0.5384, rho=0, years=seq(hy, fy))

plot(srdevs) +
  geom_vline(xintercept=dy, linetype=2)

# BUILD FLom w/ SRR
om <- FLom(stock=propagate(run, it), refpts=refpts, model='bevholtss3',
  params=params(srr), deviances=srdevs)

# HINDCAST for last 10 years, uses estimated recruitment + deviances
om <- fwd(om, catch=catch(om)[, hys], sr=rec(om))

# SETUP om future: average of last 3 years **
om <- fwdWindow(om, end=fy)

# PROJECT forward for iy assumption (TAC) **
om <- fwd(om, catch=FLQuant(3675, dimnames=list(year=2024)))

# CREATE F and SSB deviances
sdevs <- shortcut_devs(om, Fcv=0.212, Fphi=0.423, SSBcv=0.10)

# SAVE
save(om, sdevs, file="data/data.rda", compress="xz")