This is done by running the dynamics for a specified number of years while keeping recruitment and other components constant to reach steady state. Then the reproductive efficiencies are retuned to achieve that level of recruitment.

steady(
  params,
  t_max = 100,
  t_per = 7.5,
  tol = 10^(-2),
  dt = 0.1,
  return_sim = FALSE,
  progress_bar = TRUE
)

Arguments

params

A MizerParams object

t_max

The maximum number of years to run the simulation. Default is 100.

t_per

The simulation is broken up into shorter runs of t_per years, after each of which we check for convergence. Default value is 7.5. This should be chosen as an odd multiple of the timestep dt in order to be able to detect period 2 cycles.

tol

The simulation stops when the relative change in the egg production RDI over t_per years is less than tol for every background species. Default value is 1/100.

dt

The time step to use in project().

return_sim

If TRUE, the function returns the MizerSim object holding the result of the simulation run. If FALSE (default) the function returns a MizerParams object with the "initial" slots set to the steady state.

progress_bar

A shiny progress object to implement a progress bar in a shiny app. Default FALSE.

Examples

if (FALSE) { params <- newTraitParams() params@species_params$gamma[5] <- 3000 params <- setSearchVolume(params) params <- steady(params) }