Exporting parameters and settings
Below are given parameters and settings used in the Bifrost model. Tables of parameters and settings are exported to the file parametersAndSettings.xls on the directory for the selected Bifrost project. For several of the parameters and settings there are in-built explanations which are given on the produced notebook, but not on the file. The explanation is in form of mouse-over and not printed on the table since some of the explanations can be a little bulky.
Parameters are here defined as quantities that determine the value of formulas in the model. All parameters must be defined in in the initialization chapter, although definition characteristics may be modified later.
The parameters given here are the parameters that are or may be used for estimation. In the model software the parameters may be vectors over sex, where the individual components are estimable parameters. If this is the case, it is noted in the explanation column.
The names of parameters and variables used here are those that actually are used in the program. The names in the menu can differ slightly for technical reasons.
Parameters
Global parameters for historic assessment
Name | Value | Status | Explanation |
capelinP2Females | 13.8378 | Estimate | capelinP2 is the length at 50% maturity. capelinP2 is a sex-dependent vector used in the model. For estimation capelinP2Females and capelinP2Males are used. |
capelinP2Males | 13.8378 | capelinP2Females | |
capelinP2Age2Addition | 0.001 | Estimate | |
capelinMaturationM | 0.0385417 | Estimate | |
uncertaintyCapelin | 0.380689 | Estimate | |
useCapelinScaling | 1.20149 | Estimate | |
consExponentQ1 | 0.001 | Estimate | |
otherFoodCodQ1 | 0.0463545 | Estimate | |
partImmaturePreyedByCodWinter | 0.291708 | Estimate | partImmaturePreyedByCodWinter is the proportion of immature capelin biomass during quarter 1 that is subject to predation by cod together with the mature capelin biomass. The total capelin prey during quarter 1 is calculated in the routine exertPredationMortality and the effect on predation per cod (used in likelihood) of mature capelin and immature capelin is calculated in the routine predationMortalityBifrost and the effect on simulations in the routine calculatePredation, however the updating of immature and mature capelin is executed in overlying routines. The effect of partImmaturePreyedByCodWinter increases from zero in the year runLastYearNoConsumptionImmatureCapelinQ1 to its full value in the year runFirstYearFullConsumptionImmatureCapelinQ1. This is calculated in the routine partImmaturePreyed and delivered as the variable usePartImmaturePreyedByCod, which is used in the calculations. |
partMatureCodPreyingOnCapelinWinter | 0.0570797 | Estimate | partMatureCodPreyingOnCapelinWinter is the proportion of mature cod that preys on capelin (immature or mature) during quarter 1. The same linear increase in effect is applied for partMatureCodPreyingOnCapelinWinter as for partImmaturePreyedByCodWinter. |
halfConsBiomassQ1 | 0.001 | Estimate | |
maxConsCodQ1 | 0.321182 | Estimate | |
uncertaintyConsumptionQ1 | 0.956352 | Estimate | |
speciesSuitabilityCodCodQ1 | 0.236946 | Estimate | |
speciesSuitabilityCapelinCodQ1 | 6.39382 | Estimate | |
runScalingConsumptionPerCod | 0.75 | Estimate | |
runScalingCod | 0.75 | Estimate | |
speciesSuitabilityCapelinCodQ4 | 0.929167 | Constant | |
speciesSuitabilityCapelinCodQ2 | 0.929167 | speciesSuitabilityCapelinCodQ4 | |
speciesSuitabilityCapelinCodQ3 | 0.929167 | speciesSuitabilityCapelinCodQ4 | |
suitAge2CapelinQ1 | 0.1 | Estimate | |
suitAge9CapelinQ1 | 0.734864 | Estimate | |
suitAge2CapelinQ4 | 0.25927 | Estimate | |
suitAge2CapelinQ2 | 0.25927 | suitAge2CapelinQ4 | |
suitAge2CapelinQ3 | 0.25927 | suitAge2CapelinQ4 | |
suitAge9CapelinQ4 | 1. | Estimate | |
suitAge9CapelinQ2 | 1. | suitAge9CapelinQ4 | |
suitAge9CapelinQ3 | 1. | suitAge9CapelinQ4 | |
betaPar1CapelinQ1 | 1.53542 | Estimate | |
betaPar2CapelinQ1 | 5.60847 | Estimate | |
betaPar1CapelinQ4 | 13.1672 | Estimate | |
betaPar1CapelinQ2 | 13.1672 | betaPar1CapelinQ4 | |
betaPar1CapelinQ3 | 13.1672 | betaPar1CapelinQ4 | |
betaPar2CapelinQ4 | 1.5 | Estimate | |
betaPar2CapelinQ2 | 1.5 | betaPar2CapelinQ4 | |
betaPar2CapelinQ3 | 1.5 | betaPar2CapelinQ4 | |
suitAge2CodQ1 | 0.675 | Estimate | |
suitAge9CodQ1 | 0.384757 | Estimate | |
suitAge2CodQ4 | 0.250545 | Estimate | |
suitAge2CodQ2 | 0.250545 | suitAge2CodQ4 | |
suitAge2CodQ3 | 0.250545 | suitAge2CodQ4 | |
suitAge9CodQ4 | 1. | Estimate | |
suitAge9CodQ2 | 1. | suitAge9CodQ4 | |
suitAge9CodQ3 | 1. | suitAge9CodQ4 | |
betaPar1CodQ1 | 1.5 | Estimate | |
betaPar2CodQ1 | 3.24376 | Estimate | |
betaPar1CodQ4 | 12.5838 | Estimate | |
betaPar1CodQ2 | 12.5838 | betaPar1CodQ4 | |
betaPar1CodQ3 | 12.5838 | betaPar1CodQ4 | |
betaPar2CodQ4 | 1.5 | Estimate | |
betaPar2CodQ2 | 1.5 | betaPar2CodQ4 | |
betaPar2CodQ3 | 1.5 | betaPar2CodQ4 | |
suitAge1OtherFoodQ1 | 0.1 | Estimate | |
suitAge9OtherFoodQ1 | 0.1 | Estimate | |
suitAge1OtherFoodQ4 | 0.529166 | Estimate | |
suitAge1OtherFoodQ2 | 0.529166 | suitAge1OtherFoodQ4 | |
suitAge1OtherFoodQ3 | 0.529166 | suitAge1OtherFoodQ4 | |
suitAge9OtherFoodQ4 | 0.216668 | Estimate | |
suitAge9OtherFoodQ2 | 0.216668 | suitAge9OtherFoodQ4 | |
suitAge9OtherFoodQ3 | 0.216668 | suitAge9OtherFoodQ4 | |
betaPar1OtherFoodQ1 | 2.20676 | Estimate | |
betaPar2OtherFoodQ1 | 2.57326 | Estimate | |
betaPar1OtherFoodQ4 | 1.5 | Estimate | |
betaPar1OtherFoodQ2 | 1.5 | betaPar1OtherFoodQ4 | |
betaPar1OtherFoodQ3 | 1.5 | betaPar1OtherFoodQ4 | |
betaPar2OtherFoodQ4 | 13.7969 | Estimate | |
betaPar2OtherFoodQ2 | 13.7969 | betaPar2OtherFoodQ4 | |
betaPar2OtherFoodQ3 | 13.7969 | betaPar2OtherFoodQ4 | |
consExponentQ4 | 0.001 | Estimate | |
consExponentQ2 | 0.001 | consExponentQ4 | |
consExponentQ2 | 0.001 | consExponentQ4 | |
consExponentQ3 | 0.001 | consExponentQ4 | |
consExponentQ3 | 0.001 | consExponentQ4 | |
otherFoodCodQ4 | 2.81719 | Estimate | |
otherFoodCodQ2 | 2.81719 | otherFoodCodQ4 | |
otherFoodCodQ3 | 2.81719 | otherFoodCodQ4 | |
halfConsBiomassQ4 | 0.001 | Estimate | |
halfConsBiomassQ2 | 0.001 | halfConsBiomassQ4 | |
halfConsBiomassQ3 | 0.001 | halfConsBiomassQ4 | |
maxConsCodQ4 | 0.313849 | Estimate | |
maxConsCodQ2 | 0.313849 | maxConsCodQ4 | |
maxConsCodQ3 | 0.313849 | maxConsCodQ4 | |
uncertaintyConsumptionQ4 | 0.95594 | Estimate | |
speciesSuitabilityCodCodQ4 | 0.644792 | Estimate | |
speciesSuitabilityCodCodQ2 | 0.644792 | speciesSuitabilityCodCodQ4 | |
speciesSuitabilityCodCodQ3 | 0.644792 | speciesSuitabilityCodCodQ4 | |
capelinM | 0.0016667 | Constant | capelinM is monthly residual natural mortality for immature capelin. |
capelinP1Females | 3.5 | Constant | capelinP1 is the maturation as function of length at capelinP2 - the length at 50% maturity. capelinP1 is a sex-dependent vector used in the model. For estimation capelinP1Females and capelinP1Males are used. |
capelinP1Males | 3.5 | capelinP1Females | |
capelinExtensionExponentQ1 | 1. | Constant | |
capelinExtensionExponentQ2 | 1. | capelinExtensionExponentQ4 | |
capelinExtensionExponentQ3 | 1. | capelinExtensionExponentQ4 | |
capelinExtensionExponentQ4 | 1. | Constant | |
halfCapelinExtensionQ1 | 0.0001 | Constant | |
halfCapelinExtensionQ2 | 0.0001 | halfCapelinExtensionQ4 | |
halfCapelinExtensionQ3 | 0.0001 | halfCapelinExtensionQ4 | |
halfCapelinExtensionQ4 | 0.0001 | Constant | |
codExtensionExponentQ1 | 1. | Constant | |
codExtensionExponentQ2 | 1. | codExtensionExponentQ4 | |
codExtensionExponentQ3 | 1. | codExtensionExponentQ4 | |
codExtensionExponentQ4 | 1. | Constant | |
consExponent | 1.34 | Constant | |
halfCodExtensionQ1 | 0.0001 | Constant | |
halfCodExtensionQ2 | 0.0001 | halfCodExtensionQ4 | |
halfCodExtensionQ3 | 0.0001 | halfCodExtensionQ4 | |
halfCodExtensionQ4 | 0.0001 | Constant | |
codWeightExponent | 0. | Constant | |
otherFoodBottomStart | 0.5 | Constant | |
sizeCodArea | 0.8 | Constant | |
sizeCapelinArea | 0.8 | Constant | |
zerogroupCapelinCatchability | 13.56 | Constant | |
zerogroupCapelinUncertainty | 1.07 | Constant | |
zerogroupCodCatchability | 1416. | Constant | |
onegroupCodCatchability | 145. | Constant | |
uncertaintyOnegroupCod | 0.5 | Constant | |
twogroupCodCatchability | 1023. | Constant | |
uncertaintyTwogroupCod | 0.5 | Constant | |
empiricalConsumptionOtherScaling | 1. | Constant | |
uncertaintyConsumptionQ2 | 0.67 | Estimate | |
uncertaintyConsumptionQ3 | 0.67 | Estimate | |
uncertaintyRecruitsCodRec | 1. | Estimate | |
uncertaintyRecruitsCapM | 0.41 | Estimate | |
uncertaintyRecruitsCapAge1Addition | 1. | Estimate |
Yearly parameters for historic assessment
Yearly parameters are recruitment of capelin as age 1, residual mortality of capelin and recruitment of cod as age 0. For practical reasons the recruitment parameters are implemented as addition to the original quantities. The yearly parameters are tagged with the year in which the are compared to data for estimation, not the year in which they are applied.
capelinAge1Addition | capelinMCorrection | recruitmentCodCorrection | |
Year | capelinAge1Addition is a list where each etry is a year and an addition to the number of 1 year old capelin the year before. | capelinMCorrection is a list where each entry is a year and an addition in the year before to the default monthly residual mortality of 0.05. | codRecruitmentCorrection is a list where each entry is a year and an addition to the number of cod as 0 year old 3 years before. |
1973 | 458.201 | 0.004 | |
1974 | 500.291 | -0.0196 | |
1975 | 368.46 | -0.0185 | 0. |
1976 | 260.3 | -0.0118 | 10.9991 |
1977 | 166.737 | -0.027 | 4.41333 |
1978 | 389.834 | -0.026 | 5.96334 |
1979 | 601.17 | 0.0239 | 1.835 |
1980 | 394.036 | -0.05 | 1.00667 |
1981 | 106.719 | -0.0483333 | 0.568334 |
1982 | 298.858 | 0.0684333 | 0.496667 |
1983 | 662.178 | 0.218083 | 0.80875 |
1984 | 48.8158 | 0.0939167 | 1.23333 |
1985 | 219.77 | 0.20875 | 1.41833 |
1986 | 20.4086 | 0.293999 | 3.335 |
1987 | 53.8173 | 0.464166 | 1.85698 |
1988 | 32.0684 | 0.01225 | 1.32084 |
1989 | 105.847 | 0.191666 | 1.1 |
1990 | 72.9542 | -0.0483333 | 1.08375 |
1991 | 313.391 | -0.01 | 1.53667 |
1992 | 319.403 | 0.0955003 | 2.06626 |
1993 | 162.323 | 0.214582 | 2.78501 |
1994 | 35.1667 | 0.20775 | 4.72501 |
1995 | -2.73959 | -0.05 | 7.80667 |
1996 | 18.6666 | -0.05 | 7.66833 |
1997 | -1.97636 | -0.05 | 9.62413 |
1998 | 66.2829 | -0.00836673 | 9.90998 |
1999 | -23.1354 | -0.0483333 | 3.94 |
2000 | 7.42916 | -0.049 | 2.74334 |
2001 | -125.514 | -0.05 | 1.67334 |
2002 | 146.028 | 0.0425833 | 1.42167 |
2003 | 1.58764 | 0.158 | 2.85833 |
2004 | -33.2323 | -0.05 | 2.04875 |
2005 | 57.0076 | 0.189749 | 4.76541 |
2006 | 16.3153 | -0.05 | 6.05667 |
2007 | 64.3688 | -0.05 | 12.5425 |
2008 | 291.038 | -0.0483333 | 12.321 |
2009 | 78.7896 | -0.0483333 | 6.44 |
2010 | 225.62 | -0.0483333 | 3.76166 |
2011 | 267.863 | -0.0463333 | 6.25326 |
2012 | 371.062 | -0.0484954 | 6.875 |
Settings
Settings are defined as quantities that are used to administer estimations and runs by selecting options for the use of different versions of equations. Also settings can be quantities that in principle could be estimated from data, but in practice are not. Settings are defined in the initialization chapter, but some of the settings can get their values changed from the menu.
Settings for parameter estimations - historic assessment
Name | Value | Explanation |
runOverallEstimationMethod | simultaneous | The setting runOverallEstimationMethod decides how the estimation of parameters is rigged. The variable is set in the estimation menu. The values are: stepwise simultaneous |
bifrostEstimationMethod | Custom | The setting bifrostEstimationMethod originally denoted the method the inbuilt Mathematica estimation algorithm should use, e.e. "NelderMead". It can also take the value "Custom", in which case the custom made method is used. |
basicEstimationMethod | nminimize | The setting basicEstimationMethod selects the main estimation scheme if Mathematica inbuilt estimation is used. Possible values are nminimize, findfit and findminimum. This setting is not relevant if bifrostEstimationMethod is set to "Custom" |
stagesToPerform | {1} | stagesToPerform is a list of estimation stages to be performed by the routine bifrostEstimationWork. The stages and assumptions and parameters connected to each stage is defined in the lists stageSpecificSettings and estimationDeviations, defined in the uppermost project definition routine. |
componentsInLikelihoodOriginal | {{capelinAbundanceLikelihoodComponent, consumptionLikelihoodComponent, capelinRecruitsLikelihoodComponent, capelinMLikelihoodComponent, codRecruitsLikelihoodComponent}} | componentsInLikelihoodOriginal is defined at initialization of each project together with stagesToPerform. For each stage a vector of likelihood components is defined. At each stage the setting componentsInLikelihood is used, taken from componentsInLikelihoodOriginal. |
excludeFromLikelihoodCalculation | {} | excludeFromLikelihoodCalculation is a list of likelihood components not to be used when the likelihood is evaluated. |
useAge3InEstimation | False | If useAge3InEstimation is True the number of capelin at age 3 (sex divided or sex aggregated depending on the setting of sumSexForEstimation) is used for constructing the likelihood function when the maturation parameters are estimated. |
useAge4InEstimation | True | If useAge4InEstimation is True the number of capelin at age 4 (sex divided or sex aggregated depending on the setting of sumSexForEstimation) is used for constructing the likelihood function when the maturation parameters are estimated. |
sumSexForEstimation | True | If sumSexForInEstimation is True number a age for female and male capelin are summed both for the data and simulation before the likelihood is evaluated. This is relevant for the estimation of maturation parameters. If the maturation parameters do not differ by sex (are bound to each other) it would be natural to set sunSexForEstimation to True. |
useExogeneousCapelinCV | False | If useExogeneousCapelinCV is True the uncertainty in the capelin measurements used when constructing terms in the likelihood function is taken from uncertainties calculated from bootstrap replicates of the September estimate. |
exogeneousCapelinCV | 0.2 | If useExogeneousCapelinCV is False exogeneousCV is an uncertainty applied to the input September data to construct replicates of the data. |
lastCapelinEstimationYear | 1980 | lastCapelinEstimationYear is the last year where data from the September survey is used for estimating maturation parameters for capelin. |
runLastYearSuppressCodPredationCapelinAutumn | 1983 | Overlap between cod and capelin in the autumn is a rather new feature of the Barents Sea ecosystem and supposing predation of capelin by cod is small in the autumn allows for estimating maturation parameters for capelin without being bothered with the parameters being confounded with the predation parameters for cod. |
runLastYearNoConsumptionImmatureCapelinQ1 | 1983 | |
runFirstYearFullConsumptionImmatureCapelinQ1 | 1995 | After a liner increase from runLastYearNoConsumptionImmatureCapleinQ1 runFirstYearFullConsumptionImmatureCapelinQ1 is the first year when cod consumes capelin in Q1 according to parameter values. This setting is also used for conssumption of capelin by mature cod in Q1. |
capelinStockDistributionFunction | gamma | capelinStockDistributionFunction is the distributinal assumption for terms in the likelihood function used when the capelin maturation parameters are estimated. |
consumptionByCodDistributionFunction | gamma | consumptionByCodDistributionFunction is the distributional assumption for terms in the likelihood function when the predation parameters are estimated. |
capelinMDistributionFunction | gamma | capelinMDistributionFunction is the distributional assumption for terms in the likelihood function when the yearly residual natural mortality for immature capelin is estimated. |
capelinRecruitsDistributionFunction | SSQ | capelinRecruitsDistributionFunction is the distributional assumption for terms in the likelihood function when the yearly adjustment to number of capelin at age 1 is estimated. |
codRecruitsDistributionFunction | lognormal | codRecruitsDistributionFunction is the distributional assumption for terms in the likelihood function when the yearly adjustment to number of cod at age 0 is estimated. |
runExcludeLowCatchYears | No | If runExcludeLowCatchYears is Yes years with low capelin stock as defined in the hardcoded list excludeConsumptionYears are excluded from the likelihood. |
runTotalConsumptionInEstimation | No | If runTotalConsumptionInEstimation is Yes, the consumption per cod is meaned over cod age both in data and model before the likelihood is evaluated. |
minAgeCodConsumption | 3 | minAgeCodConsumption is the smallest age of cod for terms in the likelihood. |
maxAgeCodConsumption | 9 | maxAgeCodConsumption is the largest age of cod for terms in the likelihood. |
smallestAcceptableConsumptionMeasurement | 0.0001 | smallesAcceptableConsumptinMeasurement is the smallest value of consumption per cod for inclusion of terms in the likelihood. |
runConsumptionPerCodWeightingMethod | None | The variable runConsumptionPerCodWeightingMethod decides how the consumption per cod calculations have beeb weighted in order to account for the distribution of the cod stock. The values are: None ByNumberInHaul ByLengthStation ByMultspecArea |
quartersCodPredationInLikelihood | {1, 3, 4} | quartersCodPredationInLikelihood is a list of quarters for which empirical consumption per cod is used in the likelihood. |
includeCodInLikelihood | True | On True, the consumption of cod by cod is included in the likelihood. |
includeOtherFoodInLikelihood | True | On True, the consumption of other food by cod is included in the likelihood. |
runSplitConsumptionQ1OnMaturation | Yes | On Yes, runSplitConsumptionQ1OnMaturation leads to the consumption of capelin per cod being calculated in quarter 1 on mature and immature capelin separately. Mature and immature capelin are separated in the model and is separated in the empirical consumption data using capelin length, where capelin in the stomachs that is not measured for length are distributed over length groups proportionally to the capelin that has been measured for length. |
runUseCapelinCodSuitability | Adjustment | The variable runUseCapelinCodSuitability decides the source for the suitability of capelin for cod. The variable is set in the etsimation menu. The values are: Default A hardcoded vector by cod age StomachSamples An age-age matrix calculated when the coonsumption per cod calculaytions are shown Adjustment A vector over cod age is represented by a beta function with parameters that are estimated |
runUseCodCodSuitability | Adjustment | The variable runUseCodCodSuitability decides the source for the suitability of cod for cod. The variable is set in the estimation menu. The values are: AgeAgeOverlap An age-age suitability is based on ooverlaps estimated from biological samples StomachSamples An age-age matrix calculated when the coonsumption per cod calculations are shown Adjustment A vector over cod age is represented by a beta function with parameters that are estimated |
runUseCodCodSuitability | Adjustment | The variable runUseCodCodSuitability decides the source for the suitability of cod for cod. The variable is set in the estimation menu. The values are: AgeAgeOverlap An age-age suitability is based on ooverlaps estimated from biological samples StomachSamples An age-age matrix calculated when the coonsumption per cod calculations are shown Adjustment A vector over cod age is represented by a beta function with parameters that are estimated |
Settings for simulation
Name | Value | Explanation |
useCapelinGrowthMaturationModel | True | If useCapelinGrowthMaturationModel is True the weigth at age ana maturation at age for capelin is determined by a regression model on historic data, else weight at age is drawn from historic data in the year historicYearUsedPrognostic. |
applyBifrostCodWeightModel | True | If applyBifrostCodWeightModel is True the weigth at age for cod is determined by a regression model on historic data, else weight at age is drawn from historic data in the year historicYearUsedPrognostic. |
applyBifrostCodOgiveModel | True | If applyBifrostCodOgiveModel is True the maturation at age for cod is determined by a regression model on historic data, else weight at age is drawn from historic data in the year historicYearUsedPrognostic. |
adaptiveCodCatch | False | If adaptiveCodCatch is set to Yes the catch of cod is adjusted so that the total cod biomass does not exceed maximumCodBiomass. |
adaptiveHerringCatch | False | If adaptiveHerringCatch is True, the catch is adjusted so that the total biomass of the herring stock is kept below maximumHerringBiomass. This is done in the SeaStar notebook, routine oneYearProbabilistic |
mainRandomSeed | 1 | The random generator is reset with mainRandomSeed prior to any main routines, i.e. in the routine performProjectTask. It is thus possible to check whether any conclusions are dependent on the random generator seed by rerunning using another value of mainRandomSeed. |
bottomCapelinSpawning | 0.02 | bottomCapelinSpawning is the contribution to the spawning stock of capelin from the bottom dwelling part of the capelin stock (million tonnes). This stock component is not observed by the echosounder in the autumn only in bottom trawl hauls. It is entered to the spawning stock at January 1, so it is a part of the fishery. |
runSmallestCapelinSSB | 0.001 | runSmallestCapelinSSB is the lowest allowed value of the capelin spawning stock. In the routine natmortMatureCapelin the mature capelin biomass is also to not become smaller than runSmallesCapelinSSB. |
largestCapelinCatch | 3. | largestCapelinCatch is the largest allowed capelin catch (million tonnes) |
catchDivisionByMonth | {{0.2, 0.8}, {0.5, 0.5}, {0.3, 0.3, 0.4}, {1.}, {1.}, {1.}, {1.}, {1.}, {1.}, {1.}, {1.}, {1.}} | In months where the number of capelin may change fast because of high predation and high catch (mostly during January-March) the time step (month) may be sub-divided and a portion of catch is instantly in the middle of each subdivision, while predation is gradually through each half of the subdivision (before and after catch). Each entry in the list catchDivisionByMonth is a list of how the catch during that month is partioned. The number of entries are 12, one for each month, where the first entry corresponds to January. |
numberHistoricRuns | 20 | The setting numberHistoricRuns gives the number of historic replicates. |
numberHCRRuns | 20 | The setting numberHCRRuns is the number of long-term simulations for each series of reference points that is used during estimation of harvesting control rules. |
numberHistoricReplicatesForPrognostic | 20 | |
prognosticSimulationMode | randomYear | prognosticSimulationMode governs how the historic data are used in simulations and takes the values: randomYear Historic data are picked at random cyclicYear Cycling through historic data, with random starting point last using the last historic year fixedYear using a specific year |
replaceMatureWeight | True | If replaceMatureWeight is True, the weight at age for mature capelin is replaced by weight from the catches after January 1. |
doCalculateCatch | True | In some cases it may be convenient to skip the calculation of catch at the end of the spawning season of capelin. This is achieved by setting the variable doCalculateCatch to False |
overrideCatchCalculation | False | If one wants to determine the yearly catch exogeneously, i.e. it is not calculated from some haresting control rule, overrideCatchCalculation must be set to True. Then the catch which is applied by the end of the spawning season is equal to capelinCatch. capelinCatch becomes updated, however, by the routine prespawningCapeln for small values of the spawning stock. |
targetRatio | 0.5 | When recruitmentOption is recruitment the catch during prognostic simulations is set so that the ratio between the realised recruitment and the maximum recruitment for the currently estimated recruitment function is targetRatio. |
recruitmentAgeCod | 3 | For ages below recruitmentAgeCod the mortality from cannibalism is simulated. The number at zero years is estimated by fitting to the number in the VPA at age recruitmentAgeCod. |
maxCapelinRecruitmentFactor | 2. | If a continuous uncertainty model is used (i.e. lognormal) there may occasionally occur very high recruitments. To avoid this, the recruitment is constrained to be below maxCapelinRecruitmentFactor maxCapelinRecruitment. |
plotCapelinRecruitment | False | If plotCapelinRecruitment is True recruitment plots are made during each simulation during evaluation of harvesting control rules. Also, useful textual output is produced. |
plotCodRecruitment | False | If plotCodRecruitment is True recruitment plots are made during each simulation during evaluation of harvesting control rules |
useExogeneosYoungHerring | False | |
exogeneousYoungHerring | 1. | exogeneousYoungHerring is the value used for the variable youngHerring in the recruitment function for capelin in case of the setting useExogeneousYoungHerring being True. |
useExogeneosZerogroupHerring | False | |
exogeneousZerogroupHerring | 0.5 | exogeneousZerogroupHerring is the value used for the variable zerogroupHerring in the recruitment function for capelin in case of the setting useExogeneousZerogroupHerring being True. |
Settings for HCR evaluation
During harvest rule evaluation prognostic simulations are made for a number of historic replicates and a number of prognostic runs for each historic replicate. The prognostic simulations are made using a 3-fold loop over harvest rules, the outermost on herring harvest rules, the middle loop on cod harvest rules and the inner loop on capelin harvest rules.
Name | Value | Explanation |
runApplyTargetSSB | Yes | If runApplyTargetSSB is set to Yes, the HCR evaluation uses runs where a target SSB for capelin is varied from runMinTargetSSB to runMaxTargetSSB in steps of runDelTargetSSB, all accessible from the menu. |
runApplyProportionalToProduction | Yes | If runApplyProportionalToProduction is set to Yes, the HCR evaluation uses runs where a the capelin catch is set proprtional to capelin biomass in each age groups where the proportionality factor is varied from is varied from runMinProductionProportion to runMaxProductionProportion in steps of runDelProductionProportion, all accessible from the menu. |
runApplyPrecautionaryApproach | Yes | If runApplyPrecautionaryApproch is set to Yes, a prognostic simulation run applying the present harvest rule and the parameters runRiskLimitCapelin and runRiskLevelCapelin can be set from the menu. |
runApplyPrecautionaryApproach | Yes | If runApplyPrecautionaryApproch is set to Yes, a prognostic simulation run applying the present harvest rule and the parameters runRiskLimitCapelin and runRiskLevelCapelin can be set from the menu. |
adaptiveCodCatch | No | If adaptiveCodCatch is set to Yes the catch of cod is adjusted so that the total cod biomass does not exceed maximumCodBiomass. |