/[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.9 by gforget, Tue Sep 11 15:11:31 2012 UTC revision 1.10 by gforget, Tue Sep 11 15:15:19 2012 UTC
# Line 43  elseif doGcm2Faces; Line 43  elseif doGcm2Faces;
43      elseif strcmp(mygrid.fileFormat,'cube');      elseif strcmp(mygrid.fileFormat,'cube');
44          for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;          for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;
45      elseif strcmp(mygrid.fileFormat,'compact');      elseif strcmp(mygrid.fileFormat,'compact');
46                    v00=reshape(v0,[n1*n2 n3*n4*n5]);
47          nn=size(v0,1);          i0=0; i1=0;
48          pp=size(v0,2)/nn;          for iFace=1:mygrid.nFaces;
49          mm=(pp+4-mygrid.nFaces)/4*nn;              i0=i1+1;
50                        nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
51          v00=reshape(v0,[nn*nn*pp n3*n4*n5]);              i1=i1+nn*mm;
52          i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);              v1{iFace}=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]);  
53          end;          end;
           
54      end;      end;
55            
56      v1=gcmfaces(v1);      v1=gcmfaces(v1);
57            
58        if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end;
59        
60  else;  else;
61            
62      [n1,n2,n3,n4,n5]=size(v0{1});      [n1,n2,n3,n4,n5]=size(v0{1});
# Line 73  else; Line 68  else;
68          for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;          for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;
69      elseif strcmp(mygrid.fileFormat,'compact');      elseif strcmp(mygrid.fileFormat,'compact');
70                    
71          v0_faces=v0; clear v0;          %   is there a reason for this?
72          for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;          %         v0_faces=v0; clear v0;
73                    %         for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
74          nn=size(v0{1},1); mm=size(v0{1},2);          
75          pp=mm/nn*4+mygrid.nFaces-4;          if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end;
76            
77          n3=size(v0{1},3); n4=size(v0{1},4); n5=size(v0{1},5);          aa=0;bb=0;
78          v1=NaN*zeros(nn,nn*pp,n3,n4,n5);          for iFace=1:mygrid.nFaces;
79                        aa=aa+prod(mygrid.facesSize(iFace,:));
80          v11=NaN*zeros(nn*nn*pp,n3*n4*n5);              bb=bb+prod(size(v0{iFace}));
           
         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(:,:);  
81          end;          end;
82            bb=bb/aa;
83            v11=NaN*zeros(aa,bb);
84                    
85          v1=reshape(v11,[nn nn*pp n3 n4 n5]);          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            end;
92            nn=mygrid.ioSize(1); mm=mygrid.ioSize(2);        
93            v1=reshape(v11,[nn mm n3 n4 n5]);
94                    
95      end;      end;
96            
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22