Error in the computation of the theoretical solution for tracer decrease test case
In the vnv script, the computation of the theoretical solution replaces also the numerical solution given by TELEMAC (which makes theoretical solution = numerical solution):
# Computing theory values for degradation law over time for each point of extraction
theory = data2[0,:]
for i in range(len(times)):
theory[i] = 200.*math.exp(-2.3/6.*times[i])
The initialization of "theory" could be replaced by: theory=[0.]*len(times) - or other solutions.
Edited by Sara Pavan