matlab_class/%40gcmfaces/min.m

function r = min(p,varargin)
% MIN(p,varargin)
%
%overloaded gcmfaces min function :
%  1) if single gcmfaces argument, then returns the global min over all faces
%  2) if two gcmfaces arguments, then returns the min of the two at each point
%  3) otherwise calls double min function for each face, passing over the other arguments

if nargin==1;
   tmp1=[];
   for iFace=1:p.nFaces;
      iF=num2str(iFace);
      eval(['tmp1=[tmp1;p.f' iF '(:)];']);
   end;
   r=min(tmp1);
   return;
end;

if isa(varargin{1},'gcmfaces');
   r=p;
   for iFace=1:r.nFaces;
      iF=num2str(iFace);
      eval(['r.f' iF '=min(p.f' iF ',varargin{1}.f' iF ');']);
   end;
   return;
end;

if varargin{2}>0;
  r=p;
  for iFace=1:r.nFaces;
     iF=num2str(iFace); 
     eval(['r.f' iF '=min(p.f' iF ',varargin{:});']); 
  end;
else;
  tmp1=convert2gcmfaces(p);
  [n1,n2,n3,n4]=size(tmp1);
  tmp1=reshape(tmp1,n1*n2,n3,n4);
  r=min(tmp1,[],1);
end;


1	function r = min(p,varargin)
2	% MIN(p,varargin)
3	%
4	%overloaded gcmfaces min function :
5	% 1) if single gcmfaces argument, then returns the global min over all faces
6	% 2) if two gcmfaces arguments, then returns the min of the two at each point
7	% 3) otherwise calls double min function for each face, passing over the other arguments
8
9	if nargin==1;
10	tmp1=[];
11	for iFace=1:p.nFaces;
12	iF=num2str(iFace);
13	eval(['tmp1=[tmp1;p.f' iF '(:)];']);
14	end;
15	r=min(tmp1);
16	return;
17	end;
18
19	if isa(varargin{1},'gcmfaces');
20	r=p;
21	for iFace=1:r.nFaces;
22	iF=num2str(iFace);
23	eval(['r.f' iF '=min(p.f' iF ',varargin{1}.f' iF ');']);
24	end;
25	return;
26	end;
27
28	if varargin{2}>0;
29	r=p;
30	for iFace=1:r.nFaces;
31	iF=num2str(iFace);
32	eval(['r.f' iF '=min(p.f' iF ',varargin{:});']);
33	end;
34	else;
35	tmp1=convert2gcmfaces(p);
36	[n1,n2,n3,n4]=size(tmp1);
37	tmp1=reshape(tmp1,n1*n2,n3,n4);
38	r=min(tmp1,[],1);
39	end;
40
41