InitSuction rhow*g_const*Y "Initial suction"
Suction -p*(dl.Hp<0) "Current suction"
PorePressure p*(dl.Hp>=0) "Pore pressure"
k_rel Se "Relative permeability"
k K_sat*muw/gammaw "Soil permeability"
Cm phi0*(S_sat-S_res)*Se*alpha "Specific moisture capacity"
Se exp(alpha*dl.Hp)*(dl.Hp<0)+1*(dl.Hp>=0) "Effective saturation"
Hp_in 1/alpha*log((dl.u/K_sat+1)*exp(-alpha*Y)+U_in/K_sat) "Gardner steady state solution"
Hp_out 1/alpha*log((dl.u/K_sat+1)*exp(-alpha*Y)-U_out/K_sat) "Gardner steady state solution"