                # 0.0-30.0 days);
                26.5 #croute free scaling parameter [timestep^2/mm] (e.g., 0.0-50.0 days^2/mm);
                ]).reshape(15, 1)
incon = array([50, 0, 2.5, 2.5]).reshape(4, 1)


import tuwmodel as tuw

data = tuw.simulate(area=1, param=params, incon=incon, prec=prec,