function [X,Y,U,V,truedir,vel,dist]=calcdirvel(time,long,lat)
  [r c] = size(long);  
  %distance in km
  for k= 1:c-1
    [east(k),north(k)] = lonlat2km(long(k),lat(k),long(k+1),lat(k+1));
    %time elapsed in secs
    timediff(k) = (time(k+1)-time(k))*86400;
    X(k) = mean([long(k) long(k+1)]);
    Y(k) = mean([lat(k) lat(k+1)]);
    V(k) = (north(k)*100000)/timediff(k);
    U(k) = (east(k)*100000)/timediff(k);
    truedir(k) = (uv2angle(U(k),V(k))*180/pi);
    dist(k) = sqrt(north(k)^2+east(k)^2);
    vel(k)=(dist(k)*100000)/timediff(k);
end