Rev 671 | Rev 697 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 671 | Rev 680 | ||
---|---|---|---|
Line 234... | Line 234... | ||
234 | double current_psi = psiFromHeight(mContent); |
234 | double current_psi = psiFromHeight(mContent); |
235 | mPsi[doy] = current_psi; |
235 | mPsi[doy] = current_psi; |
236 | 236 | ||
237 | // (5) transpiration of the vegetation (and of water intercepted in canopy)
|
237 | // (5) transpiration of the vegetation (and of water intercepted in canopy)
|
238 | // calculate the LAI-weighted response values for soil water and vpd:
|
238 | // calculate the LAI-weighted response values for soil water and vpd:
|
- | 239 | double interception_before_transpiration = mCanopy.interception(); |
|
239 | double combined_response = calculateSoilAtmosphereResponse( current_psi, day->vpd); |
240 | double combined_response = calculateSoilAtmosphereResponse( current_psi, day->vpd); |
240 | et = mCanopy.evapotranspiration3PG(day, climate->daylength_h(doy), combined_response); |
241 | et = mCanopy.evapotranspiration3PG(day, climate->daylength_h(doy), combined_response); |
- | 242 | // if there is some flow from intercepted water to the ground -> add to "water_to_the_ground"
|
|
- | 243 | if (mCanopy.interception() < interception_before_transpiration) |
|
- | 244 | add_data.water_to_ground[doy]+= interception_before_transpiration - mCanopy.interception(); |
|
241 | 245 | ||
242 | mContent -= et; // reduce content (transpiration) |
246 | mContent -= et; // reduce content (transpiration) |
243 | // add intercepted water (that is *not* evaporated) again to the soil (or add to snow if temp too low -> call to snowpack)
|
247 | // add intercepted water (that is *not* evaporated) again to the soil (or add to snow if temp too low -> call to snowpack)
|
244 | mContent += mSnowPack.add(mCanopy.interception(),day->temperature); |
248 | mContent += mSnowPack.add(mCanopy.interception(),day->temperature); |
245 | 249 |