
[table([1:length(L)]','VariableNames',{'Number'}) table(datestr(L(:,1)))  T ]

%
% Found table indicies by examining fig2 of 20221207T164716 and the table above.
%
tableIndcies = [7 10 11 14 15];

j=1;
for( i = tableIndcies )
  fprintf('\n\nEpisode %d, Number %d.\n',j,i)  

  %Missed here.  Prop off:
  [startTimeIndx startTime]= closest(L(i,end),c.latitude.time);

  %Resumed here. Prop on:
  [stopTimeIndx  stopTime] = closest(L(i+1,1),c.latitude.time);

  [eStart nStart] = deg2utm(c.latitude.value(startTimeIndx)*180/pi, c.longitude.value(startTimeIndx)*180/pi);
  [eStop nStop]   = deg2utm(c.latitude.value(stopTimeIndx)*180/pi, c.longitude.value(stopTimeIndx)*180/pi);

  fprintf('Drift time = %s.\n', datestr(L(i,end)));

  driftDur=(stopTime-startTime)*86400;
  fprintf('Drift duration = %d seconds.\n', round(driftDur));

  driftDist = norm([eStop nStop] - [eStart nStart]);
  fprintf('Drift distance = %d meters.\n', round(driftDist));

  driftSpeed = driftDist/driftDur;
  fprintf('Drift speed = %.1f cm/s.\n', driftSpeed*100);
  driftSpeedArray(j,1) = driftSpeed;

  driftDir = atan2(eStop-eStart,nStop-nStart);
  if(driftDir < 0) driftDir = driftDir + 2*pi;end
  fprintf('Drift direction = %d degrees.\n', round(driftDir*180/pi));
  driftDirArray(j,1) = driftDir;
  j = j + 1;
end
fprintf('\n\n')
fprintf('Average Current = %.1f cm/sec.\n', mean(driftSpeedArray)*100);
fprintf('Standard Deviation of Current = %.1f cm/sec.\n', std(driftSpeedArray)*100);
fprintf('Average Direction = %d degrees.\n', round(mean(driftDirArray)*180/pi));
fprintf('Standard Deviation of Direction = %d degrees.\n', round(std(driftDirArray)*180/pi));
