Skip to contents

Solve for equilibrium numbers at length by iterating the population dynamics.

Usage

get_equilibrium_numbers(
  n_iter = 200,
  tau = 0.5,
  R0,
  F0 = 0,
  recruitment_size_sl,
  growth_ytrsll,
  M_ytrsl,
  catch_ytrf,
  selectivity_ytrsfl,
  legal_ytrsfl,
  retained_ytrsfl,
  handling_mortality_y
)

Arguments

n_iter

the number of iterations to loop over.

tau

the time step.

R0

the average recruitment.

F0

the initial exploitation rate.

recruitment_size_sl

the size at recruitment (cm) by sex.

growth_ytrsll

the growth matrix by year, season, region, and sex.

M_ytrsl

the natural mortality by year, season, region, sex, and length.

catch_ytrf

catch by year, season, region, and fishery

selectivity_ytrsfl

the selectivity by year, season, region, sex, fishery, and length.

legal status by year, season, region, sex, fishery, and length.

retained_ytrsfl

legal status by year, season, region, sex, fishery, and length.

handling_mortality_y

legal status by year, season, region, sex, fishery, and length.

Value

a array of the initial numbers in the population by region, sex, and length.