function f=dhdt(h)
A=0.8;
f=(Qin(h)-Qut(h))/A;