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
Fslide<-0.19   # Fmsy*SSB2026/MSYBtrigger
Fmsy<- 0.23     # Fmsy refeence point


set.seed(12345)
FC[[length(FC)+1]]  <- forecast(fit, fscale=c(NA,NA,NA,NA,NA,NA), fval=c(NA, NA,Fslide, 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(NA, NA,Fslide, Fmsy, Fmsy, Fmsy), catchval=c(6709,5500,NA,NA,NA,NA),  rec.years=Ry, label="Catch in 2025, 5500 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(NA, NA,Fslide, Fmsy, Fmsy, Fmsy), catchval=c(6709,6000,NA,NA,NA,NA),  rec.years=Ry, label="Catch in 2025, 6000 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(NA, NA,Fslide, Fmsy, Fmsy, Fmsy), catchval=c(6709,5707,NA,NA,NA,NA),  rec.years=Ry, label="Catch in 2025, 5707 t in 2026 (AVG 2023-25) therafter Fmsy*SSB2026/MSY Btrigger", processNoise=FALSE,addTSB=TRUE, year.base=yearbase, savesim=TRUE) 

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

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