matlab_class/%40gcmfaces/gcmfaces.m

function s=gcmfaces(varargin);
%object:    create an empty gcmfaces object
%input:     fld cell array containing each face data
%      OR   nFaces number of faces (empty faces in output)
%      OR   (none; results in empty face data, with nFaces=mygrid.nFaces)
%output:    gcmfaces object


%object:    create an empty gcmfaces object
%input:     fld cell array containing each face data
%      OR   nFaces,fld          (to add consistency check)
%      OR   nFaces number of faces (empty faces in output)
%      OR   (none)                 (empty faces in output)
%output:    gcmfaces object formatted/filled accordingly


global mygrid;

if nargin==2;
    nFaces=varargin{1};
    fld=varargin{2};
    if ~iscell(fld)|length(fld)~=nFaces;
        error('inconsistent spec. of nFaces,fld');
    end;
elseif nargin==1;
    tmp1=varargin{1};
    if iscell(tmp1);
        fld=tmp1;
        nFaces=length(fld);
    else;
        nFaces=varargin{1};
        fld=[];
    end;
elseif isfield(mygrid,'nFaces');
    nFaces=mygrid.nFaces;
    fld=[];
else;
    error('nFaces needs to be specified');
end;

if nFaces==1; gridType='ll';
elseif nFaces==4; gridType='llpc';
elseif nFaces==5; gridType='llc';
elseif nFaces==6; gridType='cube';
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;
else;
    s.nFaces=nFaces;
    s.gridType=gridType;
    for iF=1:nFacesMax;
        eval(['s.f' num2str(iF) '=[];']);
    end;
end;

s = class(s,'gcmfaces');

1	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
9	%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;
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;
49
50	if iscell(fld);
51	s.nFaces=length(fld);
52	s.gridType=gridType;
53	for iF=1:s.nFaces;
54	eval(['s.f' num2str(iF) '=fld{iF};']);
55	end;
56	for iF=s.nFaces+1:nFacesMax;
57	eval(['s.f' num2str(iF) '=[];']);
58	end;
59	else;
60	s.nFaces=nFaces;
61	s.gridType=gridType;
62	for iF=1:nFacesMax;
63	eval(['s.f' num2str(iF) '=[];']);
64	end;
65	end;
66
67	s = class(s,'gcmfaces');
68