Water dynamics

change_water_reserve(;
    WR, precipitation, LAItot,
    PET, WHC, PWP)

TBW

evaporation(; WR, WHC, PET, LAItot)

Aev(t) = WR(t) / WHC * PET(t)*[1 - min(1; LAI_tot(t)/3) ]

transpiration(;
WR,
PWP, WHC,
PET,
LAItot)

TBW