/[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.5 by gforget, Tue Jun 21 21:18:48 2011 UTC revision 1.7 by gforget, Tue Jul 5 17:51:13 2011 UTC
# 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]=size(v0);    [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      if mygrid.nFaces==1; gridType='ll';
25      elseif mygrid.nFaces==5; gridType='llc';
26      elseif mygrid.nFaces==6; gridType='cube';
27      end;
28      v1=gcmfaces(v1,gridType);
29    
30    elseif ~doGcm2Faces&mygrid.gcm2facesFast;
31    
32      [n1,n2,n3,n4,n5]=size(v0{1});
33    
34      v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5);
35      for iFace=1:mygrid.nFaces;
36        nn=size(mygrid.XC{iFace});
37        v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]);
38      end;
39      v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]);
40    
41    elseif doGcm2Faces;
42    
43      [n1,n2,n3,n4,n5]=size(v0);
44    
45    if strcmp(mygrid.fileFormat,'straight');    if strcmp(mygrid.fileFormat,'straight');
46      v1={v0};      v1={v0};
47    elseif strcmp(mygrid.fileFormat,'cube');    elseif strcmp(mygrid.fileFormat,'cube');
48      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;
49    elseif strcmp(mygrid.fileFormat,'compact');    elseif strcmp(mygrid.fileFormat,'compact');
50    
51      nn=size(v0,1);      nn=size(v0,1);
52      pp=size(v0,2)/nn;      pp=size(v0,2)/nn;
53      mm=(pp+4-mygrid.nFaces)/4*nn;      mm=(pp+4-mygrid.nFaces)/4*nn;
54    
55      v00=reshape(v0,[nn*nn*pp n3*n4]);      v00=reshape(v0,[nn*nn*pp n3*n4*n5]);
56      i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4]);      i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
57      i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4]);      i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
58      i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4]);      i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
59      i0=i1+1; i1=i1+nn*mm; v1{4}=reshape(v00(i0:i1,:),[mm nn n3 n4]);      i0=i1+1; i1=i1+nn*mm; v1{4}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
60      i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4]);      i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
61      if mygrid.nFaces==6;      if mygrid.nFaces==6;
62         i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4]);         i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
63      end;      end;
64    
65    end;    end;
# Line 46  if doGcm2Faces; Line 72  if doGcm2Faces;
72    
73  else;  else;
74    
75    [n1,n2,n3,n4]=size(v0{1});    [n1,n2,n3,n4,n5]=size(v0{1});
76    
77    if strcmp(mygrid.fileFormat,'straight');    if strcmp(mygrid.fileFormat,'straight');
78      v1=v0{1};      v1=v0{1};
79    elseif strcmp(mygrid.fileFormat,'cube');    elseif strcmp(mygrid.fileFormat,'cube');
80      v1=zeros(n2*6,n2,n3,n4);      v1=zeros(n2*6,n2,n3,n4,n5);
81      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;
82    elseif strcmp(mygrid.fileFormat,'compact');    elseif strcmp(mygrid.fileFormat,'compact');
83    
84      v0_faces=v0; clear v0;      v0_faces=v0; clear v0;
# Line 61  else; Line 87  else;
87      nn=size(v0{1},1); mm=size(v0{1},2);      nn=size(v0{1},1); mm=size(v0{1},2);
88      pp=mm/nn*4+mygrid.nFaces-4;      pp=mm/nn*4+mygrid.nFaces-4;
89    
90      n3=size(v0{1},3); n4=size(v0{1},4);      n3=size(v0{1},3); n4=size(v0{1},4); n5=size(v0{1},5);
91      v1=NaN*zeros(nn,nn*pp,n3,n4);      v1=NaN*zeros(nn,nn*pp,n3,n4,n5);
92    
93      v11=NaN*zeros(nn*nn*pp,n3*n4);      v11=NaN*zeros(nn*nn*pp,n3*n4*n5);
94        
95     i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4]); v11(i0:i1,:)=tmp1(:,:);     i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
96     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4]); v11(i0:i1,:)=tmp1(:,:);     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
97     i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);     i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
98     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
99     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);     i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
100     if mygrid.nFaces==6;     if mygrid.nFaces==6;
101        i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4]); v11(i0:i1,:)=tmp1(:,:);        i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
102     end;     end;
103    
104     v1=reshape(v11,[nn nn*pp n3 n4]);     v1=reshape(v11,[nn nn*pp n3 n4 n5]);
105    
106    end;    end;
107    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22