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

Ry<-2020:2024     # latest 5 years as recommended from ADGdeep 2022 (wgdeep2023) 
Fsq<-0.798       # latest value from SAM 
Blim<-9340		# Blim reference point
Bpa<-11627		# Bpa, MSYBtrigger reference point
yearbase<- 2025   # assessment year-1
Fmsy<-0.16   # Fmsy reference point


set.seed(12345)
FC[[length(FC)+1]]  <- forecast(fit, fscale=c(NA,NA,NA,NA,NA,NA), fval=c(NA, NA,Fmsy, Fmsy, Fmsy, Fmsy), catchval=c(6709,5000,NA,NA,NA,NA),  rec.years=Ry, label="Catch in 2025,5000 t in 2026 therafter Fmsy*SSB2026/MSY Btrigger", processNoise=FALSE,addTSB=TRUE, year.base=yearbase, savesim=TRUE) 
set.seed(12345)
FC[[length(FC)+1]]  <- forecast(fit, fscale=c(NA,NA,NA,NA,NA,NA), fval=c(Fsq, NA,Fmsy, Fmsy, Fmsy, Fmsy), catchval=c(NA,5000,NA,NA,NA,NA),  rec.years=Ry, label="Fsq in 2025, 5000 t in 2026 therafter Fmsy*SSB2026/MSY Btrigger", processNoise=FALSE,addTSB=TRUE, year.base=yearbase, savesim=TRUE) 

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