plot: Plot esd objects

In metno/esd.test: Climate analysis and empirical-statistical downscaling (ESD) package for monthly and daily data.

Description

These plot functions are S3 methods for esd objects.

Usage

nam2expr(x)
plot.station(x,plot.type="single",new=TRUE,
             lwd=3,type='l',pch=0,main=NULL,col=NULL,
             xlim=NULL,ylim=NULL,xlab="",ylab=NULL,...)
plot.eof(x,new=TRUE,xlim=NULL,ylim=NULL,pattern=1,what=c("pc","eof","var"),...)
plot.eof.field(x,new=TRUE,xlim=NULL,ylim=NULL,pattern=1,what=c("pc","eof","var"),...)
plot.eof.comb(x,new=TRUE,xlim=NULL,ylim=NULL,
                          pattern=1,col=c("red"),what=c("pc","eof","var"),...)
plot.eof.var(x,new=TRUE,xlim=NULL,ylim=NULL,pattern=20,...)
plot.ds(x,plot.type="single",what=c("map","ts"),new=TRUE,
        lwd=3,type='b',pch=0,main=NULL,col=NULL,
        xlim=NULL,ylim=NULL,xlab="",ylab=NULL,...)
plot.CCA(x,icca=1) 
plot.field(x,is=NULL,it=NULL,FUN="mean",...)
plot.spell(x)
plot.diagnose(x,...) 
plot.diagnose.comb.eof(x,...)
plot.xval(x,...)
plot.dsensemble(x,pts=FALSE,showci=TRUE,showtrend=TRUE,legend=TRUE,it=0,
                envcol=c(1,0,0,0.2),...)

Arguments

x	the object to be plotted
pattern	Which EOF pattern (mode) to plot
col	Colour
icca	Which CCA pattern to plot
is	For subsetting in space - See link{subset}, but can also be a station value and if provided, the plotting will involve an interploation to the same coordinates as defined by is.
it	For subsetting in time - See link{subset}.'it=0' returns the annual means (mean of DJF + MAM + JJA + SON)
FUN	function

Value

A field object

Author(s)

R.E. Benestad

Examples

# Example: use aggregate to compute annual mean temperature for Svalbard:
data(Svalbard)
year <- as.numeric( format(index(Svalbard), '%Y') ) 
y <- aggregate(Svalbard, by=year, mean, na.rm = FALSE) 
plot(y)

# Example with downscaling:
lon <- c(-12,37)
lat <- c(52,72)
t2m <- t2m.ERA40(lon=lon,lat=lat)
data(Oslo)
ds <- DS(Oslo,t2m)

# Plot the results for January month
plot(ds$Jan)

# Plot the residuals:
residual <- as.residual(ds)
obs <- as.anomaly(as.calibrationdata(ds))

plot.zoo(obs,lwd=2)
lines(residual,col="red")

print("Global climate model simulation NorESM")
T2m <- t2m.NorESM.M(lon=lon,lat=lat)

# Plot the global mean of the field:
plot(T2m)
# Plot area mean of a sub region
plot(T2m,is=list(lon=c(0,10),lat=c(60,70)))

# Plot interpolated results corresponding to ferder
data(ferder)
plot(T2m,ferder)

# Plot Hovmuller diagram:
plot(T2m,is=list(lon=0))

print("Extract a subset - the January month")
x <- subset(t2m,it=1)
X <- subset(T2m,it=1)

print("Combine the fields for computing common EOFs:")
XX <- combine(x,X)

print("Compute common EOFs")
eofxx <- EOF(XX)
plot(eofxx)

print("Downscale the January mean temperature") 
ds.jan <- DS(Oslo,eofxx)
plot(ds.jan)

metno/esd.test documentation built on May 22, 2019, 7:49 p.m.