traceline gets stuck

[dbrakenhoff]

The following traceline gets stuck and doesn't reach the line-sink

import timflow.steady as tfs
import numpy as np

z = np.arange(20, -0.1, -0.5)
ml = tfs.Model3D(z=z, kaq=10)
L = 50
x = np.linspace(0, 50, 2)
y = np.zeros_like(x)
tfs.HeadLineSinkString(ml, xy=list(zip(x, y)), hls=0, layers=27)
tfs.Constant(ml, 0, 100, 10)
ml.solve()
tr = tfs.timtraceline(ml, 100, 50, 6, hstepmax=5, nstepmax=1000, metadata=True)