anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine read_winkx(winkx,wxscal,iwinkx,idx,numnp,ndof,nwinkx,
  -     & nwinkxe,kr,kw,idout,ierr,wxfile,ofile)
  +      subroutine read_winkx(winkx,wxscal,iwinkx,idx,numnp,nwinkx,
  +     & nwinkxe,kr,kw,idout,wxfile,ofile,ierr,errstrng)
   c
   c...  program for reading and printing data on differential winkler
   c     restoring forces
  
positive.code: c
   20   continue
          ierr=1
          close(kr)
          return
  c
  
negative.code:             winkx(nwxtot)=wxtmp(j)
            end if
          end do
          if(idout.gt.izero.and.nnz.ne.izero) then
            if(mod(nlines,npage).eq.izero) then
              write(kw,2000,err=50)
              write(kw,3000,err=50) (labeld(il),il=1,ndof)
              write(kw,2005,err=50)
            end if
            nlines=nlines+ione
            write(kw,2020,err=50) n,(wxtmp(il),iwxtmp(il),il=1,ndof)
          end if 
        end do
        close(kr)
        if(idout.gt.izero) close(kw)
  c
  c...  normal return
  c
  