matlab_class/%40gcmfaces/gcmfaces.m

function s=gcmfaces(varargin);

if nargin==0; fld=5; gridType='llc'; 
elseif nargin==1; fld=varargin{1}; gridType='llc';
elseif nargin==2; fld=varargin{1}; gridType=varargin{2}; 
else; error('wrong gcmfaces definition'); end; 

nFacesMax=6;

if iscell(fld);
   s.nFaces=length(fld);
   s.gridType=gridType;
   for iF=1:s.nFaces;
      eval(['s.f' num2str(iF) '=fld{iF};']);
   end;
   for iF=s.nFaces+1:nFacesMax;   
      eval(['s.f' num2str(iF) '=[];']);
   end;
elseif isreal(fld);
   s.nFaces=fld;
   s.gridType=gridType;
   for iF=1:nFacesMax;
      eval(['s.f' num2str(iF) '=[];']);
   end;
else;
   error('wrong gcmfaces definition');
end;

s = class(s,'gcmfaces');

1	function s=gcmfaces(varargin);
2
3	if nargin==0; fld=5; gridType='llc';
4	elseif nargin==1; fld=varargin{1}; gridType='llc';
5	elseif nargin==2; fld=varargin{1}; gridType=varargin{2};
6	else; error('wrong gcmfaces definition'); end;
7
8	nFacesMax=6;
9
10	if iscell(fld);
11	s.nFaces=length(fld);
12	s.gridType=gridType;
13	for iF=1:s.nFaces;
14	eval(['s.f' num2str(iF) '=fld{iF};']);
15	end;
16	for iF=s.nFaces+1:nFacesMax;
17	eval(['s.f' num2str(iF) '=[];']);
18	end;
19	elseif isreal(fld);
20	s.nFaces=fld;
21	s.gridType=gridType;
22	for iF=1:nFacesMax;
23	eval(['s.f' num2str(iF) '=[];']);
24	end;
25	else;
26	error('wrong gcmfaces definition');
27	end;
28
29	s = class(s,'gcmfaces');
30