Project a population forward in time using calc_population() with constant recruitment and
seasonal dynamics (growth, movement-by-season) to obtain per recruit parameters. Note that the fishing
mortality among fleets and stocks remain linked by matrix q_fs.
calc_phi_project(
ny,
nm,
na,
nf = 1,
nr,
ns = 1,
F_mfr = array(0, c(nm, nf, nr)),
sel_mafs = array(1, c(nm, na, nf, ns)),
fwt_mafs = array(1, c(nm, na, nf, ns)),
q_fs = matrix(1, nf, ns),
M_as,
mov_marrs,
mat_as,
fec_as,
m_spawn = 1,
m_advanceage = 1,
delta_s = rep(0, ns),
natal_rs = matrix(1, nr, ns),
recdist_rs = matrix(1/nr, nr, ns)
)A named list returned by calc_population().
Integer, number of years for the projection
Integer, number of seasons
Integer, number of age classes
Integer, number of fleets
Integer, number of regions
Integer, number of stocks
Equilibrium fishing mortality (per season). Matrix [m, f, r]
Selectivity by season, age, fleet, stock. Array [m, a, f, s]
Fishery weight array by season, age, fleet, stock. Array [m, a, r, r]. Can be used
calculate yield per recruit.
Relative catchability of stock s for fleet f. Defaults to 1 if missing. Matrix [f, s]
Natural mortality. Matrix [a, s]
Movement array [m, a, r, r, s]. If missing, uses a diagonal matrix (no movement among areas).
Maturity at age. Matrix [a, s]
Fecundity at age. Matrix [a, s]
Integer, season of spawning
Integer, season at which to advance integer year age classes
Numeric vector by s. Fraction of season that elapses when spawning occurs, e.g., midseason spawning when delta_s = 0.5.
Matrix [r, s]. The fraction of the mature stock s in region r that spawns at
time of spawning. See example in Dstock.
Matrix [r, s]. The fraction of the incoming recruitment of stock s that settles in region r.
The initial population vector will be the survival at age evenly divided by the number of regions nr.
calc_phi_simple()