matlab_class/%40gcmfaces/cut_T_N.m

function r = cut_T_N(p,varargin)
%object :  remove N points at the edge of each face (opposite of exch_T_N).
%input:    b is a gcmfaces object
%optional: N (1 by default) is the halo region width
%output:   a is the reduced gcmfaces object
%

r=p;

if nargin==1; N=1; else; N=varargin{1}; end;

for iFace=1:r.nFaces;
   iF=num2str(iFace); 
   eval(['r.f' iF '=p.f' iF '(1+N:end-N,1+N:end-N,:,:);']); 
end;


1	gforget	1.1	function r = cut_T_N(p,varargin)
2	gforget	1.2	%object : remove N points at the edge of each face (opposite of exch_T_N).
3			%input: b is a gcmfaces object
4			%optional: N (1 by default) is the halo region width
5			%output: a is the reduced gcmfaces object
6			%
7	gforget	1.1
8			r=p;
9
10			if nargin==1; N=1; else; N=varargin{1}; end;
11
12			for iFace=1:r.nFaces;
13			iF=num2str(iFace);
14			eval(['r.f' iF '=p.f' iF '(1+N:end-N,1+N:end-N,:,:);']);
15			end;
16
17