| plotField {clim.pact} | R Documentation |
Produce 2D plots like maps and Hovmuller diagrams. A
poor man's version of Ferret's
http://ferret.wrc.noaa.gov/Ferret/ plot
function. plot.field is a high level command that utilises
mapField or grd.box.ts whenever appropriate.
plotField(x,lon=NULL,lat=NULL,tim=NULL,mon=NULL,val.rng=NULL,
col="black",col.coast="grey",lty=1,lwd=1,what="ano",
type="l",pch=26,my.col=NULL,add=FALSE,
main=NULL,sub=NULL,xlab=NULL,ylab=NULL,
xlim=NULL,ylim=NULL)
x |
A field object. |
lon |
Position for longitude. One of lon, lat, tim must be set. |
lat |
Position for latitude. |
tim |
Position for time as time index (1:length(tim). |
mon |
Month to extract. |
val.rng |
Valid range: interval used for colour scale. |
col |
Contour line colour for levels. |
col.coast |
Contour line colour for coast lines. |
lty |
Contour line type. |
lwd |
Contour line width. |
what |
Choose between "ano" - anomaly; "abs" absolute; "cli" climatological. |
type |
same as in plot(), only works for single time series. |
pch |
same as in plot(), only works for single time series. |
my.col |
colour palette (see link{rgb}). |
add |
TRUE adds a time series to previous (see link{grd.box.ts}). |
main |
Preset main title (see plot). |
sub |
subtitle |
xlab |
x label |
ylab |
y label |
xlim |
see plot |
ylim |
see plot |
R.E. Benestad
## Not run:
skt <- retrieve.nc("skt.mon.mean.nc",x.rng=c(-90,50),y.rng=c(0,75))
# Maps of monthly mean skin temperatures:
plotField(skt,tim=1,val.rng=c(-20,20))
dev2bitmap("ncep.skt_194801.jpg",type="jpeg")
plotField(skt,tim=100,col="blue",col.coast="darkgreen",val.rng=c(-10,10))
# For adding extra points/contours:
# From filled.contour in base
mar.orig <- (par.orig <- par(c("mar","las","mfrow")))$mar
on.exit(par(par.orig))
w <- (3 + mar.orig[2]) * par('csi') * 2.54
layout(matrix(c(2, 1), nc=2), widths=c(1, lcm(w)))
par(las = 1)
mar <- mar.orig
mar[4] <- 1
par(mar=mar)
# End of section affecting the window set up.
points(0,50,pch=21,col="red")
grid()
dev2bitmap("ncep.skt_195604.jpg",type="jpeg")
# A hovmuller diagram:
plotField(skt,lon=0,val.rng=c(-10,10))
dev2bitmap("ncep.skt_lontim.jpg",type="jpeg")
# A single time series:
plotField(skt,lon=-20,lat=50)
## End(Not run)