data("GH")
sel <- which(GH$STNS == "CE1")
YMAT <- cbind(GH$JSTR[[sel[1]]][1168:1500],
GH$JSTR[[sel[2]]][1168:1500],
GH$JSTR[[sel[3]]][1168:1500])
dt <- GH$dt[ sel[1] ]
ftime <- Zdate(GH$info, sel[1], 1)
if(interactive()){
PMOT.drive(YMAT, dt, pmolabs = c("Vertical", "North", "East"),
STAMP =ftime )
}
Run the code above in your browser using DataLab