require(grDevices)

J=matrix(c(200,200,175,150,100,100, 200,175,150,100,100,100, 150,125,100,100,100,100, 175,175,125,100,100,100, 175,200,175,125,125,100, 150,150,150,100,100,100), nrow = 6, ncol=6)
filled.contour(J, color.palette =colorRampPalette(c("blue", "gray", "white")),asp = 1)

x <- 10*1:nrow(J)
y <- 10*1:ncol(J)
filled.contour(x, y, J, color.palette =colorRampPalette(c("blue", "gray", "white")),asp = 1),
plot.title = title(main = "Cloud Contours",
                   xlab = "Meters North", ylab = "Meters West"),
plot.axes = { axis(1, seq(100, 800, by = 100))
              axis(2, seq(100, 600, by = 100)) },
key.title = title(main="Height\n(meters)"),
key.axes = axis(4, seq(90, 190, by = 10)))
mtext(paste("filled.contour(.) from", R.version.string), side = 1, line = 4, adj = 1, cex = .66)