/[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.10 by gforget, Tue Sep 11 15:15:19 2012 UTC revision 1.13 by gforget, Fri Mar 17 17:41:56 2017 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)
6  %  %
7  %note:      global mygrid parameters (nFaces, fileFormat) are used  %note:      global mygrid parameters (nFaces, fileFormat) are used
8    
9  input_list_check('convert2gcmfaces',nargin);  input_list_check('convert2gcmfaces',nargin);
10    
11  aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');  aa=whos('v0');
12    if strcmp(aa.class,'single'); v0=double(v0); end;
13  global mygrid;  aa=whos('v0');
14    doGcm2Faces=strcmp(aa.class,'double');
15  if doGcm2Faces&mygrid.gcm2facesFast;  
16        global mygrid;
17      [n1,n2,n3,n4,n5]=size(v0);  
18        if doGcm2Faces&mygrid.gcm2facesFast;
19      v0=reshape(v0,[n1*n2 n3*n4*n5]);      
20      for iFace=1:mygrid.nFaces;      [n1,n2,n3,n4,n5]=size(v0);
21          v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]);      
22      end;      v0=reshape(v0,[n1*n2 n3*n4*n5]);
23            for iFace=1:mygrid.nFaces;
24      v1=gcmfaces(v1);          v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]);
25            end;
26  elseif ~doGcm2Faces&mygrid.gcm2facesFast;      
27            v1=gcmfaces(v1);
28      [n1,n2,n3,n4,n5]=size(v0{1});      
29        elseif ~doGcm2Faces&mygrid.gcm2facesFast;
30      v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5);      
31      for iFace=1:mygrid.nFaces;      [n1,n2,n3,n4,n5]=size(v0{1});
32          nn=size(mygrid.XC{iFace});      
33          v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]);      v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5);
34      end;      for iFace=1:mygrid.nFaces;
35      v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]);          nn=size(mygrid.XC{iFace});
36                v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]);
37  elseif doGcm2Faces;      end;
38            v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]);
39      [n1,n2,n3,n4,n5]=size(v0);      
40        elseif doGcm2Faces;
41      if strcmp(mygrid.fileFormat,'straight');      
42          v1={v0};      [n1,n2,n3,n4,n5]=size(v0);
43      elseif strcmp(mygrid.fileFormat,'cube');      
44          for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;      if strcmp(mygrid.fileFormat,'straight');
45      elseif strcmp(mygrid.fileFormat,'compact');          v1={v0};
46          v00=reshape(v0,[n1*n2 n3*n4*n5]);      elseif strcmp(mygrid.fileFormat,'cube');
47          i0=0; i1=0;          for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;
48          for iFace=1:mygrid.nFaces;      elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles');
49              i0=i1+1;          v00=reshape(v0,[n1*n2 n3*n4*n5]);
50              nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);          i0=0; i1=0;
51              i1=i1+nn*mm;          for iFace=1:mygrid.nFaces;
52              v1{iFace}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);              i0=i1+1;
53          end;              nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
54      end;              i1=i1+nn*mm;
55                    v1{iFace}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
56      v1=gcmfaces(v1);          end;
57            elseif strcmp(mygrid.fileFormat,'native');
58      if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end;          v1=convert2array(v0);
59            end;
60  else;      
61            if ~strcmp(mygrid.fileFormat,'native'); v1=gcmfaces(v1); end;
62      [n1,n2,n3,n4,n5]=size(v0{1});      
63            if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end;
64      if strcmp(mygrid.fileFormat,'straight');      
65          v1=v0{1};  else;
66      elseif strcmp(mygrid.fileFormat,'cube');      
67          v1=zeros(n2*6,n2,n3,n4,n5);      [n1,n2,n3,n4,n5]=size(v0{1});
68          for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;      
69      elseif strcmp(mygrid.fileFormat,'compact');      if strcmp(mygrid.fileFormat,'straight');
70                    v1=v0{1};
71          %   is there a reason for this?      elseif strcmp(mygrid.fileFormat,'cube');
72          %         v0_faces=v0; clear v0;          v1=zeros(n2*6,n2,n3,n4,n5);
73          %         for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;          for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;
74                elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles');
75          if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end;          
76            %   is there a reason for this?
77          aa=0;bb=0;          %         v0_faces=v0; clear v0;
78          for iFace=1:mygrid.nFaces;          %         for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
79              aa=aa+prod(mygrid.facesSize(iFace,:));          
80              bb=bb+prod(size(v0{iFace}));          if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end;
81          end;  
82          bb=bb/aa;          aa=0;bb=0;
83          v11=NaN*zeros(aa,bb);          for iFace=1:mygrid.nFaces;
84                        aa=aa+prod(mygrid.facesSize(iFace,:));
85          i0=0; i1=0;              bb=bb+prod(size(v0{iFace}));
86          for iFace=1:mygrid.nFaces;          end;
87              i0=i1+1;          bb=bb/aa;
88              nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);          v11=NaN*zeros(aa,bb);
89              i1=i1+nn*mm;          
90              v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]);          i0=0; i1=0;
91          end;          for iFace=1:mygrid.nFaces;
92          nn=mygrid.ioSize(1); mm=mygrid.ioSize(2);                      i0=i1+1;
93          v1=reshape(v11,[nn mm n3 n4 n5]);              nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2);
94                        i1=i1+nn*mm;
95      end;              v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]);
96                end;
97  end;          nn=mygrid.ioSize(1); mm=mygrid.ioSize(2);        
98            v1=reshape(v11,[nn mm n3 n4 n5]);
99    
100        elseif strcmp(mygrid.fileFormat,'native');
101            v1=convert2array(v0);        
102        end;
103        
104    end;
105    
106    
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22