Running Topmodel - Revision history

Wbuytaert at 16:26, 12 February 2009

2009-02-12T16:26:24Z

Wbuytaert at 16:11, 12 February 2009

2009-02-12T16:11:47Z

Wbuytaert at 17:54, 21 November 2008

2008-11-21T17:54:05Z

Wbuytaert at 22:30, 3 August 2008

2008-08-03T22:30:39Z

Wbuytaert at 18:11, 2 August 2008

2008-08-02T18:11:01Z

Wbuytaert at 09:45, 2 August 2008

2008-08-02T09:45:57Z

New page

[[Category:Projects]]
[[Category:Hydrology in R]]

===Introduction===

===Running the model===

Topmodel needs the following parameters:

* Qs0 = initial flow at timestep 0.
* LnTe
* Sr0
* SrMax
* td
* vch
* vr
* k0
* psi
* dtheta
* dt = timestep (h)

We need to construct a vector with the parameter values (the order is important):

> parameters<-c(qs0,lnTe,m,Sr0,Srmax,td,vch,vr,k0,psi,dtheta,dt)

Then the model is run with the inputs described above:

> Qsim <- topmodel(parameters, topidx, delay, rain, ET0)

===Special features===

* The model can be run with several parameter sets at the time. Then "parameters" should be a 2D matrix with a parameter set in each row.

* The NS efficiency is calculated directly when observed discharge ("Qobs") is passed to the function:

> NS <- topmodel(parameters, topidx, delay, rain ,ET0, Qobs)

@@ Line 15: / Line 15: @@
    library(topmodel)
-   data(Huagrahuma)
+   data(huagrahuma)
 ===Running the model===

@@ Line 49: / Line 49: @@
 * The Nash-Sutcliffe efficiency is calculated directly when observed discharge ("Qobs") is passed to the function:
-    NS <- topmodel(parameters, topidx, delay, rain ,ET0, Qobs)
+    NS <- topmodel(parameters, topidx, delay, rain ,ET0, Qobs = Qobs)
 ===References===

@@ Line 21: / Line 21: @@
 Topmodel needs the following parameters:
-* qs0 = Initial subsurface flow per unit area (m)
+* qs0 = Initial subsurface flow per unit area [m]
-* lnTe = log of the areal average of T0 (m2/h)
+* lnTe = log of the areal average of T0 [m2/h]
-* m = Model parameter
+* m = Model parameter [m]
-* Sr0 = Initial root zone storage deficit (m)
+* Sr0 = Initial root zone storage deficit [m]
-* Srmax = Maximum root zone storage deficit (m)
+* Srmax = Maximum root zone storage deficit [m]
-* td = Unsaturated zone time delay per unit storage deficit
+* td = Unsaturated zone time delay per unit storage deficit [h/m]
-* vch = channel flow outside the catchment catchment (m/h)
+* vch = channel flow outside the catchment catchment [m/h]
-* vr = channel flow inside catchment (m/h)
+* vr = channel flow inside catchment [m/h]
-* k0 = Surface hydraulic conductivity (m/h)
+* k0 = Surface hydraulic conductivity [m/h]
-* CD = capillary drive, see Morel-Seytoux and Khanji (1974)
+* CD = capillary drive, see Morel-Seytoux and Khanji (1974) [m]
-* dt = The timestep (hours)
+* dt = The timestep [hours]
 We need to construct a vector with the parameter values (the order is important):

@@ Line 9: / Line 9: @@
 In its current stage, the model does very little checking of input data. Feeding the function with unrealistic parameter values, such as very small channel flows (<50 m / hour) are known to crash the underlying compiled libraries (and probably R with it). So be careful (and let me know any problems you come across)!
 ===Running the model===
@@ Line 29: / Line 36: @@
     parameters<-c(qs0,lnTe,m,Sr0,Srmax,td,vch,vr,k0,CD,dt)
 Then the model is run with the inputs described above:
@@ Line 38: / Line 47: @@
 * The model can be run with several parameter sets at the time. Then "parameters" should be a 2-dimensional matrix with a parameter set in each row.
-* The NS efficiency is calculated directly when observed discharge ("Qobs") is passed to the function:
+* The Nash-Sutcliffe efficiency is calculated directly when observed discharge ("Qobs") is passed to the function:
     NS <- topmodel(parameters, topidx, delay, rain ,ET0, Qobs)