utils:build
library(utils)
?utils
source("/tmp/r-plugin-kelley/Rsource-66034-0.R", echo=TRUE)
plot(argo,which=2,type='p')
a<-function(x) { dan <- dan + 1; cat('x=',x,'\n')}
dan<-0
a('hi')
dan
data(co2)
source('R/spectral.R');pw<-pwelch(co2); plot(pw$freq, pw$spec, log='y')
dev.new();plot(co2)
dev.set(2);source('R/spectral.R');pw<-pwelch(co2); plot(pw$freq, pw$spec, log='y')
dev.set(2);source('R/spectral.R');pw<-pwelch(co2); plot(pw$freq, pw$spec, log='y',type='l')
dev.new;spectrum(co2)
spectrum(co2); source('R/spectral.R');pw<-pwelch(co2); line(pw$freq, pw$spec, col='red')
spectrum(co2); source('R/spectral.R');pw<-pwelch(co2); lines(pw$freq, pw$spec, col='red')
spectrum(co2,log='no'); source('R/spectral.R');pw<-pwelch(co2); lines(pw$freq, pw$spec, col='red')
plot(pwelch(co2))
s<-spectrum(co2)
spectrum(co2,log='no')#; source('R/spectral.R');pw<-pwelch(co2); plot(pw);#lines(pw$freq, pw$spec, col='red')
s$method
s$series
spectrum.pgm
spectrum.prgm
spectrum.pgram
spectrum(co2,log='no'); source('R/spectral.R');pw<-pwelch(2*co2); plot(pw);#lines(pw$freq, pw$spec, col='red')
dan<-co2
spectrum(co2,log='no'); source('R/spectral.R');pw<-pwelch(dan); plot(pw);#lines(pw$freq, pw$spec, col='red')
spectrum(co2)
s$bandwidth
s$df
spectrum(co2,log='no'); source('R/spectral.R');pw<-pwelch(co2); plot(pw);#lines(pw$freq, pw$spec, col='red')
spectrum(co2,log='no'); source('R/spectral.R');pw<-pwelch(co2); #plot(pw);#lines(pw$freq, pw$spec, col='red')
names(s)
spectrum
source('R/spectral.R');pwelch(co2)
?spectrum
spec.pgram
pwelch(co2)
pwelch(co2, window=3)
pwelch(co2, window=10)
pwelch(co2, window=1)
pwelch(co2, window=2)
pwelch(co2, window=2,log='no')
par(mfrow=c(2,1)); spectrum(co2); pwelch(co2)
par(mfrow=c(2,1)); spectrum(co2,log='no'); pwelch(co2, log='no')
par(mfrow=c(2,1)); spectrum(co2,log='no'); pwelch(co2, window=4, log='no')
par(mfrow=c(2,1)); spectrum(co2,log='no'); pwelch(co2, window=3, log='no')
20000/50
20000/50/8
dir("/data/archive/sleiwex/2008/moorings/m05/")
dir("/data/archive/sleiwex/2008/moorings/m05/adv")
dir("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943")
dir("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r")
dir("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_enu.rda")
adv <- read.oce("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_enu.rda")
load("/data/archive/sleiwex/2008/moorings/m05/adv/nortek_1943/r/m05_vector_enu.rda")
ls()
plot(m05VectorEnu)
adv<-subset(m05VectorEnu, 1,10)
?subset.adv
?subset.oce
adv<-window(m05VectorEnu, as.POSIXct("2008-07-01 00:00:00", tz="UTC"), as.POSIXct("2008-07-01 00:00:10")
??size
adv<-window(m05VectorEnu, as.POSIXct("2008-07-01 00:00:00", tz="UTC"), as.POSIXct("2008-07-01 00:00:01")
adv<-window(m05VectorEnu, as.POSIXct("2008-07-01 00:00:00", tz="UTC"), as.POSIXct("2008-07-01 00:00:01"))
adv<-window(m05VectorEnu, as.POSIXct("2008-07-01 00:00:00", tz="UTC"), as.POSIXct("2008-07-01 00:00:01",tz="UTC"))
adv<-window(m05VectorEnu, as.POSIXct("2008-07-01 00:00:00", tz="UTC"), as.POSIXct("2008-07-01 00:01:00",tz="UTC"))
dev.new();spectrum(adv$data$v[,3])
spectrum(adv$data$v[,3],log='both')
spectrum(adv$data$v[,3],log='yes')
spectrum(adv$data$v[,3],log='dB')
s<-spectrum(adv$data$v[,3],log='dB')
plot(log10(s$freq), log10(s$spec))
plot(log10(s$freq), log10(s$spec),type='l')
plot(log10(s$freq), log10(s$spec),type='l');abline(a=seq(-10,10,1),b=-5/3)
plot(log10(s$freq), log10(s$spec),type='l');abline(a=seq(-10,10,.1),b=-5/3)
plot(log10(s$freq), log10(s$spec),type='l');abline(a=seq(-10,10,.1),b=rep(-5/3,20))
plot(log10(s$freq), log10(s$spec),type='l');abline(a=-5,b=-5/3)
plot(log10(s$freq), log10(s$spec),type='l');abline(a=-3,b=-5/3)
object.size(adv)
object.size(adv)/1000
save(adv, file='~/test.rda')
system("ls -lh ~/test.rda")
adv<-window(m05VectorEnu, as.POSIXct("2008-07-01 00:00:00", tz="UTC"), as.POSIXct("2008-07-01 00:02:00",tz="UTC"))
s<-spectrum(adv$data$v[,3],plot=FALSE)
plot(log10(s$freq), log10(s$spec),type='l');abline(a=-7,b=-5/3)
dev.new()
plot(adv)
summary(adv)
setwd("/Users/kelley/src/R-kelley/oce")
data(adp)#
   plot(adp, which=1, zlim=c(-0.5, 0.5))
traceback()
debug(plot.adp)
source('R/adp.R')
plot(adp, which=1, zlim=c(-0.5, 0.5))
y.look
w
zlim
zlim[1,]
zlim.given
zlim.given[1,]
source('R/adp.R');plot.adp(adp, which=1, zlim=c(-0.5, 0.5))
page(plot.adp)
?swSpice
source('R/adp.R');plot(adp, which=1)#, zlim=c(-2, 2))
source('R/adp.R');plot(adp, which=1, zlim=c(-2, 2))
source('R/imagep.R');source('R/adp.R');plot(adp, which=1, zlim=c(-2, 2))
par(mfrow=c(2,1));source('R/imagep.R');source('R/adp.R');plot(adp, which=1, zlim=c(-2, 2));plot(adp,which=1)
par(mfrow=c(2,1));source('R/imagep.R');source('R/adp.R');plot(adp,which=1);plot(adp, which=1, zlim=c(-2, 2))
?pt
?pT
?rbr
?ptr
ptr <- read.oce("/data/archive/sleiwex/2008/moorings/m08/pt/rbr_011855/manufacturer/pt_rbr_011855.dat", by=600)
readLines("/data/archive/sleiwex/2008/moorings/m08/pt/rbr_011855/manufacturer/pt_rbr_011855.dat", n=10)
setwd("/Users/kelley/Documents/literature/g/garratt1977rodc_fig3")
plot(profileMethod)
profileMethod <- read.table("profileMethod.txt", skip=1, col.names=c("V","CDN"))
eddyCorrelationMethod <- read.table("eddyCorrelationMethod.txt", skip=1, col.names=c("V","CDN"))
eddyCorrelationMethod
??merge
profileMethod <- data.frame(profileMethod, "profile")
profileMethod <- data.frame(profileMethod, method="profile")
data <- data.frame(profileMethod, eddyCorrelationMethod)
profileMethod
rbind(profileMethod, eddyCorrelationMethod)
data[1,]
lm(CDN~V+method)
lm(CDN~V+method,data=data)
aov(lm(CDN~V+method,data=data))
anova(lm(CDN~V+method,data=data))
plot(data)
names(data)
plot(data$V, data$CDN, cex=method)
plot(data$V, data$CDN, cex=as.factor(method))
?as.factor
as.factor(drag$method)
plot(data$V, data$CDN, levels(drag$method))
plot(data$V, data$CDN, pch=levels(drag$method))
?levels
factor(drag$method)
as.numeric(drag$method)
?points
points(1,1,pch=1)
plot(1,1,pch=1)
plot(1,1,pch=2)
plot(10:20,10:20,pch=10:20)
plot(1:10,1:10,pch=1:10)
drag$method
levels(drag$method)
as.level(drag$method)
as.levels(drag$method)
drag$method == "eddy"
summary(profileModel)
summary(eddyCorrelationModel)
?predict.lm
plot(data$U, data$Cd, pch=ifelse(drag$method == "profile", 1, 19))
new
profileCI
?matplot
aov(lm(Cd~U+method, data=drag))
anova(lm(Cd~U+method, data=drag))
source("/tmp/r-plugin-kelley/Rsource-8206-01.R", echo=TRUE)
setwd("/Users/kelley/src/R-kelley/oce-data/drag")
load('drag.rda')
profileModel <- lm(Cd ~ U, data=profileMethod)
subset(drag, method=="profile)
")
subset(drag, method=="profile")
names(drag)
drag
source("/tmp/r-plugin-kelley/Rsource-9871-02.R", echo=TRUE)
setwd("/Users/kelley/src/R-kelley")
()
?drag
example(drag)
example(plot.adp)
h <- seq(0, 24*365)
elevation <- 2.0 * sin(2*pi*h/12.4172)
sl <- as.sealevel(elevation)
source('oce/R/sealevel.R');sl <- as.sealevel(elevation)
summary(sl)
source('oce/R/sealevel.R');source('oce/R/misc.R');summary(sl)
class(sl)
sl
source('oce/R/sealevel.R');source('oce/R/misc.R');sl<-as.sealevel(elevation,time);summary(sl)
isS4(object)
object@metadata$stationNumber
object@metadata$version
object@latitude
object@metadata$latitude
source('oce/R/sealevel.R');source('oce/R/misc.R');sl<-as.sealevel(elevation);summary(sl)
source('oce/R/sealevel.R');source('oce/R/misc.R');sl<-as.sealevel(elevation);summary.sealevel(sl);data(ctd);summary(ctd)
source('oce/R/sealevel.R');source('oce/R/misc.R');sl<-as.sealevel(elevation);summary.sealevel(sl)
source('oce/R/sealevel.R');source('oce/R/misc.R');summary(section)
source('oce/R/section.R');source('oce/R/sealevel.R');source('oce/R/misc.R');summary(section)
plot(swN2(ctd), ctd$data$pressure,
     ylim=rev(range(ctd$data$pressure)), xlab="N2",ylab="p", type='l')
plot(swN2(ctd), ctd$data$pressure, ylim=rev(range(ctd$data$pressure)), xlab="N2",ylab="p", type='l')
source('oce/R/sw.R');source('oce/R/ctd.R');swN2(ctd)
source('oce/R/sw.R');source('oce/R/ctd.R');plot(swN2(ctd), ctd$data$pressure, ylim=rev(range(ctd$data$pressure)), xlab="N2",ylab="p", type='l')
lines(swN2(ctd, df=10), p, col="blue")
lines(swN2(ctd, df=2), p, col="blue")
p <- pressure(ctd)
plot(swN2(ctd), p, ylim=rev(range(p)), xlab="N2",ylab="p", type='l')
lines(swN2(ctd, df=3), p, col="blue")
grid()
tide <- tidem(sealevelTuktoyaktuk)
summary(tide)
t <- seq(0, 10*86400, 3600)
eta <- sin(0.080511401 * t * 2 * pi / 3600)
sl <- as.sealevel(eta)
m <- tidem(sl)
source('oce/R/tidem.R');m <- tidem(sl)
source('oce/R/ctd.R');data(ctd);write.ctd(ctd,"ctd.csv")
source('oce/R/ctd.R');data(ctd);write.ctd(ctd,"ctd.csv");swRho(ctd)
source('oce/R/sw.R');source('oce/R/ctd.R');data(ctd);write.ctd(ctd,"ctd.csv");swRho(ctd)
source('oce/R/sw.R');source('oce/R/ctd.R');data(ctd);swRho(ctd)
all.equal.numeric(swRho(ctd),          1030.818, 1e-6)
all.equal.numeric(swRho(35, 13, 1000), 1030.818, 1e-6)
source('oce/R/sw.R');source('oce/R/ctd.R');source('oce/tests/eos.R')
lr <- swLapseRate(40, 40, 10000)
stopifnot(all.equal.numeric(lr, 3.255976e-4, 1e-7))
cond <- swConductivity(35,10,100)
stopifnot(all.equal.numeric(cond, 0.618569, 1e-5))
visc <- swViscosity(30,10)
stopifnot(all.equal.numeric(visc, 0.001383779, 1e-7))
ctd <- as.ctd(35, 13, 1000)
stopifnot(all.equal.numeric(swRho(35, 13, 1000), 1030.818, 1e-6))
stopifnot(all.equal.numeric(swRho(ctd),          1030.818, 1e-6))
stopifnot(all.equal.numeric(swSigma(35, 13, 1000), 30.818, 1e-5))
stopifnot(all.equal.numeric(swSigma(ctd),          30.818, 1e-5))
stopifnot(all.equal.numeric(swTheta(35, 13, 1000), 12.858, 1e-3))
stopifnot(all.equal.numeric(swTheta(ctd),          12.858, 1e-3))
stopifnot(all.equal.numeric(swSCTp(1,   15,   0), 35.000000, 1e-6))
stopifnot(all.equal.numeric(swSCTp(1.2, 20,2000), 37.245628, 1e-6))
stopifnot(all.equal.numeric(swSCTp(0.65, 5,1500), 27.995347, 1e-6))
Cp <- swSpecificHeat(40, 40, 10000)
stopifnot(all.equal.numeric(Cp, 3849.500, 1e-3))
TT <- swTFreeze(40, 500)
stopifnot(all.equal.numeric(TT, -2.588567, 1e-6))
ctd <- as.ctd(40, 10, 4000)
stopifnot(all.equal.numeric(swTheta(40,10,4000,0,"unesco"), 9.42648, 1e-5))
stopifnot(all.equal.numeric(swTheta(ctd,0,"unesco"), 9.42648, 1e-5))
ab <- swAlphaOverBeta(40, 10, 4000, isTheta=TRUE)
ab <- swAlphaOverBeta(ctd, isTheta=TRUE)
stopifnot(all.equal.numeric(ab, 0.34763, 0.00005))
b <- swBeta(40, 10, 4000, isTheta=TRUE)
b <- swBeta(ctd, isTheta=TRUE)
stopifnot(all.equal.numeric(b, 7.2088e-4, 0.0005, scale=1e-4))
a <- swAlpha(40, 10, 4000, isTheta=TRUE)
a <- swAlpha(ctd, isTheta=TRUE)
stopifnot(all.equal.numeric(a, 2.5060e-4, 0.0005, scale=1e-4))
source('oce/R/sw.R')
v <- swSoundSpeed(40, 40, 10000)
v <- swSoundSpeed(as.ctd(40,40,10000))
stopifnot(all.equal.numeric(v, 1731.995, 0.001))
## spice (not from any trusted source, merely from the code [2008-10-02]
sp <- swSpice(35,10,100)
stopifnot(all.equal.numeric(sp, 1.131195, 0.0000015))
depth <- swDepth(10000, 30)
stopifnot(all.equal.numeric(depth, 9712.653, 0.001))
pycnocline <- subset.oce(ctd, 5<=pressure & pressure<=12)
pycnocline <- subset(ctd, 5<=pressure & pressure<=12)
?read.csv
methods('ctd')
showMethods('ctd')
showMethods(class='ctd')
plot?ctd
ctd?plot
taraM <- readLines("~/Desktop/8272011165919Dc-36Download1.txt")
taraM
n<-length(tara)
n<-length(taraM)
n
strsplit(taraM[1], ",")
taraM[1,]
taraM[2]
taraM <- readLines("~/Desktop/8272011165919Dc-36Download1.txt",skip=7)
tara<-taraM[7:n]
tara[1]
tara[2]
tara[3]
tara[4]
strsplit(tara[1])
strsplit(tara[1],',')
strsplit(tara[1],',')[[1]]
length(strsplit(tara[1],',')[[1]])
ncol <- length(strsplit(tara[1],',')[[1]])
dim(m)
m[,1]
m[,2]
m[,3]
m[1,]
m[,12]
m[,13]
m[,14]
m[,15]
m[,16]
m[1:3,]
m[,6]
m[1:10,]
source("/tmp/r-plugin-kelley/Rsource-16985-tt.R", echo=TRUE)
?a03
a03 <- read.section("http://cchdo.ucsd.edu/data/onetime/atlantic/a03/a03_hy1.csv", #
	            sectionId="a03",#
                    ship="R/V Professor Multanovskiy", #
		    scientist="Vladimir Tereschenkov",#
		    institute="SIO")
str(a03)
GS <- subset(a03, indices=124:102)
source('oce/R/oce.R');source('oce/R/section.R');GS <- subset.oce(a03, indices=124:102)
source('oce/R/oce.R');source('oce/R/section.R');GS <- subset(a03, indices=124:102)
data(a03)
source('oce/R/section.R');a03 <- read.section("http://cchdo.ucsd.edu/data/onetime/atlantic/a03/a03_hy1.csv", #
	            sectionId="a03",#
                    ship="R/V Professor Multanovskiy", #
		    scientist="Vladimir Tereschenkov",#
		    institute="SIO")
summary(a03)
processingLog(a03)
processingLogShow(a03)
source('oce/R/section.R');GSg <- sectionGrid(GS, p=seq(0,2000,100))
summary(GSg)
a03@metadata
source('oce/R/section.R');a03 <- read.section("http://cchdo.ucsd.edu/data/onetime/atlantic/a03/a03_hy1.csv", #
	            sectionId="a03",#
                    ship="R/V Professor Multanovskiy", #
		    scientist="Vladimir Tereschenkov",#
		    institute="SIO");summary(a03)
GS <- subset.oce(a03, indices=124:102)
GSg <- sectionGrid(GS, p=seq(0,2000,100))
plot(GSg, coastline=coastlineWorld, map.xlim=c(-80,-60),
     adorn=expression({abline(v=200,col='blue')}))
save(a03, file="oce/data/a03.rda")
plot(GSg, coastline=coastlineWorld, map.xlim=c(-80,-60))
plot(GSg)
source('oce/R/section.R');plot(GSg)
source('oce/R/section.R');plot.section(GSg)
class(GSg)
cl <- read.coastline("../../oce-data/coastline/coastline-sle.dat")
cl <- read.coastline("~/src/R-kelley/oce-data/coastline/coastline-sle.dat")
cl["latitude"]
cl[["latitude"]]
summary(cl)
isS4(cl)
cl
list(a=3,b=1:10)
list(a=3,b=1:10)['a']
list(a=3,b=1:10)[['a']]
source('~/src/R-kelley/oce/R/coastline.R');cl<-read.coastline('~/src/R-kelley/oce-data/coastline/coastline-sle.dat')#;cl[["latitude"]]
source('~/src/R-kelley/oce/R/coastline.R');cl<-read.coastline('~/src/R-kelley/oce-data/coastline/coastline-sle.dat');cl[["latitude"]]
source('~/src/R-kelley/oce/R/coastline.R');cl<-read.coastline('~/src/R-kelley/oce-data/coastline/coastline-sle.dat',debug=3);cl[["latitude"]]
source('~/src/R-kelley/oce/R/coastline.R');cl<-read.coastline('~/src/R-kelley/oce-data/coastline/coastline-sle.dat',debug=3);cl["latitude"]
source('~/src/R-kelley/oce/R/coastline.R');cl<-read.coastline('~/src/R-kelley/oce-data/coastline/coastline-sle.dat',debug=3);lat<-cl["latitude"];lon<-cl["lo"]
plot(cl)
source('~/src/R-kelley/oce/R/coastline.R');cl<-read.coastline('~/src/R-kelley/oce-data/coastline/coastline-sle.dat',debug=0);lat<-cl["latitude"];lon<-cl["lo"]
str(cl)
cl@processingLog
str(cl@processingLog)
processingLog(cl)
processingLogShow(cl)
data(section)
summary(section)
processingLog(section)
processingLogShow(section)
source('~/src/R-kelley/oce/R/coastline.R');cl<-read.coastline('~/src/R-kelley/oce-data/coastline/coastline-sle.dat',debug=0);lat<-cl["latitude"];lon<-cl["lo"];plot(cl);summary(cl)
lat<-coastlineHalifax[["latitude"]]
lon<-coastlineHalifax[["lon"]]
plot(lon,lat)
lat
lon<-coastlineHalifax$datalon
lon<-coastlineHalifax$data$lon
lon
lat<-coastlineHalifax$data$lat
plot(as.coastline(lat,lon))
source('oce/R/coastline.R');plot(as.coastline(lat,lon))
coastlineHalifax<-as.coastline(lat,lon)
plot(coastlineHalifax)
summary(coastlineHalifax)
save(coastlineHalifax, file='oce/data/coastlineHalifax.rda')
str(section)
str(coastlineHalifax)
plot(section, coastline=coastlineHalifax)
source('oce/R/section.R');plot(section, coastline=coastlineHalifax)
data(coastlineMaritimes)
#coastlineMaritimes <- as.coastline(coastlineMaritimes$data$lat, coastlineMaritimes$data$lon)
str(coastlineMaritimes)
coastlineMaritimes <- as.coastline(coastlineMaritimes$data$lat, coastlineMaritimes$data$lon)
plot(coastlineMaritimes)
summary(coastlineMaritimes)
save(coastlineMaritimes, file="oce/data/coastlineMaritimes.rda")
data(coastlineSLE)
str(coastlineSLE)
coastlineSLE <- as.coastline(coastlineSLE$data$lat, coastlineSLE$data$lon)
plot(coastlineSLE)
save(coastlineSLE, file='oce/data/coastlineSLE.rda')
data(coastlineWorld)
str(coastlineWorld)
coastlineWorld <- as.coastline(coastlineWorld$data$lat, coastlineWorld$data$lon, filla=TRUE,filename=coastlineWorld$metadata$filename)
coastlineWorld <- as.coastline(coastlineWorld$data$lat, coastlineWorld$data$lon, filla=TRUE)
coastlineWorld <- as.coastline(coastlineWorld$data$lat, coastlineWorld$data$lon)
coastlineWorld@metadata$fillable<-TRUE
coastlineWorld@metadata$filename<-"M_WORLD_BORDERS-0.3.shp"
plot(coastlineWorld)
summary(coastlineWorld)
save(coastlineWorld, file='oce/data/coastlineWorld.rda')
library(Matrix)
Matrix?[]
Matrix?"["
Matrix?mean
?Matrix
promptClass('coastline')
promptMethods(coastline)
promptClass(coastline)
methods::promptClass(coastline)
read.table('~/lock-exchange/brix.dat', header=TRUE)
d<-read.table('~/lock-exchange/brix.dat', header=TRUE)
lm(S~brix, data=d)
m<-lm(S~brix, data=d)
predict(m)
predict(m,4)
predict(m,new=4)
predict(m,new=list(brix=4))
predict(m,new=list(brix=4.3))
predict(m,new=list(brix=4.2))
plot(m)
plot(d$brix, d$S)
plot(d$brix, d$S, type='l')
plot(d$brix, d$S, type='b')
plot(d$brix, d$S, type='p');abline(m)
m<-lm(S~brix-1, data=d)
plot(d$brix, d$S, type='p');abline(m);title(coef(m))
plot(d$brix, d$S, type='p');abline(m);title(sprintf("1 brix = %.2f PSU", Scoef(m)))
plot(d$brix, d$S, type='p');abline(m);title(sprintf("1 brix = %.2f PSU", coef(m))
plot(d$brix, d$S, type='p');abline(m);title(sprintf("1 brix = %.2f PSU", coef(m)))
plot(d$brix, d$S, xlab="Brix value", ylab="S [PSU]", type='p');abline(m);title(sprintf("1 brix = %.2f PSU", coef(m)))
coastline?[[]]
coastline?[[
coastline?"[["
promptMethods('coastline')
data(coastlineHalifax)
coastlineHalifax[["la"]]
coastlineHalifax[["lo"]]
?switch
switch("a", "A"=3)
switch("a", "A"=3,"..."=NA)
switch("a", "A"=3,"a"=NA)
switch("a", "A"=3,".."=stop('not found'))
switch("a", "A"=3,""=stop('not found'))
source('~/src/R-kelley/oce/R/ctd.R');
source('~/src/R-kelley/oce/R/ctd.R'); ctd[["sa"]]
source('~/src/R-kelley/oce/R/ctd.R'); ctd[["te"]]
source('~/src/R-kelley/oce/R/ctd.R'); ctd[["p"]]
source('~/src/R-kelley/oce/R/ctd.R'); ctd[["f"]]
?coastline-class
?"[["
?'coastline-class'
?coastline
??coastline
ctd?"[["
?'ctd-class'
data(ctd)
show(ctd)
          definition=function(.Object,longitudesalinity,temperature,pressure,latitude, longitude, time, filename) {
)
?read.coastline
?read.drifter
d <- read.drifter("/data/OAR/6900388_prof.nc")
summary(d)
d[["latit"]]
source('oce/R/drifter.R');summary(d)
source('oce/R/drifter.R');d<-read.drifter('/data/OAR/6900388_prof.nc')
source('oce/R/drifter.R');d<-read.drifter('/data/OAR/6900388_prof.nc');summary(d)
plot(d)
show(d)
data(argo)
argo
argo<-read.drifter('/data/OAR/6900388_prof.nc')
plot(argo)
summary(argo)
save(argo, file='oce/data/argo.rda')
str(sealevelHalifax)
a<-sealevelHalifax
sealevelHalifax<-new('sealevel')
sealevelHalifax<-new('sealevel',elevation=a$data$elevation, time=a$data$time)
sealevelHalifax@metadata<-a$data$metadata
sealevelHalifax@metadata<-a$metadata
sealevelHalifax@processingLog<-unclass(a$processingLog)
save(sealevelHalifax, file='oce/data/sealevelHalifax.rda')
library(oce)
summary(sealevelHalifax)
plot(sealevelHalifax)
source('oce/R/sealevel.R');plot(sealevelHalifax)
source('oce/R/sealevel.R');plot(sealevelHalifax);summary(sealevelHalifax)
data(sealevelTuttoyaktuk)
data(sealevelTuktoyaktuk)
a<-sealevelTuktoyaktuk
sealevelTuktoyaktuk<-new('sealevel',elevation=a$data$elevation,time=a$data$time)
sealevelTuktoyaktuk@metadata<-a$metadata
sealevelTuktoyaktuk@processingLog<-unclass(a$processsingLog)
sealevelTuktoyaktuk@processingLog<-unclass(a$processingLog)
plot(sealevelTuktoyaktuk)
summary(sealevelTuktoyaktuk)
save(sealevelTuktoyaktuk, file='oce/data/sealevelTuktoyaktuk.rda')
setwd("/Users/kelley/src/R-kelley/oce/inst/doc")
data(sealevelHalifax)
spectrum(sealevelHalifax$data$elevation, spans=c(3,7))
source("/tmp/r-plugin-kelley/Rsource-70759-oce.R", echo=TRUE)
q()
