#Data Electron Microprobe from Nov3, 2005 (1103)
quartz(wi=10.8, he=7.9)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
filnam<-"BL9LL_2"
xy<-read.table(paste(filnam,"_bse.size",sep=''), skip=2)
attach(xy)
xi<-V1; yi<-V2
detach(xy)
xy<-read.table(paste(filnam,"_bse.pgm",sep=''), skip=4)
attach(xy)
XX<-V1;
detach(xy)
XX<-matrix(XX,xi,yi)
XX[1025:1080,1:790]<-XX[1025:1080,790:1]
XX[1:1024,790:770]<-XX[1:1024,770:790]
image(XX, col=gray(1:255/255),xlab='', ylab='')

quartz(wi=7, he=7)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
xy<-read.table(paste(filnam,"_2_Ca.txt",sep=''), skip=33)
attach(xy)
YY<-t(as.matrix(xy))
# print(names(xy))
image(YY, col=gray(1:255/255),xlab='', ylab='')
quartz(wi=2, he=2)
# lx = left x; ly = lower y
lx<-15; ux<-165; ly<-130; uy<-280
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
image(YY[lx:ux,ly:uy], col=gray(1:255/255),xlab='', ylab='')
yy<-YY[lx:ux,ly:uy]
detach(xy)



quartz(wi=7, he=7)
npar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
xy<-read.table(paste(filnam,"_3_P.txt",sep=''), skip=33)
attach(xy)
ZZ<-t(as.matrix(xy))
# print(names(xy))
image(ZZ, col=gray(1:255/255),xlab='', ylab='')
quartz(wi=2, he=2)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
image(ZZ[lx:ux,ly:uy], col=gray(1:255/255),xlab='', ylab='')
zz<-ZZ[lx:ux,ly:uy]
detach(xy)

quartz(wi=9, he=9)
npar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
plot(zz[1:3000], yy[1:3000], main=paste(filnam,"_2_Ca.txt",sep=''))

# Narrow down Na & K to surface examples
cat("Na and K not interesting in this example! \n")

par(opar)