Unterschiede

Hier werden die Unterschiede zwischen zwei Versionen der Seite angezeigt.

--- hydro:jahreskurve [2017/06/08 10:23] – ckuells
+++ hydro:jahreskurve [2018/08/12 23:38] – angelegt ckuells
@@ Zeile 1: / Zeile 1: @@
-# Jahreskurve
+====== Jahresgang ======
 Eine Jahreskurve der Verdunstung kann mit der folgenden Formel beschrieben werden:
-$$E_P = E_{MP}*\left[1+sin\left(\frac{360*i}{365}-90\right)\right]$$
+$$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 Jahresverdunstung $E_{MP}$ und dem Julianischen Tag $i$.
+mit der mittleren potentiellen Tagesverdunstung $E_{MP}$ und dem Julianischen Tag $i$.
 <code | Jahreskurve.py>
 from pylab import *
-from numpy import *
+def Jahresgang(t,S,Em,Eb):
-from scipy import *
+    E = (Em-Eb)/2*(1+sin((2*3.14)*t/365-S))+Eb
-from math import *
-def Jahresgang(t,S,Em):
-    E = Em*(1+sin(t/365*(2*3.14)-S))
     return E
-t = arange(0.0, 365.0, 1)
+t = arange(1.0, 365.0, 7.0)
-plot(t, Jahresgang(t,90,540))
+plot(t, Jahresgang(t,90,5.0,0.5))
 ytext = ylabel('Verdunstung (mm)')
 xtext = xlabel('Zeit (Tag)')