timml
timml copied to clipboard
traceline gets stuck
The following traceline gets stuck and doesn't reach the line-sink
z = np.arange(20, -0.1, -0.5)
ml = Model3D(z=z, kaq=10)
L = 50
x = np.linspace(0, 50, 2)
y = np.zeros_like(x)
HeadLineSinkString(ml, xy=list(zip(x, y)), hls=0, layers=27)
Constant(ml, 0, 100, 10)
ml.solve()
tr = timtraceline(ml, 100, 50, 6, hstepmax=5, nstepmax=1000, metadata=True)