HIERARCHAL INDUCTIVE PROCESS MODELING AND ANALYSIS ...

{"N.conc": "-1 * 1/( N.toCratio * 12.0107)
* P.growth_rate * P.conc"}
);
lib.add_generic_process(
"limited_growth", "limited_growth",
[("P",[pe],1,1), ("N",[no3,fe],1,100), ("E",[ee],1,1)],
[("light_lim", ["P","E"], 0), ("nutrient_lim",["P","N"], 0)],
{},
{"P.growth_rate": "(1-E.ice) * P.max_growth
* exp(0.06933 * E.TH2O) * P.growth_lim"},
{}
);
# ------ P.growth_lim --
# there are multiple factors (and formulations of factors)
# that might limit growth.
# In this library nutrient and light limitations are combined
# into P.growth_lim using a minimum function
# so that only one operates at a time (i.e., they are substitutable).
# The disadvantage
# of this encoding is that it will not be possible to determine
# which factor is operating at a given time. Temperature
# is a multiplicative control factor encoded in the P.growth_rate
# equation, and in the present library we do not consider
# alternative temperature effect functions.
# --light lim --
lib.add_generic_process(
79