rm(list=ls())
require(lattice)
# quartz(width=12, height=7.6)
# op <- par(mfrow = c(2, 3), pty = "m")       # square plotting region,
# layout(matrix(1:6, 2,3))
# skipi<-c(20, 159,298,439,578,717)   # 043
# skipi<-c(22, 161,300,451,590,729)   # 044
# skipi<-c(21, 160,299,439,578,717)   # 045
# skipi<-c(20,159,298, 439,578,717)   # 049 = 043 skips
# skipi<-c(23,217,411, 440,579,718)   # 051 skips
  skipi<-c(23,217,411, 606,800,994)   # 053 skips
URL<-"http://www.bio.umass.edu/biology/kunkel/LabWiki/images/7/7c/Chris055.txt"
# LIX<-'Chloride X-uV'
LIX<-'H mV'
Nam<-'Ch055'
filnam<-strsplit(URL,'/')
lfil<-length(filnam[[1]])
filnam<-filnam[[1]][lfil]
for(i in 1:6) { skiplin<-skipi[i]
XX<-read.csv(URL, skip=skiplin, nrows=176)
attach(XX)
ny<-16; nz<-11
quartz(width=6*ny/nz, height=6)
# print(names(XX))
mVav<-(NagelSamp0_1_mV + NagelSamp0_2_mV + NagelSamp0_3_mV)/3    # Cl
# mVav<-(NagelRead1_1_mV + NagelRead1_2_mV + NagelRead1_3_mV)/3    # H
opar=par(mar=c(5,4,4,2)+0.1)
Y<- Y_microns +75
Z<- Z_microns +50
W<- Z-Y/20
IW<-sort(W, index.return=TRUE)
IW<-IW$ix
ClmV<-matrix(mVav[IW],16,11)
# contourplot(ClXuV)
# The XuV goes from -500 to 10,000 in the 3 graphs.
cat('min/max ClmV =', min(ClmV), max(ClmV), '\n')
filled.contour(ClmV, 
              main=paste(Nam, LIX, 'skip', skiplin, 'Hrs from =', Hrs_From[1]),
              zlim=c(120,181), 
              color.palette = rainbow) 
  #           zlim=c(120,180), 
detach(XX)
              }