library(stockassessment)

setwd("run")
load("model.RData")

# not used
FC<-list()


#ref points
Fmsy<-0.35  #ne fmsy=fmsylower
Fpa<-0.51
Fsq<-0.085

Bpa<-201845
Blim<-144175

rec.years<-c(1983:2023)
base<-max(fit$data$years)-1


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



set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, ave.years=bioYears,year.base = base, rec.years=rec.years,overwriteSelYears=selYears, fval=c(Fsq,Fsq,Fmsy,Fmsy),splitLD =FALSE, label="SQ then Fmsy")

set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, ave.years=bioYears,year.base = base, rec.years=rec.years,overwriteSelYears=selYears,fval=c(Fsq,Fsq,Fsq,Fsq),splitLD =FALSE, label="FSQ, all years")


set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, ave.years=bioYears,year.base = base, rec.years=rec.years,overwriteSelYears=selYears, fval=c(Fsq,Fsq,0.000001,0.000001),splitLD =FALSE, label="SQ then zero")


set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, ave.years=bioYears, year.base = base,rec.years=rec.years,overwriteSelYears=selYears, fval=c(Fsq,Fsq,Fpa,Fpa),splitLD =FALSE, label="SQ then Fpa")



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