document forecast process

[iantaylor-NOAA]

Inspired by chat with @okenk and @chantelwetzel-noaa. Please edit this list as needed.

Documentation needed for things like

• Forecast settings
  • Adding the Pstar buffer to your forecast file using PEPtools::get_buffer()
    • default sigma input to that function is 0.5 for Category 1 stocks and 1.0 for Cat 2 but if uncertainty in model (either model$Pstar_sigma or model$OFL_sigma?) is larger, use that instead.
    • pstar input is set by the council (check with GMT rep)
  • Make sure you are using 3 # Control rule method (0: none; 1: ramp does catch=f(SSB), buffer on F; 2: ramp does F=f(SSB), buffer on F; 3: ramp does catch=f(SSB), buffer on catch; 4: ramp does F=f(SSB), buffer on catch)
  • Fixed catches for first 2 years of forecast
    • typically provided by GMT rep to STAR
    • OFL/ACL values are in GMT016 table in https://reports.psmfc.org/pacfin/f?p=501:1000:4357688221402::::: (remember to filter for your species)
• Creating the forecast table of OFL/ABC/ACL
  • r4ss::SSexecutivesummar() does most of the work.
    • You may want to specify the OFLs (and ACLs) for the next two years: SSexecutivesummary(model, forecast_ofl = c(3763, 3563))