adapted FODE function in e-RK4 to (propably) be

able to handle negative outflux velocities
-> influx from pipeline

Ausgleichsbecken.py

@@ -39,7 +39,7 @@ def get_p_halfstep(p0, p1):
     p_halfstep = (p0+p1)/2
 
 def FODE_function(x, h, alpha, p, rho=1000., g=9.81):
-    f = x**2/h*alpha+g-p/(rho*h)
+    f = x*abs(x)/h*alpha+g-p/(rho*h)
     return f
 
 

Ausgleichsbecken_class_file.py

@@ -67,6 +67,7 @@ class Ausgleichsbecken_class:
 
 
     def e_RK_4(self):
+        # Update to deal with non constant pipeline pressure!
         yn = self.outflux/self.area_outflux
         h = self.level
         dt = self.timestep