/[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.2 by gforget, Tue Jun 21 21:17:31 2011 UTC revision 1.5 by gforget, Tue Jan 12 14:37:03 2016 UTC
# Line 1  Line 1 
1  function s=gcmfaces(varargin);  function s=gcmfaces(varargin);
2  %object:    create a gcmfaces object  %object:    create an empty gcmfaces object
3  %input:     cell array (optional)  %input:     fld cell array containing each face data
4  %  %      OR   nFaces number of faces (empty faces in output)
5  %note:      - the current format is determined based on the global mygrid.  %      OR   (none; results in empty face data, with nFaces=mygrid.nFaces)
6  %           - if no input is specified, create an empty gcmfaces object.  %output:    gcmfaces object
7  %           - otherwise each record of the cell array must provide the data for one face.  
8  %           - this routine used to have several inputs of various types;  
9  %           now only the first input is used, and only if it is a cell array.  %object:    create an empty gcmfaces object
10    %input:     fld cell array containing each face data
11    %      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;  global mygrid;
18    
19  nFaces=mygrid.nFaces;  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        nFaces=1;
39        fld=[];
40        warning('nFaces set to 1 by default');
41    end;
42    
43  if nFaces==1; gridType='ll';  if nFaces==1; gridType='ll';
44  elseif nFaces==4; gridType='llpc';  elseif nFaces==4; gridType='llpc';
# Line 21  end; Line 49  end;
49    
50  nFacesMax=6;  nFacesMax=6;
51    
 if nargin>1; fld=varargin{1}; else; fld=[]; end;  
   
52  if iscell(fld);  if iscell(fld);
53     s.nFaces=length(fld);      s.nFaces=length(fld);
54     s.gridType=gridType;      s.gridType=gridType;
55     for iF=1:s.nFaces;      for iF=1:s.nFaces;
56        eval(['s.f' num2str(iF) '=fld{iF};']);          eval(['s.f' num2str(iF) '=fld{iF};']);
57     end;      end;
58     for iF=s.nFaces+1:nFacesMax;      for iF=s.nFaces+1:nFacesMax;
59        eval(['s.f' num2str(iF) '=[];']);          eval(['s.f' num2str(iF) '=[];']);
60     end;      end;
61  else;  else;
62      s.nFaces=nFaces;      s.nFaces=nFaces;
63      s.gridType=gridType;      s.gridType=gridType;
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22