anchor.diff:  c
   c...  included dimension and type statements
   c
  +      include "nforce_dim.inc"
         include "nsysdat_dim.inc"
         include "npar_dim.inc"
         include "ntimdat_dim.inc"
  
positive.code:         if(ierr.ne.izero) return
  c
          call formdf_ss(
       &   b,neq,                                                         ! force
       &   x,d,dcur,numnp,                                                ! global
       &   s,stemp,                                                       ! stiff
       &   dmat,ien,lm,lmx,infiel,iddmat,ndmatsz,numelt,nconsz,           ! elemnt
  
negative.code: clater        if(ierr.ne.izero) return
  c
  clater        call formdf_ld(
  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
  