anchor.diff:          if(ierr.ne.izero) return
   c
         else if(lgdefp.eq.izero.and.intord.eq.ithree) then
  +        write(kto,650)
           call matinit_drv(
        &   alnz,ja,nnz,neq,                                               ! sparse
        &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
positive.code: clater        if(ierr.ne.izero) return
  c
  clater      else if(lgdefp.eq.1.and.intord.eq.ithree) then
  clater        call matinit_drv(
  clater     &   alnz,ja,nnz,neq,                                               ! sparse
  clater     &   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
  