anchor.diff:  clater        if(ierr.ne.izero) return
   c
   clater        call formdf_ld(
  -clater     &   b,neq,                                                         ! force
  +clater     &   bintern,neq,                                                   ! force
   clater     &   x,d,dcur,numnp,                                                ! global
   clater     &   s,stemp,                                                       ! stiff
   clater     &   dmat,ien,lm,lmx,infiel,iddmat,ndmatsz,numelt,nconsz,           ! elemnt
  
positive.code:         if(ierr.ne.izero) return
        end if
  c
  c...transfer boundary conditions into global load vector btot(neq)
  c   and displacement increment vector deld(ndof,numnp)
  c
        call load(id,ibond,bond,dcur,deld,btot,histry,deltp,numnp,neq,
       & nhist,nstep,lastep,ierr,errstrng)
  cdebug      write(6,*) "After load:"
        if(ierr.ne.izero) return
  c
  
negative.code:       end
  c
  c version
  c $Id: elastc.f,v 1.12 2005/01/05 22:13:53 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  