data("GH")
sel<- which(GH$STNS == "CE1")
temp <- cbind(GH$JSTR[[sel[1]]][1168:1500],
GH$JSTR[[sel[2]]][1168:1500], GH$JSTR[[sel[3]]][1168:1500])
dt <- GH$dt[ sel[1] ]
STAMP <- "GH"
PMOT.drive(temp, dt,
pmolabs = c("Vertical", "North", "East"), STAMP = STAMP)
## ids <- idpoints.hodo(temp, sx, zloc$x[sn1], zloc$y[sn1])
Run the code above in your browser using DataLab