qt = 1
cp = -5.44 # permil, matrix
ct = -6.33 # permil ground water, mixing
ce = -7.91 # permil, event or macropore flow
qe = qt*(ct-cp)/(ce-cp)
qe