matlab_class/gcmfaces_convert/convert2array_llc.m

function [FLD]=convert2array_llc(fld,varargin);

if isa(fld,'gcmfaces'); do_gcmfaces2array=1; else; do_gcmfaces2array=0; end;

if do_gcmfaces2array;
   n3=max(size(fld.f1,3),1); n4=max(size(fld.f1,4),1);
   n1=size(fld.f1,1)*4; n2=size(fld.f1,2)+size(fld.f1,1);
   FLD=squeeze(zeros(n1,n2,n3,n4));
else;
   n3=max(size(fld,3),1); n4=max(size(fld,4),1);
   FLD=NaN*varargin{1};
   for k3=1:n3; for k4=1:n4; for iFace=1:FLD.nFaces; 
      iF=num2str(iFace); eval(['FLD.f' iF '(:,:,k3,k4)=FLD.f' iF '(:,:,1,1);']);
   end; end; end;
end;

for k3=1:n3; for k4=1:n4;

if do_gcmfaces2array;
   %ASSEMBLE "LATLON" FACES:
   %----------------------
   FLD0=[fld.f1(:,:,k3,k4);fld.f2(:,:,k3,k4);sym_g(fld.f4(:,:,k3,k4),7,0);sym_g(fld.f5(:,:,k3,k4),7,0)];
   %ADD POLAR CAP:
   %--------------
   pp=fld.f3(:,:,k3,k4); FLDp=[sym_g(pp,5,0);pp.*NaN;sym_g(pp.*NaN,7,0);sym_g(pp.*NaN,6,0)]; 
   FLD1=[FLD0 FLDp];
   %store:
   %------
   FLD(:,:,k3,k4)=FLD1;
else;
   n1=size(FLD.f1,1); n2=size(FLD.f1,2);
   FLD.f1(:,:,k3,k4)=fld(1:n1,1:n2,k3,k4);
   FLD.f2(:,:,k3,k4)=fld(n1+[1:n1],1:n2,k3,k4);
   FLD.f3(:,:,k3,k4)=sym_g(fld([1:n1],n2+1:n2+n1,k3,k4),7,0);
   FLD.f4(:,:,k3,k4)=sym_g(fld(n1*2+[1:n1],1:n2,k3,k4),5,0);
   FLD.f5(:,:,k3,k4)=sym_g(fld(n1*3+[1:n1],1:n2,k3,k4),5,0);
end;

end; end;

1	gforget	1.1	function [FLD]=convert2array_llc(fld,varargin);
2
3			if isa(fld,'gcmfaces'); do_gcmfaces2array=1; else; do_gcmfaces2array=0; end;
4
5			if do_gcmfaces2array;
6			n3=max(size(fld.f1,3),1); n4=max(size(fld.f1,4),1);
7			n1=size(fld.f1,1)*4; n2=size(fld.f1,2)+size(fld.f1,1);
8			FLD=squeeze(zeros(n1,n2,n3,n4));
9			else;
10			n3=max(size(fld,3),1); n4=max(size(fld,4),1);
11			FLD=NaN*varargin{1};
12			for k3=1:n3; for k4=1:n4; for iFace=1:FLD.nFaces;
13			iF=num2str(iFace); eval(['FLD.f' iF '(:,:,k3,k4)=FLD.f' iF '(:,:,1,1);']);
14			end; end; end;
15			end;
16
17			for k3=1:n3; for k4=1:n4;
18
19			if do_gcmfaces2array;
20			%ASSEMBLE "LATLON" FACES:
21			%----------------------
22			FLD0=[fld.f1(:,:,k3,k4);fld.f2(:,:,k3,k4);sym_g(fld.f4(:,:,k3,k4),7,0);sym_g(fld.f5(:,:,k3,k4),7,0)];
23			%ADD POLAR CAP:
24			%--------------
25			pp=fld.f3(:,:,k3,k4); FLDp=[sym_g(pp,5,0);pp.NaN;sym_g(pp.NaN,7,0);sym_g(pp.*NaN,6,0)];
26			FLD1=[FLD0 FLDp];
27			%store:
28			%------
29			FLD(:,:,k3,k4)=FLD1;
30			else;
31			n1=size(FLD.f1,1); n2=size(FLD.f1,2);
32			FLD.f1(:,:,k3,k4)=fld(1:n1,1:n2,k3,k4);
33			FLD.f2(:,:,k3,k4)=fld(n1+[1:n1],1:n2,k3,k4);
34			FLD.f3(:,:,k3,k4)=sym_g(fld([1:n1],n2+1:n2+n1,k3,k4),7,0);
35			FLD.f4(:,:,k3,k4)=sym_g(fld(n1*2+[1:n1],1:n2,k3,k4),5,0);
36			FLD.f5(:,:,k3,k4)=sym_g(fld(n1*3+[1:n1],1:n2,k3,k4),5,0);
37			end;
38
39			end; end;
40