Running Topmodel

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)