Rev 209 | Rev 211 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 209 | Rev 210 | ||
|---|---|---|---|
| Line 142... | Line 142... | ||
| 142 | /** calculation for the CO2 response for the ambientCO2 for the water- and nitrogen responses given.
|
142 | /** calculation for the CO2 response for the ambientCO2 for the water- and nitrogen responses given.
|
| 143 | The calculation follows Friedlingsstein 1995 (see also links to equations in code)
|
143 | The calculation follows Friedlingsstein 1995 (see also links to equations in code)
|
| 144 | */
|
144 | */
|
| 145 | double SpeciesSet::co2Response(const double ambientCO2, const double nitrogenResponse, const double soilWaterResponse) const |
145 | double SpeciesSet::co2Response(const double ambientCO2, const double nitrogenResponse, const double soilWaterResponse) const |
| 146 | {
|
146 | {
|
| - | 147 | if (nitrogenResponse==0) |
|
| - | 148 | return 0.; |
|
| - | 149 | ||
| 147 | double co2_water = 2. - soilWaterResponse; |
150 | double co2_water = 2. - soilWaterResponse; |
| 148 | double beta = mCO2beta0 * co2_water * nitrogenResponse; |
151 | double beta = mCO2beta0 * co2_water * nitrogenResponse; |
| 149 | 152 | ||
| 150 | double r =1. + M_LN2 * beta; // NPP increase for a doubling of atmospheric CO2 (Eq. 17) |
153 | double r =1. + M_LN2 * beta; // NPP increase for a doubling of atmospheric CO2 (Eq. 17) |
| 151 | 154 | ||