sandbox/jlh/30_oddCycle/oddCycle.R

In dankelley/argoFloats: Analysis of Oceanographic Argo Floats

library(argoFloats)
ID <- "1901584"
argos <- readProfiles(getProfiles(subset(getIndex(), ID=ID)))
# Determining levels in argos
for (i in 1:154){
    print(length(argos[['temperature']][[i]]))
}
# Determining max pressure
for (i in 1:154){
    print(max(argos[['pressure']][[i]]))
}

levels <- c(67 ,68 ,68 ,67 ,67 ,68 ,66 ,68 ,68 ,66 ,68 ,68 ,67 ,67 ,68 ,69 ,67 ,67 ,67 ,67 ,69 ,69 ,68 ,69 ,67 ,68 ,69 ,68 ,68 ,67 ,68 ,59 ,69 ,66 ,66 ,66 ,66 ,66 ,68 ,67 ,68 ,68 ,68 ,67 ,67 ,68 ,66 ,70 ,68 ,68 ,67 ,67 ,68 ,68 ,69 ,68 ,69 ,68 ,68 ,66 ,68 ,69 ,66 ,66 ,66 ,67 ,66 ,69 ,66 ,69 ,66 ,68 ,66 ,66 ,69 ,70 ,70 ,67 ,67 ,67 ,67 ,67 ,67 ,67 ,67 ,66 ,67 ,70 ,66 ,67 ,70 ,67 ,70 ,67 ,70 ,70 ,68 ,66 ,66 ,67 ,66 ,67 ,67 ,67 ,67 ,67 ,68 ,66 ,69 ,66 ,66 ,70 ,66 ,66 ,66 ,67 ,67 ,67 ,67 ,70 ,72 ,69 ,67 ,69 ,72 ,69 ,67 ,66 ,64 ,72 ,199 ,63 ,54 ,63 ,54 ,63 ,63 ,72 ,64 ,72 ,72 ,63 ,66 ,63 ,63 ,63 ,27 ,46 ,37 ,9 ,69 ,85 ,71 ,67)
maxPressure <- c(1250 ,1300 ,1300 ,1250 ,1250 ,1300 ,1200 ,1300 ,1300 ,1200 ,1300 ,1300 ,1250 ,1250 ,1300 ,1350 ,1250 ,1250 ,1250 ,1250 ,1350 ,1350 ,1300 ,1350 ,1250 ,1300 ,1350 ,1300 ,1300 ,1250 ,1300 ,850 ,1350 ,1200 ,1200 ,1200 ,1200 ,1200 ,1300 ,1250 ,1300 ,1300 ,1300 ,1250 ,1250 ,1300 ,1200 ,1400 ,1300 ,1300 ,1250 ,1250 ,1300 ,1300 ,1350 ,1300 ,1350 ,1300 ,1300 ,1200 ,1300 ,1350 ,1200 ,1200 ,1200 ,1250 ,1200 ,1350 ,1200 ,1350 ,1200 ,1300 ,1200 ,1200 ,1350 ,1400 ,1400 ,1250 ,1250 ,1250 ,1250 ,1250 ,1250 ,1250 ,1250 ,1200 ,1250 ,1400 ,1200 ,1250 ,1400 ,1250 ,1400 ,1250 ,1400 ,1400 ,1300 ,1200 ,1200 ,1250 ,1200 ,1250 ,1250 ,1250 ,1250 ,1250 ,1300 ,1200 ,1350 ,1200 ,1200 ,1400 ,1200 ,1200 ,1200 ,1250 ,1250 ,1250 ,1250 ,1400 ,1500 ,1350 ,1250 ,NA ,1500 ,1350 ,1250 ,1200 ,NA ,1500 ,NA ,NA ,600 ,1050 ,NA ,1050 ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,45 ,1350 ,2150 ,1450 ,NA)

latitude <- unlist(argos[["latitude"]])
longitude <- unlist(argos[["longitude"]])
cycle <- unlist(argos[["cycle"]])
dataStateIndicator <- unlist(argos[["dataStateIndicator"]])
df <- data.frame(cycle=cycle, levels=levels, maxPressure=maxPressure, latitude=latitude,longitude=longitude, dataStateIndicator=dataStateIndicator)    
head(df)