layout(matrix(1:4, 1, 1, byrow = TRUE))

filnam="L2Lgrid.TXT"
zz<-read.table(filnam, header=TRUE, skip=23, nrows=168, sep=',')
print(names(zz))
attach(zz)

LimZ<- c(0,1200)
Xmax<-max(c(Xo.1.mV,Xo.2.mV,Xo.3.mV))
Xmin<-min(c(Xo.1.mV,Xo.2.mV,Xo.3.mV))
LimXm<-c(Xmin,Xmax)
Xmax<-max(c(Xo.X.1.uV,Xo.X.2.uV,Xo.X.3.uV))
Xmin<-min(c(Xo.X.1.uV,Xo.X.2.uV,Xo.X.3.uV))
LimXu<-c(Xmin,Xmax)
Zmax<-max(c(Zo.1.mV,Zo.2.mV,Zo.3.mV))
Zmin<-min(c(Zo.1.mV,Zo.2.mV,Zo.3.mV))
LimZm<-c(Zmin,Zmax)
Zmax<-max(c(Zo.Z.1.uV,Zo.Z.2.uV,Zo.Z.3.uV))
Zmin<-min(c(Zo.Z.1.uV,Zo.Z.2.uV,Zo.Z.3.uV))
LimZu<-c(Zmin,Zmax)

ZmVmed<-(Zo.1.mV+Zo.2.mV+Zo.3.mV)/3
ZuVmed<-(Zo.Z.1.uV+Zo.Z.2.uV+Zo.Z.3.uV)/3
plot(LimZ, LimZm, typ='n', main="mV vs height", cex.axis=1.5, cex.lab=1.5)
lines(ProbeZ, Zo.1.mV, typ='p', lty=1)
lines(ProbeZ, Zo.2.mV, typ='p', lty=2)
lines(ProbeZ, Zo.3.mV, typ='p', lty=3)
for (i in seq(1,168,4)) {

lines(ProbeZ[i+0:3], ZmVmed[i+0:3], typ='l', lty=3, lwd=2, col='red')
                         }
                         
layout(matrix(1:4, 1, 2, byrow = TRUE))

zzsort<- sort(ProbeZ, index.return=TRUE)
zzsi<-zzsort$ix
Zo.Av.mV<- (Zo.1.mV+Zo.2.mV+Zo.3.mV)/3
plot(HoursFrom[zzsi], Zo.Av.mV[zzsi], typ='n', main=paste(filnam, "Zo mV vs Time"), cex.axis=1.5, cex.lab=1.5)
lines(HoursFrom[zzsi[1:42]], Zo.Av.mV[zzsi[1:42]], typ='b', col='red', lwd=2)
lines(HoursFrom[zzsi[43:84]], Zo.Av.mV[zzsi[43:84]], typ='b', col='green', lwd=2)
lines(HoursFrom[zzsi[85:126]], Zo.Av.mV[zzsi[85:126]], typ='b', col='blue', lwd=2)
lines(HoursFrom[zzsi[127:168]], Zo.Av.mV[zzsi[127:168]], typ='b', col='purple', lwd=2)

Xo.Av.mV<- (Xo.1.mV+Xo.2.mV+Xo.3.mV)/3
plot(HoursFrom[zzsi], Xo.Av.mV[zzsi], typ='n', main=paste(filnam, "Xo mV vs Time"), cex.axis=1.5, cex.lab=1.5)
lines(HoursFrom[zzsi[1:42]], Xo.Av.mV[zzsi[1:42]], typ='b', col='red', lwd=2)
lines(HoursFrom[zzsi[43:84]], Xo.Av.mV[zzsi[43:84]], typ='b', col='green', lwd=2)
lines(HoursFrom[zzsi[85:126]], Xo.Av.mV[zzsi[85:126]], typ='b', col='blue', lwd=2)
lines(HoursFrom[zzsi[127:168]], Xo.Av.mV[zzsi[127:168]], typ='b', col='purple', lwd=2)

detach(zz)