rm(list=ls())
filnam<- "jgk080729_clam_032.csv"
kline<-290
id<-1.1
XX<-read.csv(filnam, skip=213, nrows=-1)  
attach(XX)
print(names(XX))
ni<-length(X_microns)
minHr<-min(Hrs_From)
maxHr<-max(Hrs_From)
minX<-min(X_microns)
maxX<-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]<- median(c(gradH_8_uV[i],gradH_9_uV[i],gradH_10_uV[i]))
	                medCauV[i]<- median(c(gradCa_8_uV[i],gradCa_9_uV[i],gradCa_10_uV[i]))}

quartz(width=19, height=7)
layout(matrix(1:3, 1, 3, byrow = TRUE))
par(mar=c(5, 5, 4, 4.5) + 0.1)


buf0<- 13.109770
minY<- 87;  maxY<- 97
HmV<- (NagelSamp0H_10_mV+NagelSampZH_10_mV)/2
CamV<- (NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2
plot(60*(c(minHr, maxHr)-buf0),c(minY, maxY),typ='n',
main=paste(filnam,"mV vs Time at pH 7.88"), ylab='[H] mV', xlab='minutes after buffer change', cex.lab=2, cex.axis=1.5)
lines(60*(Hrs_From[Z_microns == 0]-buf0),HmV[Z_microns == 0], typ='b', col='red', cex=2, lwd=4, lty=2)
lines(60*(Hrs_From[Z_microns == 100]-buf0),HmV[Z_microns == 100], typ='b', col='green', cex=2, lwd=3, lty=3)
lines(60*(Hrs_From[Z_microns == 200]-buf0),HmV[Z_microns == 200], typ='b', col='blue', cex=2, lwd=2, lty=4)
lines(60*(Hrs_From[Z_microns == 300]-buf0),HmV[Z_microns == 300], typ='b', col='orange', cex=2, lwd=2, lty=5)
lines(60*(Hrs_From[Z_microns == 400]-buf0),HmV[Z_microns == 400], typ='b', col='purple', cex=2, lwd=2, lty=6)
par(new=TRUE, ann=FALSE)
minY<- -5;  maxY<- 18
plot(60*(c(minHr, maxHr)-buf0),c(minY, maxY),typ='n', yaxt='n', xaxt='n')
lines(60*(Hrs_From[Z_microns == 0]-buf0),CamV[Z_microns == 0], typ='b', col='black', cex=2, lwd=4, lty=2)
lines(60*(Hrs_From[Z_microns == 100]-buf0),CamV[Z_microns == 100], typ='b', col='brown', cex=2, lwd=3, lty=3)
lines(60*(Hrs_From[Z_microns == 200]-buf0),CamV[Z_microns == 200], typ='b', col='grey', cex=2, lwd=2, lty=4)
lines(60*(Hrs_From[Z_microns == 300]-buf0),CamV[Z_microns == 300], typ='b', col='purple', cex=2, lwd=1, lty=5)
lines(60*(Hrs_From[Z_microns == 400]-buf0),CamV[Z_microns == 400], typ='b', col='purple', cex=2, lwd=1, lty=6)
axis(side=4, cex.axis=1.5)
mtext("[Ca] mV", side=4, line=2.5, cex=1.25)

par(ann=TRUE)
detach(XX)

rm(list=ls())