matlab_class/gcmfaces_IO/convert2gcmfaces.m

function [v1]=convert2gcmfaces(v0,varargin);
%object:    converts model output to gcmfaces object 
%           or vice versa
%input:     v0 model output array (resp. gcmfaces object)
%output:    v1 gcmfaces object (model output array)
%
%note:      global mygrid parameters (nFaces, fileFormat) are used

input_list_check('convert2gcmfaces',nargin);

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

global mygrid;

if doGcm2Faces&mygrid.gcm2facesFast;

  [n1,n2,n3,n4,n5]=size(v0);

  v0=reshape(v0,[n1*n2 n3*n4*n5]);
  for iFace=1:mygrid.nFaces;
    v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]);
  end;

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

elseif ~doGcm2Faces&mygrid.gcm2facesFast;

  [n1,n2,n3,n4,n5]=size(v0{1});

  v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5);
  for iFace=1:mygrid.nFaces;
    nn=size(mygrid.XC{iFace});
    v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]);
  end;
  v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]);

elseif doGcm2Faces;

  [n1,n2,n3,n4,n5]=size(v0);

  if strcmp(mygrid.fileFormat,'straight');
    v1={v0};
  elseif strcmp(mygrid.fileFormat,'cube');
    for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;
  elseif strcmp(mygrid.fileFormat,'compact');

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

    v00=reshape(v0,[nn*nn*pp n3*n4*n5]);
    i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
    i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
    i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
    i0=i1+1; i1=i1+nn*mm; v1{4}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
    i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
    if mygrid.nFaces==6;
       i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
    end;

  end;

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

else;

  [n1,n2,n3,n4,n5]=size(v0{1});

  if strcmp(mygrid.fileFormat,'straight');
    v1=v0{1};
  elseif strcmp(mygrid.fileFormat,'cube');
    v1=zeros(n2*6,n2,n3,n4,n5);
    for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;
  elseif strcmp(mygrid.fileFormat,'compact');

    v0_faces=v0; clear v0; 
    for iFace=1:mygrid.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+mygrid.nFaces-4;

    n3=size(v0{1},3); n4=size(v0{1},4); n5=size(v0{1},5);
    v1=NaN*zeros(nn,nn*pp,n3,n4,n5);

    v11=NaN*zeros(nn*nn*pp,n3*n4*n5);
  
   i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
   i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
   i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
   i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
   i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
   if mygrid.nFaces==6; 
      i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
   end;

   v1=reshape(v11,[nn nn*pp n3 n4 n5]);

  end;
 
end;


1	gforget	1.1	function [v1]=convert2gcmfaces(v0,varargin);
2	gforget	1.5	%object: converts model output to gcmfaces object
3			% or vice versa
4			%input: v0 model output array (resp. gcmfaces object)
5			%output: v1 gcmfaces object (model output array)
6			%
7			%note: global mygrid parameters (nFaces, fileFormat) are used
8	gforget	1.1
9	gforget	1.5	input_list_check('convert2gcmfaces',nargin);
10	gforget	1.3
11	gforget	1.1	aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');
12
13	gforget	1.5	global mygrid;
14
15	gforget	1.7	if doGcm2Faces&mygrid.gcm2facesFast;
16
17			[n1,n2,n3,n4,n5]=size(v0);
18
19			v0=reshape(v0,[n1n2 n3n4*n5]);
20			for iFace=1:mygrid.nFaces;
21			v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]);
22			end;
23
24			if mygrid.nFaces==1; gridType='ll';
25			elseif mygrid.nFaces==5; gridType='llc';
26			elseif mygrid.nFaces==6; gridType='cube';
27			end;
28			v1=gcmfaces(v1,gridType);
29
30			elseif ~doGcm2Faces&mygrid.gcm2facesFast;
31
32			[n1,n2,n3,n4,n5]=size(v0{1});
33
34			v1=zeros(mygrid.faces2gcmSize(1)mygrid.faces2gcmSize(2),n3n4*n5);
35			for iFace=1:mygrid.nFaces;
36			nn=size(mygrid.XC{iFace});
37			v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)nn(2) n3n4*n5]);
38			end;
39			v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]);
40
41			elseif doGcm2Faces;
42	gforget	1.1
43	gforget	1.6	[n1,n2,n3,n4,n5]=size(v0);
44	gforget	1.3
45	gforget	1.5	if strcmp(mygrid.fileFormat,'straight');
46	gforget	1.3	v1={v0};
47	gforget	1.5	elseif strcmp(mygrid.fileFormat,'cube');
48	gforget	1.6	for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;
49	gforget	1.5	elseif strcmp(mygrid.fileFormat,'compact');
50	gforget	1.1
51			nn=size(v0,1);
52			pp=size(v0,2)/nn;
53	gforget	1.5	mm=(pp+4-mygrid.nFaces)/4*nn;
54	gforget	1.1
55	gforget	1.6	v00=reshape(v0,[nnnnpp n3n4n5]);
56			i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
57			i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
58			i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
59			i0=i1+1; i1=i1+nn*mm; v1{4}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
60			i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
61	gforget	1.5	if mygrid.nFaces==6;
62	gforget	1.6	i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
63	gforget	1.2	end;
64	gforget	1.1
65			end;
66
67	gforget	1.5	if mygrid.nFaces==1; gridType='ll';
68			elseif mygrid.nFaces==5; gridType='llc';
69			elseif mygrid.nFaces==6; gridType='cube';
70			end;
71	gforget	1.1	v1=gcmfaces(v1,gridType);
72
73			else;
74
75	gforget	1.6	[n1,n2,n3,n4,n5]=size(v0{1});
76	gforget	1.3
77	gforget	1.5	if strcmp(mygrid.fileFormat,'straight');
78	gforget	1.3	v1=v0{1};
79	gforget	1.5	elseif strcmp(mygrid.fileFormat,'cube');
80	gforget	1.6	v1=zeros(n2*6,n2,n3,n4,n5);
81			for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;
82	gforget	1.5	elseif strcmp(mygrid.fileFormat,'compact');
83	gforget	1.1
84			v0_faces=v0; clear v0;
85	gforget	1.5	for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
86	gforget	1.1
87			nn=size(v0{1},1); mm=size(v0{1},2);
88	gforget	1.5	pp=mm/nn*4+mygrid.nFaces-4;
89	gforget	1.1
90	gforget	1.6	n3=size(v0{1},3); n4=size(v0{1},4); n5=size(v0{1},5);
91			v1=NaNzeros(nn,nnpp,n3,n4,n5);
92	gforget	1.1
93	gforget	1.6	v11=NaNzeros(nnnnpp,n3n4*n5);
94	gforget	1.2
95	gforget	1.6	i0=1; i1=nnmm; tmp1=reshape(v0{1},[nnmm n3n4n5]); v11(i0:i1,:)=tmp1(:,:);
96			i0=i1+1; i1=i1+nnmm; tmp1=reshape(v0{2},[nnmm n3n4n5]); v11(i0:i1,:)=tmp1(:,:);
97			i0=i1+1; i1=i1+nnnn; tmp1=reshape(v0{3},[nnnn n3n4n5]); v11(i0:i1,:)=tmp1(:,:);
98			i0=i1+1; i1=i1+nnmm; tmp1=reshape(v0{4},[mmnn n3n4n5]); v11(i0:i1,:)=tmp1(:,:);
99			i0=i1+1; i1=i1+nnmm; tmp1=reshape(v0{5},[mmnn n3n4n5]); v11(i0:i1,:)=tmp1(:,:);
100	gforget	1.5	if mygrid.nFaces==6;
101	gforget	1.6	i0=i1+1; i1=i1+nnnn; tmp1=reshape(v0{6},[nnnn n3n4n5]); v11(i0:i1,:)=tmp1(:,:);
102	gforget	1.2	end;
103	gforget	1.1
104	gforget	1.6	v1=reshape(v11,[nn nn*pp n3 n4 n5]);
105	gforget	1.1
106			end;
107
108			end;
109
110