Metainformationen zur Seite
  •  

Unterschiede

Hier werden die Unterschiede zwischen zwei Versionen angezeigt.

Link zu dieser Vergleichsansicht

Beide Seiten der vorigen RevisionVorhergehende Überarbeitung
Nächste Überarbeitung
Vorhergehende Überarbeitung
hydro:speichermodell [2022/05/13 16:53] ckuellshydro:speichermodell [2024/04/10 10:02] (aktuell) – Externe Bearbeitung 127.0.0.1
Zeile 1: Zeile 1:
-=== Speichermodell ===+==== Speichermodell ==== 
 + 
 +=== Theorie ===
  
 Das Speichermodell kann zur einfachen Modellierung von Einzugsgebieten verwendet werden. Es beruht auf den Formeln Das Speichermodell kann zur einfachen Modellierung von Einzugsgebieten verwendet werden. Es beruht auf den Formeln
  
-$$ Q_t = [N \cdot c ] \cdot (1-e^{(-\frac{t}{K})})$$ für den Anstieg und +$$ Q_t = (N \cdot C ) \cdot (1-e^{(-\frac{t}{K})})$$ für den Anstieg und 
  
 $$ Q_t = Q_{(t-t_e)} \cdot e^{ (-\frac{(t-t_e)}{K})}$$  $$ Q_t = Q_{(t-t_e)} \cdot e^{ (-\frac{(t-t_e)}{K})}$$ 
  
-für den Rückgang: +für den Rückgang. Dabei sind N der Niederschlag in mm, C der Abflusskoeffizient (0<=C<=1), t die Zeit, K der Rückgangsfaktor (entspricht der Dauer, in der der Abfluss um eine logarithmische Einheit zurückgeht). Die Variable $t_e$ steht für die Zeit, bei der der Spitzenabfluss erreicht wird und bei der der Niederschlag aufhört (Ende des Niederschlages).  
 + 
 +=== Funktion in R === 
 + 
 +Das Speichermodell wurde als Funktion in R programmiert. Die Funktion benötigt die Variablen N, C, K, $t_a$=(Anfangszeit, Standard ist 0.0), $t_e$=Ende des Niederschlages (muss größer als $t_a$ und kleiner als $t_m$ sein, $t_m$=Ende der Zeitreihe (muss größer als $t_e$ sein). Sie kann als ELS() aufgerufen werden. Dabei sind alle Werte mit Standardwerten vor belegtN=50, C=0.1, K=7, $t_a$=0.0, $t_e$=7 (Tage), $t_m$= 30 (Tage). Wenn die vorbelegten Werte geändert werden sollen muss die Funktion entsprechend aufgerufen werden, z. Bsp. ELS(N=30,C=0.2,K=10), alle Variablen, die nicht aufgerufen werden, bleiben im Standardzustand. 
 + 
 +Wenn die Eingabe fehlerhaft ist, wird eine Variable //errorflag// gesetzt und eine Fehlermeldung ausgegeben, die Erstellung der Graphik mit plot() erfolgt nicht.
  
 <code rsplus |Speichermodell.r> <code rsplus |Speichermodell.r>