growthRate.H

Go to the documentation of this file.

double growthRateConst();
double growthRateCO2(double &, double &);
double growthRateBA (double &, double &);

double growthRateConst()
{
    return 0.0955;
}

double growthRateCO2(double &CO2_l, double &wCO2_Max)
{
    double G0 = 1.0e-14;
    if (CO2_l > wCO2_Max)
    {
        return (G0*(CO2_l - wCO2_Max)/max(wCO2_Max,4.4e-4));
    }
    else
    {
        return (0.0);
    }
}

double growthRateBA(double &L_l, double &LMAX)
{
    double G0 = 1.0e-14;

    if (LMAX != 0.0 && L_l > LMAX)
    {
        return (G0*(L_l - LMAX)/LMAX);
    }
    else
    {
        return (0.0);
    }
}