anchor.diff:  c
         include "rtimdat_dim.inc"
         include "ntimdat_dim.inc"
  +      include "rgiter_dim.inc"
   c
   c...  external routines
   c
  
positive.code: c
        include "rtimdat_def.inc"
        include "ntimdat_def.inc"
  c
  cdebug      write(6,*) "Hello from td_matinit_cmp_ss_f!"
  c
  
negative.code:           ngauss=infetype(1,ietype)
            inddmatg=inddmat
            indstateg=indstate
  cdebug          write(6,*) "ind,matgpt,iel,ietype,indstate,inddmat,ngauss:",
  cdebug     &     ind,matgpt,iel,ietype,indstate,inddmat,ngauss
            do l=1,ngauss
  cdebug            write(6,*) "l,inddmatg,indstateg:",l,inddmatg,indstateg
              call td_matinit(state(1,indstateg),dstate(1,indstateg),
       &       dmat(1,inddmatg),prop,rtimdat,rgiter,ntimdat,iddmat,tmax,
       &       nstate,nprop,matchg,ierr,errstrng)
              if(ierr.ne.izero) return
              tminmax=min(tminmax,tmax)
              inddmatg=inddmatg+ione
  