library(stockassessment)
setwd("run")
load("data.RData")
conf<-loadConf(dat,"../conf/model.cfg", patch=TRUE)
par<-defpar(dat,conf)

# Fixing SRR parameters:
par$rec_pars <- c(17.78250, 12.65671) # estimates from benchmark
map <- list()
map$rec_pars <- factor(c(NA,NA))

fit<-sam.fit(dat,conf,par, map=map)
if(fit$opt$convergence!=0) stop("Model did not converge.")
save(fit, file="model.RData")