#Set species name
species <- "Pacific_sardine"
#Directory with SS files
model_dir <- "Sardine_SS_files/"
#Name of SS data file
datfile_name <- "sardEM_3_3.dat"
#Efficiency
eff_case <- 0.5
#Age classes of your species
age_classes <- 1:12
#Selectivity function - a vector equivalent to length (age classes)
sel_by_age <- rep(1,length(age_classes))
timestep <- stepperyr #5
#Which atlantis timestep does the survey run in?--now from census_spec.R
# with 5 output steps per year, 0 is Jan-Feb-midMar, 1 is midMar-Apr-May,
# 2 is June-July-midAug, 3 is midAug-Sept-Oct, 4 is Nov-Dec (ish)
survey_sample_time <- midptyr #3 this is more late summer-fall, changed to 2
#The last timestep to sample
total_sample <- noutsteps-1 #495
#Vector of indices of survey times to pull
survey_sample_full <- seq(survey_sample_time,
total_sample, by=timestep)
#Effective sample size for composition data (annual total samples)
surveyEffN <- 1000
fisheryEffN <- 1000
#CVs for length at age, catch, and survey
CVs <- list("lenage"=0.1, "fishery"=0.01, "survey"=0.1)
#Number of years of data to pull
nyears <- 50
#Atlantis initialization period in years
burnin <- 30
#Maximum size of a fish in cm
maxbin <- 150
#The last timestep to sample
total_sample <- noutsteps-1 #495
#Vector of indices of catch in numbers to pull (by timestep to sum)
fish_sample_full <- c(0:total_sample)
fish_burnin <- burnin*fstepperyr+1
fish_nyears <- nyears*fstepperyr
fish_times <- fish_sample_full[fish_burnin:(fish_burnin+fish_nyears-1)]
fish_timesteps <- seq(fish_times[5], max(fish_times), by=fstepperyr) #last timestep
#fish_years <- unique(floor(fish_times/fstepperyr)+1) # my original
fish_years <- unique(floor(fish_times/fstepperyr))
#Years to survey, assuming survey is once per year
survey_years <- survey_sample_full[burnin:(burnin+nyears-1)]
#Month of survey/fishing
survey_month <- 7
fishing_month <- 1
#Efficiency of the survey
effic <- data.frame(species=species, efficiency=eff_case)
# Assume uniform selectivity, same as selex1 here
sel<- data.frame(species=rep(species,length(age_classes)),
agecl=age_classes,
selex=sel_by_age)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.