byte2binary(.vsd[24])
.vsd
.vsd[9]
.vsd[22]
.vsd[21:22]
readBin(.vsd[21:22], endian="little")
readBin(.vsd[21:22], endian="little",what="integer",size=2)
readBin(.vsd[21:22], endian="little",what="integer",size=2) * 0.01
library()
?gc
gcinfo()
gcinfo(verbose=T)
gc()
gcinfo(verbose=TRUE)
rm(list=ls())
load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_enu.rda")
enu.snippet <- subset(m05.vector.enu,#
                      as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time#
                      & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))
plot(enu.snippet)
abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC"))
plot(enu.snippet,which="u1")
abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+3*3600)
plot(d,which="u1")
abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+3*3600,col="red")
abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+0*3600,col="blue")
save(d,file="d.rda")
rm(d)
load(d.rda)
load('d.rda')
x11();plot(d,which='u1');abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+0*3600,col="blue")
load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_beam.rda")
beam.snippet <- subset(m05.vector.beam,#
                      as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time#
                      & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))
beam.snippet <- subset(m05.vector.beam,#
                      as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time#
                      & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"),debug=4)
which(c(T,T,T,F,T)==F)
which(c(T,T,T,F,T)==F)[1]
x11();plot(beam.snippet,which='u1');abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+0*3600,col="blue")
attr(beam.snippet$data$ts$time[1],"tzone")
source('~/src/R-kelley/oce/R/oce.R');beam.snippet <- subset.oce(m05.vector.beam,#
                      as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time#
                      & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"),debug=4)
plot(beam.snippet$data$ts$time, beam.snippet$data$ma$v[,1])
plot(beam.snippet$data$ts$time, beam.snippet$data$ma$v[,1],type='l')
plot(beam.snippet$data$ts$time, beam.snippet$data$ma$v[,1],type='l');abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+3*3600)
plot(beam.snippet$data$ts$time, beam.snippet$data$ma$v[,1],type='l');abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+3*3600,col="red")
plot(beam.snippet$data$ts$time, beam.snippet$data$ma$v[,1],type='l');abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+3*3600,col="red",lwd=3)
?print()
pdf("test.pdf");plot(beam.snippet,which="u1",type='l');abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+3*3600,col="red",lwd=3);dev.off()
pdf("test.pdf");plot(beam.snippet,which="u1",type='l');abline(v=as.POSIXct("2008-06-29 01:31:24",tz="UTC")+3*3600,col="red",lwd=3);mtext("2008-06-29 04:31:24 UTC",side=3,at=as.POSIXct("2008-06-29 01:31:24",tz="UTC"));dev.off()
pdf("test.pdf");plot(beam.snippet,which="u1",type='l');abline(v=as.POSIXct("2008-06-29 04:31:24",tz="UTC"),col="red",lwd=3);mtext("2008-06-29 04:31:24 UTC",side=3,at=as.POSIXct("2008-06-29 04:31:24",tz="UTC"));dev.off()
pdf("test.pdf");plot(beam.snippet,which="u1",type='l');abline(v=as.POSIXct("2008-06-29 04:31:24",tz="UTC"),col="red",lwd=3);mtext("2008-06-29 04:31:24 UTC",side=3,at=as.POSIXct("2008-06-29 04:31:24",tz="UTC"));mtext("is my line as CR's at 01:31:24 UTC?",side=3,line=-3)dev.off()
pdf("test.pdf");plot(beam.snippet,which="u1",type='l');abline(v=as.POSIXct("2008-06-29 04:31:24",tz="UTC"),col="red",lwd=3);mtext("2008-06-29 04:31:24 UTC",side=3,at=as.POSIXct("2008-06-29 04:31:24",tz="UTC"));mtext("is my line as CR's at 01:31:24 UTC?",side=3,line=-3);dev.off()
x11();source('~/src/R-kelley/oce/R/adv.R');source('~/src/R-kelley/oce/R/oce.R');plot.adv(beam.snippet,which="u1",type='l');abline(v=as.POSIXct("2008-06-29 04:31:24",tz="UTC"),col="red",lwd=3);mtext("2008-06-29 04:31:24 UTC",side=3,at=as.POSIXct("2008-06-29 04:31:24",tz="UTC"));mtext("is my line as CR's at 01:31:24 UTC?",side=3,line=-3)
x11();source('~/src/R-kelley/oce/R/adv.R');source('~/src/R-kelley/oce/R/oce.R');plot.adv(beam.snippet,which="u1",type='l',debug=4);abline(v=as.POSIXct("2008-06-29 04:31:24",tz="UTC"),col="red",lwd=3);mtext("2008-06-29 04:31:24 UTC",side=3,at=as.POSIXct("2008-06-29 04:31:24",tz="UTC"));mtext("is my line as CR's at 01:31:24 UTC?",side=3,line=-3)
x11()#
source('~/src/R-kelley/oce/R/adv.R')#
source('~/src/R-kelley/oce/R/oce.R')#
plot.adv(beam.snippet,which="u1",type='l',debug=4)#
t0 <- as.POSIXct("2008-06-29 04:31:24",tz="UTC")#
abline(v=t0,col="red",lwd=3)#
mtext(format(t0),side=3,at=t0)#
mtext("is my line as CR's at 01:31:24 UTC?",side=3,line=-3)
d <- read.adv(file="/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",from=as.POSIXct("2008-06-29 00:00:00", tz="UTC"),to=as.POSIXct("2008-06-29 06:00:00", tz="UTC"))
from <- as.POSIXct("2008-06-29 00:00:00", tz="UTC")#
to <- as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
flag <- as.POSIXct("2008-06-29 04:31:24",tz="UTC")
from <- as.POSIXct("2008-06-29 00:00:00", tz="UTC")#
to <- as.POSIXct("2008-06-29 06:00:00", tz="UTC")#
flag <- as.POSIXct("2008-06-29 04:31:24",tz="UTC")
library(oce)#
from <- as.POSIXct("2008-06-29 00:00:00", tz="UTC")#
to <- as.POSIXct("2008-06-29 06:00:00", tz="UTC")#
flag <- as.POSIXct("2008-06-29 04:31:24",tz="UTC")
d <- read.adv(file="/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",from=from,to=to)#
x11(title='focussed reading')#
plot(d, which='u1')#
abline(v=flag,col="red",lwd=3)#
mtext(format(flag),side=3,at=flag)
flag.cr <- as.POSIXct("2008-06-29 01:31:24",tz="UTC")
d <- read.adv(file="/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",from=from,to=to)#
x11(title='focussed reading')#
plot(d, which='u1')#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag)
flag.dk <- flag.cr + 3 * 3600
8/(8-7.698)
(8-7.698)/8
(8-7.698)/8 * 3600
d <- read.adv(file="/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",from=from,to=to,debug=4)#
x11(title='focussed reading')#
plot(d, which='u1')#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag)
0x30
byte2binary(0x30)
source('~/src/R-kelley/oce/R/adv.nortek.R')#
source('~/src/R-kelley/oce/R/adv.R')#
source('~/src/R-kelley/oce/R/oce.R');beam.snippet <- subset(m05.vector.beam, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet,which="u1",type='l'ug=4)#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk)
source('~/src/R-kelley/oce/R/adv.nortek.R')#
source('~/src/R-kelley/oce/R/adv.R')#
source('~/src/R-kelley/oce/R/oce.R');beam.snippet <- subset(m05.vector.beam, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet,which="u1",type='l')#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk)
source('~/src/R-kelley/oce/R/adv.nortek.R')#
source('~/src/R-kelley/oce/R/adv.R')#
source('~/src/R-kelley/oce/R/oce.R');beam.window <- window(m05.vector.beam, start=as.POSIXct("2008-06-29 00:00:00", tz="UTC"), end=as.POSIXct("2008-06-29 06:00:00", tz="UTC"),debug=4)#
plot.adv(beam.snippet,which="u1",type='l')#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk)
source('~/src/R-kelley/oce/R/adv.nortek.R')#
source('~/src/R-kelley/oce/R/adv.R')#
source('~/src/R-kelley/oce/R/oce.R');beam.snippet <- subset(m05.vector.beam, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"),debug=4)#
plot.adv(beam.snippet,which="u1",type='l')#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk)
166246*8
166246/8
166246/8 / 6
3600*8*6
3600*8*6 / 166246
hist(diff(m05.vector.beam$data$ts$time))
hist(diff(as.numeric(m05.vector.beam$data$ts$time)))
hist(diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
hist(diff(as.numeric(m05.vector.beam$data$ts.slow$time)),breaks=100)
hist(diff(as.numeric(m05.vector.beam$data$ts.slow$time)),breaks=1000)
summary(diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
stem(diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
unique(diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
which(600=diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
which(600=diff(as.numeric(m05.vector.beam$data$ts.slow$time))
which(600==diff(as.numeric(m05.vector.beam$data$ts.slow$time))
which(600==diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
diff(which(600==diff(as.numeric(m05.vector.beam$data$ts.slow$time))))
wtf <- which(600==diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
m05.vector.beam$data$ts.slow$time[wtf]
m05.vector.beam$data$ts.slow$time[sort(wtf,wtf+1)]
m05.vector.beam$data$ts.slow$time[sort(c(wtf,wtf+1))]
m05.vector.beam$data$ts.slow$time[sort(c(wtf,wtf-1,wtf+1))]
plot(m05.vector.beam$data$ts.slow$time[sort(c(wtf,wtf-1,wtf+1))])
plot(m05.vector.beam$data$ts.slow$time[sort(c(wtf,wtf-1,wtf+1))],type='l')
plot(m05.vector.beam$data$ts.slow$time[wtf],type='l')
plot(m05.vector.beam$data$ts.slow$time[wtf],rep(1,length(wtf)),type='l')
plot(m05.vector.beam$data$ts.slow$time[wtf],rep(1,length(wtf)),type='p')
plot(m05.vector.beam$data$ts.slow$time[wtf],m05.vector.beam$data$ts.slow$time[wtf]-m05.vector.beam$data$ts.slow$time[wtf-1],type='p')
plot(m05.vector.beam$data$ts.slow$time[wtf],m05.vector.beam$data$ts.slow$time[wtf]-m05.vector.beam$data$ts.slow$time[wtf+1],type='p')
plot(m05.vector.beam$data$ts.slow$time[-1],diff(as.numeric(m05.vector.beam$data$ts.slow$time))
)
plot(m05.vector.beam$data$ts.slow$time[-1],diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
plot(m05.vector.beam$data$ts.slow$time[-1],diff(as.numeric(m05.vector.beam$data$ts.slow$time)),type='l')
plot(m05.vector.beam$data$ts$time[-1],diff(as.numeric(m05.vector.beam$data$ts$time)),type='l')
plot(m05.vector.beam$data$ts$time[-1],diff(as.numeric(m05.vector.beam$data$ts$time)),type='p')
par(mfrow=c(2,1))#
oce.plot.ts(m05.vector.beam$data$ts.slow$time[-1],diff(as.numeric(m05.vector.beam$data$ts.slow$time)),type='p')#
oce.plot.ts(m05.vector.beam$data$ts$time[-1],diff(as.numeric(m05.vector.beam$data$ts$time)),type='p')
x11(title='focussed reading')#
plot(d, which='u1')#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr)
x11(title='focussed reading')#
plot(d, which='u1')#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr, line=-2)
x11(title='subsetting (moved time)')#
load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_beam.rda")#
beam.snippet <- subset(m05.vector.beam, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet,which="u1",type='l',debug=4)#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk)
1/.1299
oce.plot.ts(m05.vector.beam$data$ts.slow$time[-1],diff(as.numeric(m05.vector.beam$data$ts.slow$time)),type='p')
mean(diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
t.test(diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
boxplot(1/diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
abline(h=0)
abline(h=1)
boxplot(diff(as.numeric(m05.vector.beam$data$ts.slow$time)))
length(m05.vector.beam$data$ts.slow$time)
boxplot(diff(as.numeric(m05.vector.beam$data$ts.slow$time[1:623934])))
boxplot(diff(as.numeric(m05.vector.beam$data$ts.slow$time[1:(.5*623934)])))
unique(diff(as.numeric(m05.vector.beam$data$ts.slow$time[1:(.5*623934)])))
unique(diff(as.numeric(m05.vector.beam$data$ts.slow$time[1:(.9*623934)])))
unique(diff(as.numeric(m05.vector.beam$data$ts.slow$time[1:(.95*623934)])))
unique(diff(as.numeric(m05.vector.beam$data$ts.slow$time[1:(.99*623934)])))
unique(diff(as.numeric(m05.vector.beam$data$ts.slow$time[1:(.97*623934)])))
source('~/src/R-kelley/oce/R/adv.nortek.R')#
source('~/src/R-kelley/oce/R/adv.R')#
source('~/src/R-kelley/oce/R/oce.R')
source('~/src/R-kelley/oce/R/adp.nortek.R')
m05.vector.beam2 <- read.adv.nortek("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",#
                                   from=as.POSIXct("2008-06-26 00:00:00", tz="UTC"),#
                                   to=as.POSIXct("2008-07-02 00:00:00", tz="UTC"),#
                                   debug=4)
unique(diff(as.numeric(m05.vector.beam2$data$ts.slow$time)))
x11(title='subsetting (moved time) beam2')#
load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_beam.rda")#
beam.snippet2 <- subset(m05.vector.beam2, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet2,which="u1",type='l',debug=4)#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk)
x11(title='subsetting (moved time) beam2')#
load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_beam.rda")#
beam.snippet2 <- subset(m05.vector.beam2, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet2,which="u1",type='l',debug=4)#
abline(v=flag.dk,col="red",lwd=3);mtext(format(flag.dk),side=3,at=flag.dk);abline(v=flag.cr,col="blue",lwd=3);mtext(format(flag.dk),side=3,at=flag.ck,line=-2);
x11(title='subsetting (moved time) beam2')#
load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_beam.rda")#
beam.snippet2 <- subset(m05.vector.beam2, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet2,which="u1",type='l',debug=4)#
abline(v=flag.dk,col="red",lwd=3);mtext(format(flag.dk),side=3,at=flag.dk,line=-2);abline(v=flag.cr,col="blue",lwd=3);mtext(format(flag.dk),side=3,at=flag.cr,line=-2);
m05.vector.beam2 <- read.adv.nortek("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",#
                                   from=as.POSIXct("2008-06-26 00:00:00", tz="UTC"),#
                                   to=as.POSIXct("2008-07-02 00:00:00", tz="UTC"),#
                                   debug=4)#
x11(title='subsetting (moved time) beam2')#
beam.snippet2 <- subset(m05.vector.beam2, #
                        as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time &#
                        time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet2,which="u1",type='l',debug=4)#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk,line=-2)#
abline(v=flag.cr,col="blue",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.cr,line=-2)
m05.vector.beam2 <- read.adv.nortek("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",#
                                   from=as.POSIXct("2008-06-26 00:00:00", tz="UTC"),#
                                   to=as.POSIXct("2008-07-03 12:00:00", tz="UTC"),#
                                   debug=4)#
x11(title='subsetting (moved time) beam2 OLDOLD')#
beam.snippet2 <- subset(m05.vector.beam2, #
                        as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time &#
                        time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet2,which="u1",type='l',debug=4)#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.dk,line=-2)
summary(beam.snippet2)
library(oce)#
from <- as.POSIXct("2008-06-29 00:00:00", tz="UTC")#
to <- as.POSIXct("2008-06-29 06:00:00", tz="UTC")#
flag.cr <- as.POSIXct("2008-06-29 01:31:24",tz="UTC")#
flag.dk <- flag.cr + 3 * 3600
m05.vector.beam <- read.adv.nortek("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",#
                                   from=as.POSIXct("2008-06-26 00:00:00", tz="UTC"),#
                                   to=as.POSIXct("2008-07-03 12:00:00", tz="UTC"),#
                                   debug=4)
beam.snippet <- subset(m05.vector.beam, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet,which="u1",type='l',debug=4)#
abline(v=flag.dk,col="red",lwd=3)#
mtext(format(flag.dk),side=3,at=flag.dk)
plot.adv(beam.snippet,which="u1",type='l',debug=4)#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr)
plot.adv(beam.snippet,which="u1",type='l',debug=4)#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr, line=-3)
plot.adv(beam.snippet,which="u1",type='l',debug=4)#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr, line=-1)
Incorrect timing (prior to 2010-09-15), correct after that time.#
x11(title='demonstrate bug fix on adv timing')#
#load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_beam.rda")#
m05.vector.beam <- read.adv.nortek("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",#
                                   from=as.POSIXct("2008-06-26 00:00:00", tz="UTC"),#
                                   to=as.POSIXct("2008-07-03 12:00:00", tz="UTC"),#
                                   debug=4)#
beam.snippet <- subset(m05.vector.beam, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet,which="u1",type='l',debug=4)#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr, line=-1)
library(oce) # needs version as of 946AM, Wed 2010-09-15#
x11(title='demonstrate bug fix on adv timing')#
#load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_beam.rda")#
m05.vector.beam <- read.adv.nortek("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/raw/adv_nortek_1943.vec",#
                                   from=as.POSIXct("2008-06-26 00:00:00", tz="UTC"),#
                                   to=as.POSIXct("2008-07-03 12:00:00", tz="UTC"))#
beam.snippet <- subset(m05.vector.beam, as.POSIXct("2008-06-29 00:00:00", tz="UTC") <= time & time <= as.POSIXct("2008-06-29 06:00:00", tz="UTC"))#
plot.adv(beam.snippet,which="u1",type='l')#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr, line=-1)
plot.adv(beam.snippet,which="u1",type='l')#
flag.cr <- as.POSIXct("2008-06-29 01:31:24",tz="UTC")#
abline(v=flag.cr,col="red",lwd=3)#
mtext(format(flag.cr),side=3,at=flag.cr, line=-1)
load("/data/archive/sleiwex/2008/moorings/m10/adp/rdi_1313/r/m10_rdi_enu.rda")
plot(m10.rdi.enu, which='progressive vector')
plot(m10.rdi.enu)
plot(m10.rdi.enu, which='uv+ellipse+arrow')
source('~/src/R-kelley/oce/R/adp.R')
plot.adp(m10.rdi.enu, which='uv+ellipse+arrow')
??s4
f<-"data/archive/sleiwex/2008/moorings/m11/cm/interocean_0811786/manufacturer/cm_interocean_0811786.s4a.tab"
f<-"/data/archive/sleiwex/2008/moorings/m11/cm/interocean_0811786/manufacturer/cm_interocean_0811786.s4a.tab"
source('~/src/R-kelley/oce/R/oce.R'); magic(f)
source('~/src/R-kelley/oce/R/oce.R'); d<read.oce(f, debug=4)
source('~/src/R-kelley/oce/R/oce.R'); d<-read.oce(f, debug=4)
magick(f)
magic(f)
read.cm.s4
l[1]
spllit(l[1])
split(l[1])
split(l[1],'\t')
split('\t',l[1])
split('\\t',l[1])
?split
?strsplit
strsplit(l[1],'\t')
strsplit(l[1],'\t')[[1]]
source('~/src/R-kelley/oce/R/cm.R');source('~/src/R-kelley/oce/R/oce.R'); d<-read.oce(f, debug=4)
grep(l1.items[1], "Serial")
l1.items[1]
agrep(l1.items[1], "Serial")
?agrep
agrep("Serial", l1.items[1])
agrep("Saerial", l1.items[1])
grep("Saerial", l1.items[1])
grep("Serial", l1.items[1])
source('~/src/R-kelley/oce/R/cm.R');source('~/src/R-kelley/oce/R/oce.R'); d<-read.oce(f, debug=4);d
d
?read.table
d$data
length(d$data)
?readLines
system.time(l<-readLines("/data/archive/sleiwex/2008/moorings/m11/cm/interocean_0811786/manufacturer/cm_interocean_0811786.s4a.tab"))
length(l)
d$data[1]
d$data[2]
strsplit('\t', d$data[2])
strsplit(d$data[2],'\t')
strsplit(d$data[2],'\t')[[1]]
length(strsplit(d$data[2],'\t')[[1]])
length(strsplit(d$data[1],'\t')[[1]])
str(l)
for (i in 1:100)print(length(strsplit(d$data[i],'\t')[[1]]))
d<-read.table(f, skip=6, split='\t')
d<-read.table(f, skip=6, sep='\t')
d<-read.table(f, skip=60, sep='\t')
d<-read.table(f, skip=30, sep='\t')
d<-read.table(f, skip=20, sep='\t')
d<-read.table(f, skip=10, sep='\t')
d<-read.table(f, skip=11, sep='\t')
d<-read.table(f, skip=12, sep='\t')
d<-read.table(f, skip=13, sep='\t')
dim(d)
plot(d[,4],d[,5])
d[,1]
d[,10]
d[1:3,10]
d[1,10:15]
d[10:15,1]
d[10:15,2]
d[10:15,3]
d[10:15,4]
d[10:15,5]
d[10:15,6]
plot(d[,5],d[,6])
plot(d[,5],d[,6],type='l')
plot(d[,5],d[,6],type='l',asp=1)
d[1:20,1]
d[1,]
d<-read.table(f, skip=15, sep='\t')
d<-read.table(f, skip=16, sep='\t')
d<-read.table(f, skip=16, sep='\t',stringsAsFactors=FALSE)
plot(d[,5]/100,d[,6]/100,type='l',asp=1)
?seek()
?seek
?scan
scan(f,"character",nlines=1)
scan(f,"character",nlines=1,sep="\t")
traceback()
plot(d$data$u,d$data$v)
plot(d$data$ts$u,d$data$ts$v)
hms[10]
hms[1]
hms[2]
??as.POSIXct
strptime("6/25/2008","%m/%d/%y")
strptime("6/25/2008","%m/%d/%T")
strptime("6/25/2008","%m/%d/%Y")
strptime("6/25/2008","%m/%d/%Y",tz="UTC")
str(d)
oce.plot.ts(d$data$ts$time, d$data$ts$u)
oce.plot.ts(d$data$ts$time, d$data$ts$v)
class(d)
inherits(d,"s4")
source('~/src/R-kelley/oce/R/cm.R');source('~/src/R-kelley/oce/R/oce.R'); d<-read.oce(f, debug=4);str(d)
summary(d)
source('~/src/R-kelley/oce/R/cm.R');source('~/src/R-kelley/oce/R/oce.R'); d<-read.oce(f, debug=4);summary(d)
source('~/src/R-kelley/oce/R/cm.R');source('~/src/R-kelley/oce/R/oce.R'); d<-read.oce(f, debug=4);summary(d);plot(d)
plot(m11.cm)
load("/data/archive/sleiwex/2008/moorings/m11/cm/interocean_0811786/r/m11_cm.rda")
plot(m11.cm, which="uv+ellipse+arrow")
plot(m11.cm, which="progressive vector")
source('~/src/R-kelley/oce/R/oce.R')
source('~/src/R-kelley/oce/R/misc.R')
m11 <- apply.magnetic.declination(m11.cm, declination=-18)
plot(m11,'progressive vector')
plot(m11,'uv')
plot(m11,'uv+ellipse+arrow')
source('m11_cm.R')
summary(m11.cm)
setwd("/Users/kelley/Documents/my-papers/inprep/internal-beach/work/data/R")
f <- file('m05_sontek_adv','rb')
seek(f,0,"end")
n <- seek(f,0,"start")
n
buf <- readBin(f, "raw", n)
length(buf)
v.start.1 <- match.bytes(buf[1:min(2000, length(buf))], 0x85, 0x16)[1]
v.start.1
match.bytes(buf, 0x85, 0x16)
match.bytes(buf, 0x85, 0x16)[1]
match.bytes(buf, 0x85, 0x16)[2]
p<-match.bytes(buf, 0x85, 0x16)
diff(p)
stem(diff(p))
hist(diff(p),breaks=100)
d<-read.adv.sontek('m05_sontek_adv', from=1, to=10000,header=FALSE)
p2 <- sort(c(p, p+1))
p2[1:10]
readBin(buf[p2], what="integer", size=2, endian="little")
readBin(buf[p2], what="integer", size=2, endian="little", n=length(p))
readBin(buf[p2]+2, what="integer", size=2, endian="little", signed=FALSE, n=length(p))
readBin(buf[p2+2], what="integer", size=2, endian="little", signed=FALSE, n=length(p))
sample.number <- readBin(buf[p2+2], what="integer", size=2, endian="little", signed=FALSE, n=length(p))
plot(p, sample.number)
x11()
plot(p[-1], diff(sample.number), type='p')
ok <- diff(sample.number) == 1 # should allow for the wrap also
sum(ok)/length(ok)
1-sum(ok)/length(ok)
1/256^2
length(ok)
unique(diff(p[ok]))
stem(diff(p[ok]))
plot(p[ok],diff(p[ok]))
plot(p[ok][-1],diff(p[ok]))
plot(p[ok][-1],diff(p[ok]),type='l')
0xff
c(0xff,0xff)
readBin(c(0xff,0xff), "integer")
readBin(as.raw(c(0xff,0xff)), "integer")
readBin(as.raw(c(0xff,0xff)), "integer", size=2, endian="little", signed=FALSE)
0x85
as.raw(10)
as.raw(128)
as.raw(256)
as.raw(255)
as.raw(22)
setwd("/Users/kelley/src/R-kelley/oce/src")
.Call("ldc_sontek_adv_22", buf, 0)
readBin(as.raw(c(0xa5,0x96)), "integer", size=2, endian="little", signed=FALSE)
dyn.load("bitwise.so")
.Call("ldc_sontek_adv_22", buf, 10)
39464 -  38565
(39464 -  38565) / 20
as.raw(45)
dyn.load("bitwise.so"); .Call("ldc_sontek_adv_22", buf, 10)
f <- file('m05_sontek_adv','rb')#
  seek(f,0,"end")#
  n <- seek(f,0,"start")#
  buf <- readBin(f, "raw", n)
f
9*16+6
setwd("/Users/kelley/src/R-kelley")
source("~/a.R")
length(p)
temperature <- readBin(buf[sort(c(p, p+1))+16], "integer", signed=TRUE, size=2, endian="little")
temperature <- readBin(buf[sort(c(p, p+1))+16], "integer", signed=TRUE, size=2, endian="little",n=length(p))
length(temperature)
temperature <- 0.01 * readBin(buf[sort(c(p, p+1))+16], "integer", signed=TRUE, size=2, endian="little",n=length(p))
summary(temperature)
boxplot(temperature)
buf[p[1:100]]
f <- file('~/m05_sontek_adv','rb')#
seek(f,0,"end")#
n <- seek(f,0,"start")#
buf <- readBin(f, "raw", n)
p <- .Call("ldc_sontek_adv_22", buf, 0)
p[1:10]
buf[p]
buf[p-1][1:`10`]
buf[p-1][1:10]
ps <- p[1:20]
ps
buf[ps]
buf[ps+1]
temperature <- 0.01 * readBin(buf[sort(c(p, p+1))+17], "integer", signed=TRUE, size=2, endian="little",n=length(p))
hist(temperature)
system.time(temperature <- 0.01 * readBin(buf[sort(c(p, p+1))+17], "integer", signed=TRUE, size=2, endian="little",n=length(p)))
system.time(sample.number <- readBin(buf[sort(c(p, p+1))+2], "integer", signed=FALSE, size=2, endian="little",n=length(p)))
x11(); plot(sample.number[1:1000], temperature[1:1000], type='l')
x11(); plot(sample.number[1:10000], temperature[1:10000], type='l')
plot(temperature[1:10000], type='l')
u1 <- 1e-4*readBin(buf[sort(c(p, p+1))+3], "integer", signed=FALSE, size=2, endian="little",n=length(p)))
u1 <- 1e-4*readBin(buf[sort(c(p, p+1))+3], "integer", signed=FALSE, size=2, endian="little",n=length(p))
u1 <- 1e-4*readBin(buf[sort(c(p, p+1))+5], "integer", signed=FALSE, size=2, endian="little",n=length(p))
plot(u1[1:50000],type='l')
plot(u1[1:5000],type='l')
plot(u1[1:500],type='l')
?plot.ctd
source("/tmp/r-plugin-kelley/Rsource-40792-11-b.R", echo=TRUE)
source('~/src/R-kelley/oce/R/adp.R');source("/tmp/r-plugin-kelley/Rsource-40792-11-b.R", echo=TRUE)
source("/tmp/r-plugin-kelley/Rsource-40792-11.R", echo=TRUE)
data(gantt)
plot(gantt)
plot(read.gantt('~/plan.dat'))
setwd("/Users/kelley")
library(plan)
?gantt
d<-read.gantt('~/plan.dat')
plot(d)
60 * 60 * 25)
60 * 60 * 25
source('~/src/R-kelley/oce/R/oce.R');source('~/src/R-kelley/oce/R/adv.R');plot(m03.adv.xyz,which='roll')
source('~/src/R-kelley/oce/R/oce.R');source('~/src/R-kelley/oce/R/adv.R');plot(m03.adv.xyz,which='roll',type='p')
source('~/src/R-kelley/oce/R/oce.R');source('~/src/R-kelley/oce/R/adv.R');plot(m03.adv.xyz,which='roll',type='p',debug=10)
source('~/src/R-kelley/oce/R/oce.R');source('~/src/R-kelley/oce/R/adv.R');plot(m03.adv.xyz,which='roll',type='p',debug=0)
source("/tmp/r-plugin-kelley/Rsource-64829-02.R", echo=TRUE)
?filter()
source("/tmp/r-plugin-kelley/Rsource-22371-protractor.R", echo=TRUE)
m07.vector.beam <- read.oce("/data/archive/sleiwex/2008/moorings/m07/adv/nortek_1063/raw/adv_nortek_1063.vec",
from=as.POSIXct("2008-06-24 00:00:00", tz="UTC"),
to=as.POSIXct("2009-07-06 00:00:00", tz="UTC"),
latitude=47.87961, longitude=-69.72706)
m07.vector.beam <- subset(m07.vector.beam, as.POSIXct("2008-06-25 14:00:00", tz="UTC") <= time &
time <= as.POSIXct("2008-07-04 18:00:00", tz="UTC"))
median(m07.vector.beam$data$ts$heading)
heading <- despike(angle.remap(m07.vector.beam$data$ts.slow$heading) - 90, "trim", min=-180, max=-90)
heading <- ifelse(heading < 0, heading + 360, heading) # remap to range 0 to 360
m07.vector.beam$data$ts.slow$heading <- heading
median(m07.vector.beam$data$ts.slow$heading)
source("/tmp/r-plugin-kelley/Rsource-22610-m07_rdi.R")
plot(m07.rdi.beam, which='heading')
m07.rdi.xyz.tmp <- adp.beam2xyz(m07.rdi.beam)
m07.rdi.enu.tmp <- adp.beam2xyz(m07.rdi.xyz.tmp, declination=-18.099)
m07.rdi.enu.tmp <- adp.xyz2enu(m07.rdi.xyz.tmp, declination=-18.099)
plot(m07.rdi.enu.tmp, which='uv+ellipse+arrow')
ls()
load("/data/archive/sleiwex/2008/moorings/m07/adp/rdi_1997/r/m07_rdi_enu.rda")
length(m07.rdi.enu$data$ss$distance)
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=30))
load("/data/archive/sleiwex/2008/moorings/m06/adp/rdi_5249/r/m06_rdi_enu.rda")
length(m06.rdi.enu$data$distance)
length(m06.rdi.enu$data$ss$distance)
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=1))
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=10))
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=20))
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=29))
plot(m06.rdi.enu, which='uv+ellipse+arrow', control=list(bin=1),main="m06 RDI bin 1")
plot(m06.rdi.enu, which='uv+ellipse+arrow', control=list(bin=10),main="m06 RDI bin 10")
plot(m06.rdi.enu, which='uv+ellipse+arrow', control=list(bin=20),main="m06 RDI bin 20")
plot(m06.rdi.enu, which='uv+ellipse+arrow', control=list(bin=29),main="m06 RDI bin 29")
dev.new()
par(mfrow=c(2,2))
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=1),main="m07 RDI bin 1")
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=10),main="m07 RDI bin 10")
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=20),main="m07 RDI bin 20")
plot(m07.rdi.enu, which='uv+ellipse+arrow', control=list(bin=29),main="m07 RDI bin 29")
setwd("/Users/kelley/Documents/my-papers/inprep/internal-beach/figures/map")
?plot.coastline
source("/tmp/r-plugin-kelley/Rsource-31679-map.R", echo=TRUE)
source("/tmp/r-plugin-kelley/Rsource-20219-map-3.R", echo=TRUE)
source("/tmp/r-plugin-kelley/Rsource-74388-map-4new.R", echo=TRUE)
axis(1)
axis(2)
print(axis(2))
print(axis(1))
plotInset(-69.8,47.87,-69.74,47.90,plot(coastlineLSE,clatitude=lat0,clongitude=lon0,span=300,inset=TRUE))
plotInset(-69.8,47.87,-69.74,47.90,plot(coastlineSLE,clatitude=lat0,clongitude=lon0,span=300,inset=TRUE))
plotInset(-69.8,47.87,-69.74,47.90,plot(coastlineSLE,clatitude=lat0,clongitude=lon0,span=300,inset=TRUE,mar=c(3,3,1,1)))
plotInset(-69.8,47.87,-69.74,47.90,plot(coastlineSLE,clatitude=lat0,clongitude=lon0,span=300,inset=TRUE),mar=c(3,3,1,1))
plotInset(-69.8,47.87,-69.74,47.90,plot(coastlineSLE,clatitude=lat0,clongitude=lon0,span=300,inset=TRUE),mar=c(2,2,1,1))
plot(coastlineSLE, clatitude=lat0, clongitude=lon0, span=16, lwd=2,axes=FALSE, mar=c(0.5, 0.5, 1.25, 0.5))
plotInset(-69.8,47.87,-69.74,47.90,{plot(coastlineSLE,clatitude=lat0,clongitude=lon0,span=300,inset=TRUE);rect(-69.82,47.78,-69.70,47.91)},mar=c(2,2,1,1))
?read.coastline
library(oce)#
        d <- read.coastline.shapefile("~/Dropbox/DepthContours/DepthContours.shp")#
        isna <- is.na(d[["latitude"]])#
        idx<-1+cumsum(isna)#
        lat<-split(d[["latitude"]][!isna], idx[!isna])#
        lon<-split(d[["longitude"]][!isna], idx[!isna])#
        depths <- d[["depths"]]#
        n <- length(depths)#
        D<-200#
        plot.window(c(-180,180),c(-90,90), xlab="", ylab="", asp=1)#
        for (i in 1:n) {#
            if (depths[i]==D)#
                lines(lon[[i]],lat[[i]])#
        }
source('~/src/R-kelley/oce/R/coastline.R')
d <- read.coastline.shapefile("~/Dropbox/DepthContours/DepthContours.shp")
?plotInset
rect(xleft,ybottom,yright,ytop,col='red')
rect(xleft,ybottom,ybottom,ytop,col='red')
xleft
xright
ybottom
points(xleft,ybottom)
points(xleft,ytop)
ytop<par('usr')[3]
ytop<par('usr')[4]
ytop
points(xleft,ytop,cex=10)
usr
C
par(no.readonly=TRUE)
?text
usrSLE
source("/tmp/r-plugin-kelley/Rsource-92955-map-5.R", echo=TRUE)
par('usr');par('mar')
?rect
par('mar')
uu<-par('usr')
rect(uu[1],uu[3],uu[2],uu[4],col='pink')
c(uu[1],uu[3],uu[2],uu[4])
par('usr')
source("/tmp/r-plugin-kelley/Rsource-7264-map-5.R", echo=TRUE)
help(plotInset)
par('usr');par('mar');par('mai');par('fin');par('pin')
rect
unif
unif(3)
?rnorm
dnorm(10)
rnorm(10)
runif(10)
?runif
source("/tmp/r-plugin-kelley/Rsource-97077-test.R", echo=TRUE)
setwd("/Users/kelley/src/R-kelley/oce/inst/doc")
source("/tmp/r-plugin-kelley/Rsource-16369-oce.R", echo=TRUE)
source("/tmp/r-plugin-kelley/Rsource-16990-oce.R", echo=TRUE)
add
Slim
Tlim
plot.window(Slim, Tlim,
xlab = xlab, ylab=ylab,
xaxs = if (min(x@data$salinity,na.rm=TRUE)==0) "i" else "r", # avoid plotting S<0
cex=cex, pch=pch, col=col, cex.axis=par("cex.axis"),
...)
if (!missing(bg)) {
usr <- par('usr')
rect(usr[1], usr[3], usr[2], usr[4], col=bg)
}
connectPoints
if (grid)
grid(col=col.grid, lty=lty.grid)
?plot.window
Q
dev.off()
library(oce)
swRho(34, 10, 100)
swTheta(34, 10, 100)
swRho(34, swTheta(34, 10, 100), 0)
swRho(34, swTheta(34, 10, 100, 200), 200)
source('~/src/R-kelley/oce/R/ctd.R')
plotTS(as.ctd(c(30,40),c(-2,20),rep(0,2)), grid=TRUE, col="white")
q()
