anchor.diff:       &     tfault,numfn,                                                ! fault
        &     state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,nstatesz,     ! elemnt
        &     ndmatsz,numelt,nconsz,                                       ! elemnt
  -     &     ptmp,infmat(1,imat),nprop,matgpt,elas_strs_14,               ! materl
  +     &     ptmp,nmatel,imatvar,nstate,nprop,matgpt,elas_strs_14,        ! materl
        &     td_strs_mat_14,matchg,                                       ! materl
        &     gauss,sh,shj,infetype,                                       ! eltype
        &     rtimdat,ntimdat,rgiter,                                      ! timdat
  
positive.code:      &     tfault,numfn,                                                ! fault
       &     state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,nstatesz,     ! elemnt
       &     ndmatsz,numelt,nconsz,                                       ! elemnt
       &     ptmp,infmat(1,imat),nprop,matgpt,elas_strs_13,               ! materl
       &     td_strs_mat_13,matchg,                                       ! materl
       &     gauss,sh,shj,infetype,                                       ! eltype
       &     rtimdat,ntimdat,rgiter,                                      ! timdat
  
negative.code:       do imat=1,numat
          matmodel=infmat(1,imat)
          nmatel=infmat(2,imat)
          indprop=infmat(3,imat)
          nstate=infmatmod(2,matmodel)
          nprop=infmatmod(3,matmodel)
          imatvar=infmatmod(4,matmodel)
          matchg=.false.
          call mathist(ptmp,prop(indprop),mhist(indprop),histry,nprop,
       &   imat,nstep,nhist,lastep,matchg,ierr,errstrng)
  