plotGrid = function(fillcol="#D0D0D0",gridcol="white",
                    glwd=2,outline="black"){
   xyBnds = par()$usr
   rect(xyBnds[1],xyBnds[3],xyBnds[2],xyBnds[4],
        col=fillcol,density=-1,border=NULL,xpd=T)  # Fill
   axis(side=1,tck=1,labels=F,mgp=c(2,.1,0),col=gridcol,lwd=glwd)
   axis(side=1,tck=0,mgp=c(2,.2,0))
   axis(side=2,tck=1,labels=F,mgp=c(2,.1,0),col=gridcol,lwd=glwd)
   axis(side=2,tck=0,mgp=c(2,.3,0),las=2)
   rect(xyBnds[1],xyBnds[3], xyBnds[2],xyBnds[4],
        col=outline,density=0,border=T,xpd=T) #Outline
}

plotOutline = function(col='black',lwd=1){
   xyBnds = par()$usr
rect(xyBnds[1],xyBnds[3], xyBnds[2],xyBnds[4],
        col=col,density=0,border=T,xpd=T,lwd=lwd) #Outline
}