/[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.4 by gforget, Wed Nov 3 19:53:13 2010 UTC revision 1.6 by gforget, Fri Jul 1 15:18:18 2011 UTC
# Line 1  Line 1 
1  function [v1]=convert2gcmfaces(v0,varargin);  function [v1]=convert2gcmfaces(v0,varargin);
2    %object:    converts model output to gcmfaces object
3    %           or vice versa
4    %input:     v0 model output array (resp. gcmfaces object)
5    %output:    v1 gcmfaces object (model output array)
6    %
7    %note:      global mygrid parameters (nFaces, fileFormat) are used
8    
9  global fileFormat;  input_list_check('convert2gcmfaces',nargin);
 if isempty(fileFormat);  
  global gcmfaces_verbose;  
  if gcmfaces_verbose;  
    fprintf('\nconvert2gcmfaces.m init: there are several supported file conventions. \n');  
    fprintf('  By default gcmfaces assumes MITgcm type binary formats as follows: \n')  
    fprintf('  (1 face) straight global format; (4 or 5 faces) compact global format\n');  
    fprintf('  (6 faces) cube format with one face after the other. \n');  
    fprintf('  If this is inadequate, you can change the format below.\n\n');  
  end;  
  if nargin==1; nFaces=5; else; nFaces=varargin{1}; end;  
  if nFaces==1;          fileFormat='straight';  
  elseif nFaces==6;      fileFormat='cube';  
  else;                  fileFormat='compact';  
  end;  
 end;  
10    
11  aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');  aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');
12    
13  if doGcm2Faces;  global mygrid;
14    
15    if nargin==1; nFaces=5; else; nFaces=varargin{1}; end;  if doGcm2Faces;
16    
17    [n1,n2,n3,n4]=size(v0);    [n1,n2,n3,n4,n5]=size(v0);
18    
19    if strcmp(fileFormat,'straight');    if strcmp(mygrid.fileFormat,'straight');
20      v1={v0};      v1={v0};
21    elseif strcmp(fileFormat,'cube');    elseif strcmp(mygrid.fileFormat,'cube');
22      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;
23    elseif strcmp(fileFormat,'compact');    elseif strcmp(mygrid.fileFormat,'compact');
24    
25      nn=size(v0,1);      nn=size(v0,1);
26      pp=size(v0,2)/nn;      pp=size(v0,2)/nn;
27      mm=(pp+4-nFaces)/4*nn;      mm=(pp+4-mygrid.nFaces)/4*nn;
28    
29      v00=reshape(v0,[nn*nn*pp n3*n4]);      v00=reshape(v0,[nn*nn*pp n3*n4*n5]);
30      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]);
31      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]);
32      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]);
33      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]);
34      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]);
35      if nFaces==6;      if mygrid.nFaces==6;
36         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]);
37      end;      end;
38    
39    end;    end;
40    
41    if nFaces==1; gridType='ll'; elseif nFaces==5; gridType='llc'; elseif nFaces==6; gridType='cube'; end;    if mygrid.nFaces==1; gridType='ll';
42      elseif mygrid.nFaces==5; gridType='llc';
43      elseif mygrid.nFaces==6; gridType='cube';
44      end;
45    v1=gcmfaces(v1,gridType);    v1=gcmfaces(v1,gridType);
46    
47  else;  else;
48    
49    nFaces=get(v0,'nFaces');    [n1,n2,n3,n4,n5]=size(v0{1});
   
   [n1,n2,n3,n4]=size(v0{1});  
50    
51    if strcmp(fileFormat,'straight');    if strcmp(mygrid.fileFormat,'straight');
52      v1=v0{1};      v1=v0{1};
53    elseif strcmp(fileFormat,'cube');    elseif strcmp(mygrid.fileFormat,'cube');
54      v1=zeros(n2*6,n2,n3,n4);      v1=zeros(n2*6,n2,n3,n4,n5);
55      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;
56    elseif strcmp(fileFormat,'compact');    elseif strcmp(mygrid.fileFormat,'compact');
57    
58      v0_faces=v0; clear v0;      v0_faces=v0; clear v0;
59      for iFace=1:nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;      for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
60    
61      nn=size(v0{1},1); mm=size(v0{1},2);      nn=size(v0{1},1); mm=size(v0{1},2);
62      pp=mm/nn*4+nFaces-4;      pp=mm/nn*4+mygrid.nFaces-4;
63    
64      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);
65      v1=NaN*zeros(nn,nn*pp,n3,n4);      v1=NaN*zeros(nn,nn*pp,n3,n4,n5);
66    
67      v11=NaN*zeros(nn*nn*pp,n3*n4);      v11=NaN*zeros(nn*nn*pp,n3*n4*n5);
68        
69     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(:,:);
70     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(:,:);
71     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(:,:);
72     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(:,:);
73     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(:,:);
74     if nFaces==6;     if mygrid.nFaces==6;
75        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(:,:);
76     end;     end;
77    
78     v1=reshape(v11,[nn nn*pp n3 n4]);     v1=reshape(v11,[nn nn*pp n3 n4 n5]);
79    
80    end;    end;
81    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22