dynatopmodel (version 1.2.1)

disp_output: Display output of a Dynamic TOPMODEL run

Description

Simple output of the results of a simulation.

Usage

disp_output(qsim, rain, evap = NULL, qobs = NULL, tm = NULL,
  start = min(index(qsim)), end = max(index(qsim)),
  par = get.disp.par(lwd.rain = 3, qint = 0.1), show.maxima = FALSE,
  freq = "day", pch.qsim = "+", pch.obs = "*", ...)

Arguments

qsim

Time series of simulated discharges.

rain

Time series of rainfall, at same interval as simulated values.

evap

Time series of evapotranspiration (optional), at same interval as simulated values.

qobs

Time series of evapotranspiration (optional), at same interval as simulated values.

tm

Display a vertical line at this time in the simulation. If NULL no line will be drawn.

start

Start time for plot in a format interpretable as POSIXct date time. Defaults start of simulated discharges.

end

End time for plot in a format interpretable as POSIXct date time. Defaults to end of simulated discharges.

par

Parameters controlling display output. A default set may be obtained through a call to disp.par.

show.maxima

Boolean Whether to show the daily maxima as points

freq

character Period for which maxima are identified, day by default

pch.qsim

character Symbol for plotting maxima of simulations, if stipulated

pch.obs

character Symbol for plotting maxima of observations, if these are supplied

...

Any further named parameters will be treated as graphics parameters and applied throughout the plot.

Details

This will render the hydrograph, any observations, actual evapotranspiration, if supplied, and the rainfall hyetograph.

See Also

disp.par

Examples

# NOT RUN {
# Show the output of the storm simulation, overriding label colours and vertical axis limits.
require(dynatopmodel)

data(brompton)

x11()
with(brompton$storm.run, disp_output(evap=ae*1000,qobs=qobs*1000,
                                     qsim=qsim*1000, rain=rain*1000,
                                     max.q=2, cex.main=1, col.axis="slategrey", las.time=1))
# }