anchor.diff:            sdmin=min(sdmin,sd(i))
           end do
           sd(nd)=half*sdmin
  -        call lfit(xx,yy,sd,nd,a,ma,cov,ncvm,xl)
  +        call lfit(xx,yy,sd,nd,a,ma,cov,ncvm,xl,ierr,errstrng)
  +        if(ierr.ne.izero) return
           amag=one/sqrt(a(1)*a(1)+a(2)*a(2)+a(3)*a(3))
           a(1)=a(1)*amag
           a(2)=a(2)*amag
  
positive.code:       end
  c
  c version
  c $Id: skcomp.f,v 1.3 2004/07/13 16:40:58 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code:       double precision dnm,rat
        double precision a(3),xx(5),yy(5),sd(5),xl(3,5),cov(3,3)
  c
  cdebug      write(6,*) "Hello from skcomp_f!"
  c
        ma=ithree
        nd=ifive
  