library(stockassessment)
load("run/model.RData")
retmat<-cbind(2022-0:15, 3000, 2022-0:15)
RETRO<-retro(fit, year=retmat)
save(RETRO, file="run/retro.RData")