utils/matlab/drawedges.m

function [] = drawedges(xg,yg,varargin)
% Draws edges of cube on plot
% >> drawedges(xg,yg)

[nx ny nt]=size(xg);
hnx=ceil(nx/2);
hny=ceil(ny/2);

if ~isempty(varargin)
 col=varargin{1};
else
 col='k';
end

if ndims(xg)==2 & size(xg,1)==6*size(xg,2)
 disp('1');
 [nx ny nt]=size(xg);
 x=permute( reshape(xg,[nx/6 6 ny]),[1 3 2]);
 y=permute( reshape(yg,[nx/6 6 ny]),[1 3 2]);
 [nx ny nt]=size(xg);
elseif ndims(xg)==3 & size(xg,2)==6
 x=permute( X,[1 3 2]);
 y=permute( Y,[1 3 2]);
 [nx ny nt]=size(xg);
elseif ndims(xg)==3 & size(xg,3)==6
 [nx ny nt]=size(xg);
 x=xg;y=yg;
else
 size(xg)
 size(yg)
 error('Dimensions should be 2 or 3 dimensions: NxNx6, 6NxN or Nx6xN');
end


for sh=-180:90:180;
 for k=1:6;
  h=line( longitude(x(:,1,k))+sh, y(:,1,k));set(h,'Color',col);
  h=line( longitude(x(1,:,k))+sh, y(1,:,k));set(h,'Color',col);
 end
end
1	function [] = drawedges(xg,yg,varargin)
2	% Draws edges of cube on plot
3	% >> drawedges(xg,yg)
4
5	[nx ny nt]=size(xg);
6	hnx=ceil(nx/2);
7	hny=ceil(ny/2);
8
9	if ~isempty(varargin)
10	col=varargin{1};
11	else
12	col='k';
13	end
14
15	if ndims(xg)==2 & size(xg,1)==6*size(xg,2)
16	disp('1');
17	[nx ny nt]=size(xg);
18	x=permute( reshape(xg,[nx/6 6 ny]),[1 3 2]);
19	y=permute( reshape(yg,[nx/6 6 ny]),[1 3 2]);
20	[nx ny nt]=size(xg);
21	elseif ndims(xg)==3 & size(xg,2)==6
22	x=permute( X,[1 3 2]);
23	y=permute( Y,[1 3 2]);
24	[nx ny nt]=size(xg);
25	elseif ndims(xg)==3 & size(xg,3)==6
26	[nx ny nt]=size(xg);
27	x=xg;y=yg;
28	else
29	size(xg)
30	size(yg)
31	error('Dimensions should be 2 or 3 dimensions: NxNx6, 6NxN or Nx6xN');
32	end
33
34
35	for sh=-180:90:180;
36	for k=1:6;
37	h=line( longitude(x(:,1,k))+sh, y(:,1,k));set(h,'Color',col);
38	h=line( longitude(x(1,:,k))+sh, y(1,:,k));set(h,'Color',col);
39	end
40	end