Commit 6047af60 authored by Berges, Benoit's avatar Berges, Benoit
Browse files

update scripts

parent d8250cf9
......@@ -54,7 +54,7 @@ fullPeriod <- c(histPeriod,projPeriod)
recrPeriod <- ac(2011:2021)
selPeriod <- ac(2010:2020)
fecYears <- ac(2010:2020)
nits <- 500 # number of random samples
nits <- 20 # number of random samples
load(file.path(dataPath,paste0("data_sms",SMSkeyRuns,"_sf.RData")))
......@@ -493,7 +493,7 @@ for(idxFleet in 1:length(fleets)){
# 10) Future recruitment
#-------------------------------------------------------------------------------
recPeriod <- ac(2002:2019)
recPeriod <- ac(2002:2021)
biol.sr <- fmle(as.FLSR(biol,model='segreg')) # just to populate the structure
......
......@@ -53,7 +53,7 @@ source(file.path(functionPath,"find.FCAR.r"))
# - F sel: FAsel, FCsel, FBDsel
#-------------------------------------------------------------------------------
nits <- 500
nits <- 20
# load object
load(file.path(modelPath,paste0(runName,'_MSE_init_',ac(nits),'.RData')))
......@@ -103,7 +103,8 @@ SSBHCR <- FLQuant(NA, dimnames=list(age="all",year=ac(an(proj
#------------------------------------------------------------------------------#
start.time <- Sys.time()
for (iYr in an(projPeriod)){
#an(projPeriod)
for (iYr in 2033:2042){
cat(iYr,"\n")
cat(paste("\n Time running",round(difftime(Sys.time(),start.time,unit="mins"),0),"minutes \n"))
......@@ -227,6 +228,14 @@ for (iYr in an(projPeriod)){
# Forecast
#-------------------------------------------------------------------------------
iStocks <- stkAssessment
iFishery <- fishery
cYr <- iYr
iTAC <- TAC
iTAC_var <- TAC_var
iHistMaxYr <- histMaxYr
mpPoints <- referencePoints
projNSAS <- projectNSH_FMSYAR(stkAssessment,
fishery,
iYr,
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment