anchor.diff:        end do
         nmin=nmin+1
         nmax=nmax+1
  -      if(neq.ne.0) wavg=wavg/dble(neq)
  +      if(neq.ne.izero) wavg=wavg/dble(neq)
         return
         end
   c
   c version
  -c $Id: makemsr.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  +c $Id: makemsr.f,v 1.2 2004/07/07 18:38:30 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
        include "implicit.inc"
  c
  c...  subroutine arguments
  c
        integer neq,nnz,iwork,nmin,nmax
        integer ja(nnz),indx(neq),link(iwork),nbrs(iwork)
        double precision wavg
  c
  c...  defined constants
  c
        include "nconsts.inc"
  c
  c...  intrinsic functions
  c
  
negative.code: c
        call ifill(ja,izero,nnz)
        ja(1)=neq+2
        nmin=ione
        nmax=ione
        wavg=neq
        do i=1,neq
          loc=indx(i)
          ncol=izero
   20     continue
          ja(ja(i)+ncol)=nbrs(loc)
          ncol=ncol+ione
          loc=link(loc)
          if(loc.gt.izero) goto 20
          nmin=min(nmin,ncol)
          nmax=max(nmax,ncol)
          wavg=wavg+dble(ncol)
  