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