anchor.diff:  c
           if(ierr.ne.izero) return
   c
  -      else if(lgdefp.eq.izero.and.ibbarp.eq.ione) then
  +      else if(lgdefp.eq.izero.and.intord.eq.ithree) then
           call matinit_drv(
        &   alnz,ja,nnz,neq,                                               ! sparse
        &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
positive.code: c
          if(ierr.ne.izero) return
  c
  clater      else if(lgdefp.eq.ione.and.ibbarp.eq.izero) then
  clater        call matinit_drv(
  clater     &   alnz,ja,nnz,neq,                                               ! sparse
  clater     &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
negative.code:       ntimdat(1)=nstep
        naxstp=izero
        nittot=izero
        ntimdat(6)=nittot
        nrftot=izero
        ntimdat(7)=nrftot
        ndtot=izero
        ntimdat(8)=ndtot
        ntimdat(9)=ireform
        call const(maxstp,delt,alfa,maxit,ntdinit,lgdef,utol,ftol,
       & etol,itmax,nintg,igroup,naxstp,nfirst,rtimdat,deltp,alfap,
       & ntimdat,nstep,maxitp,ntdinitp,lgdefp,itmaxp,gtol)
        if(skc) then
          call skclear(idslp,skew,numsn,numnp)
          call skcomp(x,d,skew,idslp,ipslp,ipstrs,numsn,numnp,nstep,
  