anchor.diff:        close(kp)
   c
   600   format(//,'Beginning elastic solution:',/)
  +650   format(//,"Reforming the stiffness matrix:",/)
   700   format('STEP ',i5)
   800   format(/," Total number of equilibrium iterations        = ",i7,/,
        &         " Total number of stiffness matrix reformations = ",i7,/,
  
positive.code: c
        if(lgdefp.eq.izero.and.intord.ne.ithree) then
  cdebug        write(6,*) "Before matinit_drv (1):"
          call matinit_drv(
       &   alnz,ja,nnz,neq,                                               ! sparse
       &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
negative.code:         if(ierr.ne.izero) return
  c
  clater      else if(lgdefp.eq.ione.and.intord.ne.ithree) then
  clater        write(kto,650)
  clater        call matinit_drv(
  clater     &   alnz,ja,nnz,neq,                                               ! sparse
  clater     &   x,d,iwink,wink,numnp,nwink,                                    ! global
  