matlab_class/%40gcmfaces/nanmax.m

function r = nanmax(p,varargin)
%overloaded gcmfaces nanmax function :
%  1) if single gcmfaces argument, then returns the global nanmax over all faces
%  2) if two gcmfaces arguments, then returns the nanmax of the two at each point
%  3) otherwise calls double nanmax 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=nanmax(tmp1);
   return;
end;

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


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


1	gforget	1.1	function r = nanmax(p,varargin)
2	gforget	1.2	%overloaded gcmfaces nanmax function :
3			% 1) if single gcmfaces argument, then returns the global nanmax over all faces
4	gforget	1.3	% 2) if two gcmfaces arguments, then returns the nanmax of the two at each point
5			% 3) otherwise calls double nanmax function for each face, passing over the other arguments
6	gforget	1.1
7			if nargin==1;
8			tmp1=[];
9			for iFace=1:p.nFaces;
10			iF=num2str(iFace);
11			eval(['tmp1=[tmp1;p.f' iF '(:)];']);
12			end;
13			r=nanmax(tmp1);
14			return;
15			end;
16
17	gforget	1.3	if isa(varargin{1},'gcmfaces');
18			r=p;
19			for iFace=1:r.nFaces;
20			iF=num2str(iFace);
21			eval(['r.f' iF '=nanmax(p.f' iF ',varargin{1}.f' iF ');']);
22			end;
23			return;
24			end;
25
26
27	gforget	1.1	r=p;
28			for iFace=1:r.nFaces;
29			iF=num2str(iFace);
30			eval(['r.f' iF '=nanmax(p.f' iF ',varargin{:});']);
31			end;
32
33