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, 607,578,717)   # 052 skips
URL<-"http://www.bio.umass.edu/biology/kunkel/LabWiki/images/5/51/Chris052.txt"
# LIX<-'Chloride X-uV'
LIX<-'Cl xuV'
Nam<-'Ch052'
filnam<-strsplit(URL,'/')
lfil<-length(filnam[[1]])
filnam<-filnam[[1]][lfil]
for(i in 1:4) { 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<-(NagelSamp1_1_mV + NagelSamp1_2_mV + NagelSamp1_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
ClxuV<-matrix(Avg_xgradCl_uV[IW],ny,nz)
# contourplot(ClXuV)
# The XuV goes from -500 to 10,000 in the 3 graphs.
cat('min/max ClxuV =', min(ClxuV), max(ClxuV), '\n')
filled.contour(ClxuV, 
              main=paste(Nam, LIX, 'skip', skiplin, 'Hrs from =', Hrs_From[1]),
              zlim=c(-3540,9850), 
              color.palette = rainbow) 
  #           zlim=c(95,190), 
detach(XX)
              }