
for i = 2:length(tp_time)
    delta(i) = abs(tp_sr(i) - tp_sr(i-1));
end

delta(1) = 0;
figure('Name', 'Range Delta');
plot(tp_time, delta, 'r');
hold on;

%figure('Name', 'Range Delta vs. Lat');
%plot(tp_lat, delta, '.');
%hold on;
%
%figure('Name', 'Range Delta vs. Lon');
%plot(tp_long, delta, '.');
%hold on;

bigDelta = find(delta>500);
figure('Name', 'Lat/Lon Delta > 500');
plot(tp_long(bigDelta), tp_lat(bigDelta), '.'); 
axis([min(tp_long) max(tp_long) min(tp_lat) max(tp_lat)]);

smallDelta = find(delta<500);
figure('Name', 'Lat/Lon Delta < 500');
plot(tp_long(smallDelta), tp_lat(smallDelta), '.'); 
axis([min(tp_long) max(tp_long) min(tp_lat) max(tp_lat)])

latlon(:,1) =  tp_lat(smallDelta);
latlon(:,2) =  tp_long(smallDelta);

latlon(1,2) = latlon(2,2);
latlon(1,1) = latlon(2,1);

pwr_kml('good delta', latlon);
