function [longs, lats, heights, result] = gravity (data_matrix, range)
dial_after_corr = drift_corr(data_matrix(1:range, 2),data_matrix(1:range,11),data_matrix(1:range,3),data_matrix(1:range,4),data_matrix(1:range,12),data_matrix(1:range,13),data_matrix(1:range,14));
raw_g = dial2gravity (dial_after_corr);
lat_degrees = data_matrix (1:range, 5);
lat_minutes = data_matrix (1:range, 6)./60;
longs = -114 - data_matrix(1:range, 8)./60;
lat_radians = ((lat_degrees + lat_minutes).*pi)./180;
g_after_lat = lat_correction (raw_g, lat_radians);
g_after_car= g_after_lat - g_after_lat(1); 
heights = data_matrix(1:range, 9) - data_matrix(1:range, 10) + 1.293;
g_after_faa = g_after_car + 0.307.*heights';
%freeair (g_after_car, heights);


g_after_bouguer = bouguer_corr (g_after_faa, heights');
result = g_after_bouguer;
mean_g = mean(result);
lats = lat_degrees + lat_minutes;


figure;
plot (heights, result-mean_g, 'ks');
title ('after bouguer');
figure;
plot (heights, g_after_faa, 'ks', 'markersize', 14);;
title ('after faa');
figure;
plot (heights, g_after_lat-mean(g_after_lat), 'ks');
title ('lat corr');


[x,y] = meshgrid (  min(longs) : 0.0001 : max(longs),  min(lat_degrees + lat_minutes) : 0.0001 : max(lat_degrees + lat_minutes));

z = griddata (longs, lat_degrees + lat_minutes, g_after_bouguer', x, y,'cubic');
contourf (x,y,z,30);
colorbar;
%shading  interp;
title 'g after Bouguer , Lat and FAA corrections');
hold on;
for i=1:range
	plot(longs(i), lat_degrees(i)+lat_minutes(i), '.w', 'markersize', 20);
	text(longs(i), lat_degrees(i)+lat_minutes(i), sprintf(' %d',data_matrix(i,1)));
end;
axis tight

figure;
plot3(longs, lat_degrees+lat_minutes, heights,'.-k');