Diagnostics/DiagUtility/longitude.m

function [X] = longitude(x)
% X=longitude(x);
%
% tries to determine best range of longitude (e.g. -180<180 or 0<360)
% so that coordinate (x) doesn't span a discontinuity.
%

% also works for radians which are assumed if range of x<=2*pi

minx=min(min(min(x)));
maxx=max(max(max(x)));
%if maxx-minx < 2.2*pi
% units=180/pi;
%else
 units=1;
%end
minx=min(min(min(x*units)));
maxx=max(max(max(x*units)));

X=mod(720+x*units,360);
maxP=max(max(max(X)));
minP=min(min(min(X)));

XX=mod(X+180,360)-180;
maxM=max(max(max(XX)));
minM=min(min(min(XX)));

if maxP-minP > maxM-minM
 X=XX;
end
1	function [X] = longitude(x)
2	% X=longitude(x);
3	%
4	% tries to determine best range of longitude (e.g. -180<180 or 0<360)
5	% so that coordinate (x) doesn't span a discontinuity.
6	%
7
8	% also works for radians which are assumed if range of x<=2*pi
9
10	minx=min(min(min(x)));
11	maxx=max(max(max(x)));
12	%if maxx-minx < 2.2*pi
13	% units=180/pi;
14	%else
15	units=1;
16	%end
17	minx=min(min(min(x*units)));
18	maxx=max(max(max(x*units)));
19
20	X=mod(720+x*units,360);
21	maxP=max(max(max(X)));
22	minP=min(min(min(X)));
23
24	XX=mod(X+180,360)-180;
25	maxM=max(max(max(XX)));
26	minM=min(min(min(XX)));
27
28	if maxP-minP > maxM-minM
29	X=XX;
30	end