System info: Linux slackware R Version 1.5.1 (2002-06-17) ESS 5.1.21 Emacs 21.2.1 ------------------- Colleagues Thanks very much for advice. I was not aware of the Fields package (although I've been using Geostat and Spatial). The image.plot() function in Fields does a lovely job, with an added bonus of fixing the problem I was having of balancing the number of colour levels with the xlim parameters. image.plot(time.h, z, -log10(t5), xlab='hours', ylab='depth m', col=topo.colors(64), nlevel=32, main = 'Chatham Rise DVM 12 kHz', yaxt='n', zlim=c(1,8), cex=0.5) R and the community of users and developers certainly provide a splendid tool. Sam -- Sam McClatchie, Research scientist (fisheries acoustics)))))))))) NIWA (National Institute of Water & Atmospheric Research Ltd) PO Box 14 901, Kilbirnie, Wellington, New Zealand s.mcclatchie at niwa.cri.nz Research home page <http://www.smcc.150m.com/> /\ >><xX(&> /// \\\ //// \\\\ /// <%)Xx><< ///// \\\\\\ ><(((@> ><(((%> ..>><xX(?>O<?)Xx><< -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._