ecco-godae/climatology/plot_faces.m

function plot_faces(ng,vF,k,ccB);
% plot_faces(ng,vF,k [,ccB]);
%
% make 1 figure (fig. number = ng) composed of Nface parts
% corresponding to level k (enter k=0 if 2D field) of the structured
% array "vF": vF.f001 = face 1, vF.f002 = face 2 ...
% with color range ccB(1),ccB(2) (default = ccB=[0 0] = no call to caxis)

if nargin < 4, ccB=[0 0]; end

%- extract list of faces:
%  convention: faces identificator is "f{n}" (face number n written with 3 digits)
listV=fieldnames(vF); listF=[];
for i=1:size(listV,1)
 cvar=char(listV(i));
 if length(cvar) == 4 & strncmpi('f00',cvar,3),
  if isempty(listF), listF=cvar; else listF=char(listF,cvar); end
 end
end
listF=cellstr(listF);
Nfaces=size(listF,1);
if Nfaces > 6,
  error([' Nb of faces =',int2str(Nfaces),' > 6 !'])
end

xyP=zeros(6,4);
xyP(1,:)=[0.05  0.05  0.20  0.72];
xyP(2,:)=[0.29  0.05  0.20  0.72];
xyP(3,:)=[0.29  0.80  0.20  0.18];
xyP(4,:)=[0.53  0.05  0.20  0.72];
xyP(5,:)=[0.77  0.05  0.20  0.72];
xyP(6,:)=[0.77  0.80  0.20  0.18];

xyP(3,:)=[0.05  0.80  0.20  0.18]; %- move face 3 in top-left corner
figure(ng);clf;
for n=1:Nfaces,
 cvar=char(listF(n)); var=vF.(cvar);
 if k ~= 0, var=var(:,:,k); end
 if n > 3,  var=fliplr(var'); end %- rotate face 4,5,6
 if n == 3, var=flipud(var'); end %- rotate to put face 3 in top-left corner
 AA=axes('position',xyP(n,:));
 imagesc(var'); set(gca,'YDir','normal');
 if ccB(1) < ccB(2), caxis(ccB); end
end
H=colorbar('peer',AA,'North');
set(H,'Position',[xyP(4,1) 0.90 xyP(4,3) 0.020]);

return
1	jmc	1.1	function plot_faces(ng,vF,k,ccB);
2			% plot_faces(ng,vF,k [,ccB]);
3			%
4			% make 1 figure (fig. number = ng) composed of Nface parts
5			% corresponding to level k (enter k=0 if 2D field) of the structured
6			% array "vF": vF.f001 = face 1, vF.f002 = face 2 ...
7			% with color range ccB(1),ccB(2) (default = ccB=[0 0] = no call to caxis)
8
9			if nargin < 4, ccB=[0 0]; end
10
11			%- extract list of faces:
12			% convention: faces identificator is "f{n}" (face number n written with 3 digits)
13			listV=fieldnames(vF); listF=[];
14			for i=1:size(listV,1)
15			cvar=char(listV(i));
16			if length(cvar) == 4 & strncmpi('f00',cvar,3),
17			if isempty(listF), listF=cvar; else listF=char(listF,cvar); end
18			end
19			end
20			listF=cellstr(listF);
21			Nfaces=size(listF,1);
22			if Nfaces > 6,
23			error([' Nb of faces =',int2str(Nfaces),' > 6 !'])
24			end
25
26			xyP=zeros(6,4);
27			xyP(1,:)=[0.05 0.05 0.20 0.72];
28			xyP(2,:)=[0.29 0.05 0.20 0.72];
29			xyP(3,:)=[0.29 0.80 0.20 0.18];
30			xyP(4,:)=[0.53 0.05 0.20 0.72];
31			xyP(5,:)=[0.77 0.05 0.20 0.72];
32			xyP(6,:)=[0.77 0.80 0.20 0.18];
33
34			xyP(3,:)=[0.05 0.80 0.20 0.18]; %- move face 3 in top-left corner
35			figure(ng);clf;
36			for n=1:Nfaces,
37			cvar=char(listF(n)); var=vF.(cvar);
38			if k ~= 0, var=var(:,:,k); end
39			if n > 3, var=fliplr(var'); end %- rotate face 4,5,6
40			if n == 3, var=flipud(var'); end %- rotate to put face 3 in top-left corner
41			AA=axes('position',xyP(n,:));
42			imagesc(var'); set(gca,'YDir','normal');
43			if ccB(1) < ccB(2), caxis(ccB); end
44			end
45			H=colorbar('peer',AA,'North');
46			set(H,'Position',[xyP(4,1) 0.90 xyP(4,3) 0.020]);
47
48			return