/[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.6 by gforget, Fri Jul 1 15:18:18 2011 UTC revision 1.10 by gforget, Tue Sep 11 15:15:19 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  %object:    converts model output to gcmfaces object
3  %           or vice versa  %           or vice versa
4  %input:     v0 model output array (resp. gcmfaces object)  %input:     v0 model output array (resp. gcmfaces object)
5  %output:    v1 gcmfaces object (model output array)  %output:    v1 gcmfaces object (model output array)
# Line 12  aa=whos('v0'); doGcm2Faces=strcmp(aa.cla Line 12  aa=whos('v0'); doGcm2Faces=strcmp(aa.cla
12    
13  global mygrid;  global mygrid;
14    
15  if doGcm2Faces;  if doGcm2Faces&mygrid.gcm2facesFast;
16        
17    [n1,n2,n3,n4,n5]=size(v0);      [n1,n2,n3,n4,n5]=size(v0);
18        
19    if strcmp(mygrid.fileFormat,'straight');      v0=reshape(v0,[n1*n2 n3*n4*n5]);
20      v1={v0};      for iFace=1:mygrid.nFaces;
21    elseif strcmp(mygrid.fileFormat,'cube');          v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]);
     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]);  
22      end;      end;
23        
24    end;      v1=gcmfaces(v1);
25        
26    if mygrid.nFaces==1; gridType='ll';  elseif ~doGcm2Faces&mygrid.gcm2facesFast;
27    elseif mygrid.nFaces==5; gridType='llc';      
28    elseif mygrid.nFaces==6; gridType='cube';      [n1,n2,n3,n4,n5]=size(v0{1});
29    end;      
30    v1=gcmfaces(v1,gridType);      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            v00=reshape(v0,[n1*n2 n3*n4*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            end;
54        end;
55        
56        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});
63        
64    if strcmp(mygrid.fileFormat,'straight');      if strcmp(mygrid.fileFormat,'straight');
65      v1=v0{1};          v1=v0{1};
66    elseif strcmp(mygrid.fileFormat,'cube');      elseif strcmp(mygrid.fileFormat,'cube');
67      v1=zeros(n2*6,n2,n3,n4,n5);          v1=zeros(n2*6,n2,n3,n4,n5);
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}));
81              end;
82     i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);          bb=bb/aa;
83     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);          v11=NaN*zeros(aa,bb);
84     i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);          
85     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);          i0=0; i1=0;
86     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);          for iFace=1:mygrid.nFaces;
87     if mygrid.nFaces==6;              i0=i1+1;
88        i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);              nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
89     end;              i1=i1+nn*mm;
90                v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]);
91     v1=reshape(v11,[nn nn*pp n3 n4 n5]);          end;
92            nn=mygrid.ioSize(1); mm=mygrid.ioSize(2);        
93    end;          v1=reshape(v11,[nn mm n3 n4 n5]);
94            
95        end;
96        
97  end;  end;
98    
99    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22