Subversion Repositories public iLand

Rev

Rev 547 | Rev 553 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 547 Rev 552
Line 150... Line 150...
150
{
150
{
151
    double deficit;
151
    double deficit;
152
    double ft, fw;
152
    double ft, fw;
153
    const double top_layer_content = mRU->waterCycle()->topLayerWaterContent();
153
    const double top_layer_content = mRU->waterCycle()->topLayerWaterContent();
154
    double f_sum = 0.;
154
    double f_sum = 0.;
155
    for (const ClimateDay *day=mRU->climate()->begin(); day!=mRU->climate()->end(); ++day)
-
 
-
 
155
    int iday=0;
-
 
156
    for (const ClimateDay *day=mRU->climate()->begin(); day!=mRU->climate()->end(); ++day, ++iday)
156
    {
157
    {
157
        deficit = mRU->waterCycle()->waterDeficit_mm(day->day);
-
 
-
 
158
        deficit = mRU->waterCycle()->waterDeficit_mm(iday);
158
159
159
        ft = exp(308.56*(1./56.02-1./((273.+day->temperature)-227.13)));  // empirical variable Q10 model of Lloyd and Taylor (1994), see also Adair et al. (2008)
160
        ft = exp(308.56*(1./56.02-1./((273.+day->temperature)-227.13)));  // empirical variable Q10 model of Lloyd and Taylor (1994), see also Adair et al. (2008)
160
        fw = 1. - limit(deficit / top_layer_content, 0., 1.);
161
        fw = 1. - limit(deficit / top_layer_content, 0., 1.);
161
        // the water effect: depends on the water deficit; if the deficit is higher than the parameterized
162
        // the water effect: depends on the water deficit; if the deficit is higher than the parameterized
162
        // content of the top layer (where most microbial activity is located), than then fw gets 0.
163
        // content of the top layer (where most microbial activity is located), than then fw gets 0.