/[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.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,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);    [n1,n2,n3,n4,n5]=size(v0);
44    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22