library(stockassessment)
source("src/forcastfix.R")
setwd("run")
load("model.RData")

base.y<-max(fit$data$years)-1  # in spring

# not used
FC<-list()

#ref points
Fmsy     <-0.35  
Fmsylower<-0.29
Fpa      <-0.35

MSYBtrig <-201845
Bpa  <-201845
Blim <-144175


startR<-1983  # sampling from recruitment series, start year

# use recent 3 years of fishing selectivity
selYears<-max(fit$data$years) +(-3:-1)  
bioYears<-max(fit$data$years) +(-3:-1)
recYears<-startR:(base.y-1)

#estimated total F at age
summary_F<-summary(fit)
Fsq    <-tail(summary_F[,7],2)[[1]]

FC<-list()

set.seed(12345)
FC[["Fmsy"]] <- forecast(fit,year.base=base.y,ave.years=bioYears, rec.years=recYears,overwriteSelYears=selYears, deterministic=FALSE, splitLD =FALSE, fval=c(Fsq,Fsq,Fmsy,Fmsy),label="Fmsy")                                   


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