matlab_class/gcmfaces_IO/convert2gcmfaces.m

function [v1]=convert2gcmfaces(v0,varargin);

aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');

if doGcm2Faces;

  if nargin==1; nFaces=5; else; nFaces=varargin{1}; end;

  if nFaces>1;

    nn=size(v0,1); 
    pp=size(v0,2)/nn;
    mm=(pp+4-nFaces)/4*nn; 

    n3=size(v0,3); n4=size(v0,4);
    for i3=1:n3; for i4=1:n4;
      v00=v0(:,:,i3,i4);

      v00=v00(:);
      i0=1; i1=nn*mm; v1{1}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm);
      i0=i1+1; i1=i1+nn*mm; v1{2}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm);
      i0=i1+1; i1=i1+nn*nn; v1{3}(:,:,i3,i4)=reshape(v00(i0:i1),nn,nn);
      i0=i1+1; i1=i1+nn*mm; v1{4}(:,:,i3,i4)=reshape(v00(i0:i1),mm,nn);
      i0=i1+1; i1=i1+nn*mm; v1{5}(:,:,i3,i4)=reshape(v00(i0:i1),mm,nn);
      if nFaces==6;
         i0=i1+1; i1=i1+nn*nn;
      end;

    end; end;
  else;
    v1={v0};
  end;

  if nFaces==1; gridType='ll'; elseif nFaces==5; gridType='llc'; elseif nFaces==6; gridType='cube'; end;
  v1=gcmfaces(v1,gridType);

else;

  nFaces=get(v0,'nFaces');

  if nFaces>1;

    v0_faces=v0; clear v0; 
    for iFace=1:nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end; 

    nn=size(v0{1},1); mm=size(v0{1},2); 
    pp=mm/nn*4+nFaces-4;

    n3=size(v0{1},3); n4=size(v0{1},4);
    v1=NaN*zeros(nn,nn*pp,n3,n4);
    for i3=1:n3; for i4=1:n4;

      v11=NaN*zeros(nn*nn*pp,1);
    
      i0=1; i1=nn*mm; tmp1=v0{1}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
      i0=i1+1; i1=i1+nn*mm; tmp1=v0{2}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
      i0=i1+1; i1=i1+nn*nn; tmp1=v0{3}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
      i0=i1+1; i1=i1+nn*mm; tmp1=v0{4}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
      i0=i1+1; i1=i1+nn*mm; tmp1=v0{5}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
      if nFaces==6; 
         i0=i1+1; i1=i1+nn*nn;
      end;

      v1(:,:,i3,i4)=reshape(v11,nn,nn*pp);
    end; end;

  else;
    v1=v0{1}; 
  end;
 
end;


1	function [v1]=convert2gcmfaces(v0,varargin);
2
3	aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');
4
5	if doGcm2Faces;
6
7	if nargin==1; nFaces=5; else; nFaces=varargin{1}; end;
8
9	if nFaces>1;
10
11	nn=size(v0,1);
12	pp=size(v0,2)/nn;
13	mm=(pp+4-nFaces)/4*nn;
14
15	n3=size(v0,3); n4=size(v0,4);
16	for i3=1:n3; for i4=1:n4;
17	v00=v0(:,:,i3,i4);
18
19	v00=v00(:);
20	i0=1; i1=nn*mm; v1{1}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm);
21	i0=i1+1; i1=i1+nn*mm; v1{2}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm);
22	i0=i1+1; i1=i1+nn*nn; v1{3}(:,:,i3,i4)=reshape(v00(i0:i1),nn,nn);
23	i0=i1+1; i1=i1+nn*mm; v1{4}(:,:,i3,i4)=reshape(v00(i0:i1),mm,nn);
24	i0=i1+1; i1=i1+nn*mm; v1{5}(:,:,i3,i4)=reshape(v00(i0:i1),mm,nn);
25	if nFaces==6;
26	i0=i1+1; i1=i1+nn*nn;
27	end;
28
29	end; end;
30	else;
31	v1={v0};
32	end;
33
34	if nFaces==1; gridType='ll'; elseif nFaces==5; gridType='llc'; elseif nFaces==6; gridType='cube'; end;
35	v1=gcmfaces(v1,gridType);
36
37	else;
38
39	nFaces=get(v0,'nFaces');
40
41	if nFaces>1;
42
43	v0_faces=v0; clear v0;
44	for iFace=1:nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
45
46	nn=size(v0{1},1); mm=size(v0{1},2);
47	pp=mm/nn*4+nFaces-4;
48
49	n3=size(v0{1},3); n4=size(v0{1},4);
50	v1=NaNzeros(nn,nnpp,n3,n4);
51	for i3=1:n3; for i4=1:n4;
52
53	v11=NaNzeros(nnnn*pp,1);
54
55	i0=1; i1=nn*mm; tmp1=v0{1}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
56	i0=i1+1; i1=i1+nn*mm; tmp1=v0{2}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
57	i0=i1+1; i1=i1+nn*nn; tmp1=v0{3}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
58	i0=i1+1; i1=i1+nn*mm; tmp1=v0{4}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
59	i0=i1+1; i1=i1+nn*mm; tmp1=v0{5}(:,:,i3,i4); v11(i0:i1)=tmp1(:);
60	if nFaces==6;
61	i0=i1+1; i1=i1+nn*nn;
62	end;
63
64	v1(:,:,i3,i4)=reshape(v11,nn,nn*pp);
65	end; end;
66
67	else;
68	v1=v0{1};
69	end;
70
71	end;
72
73