anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine lnklst(lm,lmx,indx,link,nbrs,nee,neq,numel,iwork,
  -     & nsizea,nnz)
  +      subroutine lnklst(neq,lm,lmx,infiel,nconsz,numelt,infetype,indx,
  +     & link,nbrs,iwork,nsizea,nnz,numsn,ierr,errstrng)
   c
   c      program to create a linked list of nonzero row and column entries
   c
  
positive.code:       call ifill(nbrs,izero,iwork)
  c
        inext=ione
        do n=1,numel
  c
  c      check that available storage is not exceeded
  c*     Temporary output to stdout that should be replaced by an
  c*     exception.
  c
          if(inext.gt.iwork) then
            write(6,700) inext-iwork
            stop
          end if
  c
          do i=1,2*nee
            irow=lm(i,n)
            if(i.gt.nee) irow=abs(lmx(i-nee,n))
            if(irow.eq.izero) goto 20
            loc=indx(irow)
            do j=1,2*nee
              icol=lm(j,n)
              if(j.gt.nee) icol=abs(lmx(j-nee,n))
  cdebug            write(6,*) "i,j,irow,icol,loc:",i,j,irow,icol,loc
              if(icol.eq.izero) goto 30
              if(icol.eq.irow) goto 30
              if(loc.eq.izero) then
                loc=inext
                indx(irow)=inext
                nbrs(inext)=icol
                link(inext)=-irow
                inext=inext+1
              else
                loctmp=loc
    40          continue
                ival=nbrs(loctmp)
                if(icol.eq.ival) goto 30
                locsav=loctmp
                loctmp=link(loctmp)
                if(loctmp.gt.izero) goto 40
                link(locsav)=inext
                nbrs(inext)=icol
                link(inext)=-irow
                inext=inext+1
              end if
    30        continue
            end do
    20      continue
          end do
        end do
        nsizea=inext-ione
        nnz=nsizea+neq+ione
   700  format("Working storage exceeded by ",i7)
        return
        end
  c
  c version
  c $Id: lnklst.f,v 1.2 2004/07/07 17:38:21 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c
  c...  local variables
  c
        integer inext,iel,indien,ietype,nee,neeq,i,ii,irow,loc,j,jj,icol
        integer loctmp,ival,locsav
  c
  cdebug      write(6,*) "Hello from lnklst_f!"
  c
  