#Lobster UNE 2.2 newly lesioned with ball mill bit.
rm(list=ls())
filnam<- "jgk080813lob_045.csv"
kline<-290
id<-1.1
XX<-read.csv(filnam, skip=432, nrows=880)      # was 79
attach(XX)
print(names(XX))
ni<-length(X_microns)
minHr<-min(Hrs_From)
maxHr<-max(Hrs_From)
minXloc<-min(X_microns)
maxXloc<-max(X_microns)
minYloc<-min(Y_microns)
maxYloc<-max(Y_microns)
minZloc<-min(Z_microns)
maxZloc<-max(Z_microns)
leng<- length(X_microns)
medCauV<-medHuV<- rep(0,leng)
for (i in 1:leng) {medHuV[i]<-  gradH_1_uV[i]     # median(c(gradH_3_uV[i],gradH_4_uV[i],gradH_5_uV[i]))
	                medCauV[i]<- gradCa_1_uV[i]   # median(c(gradCa_3_uV[i],gradCa_4_uV[i],gradCa_5_uV[i]))}
                   }
quartz(width=19.5, height=8)
layout(matrix(1:2, 1, 2, byrow = TRUE))
par(mar=c(5, 5, 4, 4.5) + 0.1)

buf0<- 13.00
minY<- 35;  maxY<- 85
HmV<- (NagelSamp0H_1_mV+NagelSampZH_1_mV)/2
CamV<- (NagelRead0Ca_1_mV+NagelReadZCa_1_mV)/2
plot(60*(c(minHr, maxHr)-buf0),c(minY, maxY),typ='n',
main=paste(filnam,"mV vs Time at pH 7.8 Life Table Sea Water"), ylab='[H] mV', xlab='minutes after buffer change', cex.lab=2, cex.axis=1.5)
mtext("Lesion 3.1 of UNE2.2")
Z1<- -799.6; Z2<- -749.7; Z3<- -699.7; Z4<- -649.7; Z5<- -599.7; 
lines(60*(Hrs_From-buf0),HmV, typ='b', col='red', cex=2, lwd=4, lty=2)
par(new=TRUE, ann=FALSE)
minY<- 4;  maxY<- 32
plot(60*(c(minHr, maxHr)-buf0),c(minY, maxY),typ='n', yaxt='n', xaxt='n')
lines(60*(Hrs_From-buf0),CamV, typ='b', col='black', cex=2, lwd=4, lty=2)
axis(side=4, cex.axis=1.5)
mtext("[Ca] mV", side=4, line=2.5, cex=1.25)
par(new=TRUE, ann=FALSE)
minY<- 4;  maxY<- 32
plot(60*(c(minHr, maxHr)-buf0),c(minZloc, maxZloc),typ='n', yaxt='n', xaxt='n')
lines(60*(Hrs_From-buf0),Z_microns, typ='l', col='green', cex=2, lwd=4, lty=1)
axis(side=4, cex.axis=1.5)
mtext("[Ca] mV", side=4, line=2.5, cex=1.25)
par(new=TRUE, ann=FALSE)
minY<- 4;  maxY<- 32
plot(60*(c(minHr, maxHr)-buf0),c(minXloc, maxXloc),typ='n', yaxt='n', xaxt='n')
lines(60*(Hrs_From-buf0),X_microns, typ='l', col='blue', cex=2, lwd=4, lty=1)
axis(side=4, cex.axis=1.5)
mtext("[Ca] mV", side=4, line=2.5, cex=1.25)

par(ann=TRUE)
browser()
detach(XX)

rm(list=ls())