function [SurfaceDrift,TimeatSurface] = SurfaceStats(s)
% Find the surface drift and how long the spray spent at the surface
% Ben Werb, adapted from Yui Takeshita plot_spraysat.m

TimeatSurface = 24*60*(s.sdn_(1,end) - s.sdn_(2,end-1));

ellipsoid = wgs84Ellipsoid('kilometer');
SurfaceDrift = abs(distance(s.lat_(1,end), s.lon_(1,end), s.lat_(2,end-1), s.lon_(2,end-1),ellipsoid));

end