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.
|