get_tpm
Gets scaled age, length at puberty, birth for std model
Contents
Syntax
varargout = get_tpm (p, f, tel_b, tau)
Description
Obtains scaled ages, lengths at puberty, birth for the std model at constant food, temperature; Notice p-b sequence in output, due to the name of the routine. State at birth (scaled age, reserve density, length) can optionally be specified E.g. for female: [tau_b, f, l_b] with zoom factor z, but for male: [tau_b, f*z/z_m, l_b*z/z_m] with zoom factor z_m
Input
- p: 5-vector with parameters: g, k, l_T, v_H^b, v_H^p
- f: optional scalar with functional response (default f = 1)
- tel_b: optional 3-vector with scaled age, reserve density and length at birth
- tau: optional n-vector with scaled times since birth
Output
- tvel: optional (n,4)-array with time since birth, scaled maturity, reserve density and length
- tau_p: scaled age at puberty \tau_p = a_p * k_M
- tau_b: scaled age at birth \tau_b = a_b * k_M
- l_p: scaled length at puberty l_b = L_b/ L_m
- l_b: scaled length at birth l_p = L_/p L_m
- info: indicator equals 1 if successful, 0 otherwise
Remarks
If tel_b is specified and different from the DEB value for p and f, see get_tb, initial growth deviates from vBert The m in get_tpm stands for male, see get_tp for female
Example of use
get_tpm([.5, .1, 0, .01, .05]) or: tel_b = [t_b, f*z/z_m, l_b*z/z_m]; % for males assumes same absolute reserve density at birth pars_tpm = [g_m k l_T v_Hb v_Hx v_Hpm]; tau = t * k_M * TC; % -, scaled time since birth corrected for temperature [tvel, tau_p, tau_b, l_p, l_b, info] = get_tpm(pars_tpm, f, tel_b, tau);