/[MITgcm]/MITgcm_contrib/gael/matlab_class/@gcmfaces/gcmfaces.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/matlab_class/@gcmfaces/gcmfaces.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by gforget, Wed Feb 10 14:43:47 2010 UTC revision 1.4 by gforget, Tue Jul 31 19:08:43 2012 UTC
# Line 1  Line 1 
1  function s=gcmfaces(varargin);  function s=gcmfaces(varargin);
2    %object:    create an empty gcmfaces object
3    %input:     fld cell array containing each face data
4    %      OR   nFaces number of faces (empty faces in output)
5    %      OR   (none; results in empty face data, with nFaces=mygrid.nFaces)
6    %output:    gcmfaces object
7    
8  if nargin==0; fld=5; gridType='llc';  
9  elseif nargin==1; fld=varargin{1}; gridType='llc';  %object:    create an empty gcmfaces object
10  elseif nargin==2; fld=varargin{1}; gridType=varargin{2};  %input:     fld cell array containing each face data
11  else; error('wrong gcmfaces definition'); end;  %      OR   nFaces,fld          (to add consistency check)
12    %      OR   nFaces number of faces (empty faces in output)
13    %      OR   (none)                 (empty faces in output)
14    %output:    gcmfaces object formatted/filled accordingly
15    
16    
17    global mygrid;
18    
19    if nargin==2;
20        nFaces=varargin{1};
21        fld=varargin{2};
22        if ~iscell(fld)|length(fld)~=nFaces;
23            error('inconsistent spec. of nFaces,fld');
24        end;
25    elseif nargin==1;
26        tmp1=varargin{1};
27        if iscell(tmp1);
28            fld=tmp1;
29            nFaces=length(fld);
30        else;
31            nFaces=varargin{1};
32            fld=[];
33        end;
34    elseif isfield(mygrid,'nFaces');
35        nFaces=mygrid.nFaces;
36        fld=[];
37    else;
38        error('nFaces needs to be specified');
39    end;
40    
41    if nFaces==1; gridType='ll';
42    elseif nFaces==4; gridType='llpc';
43    elseif nFaces==5; gridType='llc';
44    elseif nFaces==6; gridType='cube';
45    else; error('wrong gcmfaces definition');
46    end;
47    
48  nFacesMax=6;  nFacesMax=6;
49    
50  if iscell(fld);  if iscell(fld);
51     s.nFaces=length(fld);      s.nFaces=length(fld);
52     s.gridType=gridType;      s.gridType=gridType;
53     for iF=1:s.nFaces;      for iF=1:s.nFaces;
54        eval(['s.f' num2str(iF) '=fld{iF};']);          eval(['s.f' num2str(iF) '=fld{iF};']);
55     end;      end;
56     for iF=s.nFaces+1:nFacesMax;        for iF=s.nFaces+1:nFacesMax;
57        eval(['s.f' num2str(iF) '=[];']);          eval(['s.f' num2str(iF) '=[];']);
58     end;      end;
 elseif isreal(fld);  
    s.nFaces=fld;  
    s.gridType=gridType;  
    for iF=1:nFacesMax;  
       eval(['s.f' num2str(iF) '=[];']);  
    end;  
59  else;  else;
60     error('wrong gcmfaces definition');      s.nFaces=nFaces;
61        s.gridType=gridType;
62        for iF=1:nFacesMax;
63            eval(['s.f' num2str(iF) '=[];']);
64        end;
65  end;  end;
66    
67  s = class(s,'gcmfaces');  s = class(s,'gcmfaces');
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22