Eine Jahreskurve der Verdunstung kann mit der folgenden Formel beschrieben werden:

$$E_P = 0.5*(E_{MP}-E_B)*\left[1+sin\left(\frac{2*\pi*i}{365}-90\right)\right]+E_B$$

mit der mittleren potentiellen Tagesverdunstung $E_{MP}$ und dem Julianischen Tag $i$.

Jahreskurve.py

from pylab import *
def Jahresgang(t,S,Em,Eb):
    E = (Em-Eb)/2*(1+sin((2*3.14)*t/365-S))+Eb
    return E
t = arange(1.0, 365.0, 7.0)
plot(t, Jahresgang(t,90,5.0,0.5))
ytext = ylabel('Verdunstung (mm)')
xtext = xlabel('Zeit (Tag)')
show()