library(stockassessment)
load("run/model.RData")
Fbar<-tail(fbartable(fit),1)[1]
TAC<-4514
FC<-list()

# SQ all years
set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, fval=c(Fbar, Fbar, Fbar, Fbar), processNoiseF=FALSE,rec.years=c(2022:2024), label="SQ all years (recent recruitemt)")
#F=0
set.seed(12345)
FC[[length(FC)+1]] <- forecast(fit, fval=c(Fbar,Fbar,0.000001, 0.000001), processNoiseF=FALSE,rec.years=c(2022:2024), label="2025F=2024F then zero (recent recruitemt)")

# 2025= TAC the SQ all years
#set.seed(12345)
#FC[[length(FC)+1]] <- forecast(fit, catchval=c(NA,TAC,NA, NA), fval=c(Fbar, NA, 1.33, 1.33), processNoiseF=FALSE, rec.years=c(2022:2024), label="2025F=TAC then SQ F")


# SQ then zero
#set.seed(12345)
#FC[[length(FC)+1]] <- forecast(fit, catchval=c(NA,TAC,NA, NA), fval=c(Fbar,NA,0.000001, 0.000001), rec.years=c(2022:2024), processNoiseF=FALSE, label="2025F=TACF then zero")



#Fmsy
#set.seed(12345)
#Fmsy<-0.54
#FC[[length(FC)+1]] <- forecast(fit, fval=c(Fbar,Fbar,Fmsy, Fmsy), processNoiseF=FALSE, label="2025F=2024F then Fmsy")

#Fmsy
set.seed(12345)
Fmsy<-0.54
FC[[length(FC)+1]] <- forecast(fit, fval=c(Fbar,Fbar,Fmsy, Fmsy), processNoiseF=FALSE, rec.years=c(2022:2024), label="2025F=2024F then Fmsy (recent rec)")


#Fmsyower
#set.seed(12345)
#Fmsylow<-0.27
#FC[[length(FC)+1]] <- forecast(fit, fval=c(Fbar,Fbar,Fmsylow, Fmsylow), processNoiseF=FALSE, label="2024F=2023F then Fmsylower")


# 2025 catch based on TAC
#set.seed(12345)
#Fmsy<-0.54
#TAC<-4514
#FC[[length(FC)+1]] <- forecast(fit, catchval=c(NA,TAC,NA, NA),fval=c(Fbar,NA,Fmsy, Fmsy), processNoiseF=FALSE,rec.years=c(2022:2024), label="2025F=TAC then Fmsy (recent rec)")


# 2025 catch based on TAC, only recetn rec
#set.seed(12345)
#Fmsy<-0.54
#TAC<-4514
#FC[[length(FC)+1]] <- forecast(fit, catchval=c(NA,TAC,NA, NA),fval=c(Fbar,NA,Fmsy, Fmsy), processNoiseF=FALSE, rec.years=c(2022:2024), label="2025F=TAC then Fmsy (recent rec)")



## reduced based on Btrigger #########################################################
 # MSY approach: FMSY × SSB (2026)/MSY Btrigger
 
 
# 2025 catch based on TAC
#set.seed(12345)
#Fmsy<-0.39
#TAC<-4514
#FC[[length(FC)+1]] <- forecast(fit, catchval=c(NA,TAC,NA, NA),fval=c(Fbar,NA,Fmsy, Fmsy), processNoiseF=FALSE, label="2025F=TAC then Fmsy(reduce by Btrigger) ")


# 2025 catch based on TAC, only recetn rec 
#set.seed(12345)
#Fmsy<-0.39
#TAC<-4514
#FC[[length(FC)+1]] <- forecast(fit, catchval=c(NA,TAC,NA, NA),fval=c(Fbar,NA,Fmsy, Fmsy), processNoiseF=FALSE, rec.years=c(2022:2024), label="2025F=TAC then Fmsy (reduced by Btrigger) (recent rec)")



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