#V3.30.17.01;_2021_06_15;_safe;_Stock_Synthesis_by_Richard_Methot_(NOAA)_using_ADMB_12.3 #Stock Synthesis (SS) is a work of the U.S. Government and is not subject to copyright protection in the United States. #Foreign copyrights may apply. See copyright.txt for more information. #_user_support_available_at:NMFS.Stock.Synthesis@noaa.gov #_user_info_available_at:https://vlab.noaa.gov/group/stock-synthesis #C Control file for striped mullet stock assessment 2021 #C Base run #BEGIN CONTROL FILE INPUTS 0 # 0 means do not read wtatage.ss; 1 means read and use wtatage.ss and also read and use growth parameters 1 #_N_Growth_Patterns 1 #_N_platoons_Within_GrowthPattern #_Cond 1 #_Morph_between/within_stdev_ratio (no read if N_morphs=1) #_Cond 1 #vector_Morphdist_(-1_in_first_val_gives_normal_approx) #!!END BEGIN CONTROL FILE INPUTS!! #RECRUITMENT TIMING AND DISTRIBUTION 4 # recr_dist_method for parameters: 2=main effects for GP, Area, Settle timing; 3=each Settle entity 1 # not yet implemented; Future usage: Spawner-Recruitment: 1=global; 2=by area 1 # number of recruitment settlement assignments 0 # unused option #GPattern month area age (for each settlement assignment) 1 1 1 0 #!!END RECRUITMENT TIMING AND DISTRIBUTION!! #MOVEMENT #_Cond 0 # N_movement_definitions goes here if Nareas > 1 #_Cond 1.0 # first age that moves (real age at begin of season, not integer) also cond on do_migration>0 #_Cond 1 1 1 2 4 10 # example move definition for seas=1, morph=1, source=1 dest=2, age1=4, age2=10 #!!END MOVEMENT!! #BLOCKS 0 #_Nblock_Patterns #1 #_blocks_per_pattern # begin and end years of blocks # 1970 1970 #!!END BLOCKS!! #TIME-VARYING PARAMETER CONTROLS # controls for all timevary parameters 1 #_env/block/dev_adjust_method for all time-vary parms (1=warn relative to base parm bounds; 3=no bound check) # autogen 0 0 0 0 0 # autogen: 1st element for biology, 2nd for SR, 3rd for Q, 4th reserved, 5th for selex # where: 0 = autogen all time-varying parms; 1 = read each time-varying parm line; 2 = read then autogen if parm min==-12345 # #!!END TIME-VARYING PARAMETER CONTROLS!! #BIOLOGY # setup for M, growth, maturity, fecundity, recruitment distibution, movement # #NATURAL MORTALITY 2 #_natM_type:_0=1Parm; 1=N_breakpoints;_2=Lorenzen;_3=agespecific;_4=agespec_withseasinterpolate #_no additional input for selected M option; read 1P per morph 2 #_Cond 2 #_Reference_Age_for_Lorenzen_M if natM_type = 2 #!!END NATURAL MORTALITY!! #GROWTH 1 # GrowthModel: 1=vonBert with L1&L2; 2=Richards with L1&L2; 3=age_specific_K; 4=not implemented 1 #_Age(post-settlement)_for_L1;linear growth below this 999 #_Growth_Age_for_L2 (999 to use as Linf) -999 #_exponential decay for growth above maxage (fixed at 0.2 in 3.24; value should approx initial Z; -999 replicates 3.24) 0 #_placeholder for future growth feature 0 #_SD_add_to_LAA (set to 0.1 for SS2 V1.x compatibility) 0 #_CV_Growth_Pattern: 0 CV=f(LAA); 1 CV=F(A); 2 SD=F(LAA); 3 SD=F(A); 4 logSD=F(A) #!!END GROWTH!! #MATURITY-FECUNDITY 1 #_maturity_option: 1=length logistic; 2=age logistic; 3=read age-maturity matrix by growth_pattern; 4=read age-fecundity; 5=disabled; 6=read length-maturity #_Age_Maturity by growth pattern 1 #_First_Mature_Age 1 #_fecundity option:(1)eggs=Wt*(a+b*Wt);(2)eggs=a*L^b;(3)eggs=a*Wt^b; (4)eggs=a+b*L; (5)eggs=a+b*W #!!END MATURITY-FECUNDITY!! #HERMAPHRODITISM 0 #_hermaphroditism option: 0=none; 1=female-to-male age-specific fxn; -1=male-to-female age-specific fxn #!!END HERMAPHRODITISM!! #PARAMETER OFFSET METHOD 1 #_parameter_offset_approach (1=none, 2= M, G, CV_G as offset from female-GP1, 3=like SS2 V1.x) #!!END PARAMETER OFFSET METHOD!! #CATCH MULTIPLIER #!!END CATCH MULTIPLIER!! #AGEING ERROR PARAMETERS #!!END AGEING ERROR PARAMETERS!! #READ BIOLOGY PARAMETERS #_growth_parms #_ LO HI INIT PRIOR PR_SD PR_type PHASE env_var&link dev_link dev_minyr dev_maxyr dev_PH Block Block_Fxn # Sex: 1 BioPattern: 1 NatMort 0.01 1.0 0.53 0.53 0 0.8 -3 0 0 0 0 0.5 0 0 #Natmort_at_ref_age2_fem # Sex: 1 BioPattern: 1 Growth 10 40 21.0 21.0 2 1 2 0 0 0 0 0 0 0 # L_at_Amin_Fem_GP_1 20 70 48.7 48.7 2 1 4 0 0 0 0 0 0 0 # L_at_Amax_Fem_GP_1 0.05 0.8 0.401 0.401 2 1 4 0 0 0 0 0 0 0 # VonBert_K_Fem_GP_1 0.01 0.5 0.28 0.28 2 1 3 0 0 0 0 0 0 0 # CV_young_Fem_GP_1 0.01 0.5 0.21 0.21 2 1 3 0 0 0 0 0 0 0 # CV_old_Fem_GP_1 # Sex: 1 BioPattern: 1 WtLen -3 3 1.82e-05 1.82e-05 0.8 0 -3 0 0 0 0 0 0 0 # Wtlen_1_Fem_GP_1 -3 4 2.94 2.94 0.8 0 -3 0 0 0 0 0 0 0 # Wtlen_2_Fem_GP_1 # Sex: 1 BioPattern: 1 Maturity&Fecundity 1 100 31.9 31.9 0.8 0 -3 0 0 0 0 0 0 0 # Mat50%_Fem_GP_1 -3 3 -0.375 -0.375 0.8 0 -3 0 0 0 0 0 0 0 # Mat_slope_Fem_GP_1 -3 3 1 1 0.8 0 -3 0 0 0 0 0 0 0 # Eggs_scalar_Fem_GP_1 -3 3 0 0 0.8 0 -3 0 0 0 0 0 0 0 # Eggs_exp_len_Fem_GP_1 # Sex: 2 BioPattern: 1 NatMort 0.01 1.0 0.56 0.56 0 0.8 -3 0 0 0 0 0.5 0 0 #Natmort_at_ref_age2_mal # Sex: 2 BioPattern: 1 Growth 10 40 20.7 20.7 2 1 2 0 0 0 0 0 0 0 # L_at_Amin_Mal_GP_1 20 70 42.2 42.2 2 1 5 0 0 0 0 0 0 0 # L_at_Amax_Mal_GP_1 0.05 0.8 0.430 0.430 2 1 5 0 0 0 0 0 0 0 # VonBert_K_Mal_GP_1 0.01 0.5 0.25 0.25 2 1 3 0 0 0 0 0 0 0 # CV_young_Mal_GP_1 0.01 0.5 0.14 0.14 2 1 3 0 0 0 0 0 0 0 # CV_old_Mal_GP_1 # Sex: 2 BioPattern: 1 WtLen -3 3 2.02e-05 2.02e-05 0.8 0 -3 0 0 0 0 0 0 0 # Wtlen_1_Mal_GP_1 -3 4 2.91 2.91 0.8 0 -3 0 0 0 0 0 0 0 # Wtlen_2_Mal_GP_1 # Hermaphroditism # Recruitment Distribution # Cohort growth dev base 1 1 1 1 1 0 -1 0 0 0 0 0 0 0 #CohortGrowthDev # Movement # Age Error from parameters # catch multiplier # fraction female, by GP 1E-06 0.99 0.5 0.5 0.5 0 -99 0 0 0 0 0 0 0 #FracFemale_GP_1 #!!END READ BIOLOGY PARAMETERS!! #!!END BIOLOGY!! #TIME-VARYING BIOLOGY PARAMETERS #_no timevary MG parameters #!!END TIME-VARYING BIOLOGY PARAMETERS!! #SEASONAL BIOLOGY PARAMETERS #_seasonal_effects_on_biology_parms 0 0 0 0 0 0 0 0 0 0 #_femwtlen1,femwtlen2,mat1,mat2,fec1,fec2,Malewtlen1,malewtlen2,L1,K #_ LO HI INIT PRIOR PR_SD PR_type PHASE #_Cond -2 2 0 0 -1 99 -2 #_placeholder when no seasonal MG parameters #!!END SEASONAL BIOLOGY PARAMETERS!! #SPAWNER-RECRUITMENT 3 #_Spawner-Recruitment; Options: 2=Ricker; 3=std_B-H; 4=SCAA; 5=Hockey; 6=B-H_flattop; 7=survival_3Parm; 8=Shepherd_3Parm; 9=RickerPower_3parm 0 # 0/1 to use steepness in initial equ recruitment calculation 0 # future feature: 0/1 to make realized sigmaR a function of SR curvature #_ LO HI INIT PRIOR PR_SD PR_type PHASE env-var use_dev dev_mnyr dev_mxyr dev_PH Block Blk_Fxn # parm_name 6 20 10 10 1 1 1 0 0 0 0 0 0 0 # SR_LN(R0) 0.2 1.0 0.77 0.77 0.3 6 3 0 0 0 0 0 0 0 # SR_BH_steep 0 2 0.38 0.38 0.5 0 -4 0 0 0 0 0 0 0 # SR_sigmaR -5 5 0 0 0 0 -4 0 0 0 0 0 0 0 # SR_regime 0 0 0 0 0 0 -99 0 0 0 0 0 0 0 # SR_autocorr #_no timevary SR parameters 1 #do_recdev: 0=none; 1=devvector (R=F(SSB)+dev); 2=deviations (R=F(SSB)+dev); 3=deviations (R=R0*dev; dev2=R-f(SSB)); 4=like 3 with sum(dev2) adding penalty 1990 # first year of main recr_devs; early devs can preceed this era 2019 # last year of main recr_devs; forecast devs start in following year 2 #_recdev phase 1 # (0/1) to read 13 advanced options 0 #_recdev_early_start (0=none; neg value makes relative to recdev_start) -4 #_recdev_early_phase 0 #_forecast_recruitment phase (incl. late recr) (0 value resets to maxphase+1) 1 #_lambda for Fcast_recr_like occurring before endyr+1 1986.5 #_last_early_yr_nobias_adj_in_MPD 1996.3 #_first_yr_fullbias_adj_in_MPD 2018.0 #_last_yr_fullbias_adj_in_MPD 2019.2 #_first_recent_yr_nobias_adj_in_MPD 0.8535 #_max_bias_adj_in_MPD (1.0 to mimic pre-2009 models) #1900 #_last_early_yr_nobias_adj_in_MPD #1900 #_first_yr_fullbias_adj_in_MPD #2019 #_last_yr_fullbias_adj_in_MPD #2019 #_first_recent_yr_nobias_adj_in_MPD #1.0 #_max_bias_adj_in_MPD (1.0 to mimic pre-2009 models) 0 #_period of cycles in recruitment (N parms read below) -5 #min rec_dev 5 #max rec_dev 0 #_read_recdevs #_end of advanced SR options # #_placeholder for full parameter lines for recruitment cycles # read specified recr devs #_Yr Input_value # # all recruitment deviations # 1980R 1981R 1982R 1983R 1984R 1985R 1986R 1987R 1988R 1989R 1990R 1991R 1992R 1993R 1994R 1995R 1996R 1997R 1998R 1999R 2000R 2001R 2002R 2003R 2004R 2005R 2006R 2007R 2008R 2009F 2010F 2011F 2012F 2013F 2014F 2015F 2016F 2017F 2018F 2019F 2020F 2021F 2022F 2023F 2024F 2025F 2026F 2027F 2028F 2029F # 1.21353 -3.10888 -3.2481 -1.4731 0.820105 1.21806 0.839873 -0.583834 1.09464 0.684218 1.25871 1.19027 -0.701271 -0.39778 -0.268435 0.00841725 0.514882 0.555883 -0.200102 0.0905353 0.213104 0.0378078 -0.0315715 -0.138329 0.237996 0.0140324 0.163954 0.121635 -0.126248 -0.421511 0.261375 -0.149551 -0.187187 0.252072 -0.238409 -0.308064 -0.535089 -1.91311 -0.237954 -1.48106 0 0 0 0 0 0 0 0 0 0 # implementation error by year in forecast: 0 0 0 0 0 0 0 0 0 0 # #!!END SPAWNER-RECRUITMENT!! #FISHING MORTALITY METHOD #Fishing Mortality info 0.3 # F ballpark -2001 # F ballpark year (neg value to disable) 3 # F_Method: 1=Pope; 2=instan. F; 3=hybrid (hybrid is recommended) 1.0 # max F or harvest rate, depends on F_Method # no additional F input needed for Fmethod 1 # if Fmethod=2; read overall start F value; overall phase; N detailed inputs to read # if Fmethod=3; read N iterations for tuning for Fmethod 3 4 # N iterations for tuning F in hybrid method (recommend 3 to 7) #!!END FISHING MORTALITY METHOD!! #INITIAL FISHING MORTALITY #_initial_F_parms; count = 1 #_ LO HI INIT PRIOR PR_SD PR_type PHASE 0 1 0.4 0.1 0.05 0 1 # InitF_Comm 0 1 0.4 0.1 0.05 0 1 # InitF_Rec #!!END INITIAL FISHING MORTALITY!! #????? # F rates by fleet # Yr: 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 # seas: 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 # MRIP 0.188552 0.215992 0.160011 0.152514 0.132554 0.199041 0.245344 0.248651 0.22926 0.260397 0.173963 0.261557 0.203192 0.257605 0.287578 0.272366 0.407635 0.457205 0.347377 0.412134 0.654667 0.407542 0.597972 0.573779 1.19488 2.4629 1.46088e-005 0.0679173 0.159332 0.21941 0.241833 0.251529 0.256297 0.258815 0.259835 0.259835 # NC_Comm_Other 0.00122898 0.00172276 0.00145903 0.00102129 0.00122861 0.00156665 0.00218629 0.0015919 0.00131744 0.00200175 0.00367859 0.00218588 0.00244495 0.00348454 0.00291455 0.00519237 0.00743941 0.00550473 0.0062239 0.00860971 0.0097198 0.00709691 0.00780524 0.00880144 0.00980297 0.0392972 2.16631e-007 0.00100713 0.00236271 0.00325359 0.0035861 0.00372989 0.00380059 0.00383793 0.00385304 0.00385304 # NC_Comm_Pound 0.000574441 0.000464023 0.000797989 0.000516214 0.000694653 0.000439042 0.000917392 0.000928647 0.0011053 0.000501333 0.000405581 0.000574933 0.000667604 0.000787523 0.00210877 0.0026215 0.00292854 0.00412782 0.00360727 0.00777179 0.00923072 0.00748974 0.00403148 0.0095389 0.00973415 0.0234587 1.58572e-007 0.000737217 0.00172949 0.00238161 0.002625 0.00273025 0.00278201 0.00280934 0.0028204 0.0028204 # FL_Comm 0.00691357 0.00584283 0.00420521 0.0049168 0.00478212 0.00414439 0.00729197 0.00613325 0.00639235 0.00649894 0.00602002 0.00882404 0.00892915 0.00873269 0.00999761 0.010424 0.00955567 0.00915407 0.0123916 0.0170562 0.0193082 0.0172877 0.0192494 0.0295851 0.0369239 0.0994613 5.25939e-007 0.00244514 0.00573623 0.00789911 0.00870638 0.00905547 0.00922713 0.00931778 0.00935447 0.00935447 #????? #CATCHABILITY #_Q_setup for fleets with cpue or survey data #_1: fleet number #_2: link type: (1=simple q, 1 parm; 2=mirror simple q, 1 mirrored parm; 3=q and power, 2 parm; 4=mirror with offset, 2 parm) #_3: extra input for link, i.e. mirror fleet# or dev index number #_4: 0/1 to select extra sd parameter #_5: 0/1 for biasadj or not #_6: 0/1 to float #_ fleet link link_info extra_se biasadj float # fleetname 3 3 0 0 0 1 # P915 -9999 0 0 0 0 0 #_Q_parms(if_any);Qunits_are_ln(q) #_LO HI INIT PRIOR PR_SD PR_type PHASE env-var use_dev dev_mnyr dev_mxyr dev_PH Block Blk_Fxn # parm_name -25 25 0 0 10 0 -8 0 0 0 0 0 0 0 # LnQ_base_P915(3) -25 25 0 0 10 6 9 0 0 0 0 0 0 0 # Q_power_P915(3) #_no timevary Q parameters # #!!END CATCHABILITY!! #SELECTIVITY AND DISCARD #SPECIFY SIZE SELECTIVITY #_size_selex_patterns #Pattern:_0; parm=0; selex=1.0 for all sizes #Pattern:_1; parm=2; logistic; with 95% width specification #Pattern:_5; parm=2; mirror another size selex; PARMS pick the min-max bin to mirror #Pattern:_15; parm=0; mirror another age or length selex #Pattern:_6; parm=2+special; non-parm len selex #Pattern:_43; parm=2+special+2; like 6, with 2 additional param for scaling (average over bin range) #Pattern:_8; parm=8; New doublelogistic with smooth transitions and constant above Linf option #Pattern:_9; parm=6; simple 4-parm double logistic with starting length; parm 5 is first length; parm 6=1 does desc as offset #Pattern:_21; parm=2+special; non-parm len selex, read as pairs of size, then selex #Pattern:_22; parm=4; double_normal as in CASAL #Pattern:_23; parm=6; double_normal where final value is directly equal to sp(6) so can be >1.0 #Pattern:_24; parm=6; double_normal with sel(minL) and sel(maxL), using joiners #Pattern:_25; parm=3; exponential-logistic in size #Pattern:_27; parm=3+special; cubic spline #Pattern:_42; parm=2+special+3; // like 27, with 2 additional param for scaling (average over bin range) #_discard_options:_0=none;_1=define_retention;_2=retention&mortality;_3=all_discarded_dead;_4=define_dome-shaped_retention #_Pattern Discard Male Special 1 0 0 0 # 1 Comm 24 0 0 0 # 2 Rec 1 0 0 0 # 3 P915 #!!END SPECIFY SIZE SELECTIVITY!! #SPECIFY AGE SELECTIVITY # #_age_selex_patterns #Pattern:_0; parm=0; selex=1.0 for ages 0 to maxage #Pattern:_10; parm=0; selex=1.0 for ages 1 to maxage #Pattern:_11; parm=2; selex=1.0 for specified min-max age #Pattern:_12; parm=2; age logistic #Pattern:_13; parm=8; age double logistic #Pattern:_14; parm=nages+1; age empirical #Pattern:_15; parm=0; mirror another age or length selex #Pattern:_16; parm=2; Coleraine - Gaussian #Pattern:_17; parm=nages+1; empirical as random walk N parameters to read can be overridden by setting special to non-zero #Pattern:_41; parm=2+nages+1; // like 17, with 2 additional param for scaling (average over bin range) #Pattern:_18; parm=8; double logistic - smooth transition #Pattern:_19; parm=6; simple 4-parm double logistic with starting age #Pattern:_20; parm=6; double_normal,using joiners #Pattern:_26; parm=3; exponential-logistic in age #Pattern:_27; parm=3+special; cubic spline in age #Pattern:_42; parm=2+special+3; // cubic spline; with 2 additional param for scaling (average over bin range) #_Pattern Discard Male Special 0 0 0 0 # 1 Comm 0 0 0 0 # 2 Rec 0 0 0 0 # 3 P915 #!!END SPECIFY AGE SELECTIVITY!! #SIZE SELECTIVITY PARAMETERS #_ LO HI INIT PRIOR PR_SD PR_type PHASE env-var use_dev dev_mnyr dev_mxyr dev_PH Block Blk_Fxn # parm_name # 1 Comm LenSelex 0 60 31.6845 33.5 1 0 5 0 0 0 0 0 0 0 # Size_inflection_Comm(1) 0.01 40 6.21338 2.9 1 0 6 0 0 0 0 0 0 0 # Size_95%width_Comm(1) # 2 Rec LenSelex 11 53 12 12 1 0 -3 0 0 0 0 0 0 0 # Size_DblN_peak_Rec(2) -15 10 -3.8 -3.8 1 0 -3 0 0 0 0 0 0 0 # Size_DblN_top_logit_Rec(2) -2 9 2 2 1 0 -4 0 0 0 0 0 0 0 # Size_DblN_ascend_se_Rec(2) -2 9 4.2 4.2 1 0 -4 0 0 0 0 0 0 0 # Size_DblN_descend_se_Rec(2) -1000 15 -999 -999 1 0 -2 0 0 0 0 0 0 0 # Size_DblN_start_logit_Rec(2) -1000 15 -999 -999 1 0 -2 0 0 0 0 0 0 0 # Size_DblN_end_logit_Rec(2) # 3 P915 LenSelex 0 60 28.5181 33.5 1 0 3 0 0 0 0 0 0 0 # Size_inflection_P915(3) 0.01 40 2.60311 2.9 1 0 4 0 0 0 0 0 0 0 # Size_95%width_P915(3) #!!END SIZE SELECTIVITY PARAMETERS!! # timevary selex parameters #_ LO HI INIT PRIOR PR_SD PR_type PHASE # parm_name # info on dev vectors created for selex parms are reported with other devs after tag parameter section #!!END SELECTIVITY AND DISCARD!! #TWO-DIMENSIONAL AUTO-REGRESSIVE SELECTIVITY 0 # use 2D_AR1 selectivity(0/1): experimental feature #_no 2D_AR1 selex offset used #!!END TWO-DIMENSIONAL AUTO-REGRESSIVE SELECTIVITY!! #TAG-RECAPTURE PARAMETERS # Tag loss and Tag reporting parameters go next 0 # TG_custom: 0=no read; 1=read if tags exist #_Cond -6 6 1 1 2 0.01 -4 0 0 0 0 0 0 0 #_placeholder if no parameters #!!END TAG-RECAPTURE PARAMETERS!! # no timevary parameters # #VARIANCE ADJUSTMENT FACTORS # Input variance adjustments factors: #_1=add_to_survey_CV #_2=add_to_discard_stddev #_3=add_to_bodywt_CV #_4=mult_by_lencomp_N #_5=mult_by_agecomp_N #_6=mult_by_size-at-age_N #_7=mult_by_generalized_sizecomp #_Factor Fleet Value 4 1 2.0 4 2 1.6 4 3 3.7 5 1 1.4 5 3 2.0 -9999 1 0 # terminator #!!END VARIANCE ADJUSTMENT FACTORS!! #LAMBDAS 12 #_maxlambdaphase 1 #_sd_offset; must be 1 if any growthCV, sigmaR, or survey extraSD is an estimated parameter # read 0 changes to default Lambdas (default value is 1.0) # Like_comp codes: 1=surv; 2=disc; 3=mnwt; 4=length; 5=age; 6=SizeFreq; 7=sizeage; 8=catch; 9=init_equ_catch; # 10=recrdev; 11=parm_prior; 12=parm_dev; 13=CrashPen; 14=Morphcomp; 15=Tag-comp; 16=Tag-negbin; 17=F_ballpark; 18=initEQregime #like_comp fleet phase value sizefreq_method 5 1 1 0.25 1 5 3 1 0.25 1 -9999 1 1 1 1 # terminator #!!END LAMBDAS!! # lambdas (for info only; columns are phases) #CONTROLS FOR VARIANCE OF DERIVED QUANTITIES 0 # (0/1) read specs for more stddev reporting # 1 1 -1 5 1 5 1 -1 5 # selex type, len/age, year, N selex bins, Growth pattern, N growth ages, NatAge_area(-1 for all), NatAge_yr, N Natages # 5 15 25 35 43 # vector with selex std bin picks (-1 in first bin to self-generate) # 1 2 14 26 40 # vector with growth std bin picks (-1 in first bin to self-generate) # 1 2 14 26 40 # vector with NatAge std bin picks (-1 in first bin to self-generate) #!!END CONTROLS FOR VARIANCE OF DERIVED QUANTITIES!! 999