read.results: read.results

In wasim: Visualisation and analysis of output files of the hydrological model WASIM

Description

Import WASIM simulation results

Usage

read.results(record, path = "wasim/output", 
                 ts.length, subcatchments, data.types.prefix =
                 data.types$prefix, data.types.has_stat =
                 data.types$has_stat, ending, endings = rep(ending,
                 NROW(data.types.prefix)), na.values = c("999",
                 "999.00"), read.files.nr = 1:nrow(data.types),
                 remove.missing = FALSE)

read.dates(file=NULL,na.strings=c("999", "999.00","-9999.00"),sep="\t", skip=3)

read.observations(filelist, path = "wasim/input", gauge_names = NULL,
                 date_modelled = NULL, na.strings = c("999", "999.00",
                 "9999", "9999.00", "9999.000000", "-9999", "-9999.0",
                 "-9999.00"), sep = "")

Arguments

file	WASIM output file the date information is read from. If not specified, the function tries to open the first available file from the globaldata.types\$filename.
record	subdirectory/ string to be appended to path (eg. for multiple runs, etc.)
path	path the output files reside in (no ending slash)
ts.length	number of records to be read
subcatchments	number of subcatchments (columns) in WASIM-output file (without statistic column)
data.types.prefix	prefixes for filenames (better use default)
data.types.has_stat	indicator whether files contain an extra column for summary statistics
ending	suffix for filenames
endings	suffixes for filenames. Use this if endings differ for various files.
read.files.nr	numeric vector indicating which entries from data.types.prefix to process
na.values	list of strings for indicating nodata values in the WASIM-files
na.strings	Strings in input file to be interpreted as nodata (NA)
sep	field separator in input file
remove.missing	boolean indication whether to remove an entry from the default data.types if a file was not found
skip	number of lines to skip from the header
filelist	list of input files to be considered
gauge_names	names of gauges that should be read (NULL:read all)
date_modelled	timespan that will be extracted

Details

read.results causes error messages: "...line n hasn't m elements..." The global variable data.types\$has\_stat specifies an statistics column in one of the output files where there is none.

"...has length 14 expecting length 15 - filling ..." The global variable data.types\$has\_stat specifies no statistics column in one of the output files where there seems to be one.

Value

date_vec

vector of POSIX-dates

Author(s)

Dominik Reusser, Till Francke

See Also

p.storage

Examples

  unzip(system.file("weisseritz.zip", package="wasim"))
d.datum <- read.dates(file= "weisseritz/qgesww100.txt")  # Datumsreihe heisst jetzt d.datum
str(d.datum)

d.wasim.out <- read.results("", "weisseritz",ts.length=21937,subcatchments=2, ending="ww100.txt")    # Datensatz heisst jetzt d.wasim.out 

str(d.wasim.out)

wasim documentation built on May 1, 2019, 6:50 p.m.