%%%%%%%%%%%%%%%%%%%%%%%%%%%
function setGeoAxis
lat = mean(get(gca,'YLim'));
if abs(lat) >= 90,
	errorbox(['Bad latitude (' mynum2str(lat) ')']),
	return
end
set(gca,'DataAspectRatio',[1/cos(deg2rad(lat)) 1 1])