It apparently uses "raw" drifter data rather then the webserved ORACLE data.
"""
from matplotlib.dates import date2num, num2date
import datetime
import pylab
import sys
import numpy as np
import matplotlib.pyplot as plt
pydir='../'
sys.path.append(pydir)
from hx import getcodar_ctl_file,getdrift_raw_range_latlon,getdrift_raw
###############################################
inputfilename='./getcodar_bydrifter_ctl.txt'
(datetime_wanted,filename,driftnumber,url,model_option,num,interval_dtime,interval,step_size)=getcodar_ctl_file(inputfilename)
id3=int(driftnumber) #change format
(maxlon,minlon,maxlat,minlat)=getdrift_raw_range_latlon(filename,id3,interval,datetime_wanted,num,step_size)
for i in range(5):
if maxlat-minlat<=0.1:
maxlat=maxlat+0.01
minlat=minlat-0.01
if maxlon-minlon<=0.1:
maxlon=maxlon+0.01
minlon=minlon-0.01
gbox=[minlon-0.03,maxlon+0.03, minlat-0.03, maxlat+0.03]
# get the edge of the picture
for x in range(num):
(lat,lon,time1)=getdrift_raw(filename,id3,interval,datetime_wanted)
"""
import sys
from matplotlib.dates import num2date
import pylab
import matplotlib.pyplot as plt
# add the path of ocean.py
#pydir='/net/home3/ocn/jmanning/py/huanxin/work'
pydir = '../'
sys.path.append(pydir)
from hx import getcodar_ctl_file, getdrift_raw_range_latlon, getcodar_ctl_id
from hx import getdrift_raw, getcodar, gettrack_codar
#### input #########################################################
inputfilename = './getcodar_bydrifter_ctl.txt'
(datetime_wanted, filename, driftnumber, url, model_option, num,
interval_dtime, interval, step_size) = getcodar_ctl_file(inputfilename)
startdate = datetime_wanted
id3 = int(driftnumber) #change format
datetime_wanted_1 = datetime_wanted
(maxlon, minlon, maxlat,
minlat) = getdrift_raw_range_latlon(filename, id3, interval,
datetime_wanted_1, num, step_size)
id = getcodar_ctl_id(model_option, url, startdate)
(lat, lon, time) = getdrift_raw(filename, id3, interval,
datetime_wanted) #drifter lat lon time
la = lat[0]
lo = lon[0]
lat_k, lon_k = [la], [lo]
uu, vv = [0], [0]
from hx import getcodar_ctl_file, getdrift_raw_range_latlon, getcodar_ctl_id
from hx import getdrift_raw, getcodar, gettrack_codar
#### input #########################################################
inputfilename = "./getcodar_bydrifter_ctl.txt"
(
datetime_wanted,
filename,
driftnumber,
url,
model_option,
num,
interval_dtime,
interval,
step_size,
) = getcodar_ctl_file(inputfilename)
startdate = datetime_wanted
id3 = int(driftnumber) # change format
datetime_wanted_1 = datetime_wanted
(maxlon, minlon, maxlat, minlat) = getdrift_raw_range_latlon(filename, id3, interval, datetime_wanted_1, num, step_size)
id = getcodar_ctl_id(model_option, url, startdate)
(lat, lon, time) = getdrift_raw(filename, id3, interval, datetime_wanted) # drifter lat lon time
la = lat[0]
lo = lon[0]
lat_k, lon_k = [la], [lo]
uu, vv = [0], [0]
for q in range(interval):
lat_vel, lon_vel, u, v = getcodar(url, id)
lon_vel_list = lon_vel[0]