Die Infiltration kann mit dem Modell nach Horton berechnet werden mit:

in mm/h mit: Infiltrationsrate, konstante Rate nach der Anfangszeit, der anfänglichen Infiltrationsrate und einem Abklingfaktor.

| Horton.py

from pylab import *
def Horton(t):
    fu = 3.0
    f0 = 12.0
    k = 0.5
    H = fu+(f0-fu)*exp(-k*t)
    return H
t = arange(0.0, 60.0, 0.5)
plot(t, Horton(t))
ytext = ylabel('Infiltrationsrate (mm/min)')
xtext = xlabel('Zeit (Min.)')
show()

Die kumulative Infiltration beträgt:

in mm.

| HortonInt.py

from pylab import *
def HortonInt(t):
    fu = 3.0
    f0 = 12.0
    k = 0.05
    Hi = fu*t+(f0-fu)/k*(1-exp(-k*t))
    return Hi
t = arange(0.0, 60.0, 0.5)
plot(t, HortonInt(t))
ytext = ylabel('Infiltration (mm)')
xtext = xlabel('Zeit (Min.)')
show()