library(stockassessment)
load("run/model.RData")

#RETRO<-retro(fit, year=5)

retmat<-cbind(2022-0:4, 3000, 2022-0:4)
RETRO<-retro(fit, year=retmat)


save(RETRO, file="run/retro.RData")