R/sl.plot.field.R

In helgegoessling/spheRlab: Spherical Geometry, Analysis, and Plotting of Geoscientific Data on Arbitrary Grids

sl.plot.field <-
function (plot.init.res,num,lon.v,lat.v,fill=TRUE,col.fill="colbar",border=TRUE,col.border="colbar",colbar=sl.colbar.redgreyblue_256,colbar.breaks=NA,colbar.breaks.log=FALSE,border.lwd=1,border.lty=1,na.col=NA) {
	
	Npoly = nrow(lon.v)
	if (col.fill == "colbar" || col.border == "colbar") {
		colbar.res = sl.num2colbar(num,colbar,colbar.breaks,colbar.breaks.log)
		col.ind = colbar.res$colour.index
		na.col.ind = length(colbar) + 1
		colbar[[na.col.ind]] = na.col
		col.ind[is.na(col.ind)] = na.col.ind
	}
	for (np in 1:Npoly) {
		cb.fill = col.fill
		cb.border = col.border
		if (col.fill == "colbar") {cb.fill = colbar[[col.ind[np]]]}
		if (col.border == "colbar") {cb.border = colbar[[col.ind[np]]]}
		sl.plot.polygon(plot.init.res,lon.v[np,],lat.v[np,],fill=fill,col.fill=cb.fill,border=border,col.border=cb.border,border.lwd=border.lwd,border.lty=border.lty)
	}
	
	if (exists("colbar.res")) {
		return(colbar.res)
	} else {
		return(NULL)
	}
	
}