sandbox/jlh/18_tweets/useAdjusted.R

In dankelley/argoFloats: Analysis of Oceanographic Argo Floats

if (!interactive()) png("useAdjusted.png", unit="in", width=7, height=4, pointsize=11, res=100)
library(argoFloats)
ai <- getIndex(filename = 'synthetic')
index1 <- subset(ai, ID='5903586')
index2 <- subset(index1, cycle='001')
profiles <- getProfiles(index2)
raw <- readProfiles(profiles)
adjusted <- useAdjusted(raw)
rawOxygen <- unlist(raw[['oxygen']])
rawPressure <- unlist(raw[['pressure']])
adjustedOxygen <- unlist(adjusted[['oxygen']])
adjustedPressure <- unlist(adjusted[['pressure']])
xlim <- range(c(rawOxygen, adjustedOxygen), na.rm=TRUE)
plot(rawOxygen, rawPressure,
     pch=1, col=1,
     xlim=xlim, ylim=rev(range(rawPressure, na.rm=TRUE)),
     xlab=expression("Raw Oxygen ["*mu*mol/kg*"]"),
     ylab='Pressure (dbar)')
points(adjustedOxygen, adjustedPressure,
       pch=3, col=2)
legend("bottomright", pch=c(1,3), col=c(1,2), c("Raw", "Adjusted"))
if (!interactive()) dev.off()