read_95ppu {SWAT2R}R Documentation

read_95ppu

Description

This function reads the SWAT-CUP file ‘95ppu.sf2’, which contains the 95% of predictive uncertainty corresponding to the behavioural parameter sets, and optionally, plots the corresponding uncertainty bounds. Fields in '95ppu.sf2' are : observed L95PPU U95PPU Best_Sim M95PPU

Usage

read_95ppu(file, ua.algorithm, plot=TRUE, dates, date.fmt="%Y-%m-%d", 
           gof.leg= TRUE, gof.digits=2, main="95PPU Plots", ylab="Q, [m3/s]", 
           verbose=TRUE,...) 

Arguments

file

character with the name of the file which the data are to be read from. Valid values are in c('95ppu.sf2', '95ppu.out'). If it does not contain an absolute path, the file name is relative to the current working directory, getwd()

ua.algorithm

OPTIONAL, only required when file is missing.
character, indicating the SWAT-CUP algorithm that created file. Valid values are in c('GLUE', 'SUFI2')

plot

logical, indicating if the 95PPU and the observed values have to be plotted or not

dates

Dates for the correponding values in the 'sim' and 'obs' time series.
If dates is a factor, it have to be converted into ‘Date’ class, using the date format specified by date.fmt.
If dates is already of Date class, the number of dates must be equal to the number of elements in the simulated and obsserved values

date.fmt

character, with the format used to define Date.Ini, Date.Fin. See 'format' in 'as.Date'. Default value is "%Y-%m-%d"

gof.leg

logical, indicating if the p-factor and r-factor have to be computed and ploted as legends on the graph.

gof.digits

OPTIONAL, numeric. Only used when gof.leg=TRUE. Decimal places used for rounding p-factor and r-factor

main

character with an overall title for the plot: see 'title'.

ylab

character with a title for the y axis: see 'title'.

verbose

logical; if TRUE, progress messages are printed

...

additional argument to be passed to the plotbands function

Author(s)

Mauricio Zambrano-Bigiarini, mzb.devel@gmail.com

References

SWAT-CUP at http://www.eawag.ch/organisation/abteilungen/siam/software/swat/index_EN

See Also

read_sub, read_hru, read_rsv, plotbands


[Package SWAT2R version 0.1-11 Index]