get_lj
Gets scaled length at metamorphosis
Contents
Syntax
[lj, lp, lb, info] = get_lj(p, f, lb0)
Description
Type M-acceleration: Isomorph, but V1-morph between vHb and vHj This routine obtains scaled length at metamorphosis lj given scaled maturity at metamorphosis vHj. The theory behind get_lj, is discussed in the comments to DEB3. If scaled length at birth (third input) is not specified, it is computed (using automatic initial estimate); if it is specified. however, is it just copied to the (third) output. The code assumes vHb < vHj < vHp (see first input).
Input
- p: 6-vector with parameters: g, k, l_T, v_H^b, v_H^j, v_H^p
if p is a 5-vector, output lp is empty
- f: optional scalar with scaled functional responses (default 1)
- lb0: optional scalar with scaled length at birth
Output
- lj: scalar with scaled length at metamorphosis
- lp: scalar with scaled length at puberty
- lb: scalar with scaled length at birth
- info: indicator equals 1 if successful
Remarks
Similar to get_lj1, which uses root finding, rather than integration Scaled length l = L/ L_m with L_m = kap {p_Am}/ [p_M] where {p_Am} is of embryo, because the amount of acceleration is food-dependent so the value after metamorphosis is not a parameter. {p_Am} and v increase between maturities v_H^b and v_H^j till {p_Am} l_j/ l_b and v l_j/ l_b at metamorphosis. After metamorphosis l increases from l_j till l_i = f l_j/ l_b - l_T. Scaled length l can thus be larger than 1. See get_lj_foetus for foetal development.
Example of use
get_lj([.5, .1, .1, .01, .2, .3])