quartz(width=18, height=6)
Hrs_Start<- 12.73
YY<-readLines("jgk080731apatite_038.csv",  n=22)
ilog<-grep("LogEntry",YY)
YY<-read.csv("jgk080731apatite_038.csv", header=FALSE,  nrows=22)
YY<- YY[ilog,c(1,7)]
YY[,1]<- YY[,1]-Hrs_Start
print(YY)

XX<-read.csv("jgk080731apatite_038.csv", skip=59, nrows=-1)
attach(XX)

Tmin<- 0; Tmax<- 5
uVmin<- -400; uVmax<- 100
Zmin<- 0    ; Zmax<- 2000 
Z_add<- c(rep(3150,592),rep(2500,292),rep(500,169),rep(850,404),rep(1150,492))
Zamin<- min(Z_add); Zamax<- max(Z_add)
Hmvmin<- min(NagelSamp0H_1_mV); Hmvmax<- max(NagelSamp0H_1_mV)
CamVmin<- min(NagelRead0Ca_1_mV); CamVmax<- max(NagelRead0Ca_1_mV)

print(names(XX))
plot(c(Tmin,Tmax), c(uVmin, uVmax), typ='n', main="Apatite crystal uV measurements for [H+] (red) and [Ca++] (green) in the Z-axis (blue), HmV (black), CamV (purple) vs Time",
       ylab="H uV", xlab="time (hrs)")
lines(Hrs_From - Hrs_Start, AvggradHuV, typ='p', col='red')
lines(Hrs_From - Hrs_Start, AvggradCauV, typ='p', col='green')
par(new=TRUE, ann=FALSE)
plot(c(Tmin,Tmax), c(Zmin, Zmax), typ='n', main="Apatite crystal uV measurements for [H+] and [Ca++] vs Time", xaxt='n', yaxt='n')
lines(Hrs_From - Hrs_Start, Z_microns+Z_add, typ='p', col='blue')
axis(side=4, cex.axis=1.5)
par(new=TRUE, ann=FALSE)
plot(c(Tmin,Tmax), c(Zamin, Zamax), typ='n', main="Apatite crystal uV measurements for [H+] and [Ca++] vs Time", xaxt='n', yaxt='n')
lines(Hrs_From - Hrs_Start, Z_add, typ='l', col='black')
par(new=TRUE, ann=FALSE)
plot(c(Tmin,Tmax), c(Hmvmin, Hmvmax), typ='n', main="Apatite crystal uV measurements for [H+] and [Ca++] vs Time", xaxt='n', yaxt='n')
lines(Hrs_From - Hrs_Start, NagelSamp0H_1_mV, typ='p', col='black', pch=18)
par(new=TRUE, ann=FALSE)
plot(c(Tmin,Tmax), c(CamVmin, CamVmax), typ='n', main="Apatite crystal uV measurements for [H+] and [Ca++] vs Time", xaxt='n', yaxt='n')
lines(Hrs_From - Hrs_Start, NagelRead0Ca_1_mV, typ='p', col='purple', pch=18)

browser()
detach(XX)