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;
    
    v1=gcmfaces(v1);
    
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')|strcmp(mygrid.fileFormat,'nctiles');
        v00=reshape(v0,[n1*n2 n3*n4*n5]);
        i0=0; i1=0;
        for iFace=1:mygrid.nFaces;
            i0=i1+1;
            nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
            i1=i1+nn*mm;
            v1{iFace}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
        end;
    end;
    
    v1=gcmfaces(v1);
    
    if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end;
    
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')|strcmp(mygrid.fileFormat,'nctiles');
        
        %   is there a reason for this?
        %         v0_faces=v0; clear v0;
        %         for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
        
        if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end;

        aa=0;bb=0;
        for iFace=1:mygrid.nFaces;
            aa=aa+prod(mygrid.facesSize(iFace,:));
            bb=bb+prod(size(v0{iFace}));
        end;
        bb=bb/aa;
        v11=NaN*zeros(aa,bb);
        
        i0=0; i1=0;
        for iFace=1:mygrid.nFaces;
            i0=i1+1;
            nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
            i1=i1+nn*mm;
            v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]);
        end;
        nn=mygrid.ioSize(1); mm=mygrid.ioSize(2);        
        v1=reshape(v11,[nn mm n3 n4 n5]);
        
    end;
    
end;


1	gforget	1.1	function [v1]=convert2gcmfaces(v0,varargin);
2	gforget	1.9	%object: converts model output to gcmfaces object
3	gforget	1.5	% 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	gforget	1.9
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			v1=gcmfaces(v1);
25
26	gforget	1.7	elseif ~doGcm2Faces&mygrid.gcm2facesFast;
27	gforget	1.9
28			[n1,n2,n3,n4,n5]=size(v0{1});
29
30			v1=zeros(mygrid.faces2gcmSize(1)mygrid.faces2gcmSize(2),n3n4*n5);
31			for iFace=1:mygrid.nFaces;
32			nn=size(mygrid.XC{iFace});
33			v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)nn(2) n3n4*n5]);
34			end;
35			v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]);
36
37	gforget	1.7	elseif doGcm2Faces;
38	gforget	1.9
39			[n1,n2,n3,n4,n5]=size(v0);
40
41			if strcmp(mygrid.fileFormat,'straight');
42			v1={v0};
43			elseif strcmp(mygrid.fileFormat,'cube');
44			for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;
45	gforget	1.11	elseif strcmp(mygrid.fileFormat,'compact')\|strcmp(mygrid.fileFormat,'nctiles');
46	gforget	1.10	v00=reshape(v0,[n1n2 n3n4*n5]);
47			i0=0; i1=0;
48			for iFace=1:mygrid.nFaces;
49			i0=i1+1;
50			nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
51			i1=i1+nn*mm;
52			v1{iFace}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
53	gforget	1.9	end;
54	gforget	1.2	end;
55	gforget	1.9
56			v1=gcmfaces(v1);
57
58	gforget	1.10	if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end;
59
60	gforget	1.1	else;
61	gforget	1.9
62			[n1,n2,n3,n4,n5]=size(v0{1});
63
64			if strcmp(mygrid.fileFormat,'straight');
65			v1=v0{1};
66			elseif strcmp(mygrid.fileFormat,'cube');
67			v1=zeros(n2*6,n2,n3,n4,n5);
68			for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;
69	gforget	1.11	elseif strcmp(mygrid.fileFormat,'compact')\|strcmp(mygrid.fileFormat,'nctiles');
70	gforget	1.9
71	gforget	1.10	% is there a reason for this?
72			% v0_faces=v0; clear v0;
73			% for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
74	gforget	1.9
75	gforget	1.10	if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end;
76
77			aa=0;bb=0;
78			for iFace=1:mygrid.nFaces;
79			aa=aa+prod(mygrid.facesSize(iFace,:));
80			bb=bb+prod(size(v0{iFace}));
81			end;
82			bb=bb/aa;
83			v11=NaN*zeros(aa,bb);
84	gforget	1.9
85	gforget	1.10	i0=0; i1=0;
86			for iFace=1:mygrid.nFaces;
87			i0=i1+1;
88			nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
89			i1=i1+nn*mm;
90			v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]);
91	gforget	1.9	end;
92	gforget	1.10	nn=mygrid.ioSize(1); mm=mygrid.ioSize(2);
93			v1=reshape(v11,[nn mm n3 n4 n5]);
94	gforget	1.9
95			end;
96
97	gforget	1.1	end;
98
99