/[MITgcm]/MITgcm_contrib/gael/matlab_class/gcmfaces_IO/convert2gcmfaces.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/matlab_class/gcmfaces_IO/convert2gcmfaces.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by gforget, Wed Feb 10 14:46:48 2010 UTC revision 1.9 by gforget, Tue Sep 11 15:11:31 2012 UTC
# Line 1  Line 1 
1  function [v1]=convert2gcmfaces(v0,varargin);  function [v1]=convert2gcmfaces(v0,varargin);
2    %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    
9  aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');  input_list_check('convert2gcmfaces',nargin);
   
 if doGcm2Faces;  
   
   if nargin==1; nFaces=5; else; nFaces=varargin{1}; end;  
10    
11    if nFaces>1;  aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');
   
     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;  
12    
13    if nFaces==1; gridType='ll'; elseif nFaces==5; gridType='llc'; elseif nFaces==6; gridType='cube'; end;  global mygrid;
   v1=gcmfaces(v1,gridType);  
14    
15    if doGcm2Faces&mygrid.gcm2facesFast;
16        
17        [n1,n2,n3,n4,n5]=size(v0);
18        
19        v0=reshape(v0,[n1*n2 n3*n4*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    elseif ~doGcm2Faces&mygrid.gcm2facesFast;
27        
28        [n1,n2,n3,n4,n5]=size(v0{1});
29        
30        v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*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) n3*n4*n5]);
34        end;
35        v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]);
36        
37    elseif doGcm2Faces;
38        
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        elseif strcmp(mygrid.fileFormat,'compact');
46            
47            nn=size(v0,1);
48            pp=size(v0,2)/nn;
49            mm=(pp+4-mygrid.nFaces)/4*nn;
50            
51            v00=reshape(v0,[nn*nn*pp n3*n4*n5]);
52            i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
53            i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
54            i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
55            i0=i1+1; i1=i1+nn*mm; v1{4}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
56            i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
57            if mygrid.nFaces==6;
58                i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
59            end;
60            
61        end;
62        
63        v1=gcmfaces(v1);
64        
65  else;  else;
66        
67    nFaces=get(v0,'nFaces');      [n1,n2,n3,n4,n5]=size(v0{1});
68        
69    if nFaces>1;      if strcmp(mygrid.fileFormat,'straight');
70            v1=v0{1};
71      v0_faces=v0; clear v0;      elseif strcmp(mygrid.fileFormat,'cube');
72      for iFace=1:nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;          v1=zeros(n2*6,n2,n3,n4,n5);
73            for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;
74      nn=size(v0{1},1); mm=size(v0{1},2);      elseif strcmp(mygrid.fileFormat,'compact');
75      pp=mm/nn*4+nFaces-4;          
76            v0_faces=v0; clear v0;
77      n3=size(v0{1},3); n4=size(v0{1},4);          for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
78      v1=NaN*zeros(nn,nn*pp,n3,n4);          
79      for i3=1:n3; for i4=1:n4;          nn=size(v0{1},1); mm=size(v0{1},2);
80            pp=mm/nn*4+mygrid.nFaces-4;
81        v11=NaN*zeros(nn*nn*pp,1);          
82                n3=size(v0{1},3); n4=size(v0{1},4); n5=size(v0{1},5);
83        i0=1; i1=nn*mm; tmp1=v0{1}(:,:,i3,i4); v11(i0:i1)=tmp1(:);          v1=NaN*zeros(nn,nn*pp,n3,n4,n5);
84        i0=i1+1; i1=i1+nn*mm; tmp1=v0{2}(:,:,i3,i4); v11(i0:i1)=tmp1(:);          
85        i0=i1+1; i1=i1+nn*nn; tmp1=v0{3}(:,:,i3,i4); v11(i0:i1)=tmp1(:);          v11=NaN*zeros(nn*nn*pp,n3*n4*n5);
86        i0=i1+1; i1=i1+nn*mm; tmp1=v0{4}(:,:,i3,i4); v11(i0:i1)=tmp1(:);          
87        i0=i1+1; i1=i1+nn*mm; tmp1=v0{5}(:,:,i3,i4); v11(i0:i1)=tmp1(:);          i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
88        if nFaces==6;          i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
89           i0=i1+1; i1=i1+nn*nn;          i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
90        end;          i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
91            i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
92        v1(:,:,i3,i4)=reshape(v11,nn,nn*pp);          if mygrid.nFaces==6;
93      end; end;              i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
94            end;
95    else;          
96      v1=v0{1};          v1=reshape(v11,[nn nn*pp n3 n4 n5]);
97    end;          
98        end;
99        
100  end;  end;
101    
102    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22