#GeOxy60a.R script

quartz(width=11, height=9)
layout(matrix(1:1, 1, 1, byrow = TRUE))

XX<-read.csv("jgk080825lobO_060.csv", skip=1645, nrows=-1)
attach(XX)
OxuV<-AvgGradO2uV
Start<- Hrs_From[1]
Minutes<- round((Hrs_From -Start)*60,1)

# plot(c(min(Minutes),max(Minutes)),c(min(OxuV),max(OxuV)),typ='n',
minMin<-min(Minutes); maxMin<- max(Minutes)
# plot(c(minMin, maxMin),c(min(OxuV),max(OxuV)),typ='n',
plot(c(minMin, maxMin),c(-10000,3000),typ='n',
     ylab="uV",
     xlab="Time (Minutes)",
     lwd=2,
     col='red',
     cex.axis=2,
     cex.lab=1.5,
     cex.main=1.5,
     main="pH 8.0: lesion studies: Xpos(green), Ypos(blue), Zpos(purple)")

lines(Minutes, OxuV, col='red', lwd=2, typ='p')

par(new=TRUE, ann=FALSE)
minX<- min(X_microns);  maxX<-max(X_microns)
plot(c(minMin, maxMin),c(minX, maxX),typ='n', yaxt='n', xaxt='n')
lines(Minutes,X_microns, typ='p', col='green', cex=2, lwd=3, lty=2)

par(new=TRUE, ann=FALSE)
minY<- min(Y_microns);  maxY<-max(Y_microns)
plot(c(minMin, maxMin),c(minY, maxY),typ='n', yaxt='n', xaxt='n')
lines(Minutes,Y_microns, typ='p', col='blue', cex=2, lwd=3, lty=2)

par(new=TRUE, ann=FALSE)
minZ<- min(Z_microns);  maxZ<-max(Z_microns)
plot(c(minMin, maxMin), c(minZ, maxZ), typ='n', yaxt='n', xaxt='n')
lines(Minutes,Z_microns, typ='p', col='purple', cex=2, lwd=3, lty=2)

detach(XX)

mtext(" Biddeford, VIII-25-2008", line=-2, cex=2, adj=0)
mtext("Electrode 2", line=-2, cex=2, adj=1)