Calculate Reference Yield
getFref3(x, Asize, nareas, maxage, N, pyears, M_ageArray, Mat_age, Wt_age,
V, retA, Perr, mov, SRrel, Find, Spat_targ, hs, R0a, SSBpR, aR, bR, MPA,
maxF, SSB0, plusgroup = 0)
Integer, the simulation number
A matrix (dimensions nsim by nareas) with relative size of areas
The number of spatial areas
The maximum age
Array of the numbers-at-age in population. Dimensions are nsim, maxage, nyears, nareas. Only values from the first year are used, which is the current N-at-age.
The number of years to project forward. Equal to 'nyears' for optimizing for q.
An array (dimensions nsim, maxage, nyears+proyears) with the natural mortality-at-age and year
An array (dimensions nsim, maxage, nyears+proyears) with the proportion mature for each age-class
An array (dimensions nsim, maxage, nyears+proyears) with the weight-at-age and year
An array (dimensions nsim, maxage, nyears+proyears) with the vulnerability-at-age and year
An array (dimensions nsim, maxage, nyears+proyears) with the probability retained-at-age and year
A matrix (dimensions nsim, nyears+proyears) with the recruitment deviations
An array (dimensions nsim, nareas, nareas) with the movement matrix
A numeric vector nsim long specifying the recruitment curve to use
A matrix (dimensions nsim, nyears) with the historical fishing effort
A numeric vector nsim long with the spatial targeting
A numeric vector nsim long with the steepness values for each simulation
A matrix (dimensions nsim, nareas) with the unfished recruitment by area
A matrix (dimensions nsim, nareas) with the unfished spawning-per-recruit by area
A numeric vector nareas long with the Ricker SRR a values
A numeric vector nareas long with the Ricker SRR b values
A matrix of spatial closures by year
A numeric value specifying the maximum fishing mortality for any single age class
SSB0
Integer. Default = 0 = no plus-group. Use 1 to include a plus-group
logical - use the CPP code? For testing purposes only