getq3: optimize for catchability (q)
In DLMtool/DLMtool: Data-Limited Methods Toolkit

Function optimizes catchability (q, where F=qE) required to get to user-specified stock depletion

getq3(
  x,
  D,
  SSB0,
  nareas,
  maxage,
  N,
  pyears,
  M_ageArray,
  Mat_age,
  Asize,
  Wt_age,
  V,
  retA,
  Perr,
  mov,
  SRrel,
  Find,
  Spat_targ,
  hs,
  R0a,
  SSBpR,
  aR,
  bR,
  bounds = c(1e-05, 15),
  maxF,
  MPA,
  plusgroup,
  VB0,
  optVB
)

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