--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/lanarrow.m	Tue Oct 20 16:23:49 2009 -0400
@@ -0,0 +1,35 @@
+%======================================================================
+%                    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
+