FUSE-RHydro tutorial 1

RHydro01.R

# Step A: install RHydro

if(!require(RHydro)) install.packages("RHydro",repos="http://R-Forge.R-project.org")
 

# Step B: prepare your data

library(RHydro)

temp <- read.csv("dummyData.csv")
DATA <-  zooreg(temp[,2:4], order.by=temp[,1])
 

# Step C: run a simulation

myMID <- 60

myParameterSet <- list(rferr_add = 0, 
                       rferr_mlt = 1, 
                       frchzne = 0.5, 
                       fracten = 0.6828, 
                       maxwatr_1 = 106.6406, 
                       percfrac = 0.5, 
                       fprimqb = 0.5, 
                       qbrate_2a = 0.025, 
                       qbrate_2b = 0.01, 
                       qb_prms = 0.1294, 
                       maxwatr_2 = 3839.8, 
                       baserte = 50, 
                       rtfrac1 = 0.75, 
                       percrte = 703.128, 
                       percexp = 4.8594, 
                       sacpmlt = 10, 
                       sacpexp = 5, 
                       iflwrte = 500, 
                       axv_bexp = 0.7039, 
                       sareamax = 0.25, 
                       loglamb = 7.5, 
                       tishape = 3, 
                       qb_powr = 5, 
                       timedelay = 1.3355)

myDELTIM <- 1

# You can calculate the effective rainfall as follows:

Qinst <- fusesma.sim (DATA, 
                      mid=myMID, 
                      modlist, 
                      deltim=myDELTIM, 
                      states=FALSE, 
                      fluxes=FALSE, 
                      fracstate0=0.25, 
                      myParameterSet$rferr_add, 
                      myParameterSet$rferr_mlt, 
                      myParameterSet$frchzne, 
                      myParameterSet$fracten, 
                      myParameterSet$maxwatr_1,
                      myParameterSet$percfrac, 
                      myParameterSet$fprimqb, 
                      myParameterSet$qbrate_2a, 
                      myParameterSet$qbrate_2b, 
                      myParameterSet$qb_prms, 
                      myParameterSet$maxwatr_2, 
                      myParameterSet$baserte, 
                      myParameterSet$rtfrac1, 
                      myParameterSet$percrte, 
                      myParameterSet$percexp, 
                      myParameterSet$sacpmlt, 
                      myParameterSet$sacpexp, 
                      myParameterSet$iflwrte, 
                      myParameterSet$axv_bexp, 
                      myParameterSet$sareamax, 
                      myParameterSet$loglamb,
                      myParameterSet$tishape,
                      myParameterSet$qb_powr)

# Now calculate the routed runoff:

Qrout <- fuserouting.sim(Qinst, 
                         mid=myMID, 
                         modlist=modlist, 
                         timedelay=myParameterSet$timedelay, 
                         deltim=myDELTIM)

# Plot the routed discharge:

plot(Qrout,type="l")