## Read in temp profiles and plot.
hdr <- scan(file="/home/hcp/wrk/invth_pract/tprofs.dat",nlines=1)
nprofs<-hdr[1]
nlevs<-hdr[2]
zspacing<-hdr[3]

tmpdat <- scan(file="/home/hcp/wrk/invth_pract/tprofs.dat",nlines=nlevs,skip=1)
dim(tmpdat) <- c(nprofs+1,nlevs)
zeta <- tmpdat[1,]
tprofs <- tmpdat[2:(nprofs+1),]