The steady state is found by running the dynamics while keeping reproduction
and other components constant until the size spectra no longer change much (or
until time t_max
is reached, if earlier). Then the reproduction parameters
are set to values that give the level of reproduction observed in that
steady state.
steady( params, t_max = 100, t_per = 1.5, dt = 0.1, tol = 0.1 * dt, return_sim = FALSE, preserve = c("reproduction_level", "erepro", "R_max"), progress_bar = TRUE )
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 |
dt | The time step to use in |
tol | The simulation stops when the relative change in the egg
production RDI over |
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. |
preserve |
|
progress_bar | A shiny progress object to implement a progress bar in a shiny app. Default FALSE. |
if (FALSE) { params <- newTraitParams() species_params(params)$gamma[5] <- 3000 params <- steady(params) plotSpectra(params) }