Color scale bar function

color.scaleBAR <- function(colors, min, max=-min, nticks=11, ticks=seq(min, max, len=nticks), title='',digits=2){

  scale = (length(colors)-1)/(max-min)

  plot(c(0,10), c(min,max), type='n', bty='n', xaxt='n', xlab='', yaxt='n', ylab='', main=title)

  axis(2, round(ticks,digits), las=1)

  for (i in 1:(length(colors)-1)) {

    y = (i-1)/scale + min

    rect(0,y,10,y+1/scale, col=colors[i], border=NA)

  }

}