Rev 723 | Rev 779 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 723 | Rev 777 | ||
---|---|---|---|
Line 229... | Line 229... | ||
229 | double total_excess = 0.; |
229 | double total_excess = 0.; |
230 | for (; day<end; ++day, ++doy) { |
230 | for (; day<end; ++day, ++doy) { |
231 | // (1) precipitation of the day
|
231 | // (1) precipitation of the day
|
232 | prec_mm = day->preciptitation; |
232 | prec_mm = day->preciptitation; |
233 | // (2) interception by the crown
|
233 | // (2) interception by the crown
|
234 | prec_after_interception = mCanopy.flow(prec_mm, day->temperature); |
- | |
- | 234 | prec_after_interception = mCanopy.flow(prec_mm); |
|
235 | // (3) storage in the snow pack
|
235 | // (3) storage in the snow pack
|
236 | prec_to_soil = mSnowPack.flow(prec_after_interception, day->temperature); |
236 | prec_to_soil = mSnowPack.flow(prec_after_interception, day->temperature); |
237 | // save extra data (used by e.g. fire module)
|
237 | // save extra data (used by e.g. fire module)
|
238 | add_data.water_to_ground[doy] = prec_to_soil; |
238 | add_data.water_to_ground[doy] = prec_to_soil; |
239 | add_data.snow_cover[doy] = mSnowPack.snowPack(); |
239 | add_data.snow_cover[doy] = mSnowPack.snowPack(); |
Line 338... | Line 338... | ||
338 | /** Interception in crown canopy.
|
338 | /** Interception in crown canopy.
|
339 | Calculates the amount of preciptitation that does not reach the ground and
|
339 | Calculates the amount of preciptitation that does not reach the ground and
|
340 | is stored in the canopy. The approach is adopted from Picus 1.3.
|
340 | is stored in the canopy. The approach is adopted from Picus 1.3.
|
341 | Returns the amount of precipitation (mm) that surpasses the canopy layer.
|
341 | Returns the amount of precipitation (mm) that surpasses the canopy layer.
|
342 | @sa http://iland.boku.ac.at/water+cycle#precipitation_and_interception */
|
342 | @sa http://iland.boku.ac.at/water+cycle#precipitation_and_interception */
|
343 | double Canopy::flow(const double &preciptitation_mm, const double &temperature) |
- | |
- | 343 | double Canopy::flow(const double &preciptitation_mm) |
|
344 | {
|
344 | {
|
345 | // sanity checks
|
345 | // sanity checks
|
346 | mInterception = 0.; |
346 | mInterception = 0.; |
347 | mEvaporation = 0.; |
347 | mEvaporation = 0.; |
348 | if (!mLAI) |
348 | if (!mLAI) |