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]