Notebook

Looking at the drifter data from Mark. Trying to determine dates and locations of releases.

In [53]:

%matplotlib inline
from matplotlib import pylab
import matplotlib.pyplot as plt
import netCDF4 as nc
import numpy as np
import scipy.io
import datetime as dt

from salishsea_tools import nc_tools, viz_tools, tidetools, stormtools, bathy_tools

from __future__ import division

In [2]:

drifters = scipy.io.loadmat('/ocean/mhalvers/research/drifters/SoG_drifters.mat',squeeze_me=True)
ubc = drifters['ubc']
ios=drifters['ios']

In [3]:

grid = nc.Dataset('/data/nsoontie/MEOPAR/NEMO-forcing/grid/bathy_meter_SalishSea2.nc','r')
bathy = grid.variables['Bathymetry'][:, :]
X = grid.variables['nav_lon'][:, :]
Y = grid.variables['nav_lat'][:, :]

Retrieving deploy times and positions¶

In [4]:

def convert_time(matlab_time_array):
    "converts a matlab time array to python format"
    python_time_array=[]
    for t in matlab_time_array:
        python_datetime = dt.datetime.fromordinal(int(t)) + dt.timedelta(days=t%1) - dt.timedelta(days = 366)
        python_time_array.append(python_datetime)
    
    python_time_array = np.array(python_time_array)
    return python_time_array

In [5]:

def get_tracks(switch,lats,lons,ptime,in_water):
    """returns a list of tracks of each buoy, ie a trajectory for each time the buoy was released into the water"""
    all_tracks=[]
    for ind in switch:
        track_on = 1
        i = ind
        track ={'time':[], 'lat':[],'lon':[]}
        while(track_on):
            if in_water[i]!=1:
                track_on=0
            elif i==np.shape(in_water)[0]-1:    
                track['time'].append(ptime[i])
                track['lat'].append(lats[i])
                track['lon'].append(lons[i])
                track_on=0
            else:
                track['time'].append(ptime[i])
                track['lat'].append(lats[i])
                track['lon'].append(lons[i])
            i=i+1
        all_tracks.append(track)
        
    return all_tracks

In [6]:

def organize_info(buoy,btype):
    """ organizes the buoy info. Groups the buoy data into tracks for when it was released into the water. """
    #creat arrays for easier access
    buoy_name = btype[buoy][0]
    lats = btype[buoy]['lat'].flatten()
    lons = btype[buoy]['lon'].flatten()
    mtime = btype[buoy]['mtime']
    in_water = btype[buoy]['isSub'].flatten()
    #convert mtime to python datetimes
    ptime = convert_time(mtime)
    
    #loop through in_water flag to find when buoy switched from being out of water to being in water. 
    switch = []; 
    for ind in np.arange(1,in_water.shape[0]):
        if int(in_water[ind]) != int(in_water[ind-1]):
            if int(in_water[ind])==1:
                switch.append(ind)
    
    all_tracks=get_tracks(switch,lats,lons,ptime.flatten(),in_water)
    
    return buoy_name, all_tracks

In [7]:

def print_info(buoy,btype):
    """ prints the release time, lat, lon, and duration of a buoy track"""
    name, tracks=organize_info(buoy,btype)
    print name
    print 'Release times, positions and duration in hours'
    for t in tracks:
        print t['time'][0], t['lat'][0], t['lon'][0], (t['time'][-1]-t['time'][0]).total_seconds()/3600

In [8]:

def find_start(tracks, start_date):
    """returns the a list of indices for a track released on start date.
    Only checks the month and day of the start day"""
    i=0
    ind=[]
    starttimes=[]
    for t in tracks:   
        if int(t['time'][0].month) == start_date.month:
                if int(t['time'][0].day) == start_date.day:
                    ind.append(i)
        i=i+1
        
    return ind

In [72]:

def plot_buoy(ax, tracks, startdate, i=0, fancy=False):
    """ plots a buoy trajectory at the given startdate in an axis, ax.
    returns the trajectory that was plotted.
    The first track released on the startdate is plotted.
    For trajectories that were released mulitples times a day, i selects which release is plotted.
    """
    ind =find_start(tracks,startdate)
    traj=tracks[ind[i]]
    duration = (traj['time'][-1]-traj['time'][0]).total_seconds()/3600
    print 'Released', traj['time'][0], 'at', traj['lat'][0], ',' , traj['lon'][0], 'for' , duration, 'hours'
    ax.plot(traj['lon'],traj['lat'],'ob')
    ax.plot(traj['lon'][0],traj['lat'][0],'sr')
    
    [j,i]=tidetools.find_closest_model_point(float(traj['lon'][0]),float(traj['lat'][0]),X,Y,bathy)
    ax.plot(-123-np.array([18.2, 13.7, 12])/60.,49+np.array([6.4, 8, 7.6])/60.,'-k',lw=2); 
    if fancy:
        cmap = plt.get_cmap('winter_r')
        cmap.set_bad('burlywood')
        ax.pcolormesh(X, Y, bathy, cmap=cmap)
        ax.set_title('Observed Drift Track')
        ax.set_xlabel('Longitude')
        ax.set_ylabel('Latitude')
        ax.text(-123.15,49.13, "Fraser River", fontsize=12)
    else:
        viz_tools.plot_coastline(ax, grid, coords='map')
        viz_tools.plot_coastline(ax, grid, coords='map',isobath=4)
        viz_tools.plot_coastline(ax, grid, coords='map',isobath=20)
        print 'NEMO coords:', j,i
    return traj

In [10]:

def plot_tides(ax,relase_times,start,end):
    """Plots the the tides at Point Atkinson between start and end , highlighting the release times 
    listed in relase_times."""
    wlev = stormtools.load_observations(start, end,'PointAtkinson')
    ax.plot(wlev.time,wlev.slev,'-k')
    for t in release_times:
        ax.plot([t,t],[0,5],'-')
    ax.set_xlim(dt.datetime.strptime(start,'%d-%b-%Y'),dt.datetime.strptime(end,'%d-%b-%Y'))
    ax.set_ylim([0,5])
    ax.set_ylabel('Water level (m)')
    ax.set_xlabel('time (UTC)')
    ax.set_title('Observed water level at Point Atkinson ' + start +'-'+ end)

In [11]:

def plot_winds(ax_speed,ax_dir,release_times,start,end):
    """Plots Sandheads winds between start and end. Highlights release times listed in release_times"""
    [wind_speed,wind_dir,temp,time, lat, lon] = stormtools.get_EC_observations('Sandheads',start,end)
    ax=ax_speed
    ax.plot(time,wind_speed,'-k')
    for t in release_times:
        ax.plot([t,t],[0,10],'-')
    ax.set_xlim(dt.datetime.strptime(start,'%d-%b-%Y'),dt.datetime.strptime(end,'%d-%b-%Y'))
    ax.set_ylim([0,10])
    ax.set_ylabel('Wind Speed (m/s)')
    ax.set_title('Observed wind speed at Sandheads ' + start +'-'+ end)
    ax=ax_dir
    ax.plot(time,wind_dir,'-k')
    for t in release_times:
        ax.plot([t,t],[0,360],'-')
    ax.set_xlim(dt.datetime.strptime(start,'%d-%b-%Y'),dt.datetime.strptime(end,'%d-%b-%Y'))
    ax.set_ylim([0,350])
    ax.set_ylabel('Wind Direction (deg CCW of E)')
    ax.set_title('Observed wind direction at Sandheads ' + start +'-'+ end)

UBC Buoys¶

First, print some information about each buoy.

In [12]:

for num in np.arange(0,ubc.shape[0]):
    print_info(num,ubc)

UBC-I-0001
Release times, positions and duration in hours
2014-08-21 16:40:04.999997 49.1632016667 -123.358991667 0.667500000556
2014-08-21 17:50:08.000002 49.1600033333 -123.365111667 0.833333334167
2014-08-21 19:05:07.999996 49.1571733333 -123.359418333 0.999722222778
2014-08-21 20:45:06.999996 49.1571983333 -123.360888333 1.00138888861
2014-08-21 22:00:08.999998 49.1548316667 -123.366338333 1.00000000194
2014-08-28 16:30:05.000003 49.2163283333 -123.355693333 0.917499999444
2014-08-28 18:10:08.999996 49.2151716667 -123.364878333 1.17055555722
2014-08-28 19:40:08.000001 49.2150016667 -123.363986667 1.16638888806
2014-08-28 21:15:08.000004 49.2174433333 -123.365893333 0.583333331944
2014-09-04 16:30:07.999997 49.1352766667 -123.333071667 0.666666668889
2014-09-04 17:40:07.999998 49.1345266667 -123.337721667 0.75
2014-09-04 19:00:10.999997 49.1348083333 -123.338661667 0.832499999722
2014-09-04 20:10:08.000005 49.1353116667 -123.337346667 0.583333331944
2014-09-04 21:05:07.000005 49.1409816667 -123.326211667 0.167499999444
2014-09-04 22:00:08.999998 49.1339483333 -123.336955 0.416666666944
2014-09-04 22:50:07.999996 49.1376616667 -123.334643333 1.00000000167
2014-09-11 17:15:07.999997 49.1259933333 -123.432688333 0.666944445
2014-09-11 18:05:07.999999 49.1229316667 -123.429991667 0.416666666944
2014-09-11 18:40:08.000005 49.1230366667 -123.430511667 0.924722221667
2014-09-11 19:50:07.999996 49.122965 -123.430361667 1.08500000167
2014-09-11 21:10:15.999998 49.1307516667 -123.43125 0.0
2014-09-11 21:30:11.000004 49.1303116667 -123.423141667 0.832499999444
2014-09-19 10:20:08.999997 49.079315 -123.442168333 55.8333333331
2014-10-03 17:00:06.999996 49.2213966667 -123.360713333 1.58361111278
2014-10-03 19:05:07.999996 49.1345333333 -123.333225 0.50027778
2014-10-03 19:45:12.000002 49.1331933333 -123.341221667 0.915277778056
2014-10-03 20:55:07 49.132715 -123.342146667 1.75027777806
2014-10-08 16:00:06.999999 49.0862766667 -123.320955 26.5002777781
UBC-I-0002
Release times, positions and duration in hours
2014-08-21 16:50:05.999997 49.1624266667 -123.360091667 0.667222224167
2014-08-21 17:45:09.000005 49.1591366667 -123.364311667 0.833611110278
2014-08-21 19:05:07.999996 49.1566833333 -123.364288333 1.08305555667
2014-08-21 20:40:07.000003 49.1589366667 -123.364911667 0.916666665278
2014-08-21 22:10:07.000003 49.1599833333 -123.368518333 0.916666665278
2014-08-28 16:35:07.999999 49.2204116667 -123.35495 0.916388888889
2014-08-28 17:50:10.000002 49.215015 -123.364911667 1.3358333325
2014-08-28 19:40:06.999997 49.2149183333 -123.36483 1.16694444583
2014-08-28 21:15:08.000004 49.2166983333 -123.365858333 0.583333331944
2014-09-04 16:35:07.000005 49.1348366667 -123.331138333 0.499999998056
2014-09-04 17:40:07.999998 49.1329933333 -123.338713333 0.666666666111
2014-09-04 19:05:07.999996 49.1372083333 -123.3351 0.75
2014-09-04 20:15:07.999997 49.1362216667 -123.334896667 0.500000000833
2014-09-04 21:10:10.000001 49.13669 -123.335393333 0.583055555556
2014-09-04 22:00:36.999997 49.1347366667 -123.336283333 0.408888889444
2014-09-04 22:50:09 49.136615 -123.335816667 0.66638889
2014-09-11 17:15:07.999997 49.1268333333 -123.43173 0.666666668889
2014-09-11 18:05:07.999999 49.1216683333 -123.431543333 0.416666666944
2014-09-11 18:40:08.000005 49.1242266667 -123.429055 0.924722221667
2014-09-11 19:50:09 49.1238733333 -123.429336667 1.08305555694
2014-09-11 21:35:06.999998 49.1295483333 -123.424953333 0.750277778889
2014-09-19 10:20:08.999997 49.0794016667 -123.44245 55.6663888889
2014-10-03 17:00:06.999996 49.2223983333 -123.360571667 0.500000000833
2014-10-03 17:40:11.000002 49.2170533333 -123.366633333 0.9991666675
2014-10-03 18:50:07.000005 49.1879816667 -123.346305 0.0
2014-10-03 19:10:07.999998 49.1353483333 -123.330375 0.250000001944
2014-10-03 19:45:12.000002 49.1333333333 -123.340268333 0.832222220556
2014-10-03 20:55:08.000005 49.1335616667 -123.340743333 1.58333333111
2014-10-08 15:55:08.999996 49.0881183333 -123.320238333 0.0
2014-10-08 16:05:07.999996 49.0956433333 -123.320643333 26.4997222228
UBC-I-0003
Release times, positions and duration in hours
2014-08-21 17:00:09.999999 49.1622083333 -123.359383333 0.499166666667
2014-08-21 17:45:06.999996 49.1588283333 -123.362636667 0.833611112778
2014-08-21 19:10:07.000003 49.1571216667 -123.364928333 1.08388888833
2014-08-21 20:40:09.000003 49.1595783333 -123.365821667 0.916388888611
2014-08-21 22:10:07.999998 49.1610016667 -123.370008333 0.917222223333
2014-08-28 16:35:07.999999 49.2219966667 -123.354388333 0.999999999167
2014-08-28 17:50:10.000002 49.2151733333 -123.361788333 0.749444444444
2014-08-28 18:45:07.999997 49.2148183333 -123.363811667 1.41638888972
2014-08-28 20:30:08.999998 49.2153583333 -123.366746667 0.916388888889
2014-08-28 21:35:10.999996 49.2190833333 -123.36525 0.499166666667
2014-09-04 16:40:08.000001 49.1340483333 -123.331875 0.333333333056
2014-09-04 17:40:09.000003 49.1337966667 -123.338205 0.666388886944
2014-09-04 19:05:07.999996 49.1365116667 -123.335713333 0.749722223611
2014-09-04 20:10:08.000005 49.13389 -123.338846667 0.608333333333
2014-09-04 21:05:07.000005 49.1403983333 -123.327588333 0.250555554167
2014-09-04 22:00:08.000004 49.1329283333 -123.33708 0.333611109167
2014-09-04 22:50:09.999995 49.13869 -123.33338 0.9161111125
2014-09-11 17:10:08.000005 49.1304433333 -123.429628333 0.583333331944
2014-09-11 18:10:08.000001 49.1275983333 -123.42535 1.25388888833
2014-09-11 19:55:07.999998 49.1273083333 -123.425863333 1.25027777722
2014-09-11 21:35:06.999998 49.1285616667 -123.425955 0.833611113056
2014-09-19 10:20:08.999997 49.0793866667 -123.442943333 11.1675000006
2014-09-20 04:40:08.999996 49.0754766667 -123.439433333 12.5005555561
2014-09-21 10:00:08.000004 49.0190616667 -123.490253333 9.16749999944
2014-10-03 17:00:06.999996 49.2192716667 -123.360443333 0.50027778
2014-10-03 17:45:08 49.2188683333 -123.366166667 0.916666668056
2014-10-03 19:10:09.000003 49.1360583333 -123.329678333 0.249999998889
2014-10-03 19:45:12.999997 49.13374 -123.339655 0.831944444167
2014-10-03 20:55:07 49.1342816667 -123.339538333 1.0011111125
2014-10-03 22:05:12.999999 49.13534 -123.330628333 0.415277777222
2014-10-08 16:10:07.000003 49.1054983333 -123.319833333 28.7502777783
UBC-I-0004
Release times, positions and duration in hours
2014-08-21 16:50:07.999996 49.1632983333 -123.357691667 0.500000000833
2014-08-21 17:45:08 49.1591483333 -123.365163333 0.916666668056
2014-08-21 19:05:07.999996 49.1568983333 -123.362653333 1.08861111278
2014-08-21 20:40:09.000003 49.1582316667 -123.364028333 1.00027777806
2014-08-21 22:05:09 49.1585783333 -123.369058333 0.999722222778
2014-08-28 16:30:09.000001 49.2173583333 -123.355303333 1.33305555583
2014-08-28 18:10:06.999997 49.2152433333 -123.363493333 2.16750000028
2014-08-28 20:30:08.000004 49.2150966667 -123.363843333 0.833055554722
2014-08-28 21:35:08.000002 49.2209366667 -123.365721667 0.499722221944
2014-09-04 16:40:06.999997 49.1333083333 -123.332743333 0.333611111944
2014-09-04 17:40:09.000003 49.1352933333 -123.337271667 0.749722220833
2014-09-04 19:05:07.000001 49.1357216667 -123.336458333 0.750277776389
2014-09-04 20:10:08.000005 49.1316533333 -123.334875 0.583611111111
2014-09-04 21:05:07.999999 49.1400516667 -123.3292 0.333333333056
2014-09-04 22:05:07.999996 49.1376566667 -123.333716667 0.50027778
2014-09-04 22:45:10.999997 49.1338033333 -123.339186667 0.666111110833
2014-09-11 17:10:10.000004 49.1285016667 -123.430555 0.666111110556
2014-09-11 18:05:07.000005 49.1241466667 -123.428625 1.08361110917
2014-09-11 19:55:07.999998 49.1306816667 -123.421511667 1.33361111111
2014-09-11 21:35:08.999997 49.1274366667 -123.426991667 0.749722223611
2014-09-19 11:40:07.999998 49.0783233333 -123.441393333 9.83361111306
2014-09-20 04:40:08.000001 49.0755083333 -123.439546667 12.5000000011
2014-09-21 09:30:08 49.0316116667 -123.448503333 9.16722222333
2014-10-03 17:00:08 49.2203283333 -123.36073 0.499722221944
2014-10-03 17:45:06.999996 49.2204866667 -123.365993333 0.916666667778
2014-10-03 19:10:07.000003 49.1369116667 -123.329105 0.333333333056
2014-10-03 19:45:06.999999 49.1343766667 -123.338953333 0.833611110278
2014-10-03 20:55:08.000005 49.1348533333 -123.338363333 0.920833331389
2014-10-03 22:05:07.000001 49.13484 -123.332016667 0.500555556389
2014-10-08 17:30:08.000004 49.1064116667 -123.320053333 28.2499999981
UBC-I-0005
Release times, positions and duration in hours
2014-08-21 17:45:08 49.1589283333 -123.363593333 0.833055555
2014-08-21 19:05:07.000001 49.1569216667 -123.361086667 1.08333333306
2014-08-21 20:45:06.999996 49.1577733333 -123.361813333 0.916666667778
2014-08-21 22:05:07.000001 49.1574766667 -123.367741667 1.00194444417
2014-08-28 16:35:07.000005 49.2212016667 -123.3547 0.99972222
2014-08-28 17:50:08.999997 49.2151183333 -123.363361667 1.24972222194
2014-08-28 19:40:06.999997 49.2148216667 -123.365558333 1.16666666694
2014-08-28 21:15:06.999999 49.2158466667 -123.365718333 0.500000000833
2014-09-04 16:40:08.999996 49.13559 -123.330558333 0.416666666944
2014-09-04 17:40:07.999998 49.1323933333 -123.339593333 0.666666666111
2014-09-04 19:05:09 49.1379966667 -123.33448 1.583055555
2014-09-04 21:05:09.000004 49.1396683333 -123.330903333 0.416666666944
2014-09-04 22:00:08.000004 49.1373533333 -123.335253333 0.583611111111
2014-09-04 22:45:08.000004 49.13484 -123.338383333 0.666944442222
2014-09-11 17:10:08.999999 49.1276283333 -123.431205 0.833333333889
2014-09-11 18:10:08.000001 49.1263283333 -123.426596667 1.1672222225
2014-09-11 19:55:07.999998 49.1296266667 -123.422841667 1.333333335
2014-09-11 21:35:08.000002 49.1262216667 -123.428261667 0.833333334167
2014-09-19 11:40:07.999998 49.0782866667 -123.441641667 53.6747222236
2014-10-03 17:15:07.999997 49.2156316667 -123.365708333 1.33416666639
2014-10-03 19:45:06.999999 49.1349333333 -123.338005 0.833333333889
2014-10-03 20:55:07 49.135515 -123.336993333 1.91666666694
2014-10-08 17:30:06.999999 49.0965066667 -123.324513333 26.166944445
UBC-I-0006
Release times, positions and duration in hours
2014-08-21 17:55:05.000001 49.16164 -123.364875 0.834166665556
2014-08-21 19:00:08.999998 49.1567816667 -123.358705 1.08305555667
2014-08-21 20:45:09.999999 49.1566783333 -123.359868333 0.9991666675
2014-08-21 22:05:07.000001 49.1564266667 -123.366736667 0.916944444167
2014-08-28 16:35:07.000005 49.2196466667 -123.355246667 0.916944444167
2014-08-28 18:05:07.999999 49.2147783333 -123.360738333 2.083333335
2014-08-28 20:30:10.000003 49.2163833333 -123.366813333 0.93222222
2014-08-28 21:35:08.000002 49.2199966667 -123.365391667 0.499999998333
2014-09-04 16:35:07.999999 49.1341433333 -123.328878333 0.500555556389
2014-09-04 17:40:07.999998 49.1315083333 -123.339128333 2.25111111056
2014-09-04 20:10:08.999999 49.1321233333 -123.340621667 0.500000001111
2014-09-04 21:10:08.999996 49.1390016667 -123.3307 0.507777778333
2014-09-04 22:00:06.999999 49.1321733333 -123.337371667 0.333333333056
2014-09-04 22:55:07.999998 49.1385833333 -123.332368333 0.833055555
2014-09-11 17:10:07 49.1321216667 -123.428711667 0.583611111111
2014-09-11 18:10:06.999997 49.1299783333 -123.422863333 1.33472222167
2014-09-11 19:50:09 49.1263466667 -123.426193333 1.24972222194
2014-09-11 21:35:08.000002 49.12504 -123.42943 0.916666665278
2014-09-19 11:40:07.999998 49.0780116667 -123.441743333 9.666944445
2014-09-20 11:20:07.999999 49.0793333333 -123.4434 6.83361111028
2014-09-20 19:31:06.999999 49.0429283333 -123.316041667 0.150833333889
2014-09-20 20:20:07.999999 49.0969733333 -123.312955 21.5080555547
2014-10-03 17:10:05.999995 49.2169466667 -123.365928333 1.41777777972
2014-10-03 19:45:12.000002 49.1357733333 -123.337175 0.832222220556
2014-10-03 20:50:08.999997 49.1359183333 -123.337066667 0.666388889722
2014-10-03 21:40:08.999999 49.1316266667 -123.340063333 1.16611111167
2014-10-08 17:35:07.000001 49.08599 -123.323693333 23.8333333322
UBC-I-0007
Release times, positions and duration in hours
2014-08-21 19:15:06.000001 49.1581933333 -123.365513333 1.08361111222
2014-08-21 20:35:07.000001 49.1607333333 -123.3686 0.917222223333
2014-08-21 22:15:09.999999 49.1635683333 -123.370778333 0.832777778611
2014-08-28 16:40:06.999997 49.2233083333 -123.354386667 0.916666667778
2014-08-28 17:45:06.999996 49.2145266667 -123.360768333 0.58388889
2014-08-28 18:50:09.000004 49.2151983333 -123.36498 1.41638888694
2014-08-28 20:30:08.000004 49.2149683333 -123.365443333 0.835277776111
2014-08-28 21:35:34.000002 49.2218966667 -123.36585 0.493055555
2014-09-04 16:35:07.999999 49.1345116667 -123.330011667 0.583611111111
2014-09-04 17:40:11.000002 49.1367933333 -123.336155 0.832777778611
2014-09-04 19:00:06.999999 49.1324983333 -123.340858333 0.750277779167
2014-09-04 20:15:07.999997 49.1384466667 -123.332941667 0.583055555833
2014-09-04 21:10:08.000001 49.13849 -123.332313333 0.500000001111
2014-09-04 22:00:06.999999 49.1365816667 -123.335455 0.416944446111
2014-09-04 22:50:07.999996 49.1348233333 -123.337636667 0.583333334722
2014-09-11 17:10:07 49.1313433333 -123.429128333 0.583611111111
2014-09-11 18:10:06.999997 49.1287316667 -123.42408 1.25027778
2014-09-11 19:55:07.999998 49.1284333333 -123.424493333 1.33361111111
2014-09-11 21:40:08.000005 49.1243566667 -123.431325 0.833333331111
2014-09-19 13:10:07.999998 49.0804766667 -123.44053 8.00000000083
2014-09-20 11:10:08.000005 49.0787216667 -123.441358333 6.83361111
2014-09-21 09:10:06.999997 49.0447116667 -123.403411667 9.00027777889
2014-10-03 17:10:05.999995 49.2181716667 -123.366143333 1.33361111139
2014-10-03 19:45:10.000003 49.1363483333 -123.336246667 0.749444444722
2014-10-03 20:50:08.999997 49.1366116667 -123.335778333 0.583055555833
2014-10-03 21:40:08.000005 49.1324233333 -123.341191667 1.16694444306
2014-10-08 19:10:07.999998 49.0856266667 -123.325655 23.0836111111
UBC-I-0008
Release times, positions and duration in hours
2014-08-21 20:45:04.999996 49.1627533333 -123.366086667 0.750555555556
2014-08-21 22:15:08 49.1649466667 -123.372203333 0.916388888889
2014-08-28 16:35:07.999999 49.2188183333 -123.355413333 0.916944444167
2014-08-28 18:05:09.000004 49.2147566667 -123.36228 1.16694444306
2014-08-28 19:40:08.000001 49.21509 -123.363311667 1.16666666694
2014-08-28 21:15:08.000004 49.2182566667 -123.365993333 0.583333331944
2014-09-04 16:30:07.000002 49.13564 -123.334116667 0.75
2014-09-04 17:45:06.999996 49.1360516667 -123.334478333 0.666944445
2014-09-04 19:00:06.999999 49.1333016667 -123.340025 0.75
2014-09-04 20:15:07.999997 49.1377016667 -123.333703333 0.586111111389
2014-09-04 21:10:06.999997 49.1375366667 -123.334003333 0.583333334722
2014-09-04 22:00:55 49.1356816667 -123.335641667 0.4036111125
2014-09-04 22:50:36.999999 49.1355866667 -123.336808333 0.658333333333
2014-09-11 17:10:08.999999 49.1294933333 -123.430183333 0.666388889722
2014-09-11 18:05:07.000005 49.1253516667 -123.427353333 1.25055555333
2014-09-11 19:55:07.000003 49.1317516667 -123.420008333 1.33361111139
2014-09-11 21:40:07 49.1231466667 -123.432493333 0.833611110278
2014-09-19 13:10:07.999998 49.0803516667 -123.440405 8.00000000083
2014-09-20 13:20:08.999997 49.0768533333 -123.44025 5.33305555778
2014-09-21 08:50:08.999997 49.0574933333 -123.354403333 9.83361111306
2014-10-03 17:10:07 49.2190766667 -123.366333333 1.25000000083
2014-10-03 19:45:10.999997 49.1369833333 -123.335355 0.749166668611
2014-10-03 20:50:08.000002 49.1372966667 -123.334508333 0.583055555833
2014-10-03 21:40:07 49.1331266667 -123.340413333 1.00027777806
2014-10-08 19:20:08.999997 49.0943266667 -123.330828333 24.9163888889
UBC-I-0009
Release times, positions and duration in hours
2014-08-21 19:15:06.000001 49.15774 -123.367186667 1.08388888833
2014-08-21 20:40:07.000003 49.1604783333 -123.366358333 0.916666665278
2014-08-21 22:10:07.999998 49.1622216667 -123.371421667 0.916666667778
2014-08-28 16:30:05.000003 49.2155716667 -123.356013333 0.917499999444
2014-08-28 18:10:06.999997 49.2150966667 -123.365536667 1.26777778
2014-08-28 19:35:07.999999 49.2147683333 -123.362918333 1.16694444611
2014-08-28 21:15:08.000004 49.2189766667 -123.366166667 0.583333331944
2014-08-28 22:15:11.000004 49.2421083333 -123.31863 0.0
2014-09-04 16:40:06.999997 49.1363566667 -123.329793333 0.583333334722
2014-09-04 17:45:08 49.137565 -123.333421667 0.75
2014-09-04 19:01:00.999998 49.1339983333 -123.339113333 0.735277779444
2014-09-04 20:15:07.999997 49.1369683333 -123.334278333 0.583888890278
2014-09-04 21:15:08.000004 49.1356183333 -123.334686667 0.502499998333
2014-09-04 21:55:10.000001 49.1319583333 -123.340668333 0.333055556667
2014-09-04 22:55:07.999998 49.1392466667 -123.33133 0.833611113056
2014-09-11 17:15:07.999997 49.1252316667 -123.432993333 1.16666666694
2014-09-11 18:40:10.000004 49.1254533333 -123.427733333 0.924166666389
2014-09-11 19:50:07.999996 49.1251833333 -123.427653333 1.25055555611
2014-09-11 21:40:10.000004 49.1219566667 -123.433483333 0.832499999444
2014-09-19 13:10:09.000003 49.0802733333 -123.440366667 53.6663888878
UBC-I-0010
Release times, positions and duration in hours
2014-09-04 23:35:09 49.1264783333 -123.340441667 0.749722221111
2014-09-20 13:20:08.000002 49.0767216667 -123.440253333 5.33333333417
2014-09-21 06:30:09.000005 49.0710983333 -123.311388333 12.3338888883
2014-10-03 17:10:05.000001 49.2219533333 -123.366071667 1.33444444556
2014-10-03 19:50:07.999996 49.1378466667 -123.332918333 0.666666668889
2014-10-03 20:50:06.999998 49.1381583333 -123.333038333 0.584444445556
2014-10-03 21:40:08.999999 49.13374 -123.339453333 0.999722222778
2014-10-08 19:20:06.999998 49.1071733333 -123.319563333 6.66666666611

IOS Buoys¶

In [13]:

for num in np.arange(0,ios.shape[0]):
    print_info(num,ios)

51
Release times, positions and duration in hours
2014-09-11 18:01:32.999998 49.12199 -123.43108 0.250277778333
2014-09-11 20:29:03.999999 49.12461 -123.43199 1.24638888694
2014-10-03 16:46:53.999992 49.22194 -123.3623 0.481944445556
2014-10-03 19:05:22.999995 49.13684 -123.33138 3.77611111194
2014-10-08 15:53:47.999993 49.08842 -123.32025 26.1491666686
52
Release times, positions and duration in hours
2014-09-11 18:03:46.999996 49.12542 -123.42722 0.0833333338889
2014-09-11 19:50:01.999998 49.12751 -123.42473 0.0827777786111
2014-09-11 21:34:42.999997 49.12605 -123.42802 0.830555554444
2014-10-03 17:02:04.999999 49.21673 -123.3669 2.34111111139
2014-10-03 19:40:34.000001 49.13378 -123.34093 3.205
2014-10-08 15:59:00 49.0972 -123.32001 24.2474999986
53
Release times, positions and duration in hours
2014-09-11 18:08:13.000001 49.1288 -123.42386 531.681111111
2014-10-08 16:08:57.999995 49.10592 -123.31999 48.7591666672
54
Release times, positions and duration in hours
2014-10-05 19:52:30.999995 48.78801 -122.71062 85.4858333347
55
Release times, positions and duration in hours
2014-09-19 12:28:04.999993 49.08904 -123.48761 9.46527778
2014-09-20 13:11:35.999993 49.07387 -123.46614 483.087500001
56
Release times, positions and duration in hours
2014-09-19 12:16:58.999995 49.08163 -123.45322 425.832222221
57
Release times, positions and duration in hours
58
Release times, positions and duration in hours
2014-10-08 19:11:23.999998 49.09666 -123.32301 23.9202777772
59
Release times, positions and duration in hours
2014-09-19 13:02:20.000002 49.08028 -123.43953 9.53972222
2014-09-20 04:30:43.999993 49.07602 -123.439 399.183888889
60
Release times, positions and duration in hours
2014-09-19 13:26:04.999997 49.08292 -123.44423 74.8572222228
61
Release times, positions and duration in hours
2014-09-19 10:56:11.999993 49.08078 -123.45824 152.501111113
62
Release times, positions and duration in hours
2014-09-19 19:44:05.999998 49.18535 -123.53365 178.849166667
2014-10-02 23:47:43.999999 49.33023 -123.60278 186.079722223

Separating into tracks based on the "isSub" flag is not going to work for the IOS drifters. I will analyze these tracks in a separate notebook later

September 19¶

Buoy 1 - Drop 1¶

In [14]:

buoy=0
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0001
Released 2014-09-19 10:20:08.999997 at 49.079315 , -123.442168333 for 55.8333333331 hours
NEMO coords: 432 269

Out[14]:

(48.8, 49.4)

Buoy 2 - Drop 1¶

On Sept 19, Buoy 2 was released with Buoy 1. Same starting time and position and duration. So it should do the same thing as first buoy.

In [15]:

buoy=1
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0002
Released 2014-09-19 10:20:08.999997 at 49.0794016667 , -123.44245 for 55.6663888889 hours
NEMO coords: 432 269

Out[15]:

(48.8, 49.4)

Buoy 3 - Drop 1¶

On Sept 19, Buoy 3 was released with Buoy 1 and 2 but then was removed after 11 hours. It was also released on Sept 20 and 21.

In [16]:

buoy=2
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0003
Released 2014-09-19 10:20:08.999997 at 49.0793866667 , -123.442943333 for 11.1675000006 hours
NEMO coords: 432 269

Out[16]:

(48.8, 49.4)

Buoy 4 - Drop 2¶

Buoy 4: was released on Sept 19 at 11:40am for 10 hours at pretty much the same location as Buoy 1,2,3

In [17]:

buoy=3
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0004
Released 2014-09-19 11:40:07.999998 at 49.0783233333 , -123.441393333 for 9.83361111306 hours
NEMO coords: 432 269

Out[17]:

(48.8, 49.4)

Buoy 5 - Drop 2¶

Released on Sept 19 at 11:40 for 54 hours. Same drop point as above

In [18]:

buoy=4
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0005
Released 2014-09-19 11:40:07.999998 at 49.0782866667 , -123.441641667 for 53.6747222236 hours
NEMO coords: 432 269

Out[18]:

(48.8, 49.4)

Buoy 6 - Drop 2¶

Sept 19 at 11:40am, same location for about 10 hours

In [19]:

buoy=5
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0006
Released 2014-09-19 11:40:07.999998 at 49.0780116667 , -123.441743333 for 9.666944445 hours
NEMO coords: 432 269

Out[19]:

(48.8, 49.4)

Buoy 7 - Drop 3¶

Sept 19 at 13:10 for 8 hours. Same drop location

In [20]:

buoy=6
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0007
Released 2014-09-19 13:10:07.999998 at 49.0804766667 , -123.44053 for 8.00000000083 hours
NEMO coords: 433 269

Out[20]:

(48.8, 49.4)

Buoy 8 - Drop 3¶

In [21]:

buoy=7
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0008
Released 2014-09-19 13:10:07.999998 at 49.0803516667 , -123.440405 for 8.00000000083 hours
NEMO coords: 433 269

Out[21]:

(48.8, 49.4)

Buoy 9 - Drop 3¶

Sept 19 at 13:10 for 54 hours. Similar drop location

In [22]:

buoy=8
name, tracks=organize_info(buoy,ubc)
fig,ax = plt.subplots(1,1,figsize=(5,5))
print name
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,19))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0009
Released 2014-09-19 13:10:09.000003 at 49.0802733333 , -123.440366667 for 53.6663888878 hours
NEMO coords: 433 269

Out[22]:

(48.8, 49.4)

Tides¶

Look quickly at the tides on Sept 19.

In [23]:

fig,ax=plt.subplots(1,1,figsize=(10,4))

t1=dt.datetime(2014,9,19,10,20)
t2=dt.datetime(2014,9,19,11,40)
t3=dt.datetime(2014,9,19,13,10)
release_times=[t1,t2,t3];
plot_tides(ax,release_times,'18-Sep-2014', '22-Sep-2014')

So drifters where relased during an ebb tide? Check with Rich. Why did they initially go north west then? Winds?

Winds¶

In [24]:

fig,axs = plt.subplots(2,1,figsize=(10,7))
plot_winds(axs[0],axs[1],release_times,'18-Sep-2014', '22-Sep-2014')

So winds were towards the north west.

September 20¶

Drop 1¶

Buoy 3¶

In [25]:

buoy = 2
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,20))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0003
Released 2014-09-20 04:40:08.999996 at 49.0754766667 , -123.439433333 for 12.5005555561 hours
NEMO coords: 431 269

Out[25]:

(48.8, 49.4)

Buoy 4¶

In [26]:

buoy = 3
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,20))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0004
Released 2014-09-20 04:40:08.000001 at 49.0755083333 , -123.439546667 for 12.5000000011 hours
NEMO coords: 431 269

Out[26]:

(48.8, 49.4)

Drop 2¶

Buoy 7¶

In [27]:

buoy = 6
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj =plot_buoy(ax,tracks,dt.datetime(2014,9,20))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0007
Released 2014-09-20 11:10:08.000005 at 49.0787216667 , -123.441358333 for 6.83361111 hours
NEMO coords: 432 269

Out[27]:

(48.8, 49.4)

Buoy 6¶

In [28]:

buoy = 5
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj =plot_buoy(ax,tracks,dt.datetime(2014,9,20),i=0)
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0006
Released 2014-09-20 11:20:07.999999 at 49.0793333333 , -123.4434 for 6.83361111028 hours
NEMO coords: 432 269

Out[28]:

(48.8, 49.4)

Drop 3¶

Buoy 8¶

In [29]:

buoy = 7
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,20))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0008
Released 2014-09-20 13:20:08.999997 at 49.0768533333 , -123.44025 for 5.33305555778 hours
NEMO coords: 432 268

Out[29]:

(48.8, 49.4)

Buoy 10¶

In [30]:

buoy = 9
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,20))
ax.set_xlim([-123.6,-123.1])
ax.set_ylim([48.8,49.4])

UBC-I-0010
Released 2014-09-20 13:20:08.000002 at 49.0767216667 , -123.440253333 for 5.33333333417 hours
NEMO coords: 432 268

Out[30]:

(48.8, 49.4)

Tides and Winds¶

In [31]:

fig,ax=plt.subplots(1,1,figsize=(10,4))
t1=dt.datetime(2014,9,20,4,40)
t2=dt.datetime(2014,9,20,11,10)
t3=dt.datetime(2014,9,20,13,20)
release_times=[t1,t2,t3];
#tides
plot_tides(ax,release_times,'18-Sep-2014', '22-Sep-2014')
#winds
fig,axs = plt.subplots(2,1,figsize=(10,7))
plot_winds(axs[0],axs[1],release_times,'18-Sep-2014', '22-Sep-2014')

Late Sept 20 to Sept 21¶

In [32]:

buoy = 5
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,20),i=-1)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0006
Released 2014-09-20 20:20:07.999999 at 49.0969733333 , -123.312955 for 21.5080555547 hours
NEMO coords: 427 290

Out[32]:

(48.8, 49.4)

In [33]:

fig,ax=plt.subplots(1,1,figsize=(10,4))
t1=dt.datetime(2014,9,20,20,20)
release_times=[t1];
#tides
plot_tides(ax,release_times,'18-Sep-2014', '22-Sep-2014')
#winds
fig,axs = plt.subplots(2,1,figsize=(10,7))
plot_winds(axs[0],axs[1],release_times,'18-Sep-2014', '22-Sep-2014')

Winds were strongish and southeast.

September 21¶

Buoy 10¶

In [34]:

buoy = 9
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,21))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0010
Released 2014-09-21 06:30:09.000005 at 49.0710983333 , -123.311388333 for 12.3338888883 hours
NEMO coords: 422 287

Out[34]:

(48.8, 49.4)

Buoy 8¶

In [35]:

buoy = 7
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,21))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0008
Released 2014-09-21 08:50:08.999997 at 49.0574933333 , -123.354403333 for 9.83361111306 hours
NEMO coords: 422 279

Out[35]:

(48.8, 49.4)

Buoy 7¶

In [36]:

buoy = 6
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,21))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0007
Released 2014-09-21 09:10:06.999997 at 49.0447116667 , -123.403411667 for 9.00027777889 hours
NEMO coords: 423 270

Out[36]:

(48.8, 49.4)

Buoy 4¶

In [37]:

buoy = 3
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,21))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0004
Released 2014-09-21 09:30:08 at 49.0316116667 , -123.448503333 for 9.16722222333 hours
NEMO coords: 424 262

Out[37]:

(48.8, 49.4)

Buoy 3¶

In [38]:

buoy = 2
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,9,21))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0003
Released 2014-09-21 10:00:08.000004 at 49.0190616667 , -123.490253333 for 9.16749999944 hours
NEMO coords: 424 254

Out[38]:

(48.8, 49.4)

In [39]:

fig,ax=plt.subplots(1,1,figsize=(10,4))
t1=dt.datetime(2014,9,21,6,30)
t2=dt.datetime(2014,9,21,8,50)
t3=dt.datetime(2014,9,21,9,10)
t4=dt.datetime(2014,9,21,9,30)
t5=dt.datetime(2014,9,21,10,0)
release_times=[t1,t2,t3,t4,t5];
#tides
plot_tides(ax,release_times,'20-Sep-2014', '22-Sep-2014')
#winds
fig,axs = plt.subplots(2,1,figsize=(10,7))
plot_winds(axs[0],axs[1],release_times,'20-Sep-2014', '22-Sep-2014')

Vertical lines are the start times of the various drops.

October 8¶

Drop 1¶

In [73]:

buoy = 2
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(7,7))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8), fancy=True)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0003
Released 2014-10-08 16:10:07.000003 at 49.1054983333 , -123.319833333 for 28.7502777783 hours

Out[73]:

(48.8, 49.4)

In [ ]:

buoy = 2
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

In [41]:

buoy = 3
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0004
Released 2014-10-08 17:30:08.000004 at 49.1064116667 , -123.320053333 for 28.2499999981 hours
NEMO coords: 429 290

Out[41]:

(48.8, 49.4)

In [42]:

buoy = 9
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8))
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0010
Released 2014-10-08 19:20:06.999998 at 49.1071733333 , -123.319563333 for 6.66666666611 hours
NEMO coords: 429 290

Out[42]:

(48.8, 49.4)

Drop 2¶

In [43]:

buoy = 1
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8),i=-1)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

UBC-I-0002
Released 2014-10-08 16:05:07.999996 at 49.0956433333 , -123.320643333 for 26.4997222228 hours
NEMO coords: 427 288

Out[43]:

(48.8, 49.4)

In [ ]:

buoy = 4
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8),i=-1)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

In [ ]:

buoy = 7
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8),i=-1)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

Drop 3¶

In [ ]:

buoy = 0
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8),i=-1)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

In [ ]:

buoy = 5
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8),i=-1)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

In [ ]:

buoy = 6
name, tracks=organize_info(buoy,ubc)
print name
fig,ax = plt.subplots(1,1,figsize=(5,5))
traj=plot_buoy(ax,tracks,dt.datetime(2014,10,8),i=-1)
ax.set_xlim([-123.6,-123])
ax.set_ylim([48.8,49.4])

In [ ]:

fig,ax=plt.subplots(1,1,figsize=(10,4))
t1=dt.datetime(2014,10,8,16,0)
t2=dt.datetime(2014,10,8,17,30)
t3=dt.datetime(2014,10,8,19,0)
release_times=[t1,t2,t3];
#tides
plot_tides(ax,release_times,'6-Oct-2014', '10-Oct-2014')
#winds
fig,axs = plt.subplots(2,1,figsize=(10,7))
plot_winds(axs[0],axs[1],release_times,'6-Oct-2014', '10-Oct-2014')

Summary¶

** September 19**

Release 3 particles with the following paramters

NEMO coords (j,i): 432,269 at 10:20 for 11 hours
NEMO coords (j,i): 432,269 at 11:40 for 10 hours
NEMO coords (j,i): 433,269 at 13:10 for 8 hours

** September 20**

Release 3 particles with the following parameters

NEMO coords (j,i): 431,269 at 4:40 for 13 hours
NEMO coords (j,i): 432,269 at 11:10 for 7 hours
NEMO coords (j,i): 432,268 at 13:20 for 6 hours

** Late September 20**

NEMO coords (j,i): 427,290 at 20:20 for 22 hours

** September 21**

NEMO coords (j,i): 422,287 at 06:30 for 12 hours (UBC 10)
NEMO coords (j,i): 422,279 at 08:50 for 10 hours (UBC 8)
NEMO coords (j,i): 423,270 at 09:10 for 9 hours (UBC 7)
NEMO coords (j,i): 424,262 at 09:30 for 9 hours (UBC 4)
NEMO coords (j,i): 424,254 at 10:00 for 9 hours (UBC 3)

** October 8**

Drop 1

UBC-I-0003 - NEMO coods (j,i): 429, 290 at 16:00 hours for 29 hours
UBC-I-0004 - NEMO coods (j,i): 429, 290 at 17:30 hours for 29 hours
UBC-I-0010 - NEMO coods (j,i): 429, 290 at 19:30 hours for 7 hours

Drop 2

UBC-I-0002 - NEMO coods (j,i): 427, 288 at 16:00 hours for 27 hours
UBC-I-0005 - NEMO coods (j,i): 427, 289 at 17:30 hours for 27 hours
UBC-I-0008 - NEMO coods (j,i): 427, 289 at 19:30 hours for 25 hours (Note: find_model_point returned empty, but the lat/lon are similar to UBC-I-0005 so I used that one. perhaps tolerance is to restricitve on find_model_point)

Drop 3

UBC-I-0001 - NEMO coods (j,i): 425, 287 at 16:00 hours for 27 hours
UBC-I-0006 - NEMO coods (j,i): 425, 287 at 17:30 hours for 24 hours
UBC-I-0007 - NEMO coods (j,i): 425, 287 at 19:00 hours for 24 hours

Note: Mark's email indicated UBC-I-0006 was released in the "Drop 1" group and UBC-I-0004 was released in the "Drop 3" group. I have it the other way around..

Next Steps¶

1.Possibly create a drifters.py module to streamline analysis

In [ ]: