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