Subversion Repositories public iLand

    double total_excess = 0.;

    double total_excess = 0.;

    for (; day<end; ++day, ++doy) {

    for (; day<end; ++day, ++doy) {

        // (1) precipitation of the day

        // (1) precipitation of the day

        prec_mm = day->preciptitation;

        prec_mm = day->preciptitation;

        // (2) interception by the crown

        // (2) interception by the crown

        // (3) storage in the snow pack

        // (3) storage in the snow pack

        prec_to_soil = mSnowPack.flow(prec_after_interception, day->temperature);

        prec_to_soil = mSnowPack.flow(prec_after_interception, day->temperature);

        // save extra data (used by e.g. fire module)

        // save extra data (used by e.g. fire module)

        add_data.water_to_ground[doy] = prec_to_soil;

        add_data.water_to_ground[doy] = prec_to_soil;

        add_data.snow_cover[doy] = mSnowPack.snowPack();

        add_data.snow_cover[doy] = mSnowPack.snowPack();

/** Interception in crown canopy.

/** Interception in crown canopy.

    Calculates the amount of preciptitation that does not reach the ground and

    Calculates the amount of preciptitation that does not reach the ground and

    is stored in the canopy. The approach is adopted from Picus 1.3.

    is stored in the canopy. The approach is adopted from Picus 1.3.

    Returns the amount of precipitation (mm) that surpasses the canopy layer.

    Returns the amount of precipitation (mm) that surpasses the canopy layer.

    @sa http://iland.boku.ac.at/water+cycle#precipitation_and_interception */

    @sa http://iland.boku.ac.at/water+cycle#precipitation_and_interception */

{

{

    // sanity checks

    // sanity checks

    mInterception = 0.;

    mInterception = 0.;

    mEvaporation = 0.;

    mEvaporation = 0.;

    if (!mLAI)

    if (!mLAI)