Plot a Walter Lieth climate diagram

Share:

Description

This produces a Walter Lieth climate diagram from annual cycles of temperature and precipitation.

Usage

1
2
3
plotWalterLieth(X, station="", alt=NA, per, margin=c(4, 4, 5, 4),
pcol="#005ac8", tcol="#e81800", pfcol="#79e6e8", sfcol =
  "#09a0d1", shem=FALSE, ...)

Arguments

X

A dataclim object or a matrix.

station

A string describing the station. Defaults to the empty string "".

alt

The altitude of the station. Defaults to NA.

per

A string for the period over which the monthly averages are computed. Defaults to the basePeriod slot if X is of class dataclim or to the empty string "" if X is a matrix

.

margin

Vector with four numeric elements describing the margins of the plot. Defaults to c(4, 4, 5,4). See par().

pcol

Color for (p)recipitation. Defaults to "#005ac8".

tcol

Color for (t)emperature. Defaults to "#e81800".

pfcol

Color (p)robable (f)rost. Defaults to "#79e6e8".

sfcol

Color (s)ure (f)rost. Defaults to "#09a0d1".

shem

Logical. Is the station on the Southern hemisphere? Defaults to FALSE.

...

Further arguments passed to plot().

Details

If X is an object of dataclim, the extracts the monthly averages from a call to summary(X) and produces the Walter Lieth climate diagram.

If X is a matrix, X is expected to be a 4x12 matrix with no NAs, with the first row: monthly precipitation (mm), second row: monthly average maximum daily temperature (degC), third row: monthly average minimum daily temperature (degC), forth row: monthly absolute minimum daily temperature.

Value

Nothing. A plot.

Note

The code is based on code from the former CRAN package climatol by Jose A. Guijarro <jguijarrop@aemet.es>.

Author(s)

Boris Orlowsky <boris@climate-babel.org>

References

Walter, H. & Lieth, H (1960): Klimadiagramm Weltatlas. G. Fischer, Jena.

See Also

dataclim

Examples

1
2
3
4
5
6
data(potsdam)
date <- as.Date(potsdam$date)
myDataclim <- createDataclim(date=date, tmin=potsdam$tmin, tmax=potsdam$tmax,
                          prec=potsdam$prec, basePeriod=1981:2010)
plotWalterLieth(myDataclim)
plotWalterLieth(myDataclim, station="Potsdam", alt=80)