lanarrow.m

%======================================================================
%                    L A N A R R O W . M 
%                    doc: Sat Jul 17 21:21:54 2004
%                    dlm: Fri Mar  5 15:48:43 2010
%                    (c) 2004 ladcp@
%                    uE-Info: 26 22 NIL 0 0 72 0 2 4 NIL ofnI
%======================================================================

% CHANGES BY ANT:
%  Jul 17, 2004: - adapted to new return params of [getinv.m]

% ========

   ps1=ps;
   ps1.down_up=0;
   ps1.solve=0;
   for i=1:ceil(ps.outlier)
    if exist('dr')==1
     [p,dr,ps1,de,der]=getinv(di,p,ps1,dr);
    else
     [p,dr,ps1,de,der]=getinv(di,p,ps1);
    end
    dif=(di.ru-der.ru_oce-der.ru_ctd).^2+...
         (di.rv-der.rv_oce-der.rv_ctd).^2;
    [es,ii]=sort(dif(:));
    iok=find(isfinite(es));
    ln=length(iok)*0.01;
    if ln>0
     di.weight(ii(iok(end-[0:ln])))=NaN;
     disp([' give low weight to 1% of data scan#:',int2str(i)])
    end
    pause(0.01)
   end
   clear ps1