wasim (version 1.1.2)

p.grid: Read and plot a grid-file output by WaSiM

Description

These two functions help to read and plot grid-files created by WaSiM. You need the gridasci-executable if to convert binary grid data to ascii grid data.

Usage

p.grid(grid, ...)
read.grid(file, record_size=4, endian="little")

Arguments

grid

A grid read with the read.grid function

...

Additional arguments passed to plot

file

WASIM grid output file in ASCII or binary format.

record_size

number of bytes used to represent a float value in the binary grid. Usually 4, but may be platform-dependent.

endian

see readBin

Value

read.grid returns a list:

head

dataframe of header information containing the fields ncols, nrows, xllcorner, yllcorner, cellsize, nodata\_value.

tab

The actual grid

Details

The file is opened as a binary grid unless the infile has one of the following extensions: *.asc, *.ascii, *.txt. Nodata values are represented by NAs.

References

WaSiM-ETH manual

See Also

read.results, p.flow_comp for more WaSiM related function. filled.contour for plotting options

Examples

  unzip(system.file("weisseritz.zip", package="wasim"))
  sd.grid <- read.grid( file= "weisseritz/sd__ww100.grd")  # Datumsreihe heisst jetzt d.datum
  p.grid(sd.grid, 
        color.palette=topo.colors,
	      zlim=c(-20,20),
        main="Saturation deficit at the Weisseritz")