anchor.diff:  c
    20   continue
           ierr=1
  +        errstrng="read_winkx"
           close(kr)
           return
   c
  
positive.code: c
   30   continue
          ierr=3
          close(kr)
          return
  c
   2005 format(/)
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        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
  