update URL

Author: Jeffrey Whitaker

doc/users/figures/plotargo.py

@@ -1,28 +1,24 @@
-from netCDF4 import Dataset
+from netCDF4 import Dataset, num2date
 import time, calendar, datetime, numpy
 from mpl_toolkits.basemap import Basemap
 import matplotlib.pyplot as plt
-def datetomsecs(d):
-    """convert from datetime to msecs since the unix epoch began"""
-    return int(calendar.timegm(time.struct_time(d.timetuple()))*1000)
-# set date range
-date1 = datetime.datetime(2010,1,1,0)
-date2 = datetime.datetime(2010,1,8,0)
-t1 = datetomsecs(date1); t2 = datetomsecs(date2)
-# build constraint expression to get locations of floats in specified time
-# range.
-urlbase='http://dapper.pmel.noaa.gov/dapper/argo/argo_all.cdp'
-sel="?location.JULD,location.LATITUDE,location.LONGITUDE&location.JULD>%s&location.JULD<%s"%(t1,t2)
-# retrieve data
-dset = Dataset(urlbase+sel)
-lats = dset.variables['location.LATITUDE'][:]
-lons = dset.variables['location.LONGITUDE'][:]
+# data downloaded from
+# http://coastwatch.pfeg.noaa.gov/erddap/tabledap/apdrcArgoAll.html via
+# http://coastwatch.pfeg.noaa.gov/erddap/tabledap/apdrcArgoAll.nc?longitude,latitude,time&longitude>=0&longitude<=360&latitude>=-90&latitude<=90&time>=2010-01-01&time<=2010-01-08&distinct()
+dset = Dataset('apdrcArgoAll_af28_ee2b_1d10.nc')
+lats = dset.variables['latitude'][:]
+lons = dset.variables['longitude'][:]
+time = dset.variables['time']
+times = time[:]
+t1 = times.min(); t2 = times.max()
+date1 = num2date(t1, units=time.units)
+date2 = num2date(t2, units=time.units)
 # draw map with markers for float locations
 m = Basemap(projection='hammer',lon_0=180)
 x, y = m(lons,lats)
 m.drawmapboundary(fill_color='#99ffff')
 m.fillcontinents(color='#cc9966',lake_color='#99ffff')
 m.scatter(x,y,3,marker='o',color='k')
 plt.title('Locations of %s ARGO floats active between %s and %s' %\
-        (len(lats),date1.strftime('%Y%m%d'),date2.strftime('%Y%m%d')))
+        (len(lats),date1,date2),fontsize=12)
 plt.show()