filnam<-"JK2.xyz"
XX<- read.table(filnam)
print(names(XX))
Xloc<- XX[,2]
Yloc<- XX[,3]
print(summary(XX))
Ilow<-c(1, 102,202,302,402,502)
Ihi<-c(101,201,301,401,501,601)
filcol<-c("red", "green", "blue", "orange", "purple", "grey", "darkgreen")
fillwd<-c(1    ,  2     ,  2    ,  3      ,  2     ,  3      ,  2    )
fillty<-c("solid","solid","solid","solid","solid","solid", "solid")
filnlev<-c(8    , 8     , 8    ,  6      ,  6     ,  6      ,  10    )
plot(Xloc,Yloc, typ='n', main=paste(filnam,'has',6,'groups'))
for (i in 1:6){
     points(Xloc[Ilow[i]:Ihi[i]],Yloc[Ilow[i]:Ihi[i]], col=filcol[i])
	           }
 minX<- min(Xloc); maxX<- max(Xloc); minY<- min(Yloc); maxY<- max(Yloc)
 minX2<- floor(minX/200); maxX2<- floor(maxX/200)
 minY2<- floor(minY/200); maxY2<- floor(maxY/200)
 
 Xser<-200*(minX2:maxX2); Yser<-200*(minY2:maxY2)
 n<- length(Xser)
 for (i in 1:n) {lines(c(Xser[i],Xser[i]), c(minY, maxY), lty='dashed')}
 n<- length(Yser)
 for (i in 1:n) {lines(c(minX, maxX), c(Yser[i],Yser[i]),  lty='dashed')}
# Try now to extract data
qw=14; qh=6
quartz(width = qw, height = qh)
filnam<- "JKsample2compound.asc"
YY<- read.table(filnam, skip=9)
Na2CO3<- YY[,2]
MgCO3<- YY[,3]
SrCO3<- YY[,4]
P2CO35<- YY[,5]
K2CO3<- YY[,6]
Cl<- YY[,7]
CaCO3<- YY[,8]
BaCO3<- YY[,9]
FeCO3<- YY[,10]
MnCO3<- YY[,11]
CCO32<- YY[,12]
n<-length(CaCO3)
plot(1:n, CaCO3, typ='l', main=paste(filnam,'CaCO3 and P2CO35'))
lines(1:n,15*P2CO35, col='red')
lines(1:n,15*Cl, col='blue')
lines(1:n,30*Na2CO3, col='green')
lines(1:n,30*K2CO3, col='purple')
mtext("CaCO3")
mtext("15*P2(CO3)5", col='red', line=-1)
mtext("15*Cl", col='blue', line=-2)
mtext("30*Na2CO3", col='green', line=-3)
mtext("30*K2CO3", col='purple', line=-4)
quartz()
plot(P2CO35, CaCO3, main=paste(filnam,'CaCO3 vs P2CO35'))

quartz(width = qw, height = qh)
filnam<- "JKsample2wt.asc"
YY<- read.table(filnam, skip=9)
M_Na<- (Na<- YY[,2])/22.98976928
M_Mg<- (Mg<- YY[,3])/24.3050
M_Sr<- (Sr<- YY[,4])/87.62
M_P<- (P<- YY[,5])/30.973762
M_K<- (K<- YY[,6])/39.0983
M_Cl<- (Cl<- YY[,7])/35.453
M_Ca<- (Ca<- YY[,8])/40.078
M_Ba<- (Ba<- YY[,9])/137.327
M_Fe<- (Fe<- YY[,10])/55.845
M_Mn<- (Mn<- YY[,11])/54.938045
CC<- YY[,12]

n<-length(Ca)
plot(1:n, M_Ca, typ='l', main=paste(filnam,'Ca P Cl Na K Molar'))
lines(1:n,10*M_P, col='red')
lines(1:n,15*M_Cl, col='blue')
lines(1:n,20*M_Na, col='green')
lines(1:n,40*M_K, col='purple')
mtext("M_Ca")
mtext("15*M_P", col='red', line=-1)
mtext("15*M_Cl", col='blue', line=-2)
mtext("30*M_Na", col='green', line=-3)
mtext("30*M_K", col='purple', line=-4)
quartz(width = qw, height = qh)
plot(1:n, M_Ca, typ='l', main=paste(filnam,'Ca Mg Sr Ba Mn Molar'))
lines(1:n,10*M_Mg, col='red')
lines(1:n,50*M_Sr, col='blue')
lines(1:n,100*M_Ba, col='green')
lines(1:n,100*M_Mn, col='purple')
mtext("M_Ca")
mtext("10*M_Mg", col='red', line=-1)
mtext("50*M_Sr", col='blue', line=-2)
mtext("100*M_Ba", col='green', line=-3)
mtext("100*M_Mn", col='purple', line=-4)

quartz()
plot(P, Ca, main=paste(filnam,'Ca vs P'))
quartz(width = qw, height = qh)
plot(1:n, Ca, typ='l', main=paste(filnam,'CaCO3 and P2CO35'))
lines(1:n,15*P, col='red')
lines(1:n,15*Cl, col='blue')
lines(1:n,30*Na, col='green')
lines(1:n,30*K, col='purple')
mtext("Ca")
mtext("15*P", col='red', line=-1)
mtext("15*Cl", col='blue', line=-2)
mtext("30*Na", col='green', line=-3)
mtext("30*K", col='purple', line=-4)

quartz(width = qw, height = qh)
filnam<- "JKvsampTrav3compound.asc"
YY<- read.table(filnam, skip=9)
Na2CO3<- YY[,2]
MgCO3<- YY[,3]
SrCO3<- YY[,4]
P2CO35<- YY[,5]
K2CO3<- YY[,6]
Cl<- YY[,7]
CaCO3<- YY[,8]
BaCO3<- YY[,9]
FeCO3<- YY[,10]
MnCO3<- YY[,11]
CCO32<- YY[,12]
n<-length(CaCO3)
plot(1:n, CaCO3, typ='l', main=paste(filnam,'CaCO3 and P2CO35'))
lines(1:n,15*P2CO35, col='red')
lines(1:n,15*Cl, col='blue')
lines(1:n,30*Na2CO3, col='green')
lines(1:n,30*K2CO3, col='purple')
mtext("CaCO3")
mtext("15*P2(CO3)5", col='red', line=-1)
mtext("15*Cl", col='blue', line=-2)
mtext("30*Na2CO3", col='green', line=-3)
mtext("30*K2CO3", col='purple', line=-4)
quartz()
plot(P2CO35, CaCO3, main=paste(filnam,'CaCO3 vs P2CO35'))