function lat_g = lat_correction (g, rad_latitudes)
	theory = 978032.* (1 + (5.2789e-3).*((sin(rad_latitudes)).^2) - (2.35e-6).*((sin(rad_latitudes)).^4));
	temp = g - theory';
	temp = temp + 975000;
	lat_g = temp;