/[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.2 by gforget, Thu Apr 15 23:15:11 2010 UTC
# Line 13  if doGcm2Faces; Line 13  if doGcm2Faces;
13      mm=(pp+4-nFaces)/4*nn;      mm=(pp+4-nFaces)/4*nn;
14    
15      n3=size(v0,3); n4=size(v0,4);      n3=size(v0,3); n4=size(v0,4);
     for i3=1:n3; for i4=1:n4;  
       v00=v0(:,:,i3,i4);  
16    
17        v00=v00(:);      v00=reshape(v0,[nn*nn*pp n3*n4]);
18        i0=1; i1=nn*mm; v1{1}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm);      i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4]);
19        i0=i1+1; i1=i1+nn*mm; v1{2}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm);      i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4]);
20        i0=i1+1; i1=i1+nn*nn; v1{3}(:,:,i3,i4)=reshape(v00(i0:i1),nn,nn);      i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4]);
21        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{4}=reshape(v00(i0:i1,:),[mm nn n3 n4]);
22        i0=i1+1; i1=i1+nn*mm; v1{5}(:,:,i3,i4)=reshape(v00(i0:i1),mm,nn);      i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4]);
23        if nFaces==6;      if nFaces==6;
24           i0=i1+1; i1=i1+nn*nn;         i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4]);
25        end;      end;
26    
     end; end;  
27    else;    else;
28      v1={v0};      v1={v0};
29    end;    end;
# Line 48  else; Line 45  else;
45    
46      n3=size(v0{1},3); n4=size(v0{1},4);      n3=size(v0{1},3); n4=size(v0{1},4);
47      v1=NaN*zeros(nn,nn*pp,n3,n4);      v1=NaN*zeros(nn,nn*pp,n3,n4);
     for i3=1:n3; for i4=1:n4;  
48    
49        v11=NaN*zeros(nn*nn*pp,1);      v11=NaN*zeros(nn*nn*pp,n3*n4);
50          
51        i0=1; i1=nn*mm; tmp1=v0{1}(:,:,i3,i4); v11(i0:i1)=tmp1(:);     i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4]); v11(i0:i1,:)=tmp1(:,:);
52        i0=i1+1; i1=i1+nn*mm; tmp1=v0{2}(:,:,i3,i4); v11(i0:i1)=tmp1(:);     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4]); v11(i0:i1,:)=tmp1(:,:);
53        i0=i1+1; i1=i1+nn*nn; tmp1=v0{3}(:,:,i3,i4); v11(i0:i1)=tmp1(:);     i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);
54        i0=i1+1; i1=i1+nn*mm; tmp1=v0{4}(:,:,i3,i4); v11(i0:i1)=tmp1(:);     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);
55        i0=i1+1; i1=i1+nn*mm; tmp1=v0{5}(:,:,i3,i4); v11(i0:i1)=tmp1(:);     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);
56        if nFaces==6;     if nFaces==6;
57           i0=i1+1; i1=i1+nn*nn;        i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);
58        end;     end;
59    
60        v1(:,:,i3,i4)=reshape(v11,nn,nn*pp);     v1=reshape(v11,[nn nn*pp n3 n4]);
     end; end;  
61    
62    else;    else;
63      v1=v0{1};      v1=v0{1};
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22