#library(stockassessment)

#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")