function r = ones(p,varargin)
%overloaded gcmfaces ones function :
%  if the first arguments is a gcmfaces object
%  we call the double ones function for each face data
%  using face dimensions for the first two arguments and
%  passing over the other arguments

r=p;

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