composite.field {clim.pact}R Documentation

Composite maps

Description

Produce composites of maps based on station observations, weights, or years. For a station object, assign +/- 1 weights according to whether the values stored in val (monthly station objects) are greater than mean + sd or less than mean - sd.

Usage

composite.field(x,y,lsig.mask=TRUE,sig.lev=0.05,s=0.42,mon=NULL,
                      lty=1,col="black",lwd=1,main=NULL,sub=NULL)
compositeField(x,y,lsig.mask=TRUE,sig.lev=0.05,s=0.42,mon=NULL,
                      lty=1,col="black",lwd=1,main=NULL,sub=NULL)

Arguments

x A field object.
y A station object, a vector consisting of [-1,0,+1] or a vector consisting of the years (negative values are used for negative phase, eg. c(1991,1993,1998,-1961,-1963,-1994).
lsig.mask FLAG: mask the regions not statistically significant.
sig.lev Level of significance.
s Threshold for defining whether y is high or low = s * sd(y).
mon Month to analyse.
lty Contour line type.
col Contour line colour.
lwd Contour line width.
main Main title.
sub subtitle.

Value

A map object.

Author(s)

R.E. Benestad

Examples

## Not run: 
slp <- retrieve.nc("ncep_slp.nc")
data(oslo.t2m)
composite.field(slp,oslo.t2m)
## End(Not run)

[Package clim.pact version 2.2-15 Index]