Example 001

setwd("/Users/rockel/work/R/examples/example001")

#
# vector plot
#
levs<-seq(-23.,23.,2.)
pdffile<-"exp001.pdf"
title<-"Vectorplot"
col<-rainbow(length(levs)-1)
varname<-"HSURF"
ncfile<-"cclm_const.nc"
ncufile<-"U_10M.nc"
ncvfile<-"V_10M.nc"
nctfile<-"T_2M.nc"

ncu <- nc_open(ncufile)
ncv <- nc_open(ncvfile)
nct <- nc_open(nctfile)
u<-ncvar_get(ncu,"U_10M")
legscale<-10
units <- paste(legscale," ",ncatt_get(ncu, "U_10M", "units")$value,sep="")
v<-ncvar_get(ncv,"V_10M")
t2m<-ncvar_get(nct,"T_2M")
rlon<-ncvar_get(ncu,"rlon")
rlat<-ncvar_get(ncu,"rlat")
nc_close(ncu)
nc_close(ncv)
nc_close(nct)

t2m1<-t2m-273.15

# decide whether landscape or portrait looks prettier
pdf(pdffile, width=8.27,height=11.69,paper="special")
par(oma=c(3,2,5,2))
layout(matrix(c(1,2,3),nrow=3,ncol=1), heights=c(1,lcm(2.5),lcm(2.5)),respect=T)

par(mar=c(0.5,0.5,0.5,0.5))
ptm <- proc.time()
tmp <- plotmap(file=ncfile, varname=varname, rivers=F, lakes=F, cities=F,sponge=0,grid.txt=F, hires=F,map.lwd=1,interior=T)
windvec(u, v, rlon, rlat,units=units,colvar=t2m1,col=col,levels=levs,legend.scale=legscale)
tmpvec<-list(col=col, lev=levs, sea.col=NULL, flag_values=NULL,
        flag_meanings=NULL, longname="2m temperature", units="C",gridmapping=tmp$gridmapping)
class(tmpvec)  <- "plotmap"
par(mai=c(0.5,0.5,0.35,0.5))
plot_colourbar(tmp)
plot_colourbar(tmpvec,triangle_ends=c(TRUE,TRUE))
mtext(title, 3, line=1,outer=T, font=1, cex=1.4)
dev.off()