Subversion Repositories public iLand

Rev

Rev 671 | Rev 697 | Go to most recent revision | Show entire file | Ignore 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