library(stockassessment)
load("run/model.RData")
FC<-list()

#2020 forecast
#Fmsy
set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, fscale=c(NA,NA,NA,NA), fval=c(0.527,0.527,0.59,0.59), label="2018F=2019F then Fmsy", rec.years=1973:2018)

# 2020  TAC2019
# Fmsy
set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, fval=c(.527,NA,.59,.59), catchval=c(NA,20000,NA,NA),label="2018+TAC2019 then Fmsy", rec.years=1973:2018)

save(FC, file="run/forecast.RData")