%% from dick limebrener

function[sc]= aspect(theta);

% use with matlab 5 to scale longitude
% scales figure with latitude theta(degrees)

sc=cos(theta*pi/180);
set(gca, 'dataaspectratio', [1 sc 1]);