Diagnostics/DiagUtility/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)
 [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	[nx ny nt]=size(xg);
17	x=permute( reshape(xg,[nx/6 6 ny]),[1 3 2]);
18	y=permute( reshape(yg,[nx/6 6 ny]),[1 3 2]);
19	[nx ny nt]=size(xg);
20	elseif ndims(xg)==3 & size(xg,2)==6
21	x=permute( X,[1 3 2]);
22	y=permute( Y,[1 3 2]);
23	[nx ny nt]=size(xg);
24	elseif ndims(xg)==3 & size(xg,3)==6
25	[nx ny nt]=size(xg);
26	x=xg;y=yg;
27	else
28	size(xg)
29	size(yg)
30	error('Dimensions should be 2 or 3 dimensions: NxNx6, 6NxN or Nx6xN');
31	end
32
33
34	for sh=-180:90:180;
35	for k=1:6;
36	h=line( longitude(x(:,1,k))+sh, y(:,1,k));set(h,'Color',col);
37	h=line( longitude(x(1,:,k))+sh, y(1,:,k));set(h,'Color',col);
38	end
39	end