Commit d496ef2b authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

simple parameter for particle mortality

parent deb4f339
......@@ -367,12 +367,12 @@ public:
* 1 : constant mortality rate given by [1] mortRateMax (by day)
* 2 : weibullMortalityRate with [2] lambda and [3] nu;
*/
void addMortality(fullMatrix<double>* mortParams) {
void addMortality(int type, double rateMax, double weibullLambda = 0, double weibullNu = 0) {
_mortFlag = true;
_mortType = static_cast<int>(mortParams->get(0,0)+0.1);
_mortRate_max = mortParams->get(0,1);
_mortLambda = mortParams->get(0,2); //(Connolly & Baird 2010)lambda = 1/b(Pinder et al. 1978)
_mortNu = mortParams->get(0,3); //(Connolly & Baird 2010)nu = c(Pinder et al. 1978)
_mortType = type;
_mortRate_max = rateMax;
_mortLambda = weibullLambda; //(Connolly & Baird 2010)lambda = 1/b(Pinder et al. 1978)
_mortNu = weibullNu; //(Connolly & Baird 2010)nu = c(Pinder et al. 1978)
}
/** Activate particle settling and load parameters */
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment