anchor.diff:            IF( mtype == USERMATERIAL ) THEN
               CALL StressUpdate( gausses(LX), D3, dstrain, gausses(LX)%stress )
             ELSE IF( mtype == NORTON ) THEN
  -            gausses(LX)%pMaterial%mtype = mtype
  +            !gausses(LX)%pMaterial%mtype = mtype
               IF( tincr /= 0.0D0 .AND. ANY( gausses(LX)%stress /= 0.0D0 ) ) THEN
  -              gausses(LX)%pMaterial%mtype = mtype
  +              !gausses(LX)%pMaterial%mtype = mtype
                 IF( PRESENT(TT) .AND. PRESENT(T0) ) THEN
                   CALL StressUpdate( gausses(LX), D3, gausses(LX)%strain, gausses(LX)%stress, tincr, ttc, ttn )
                 ELSE
  
positive.code: 
        ELSE
  
          CALL MatlMatrix( gausses(LX), D3, D, tincr , coordsys)
  
        END IF
  
  
negative.code:     real(kind=kreal)   :: D(6,6), B(6,ndof*nn), B1(6,ndof*nn), spfunc(nn), ina(1)
      real(kind=kreal)   :: gderiv(nn,3), gdispderiv(3,3), det, WG, ttc,tt0, ttn,outa(1)
      integer(kind=kint) :: i, j, k, LX, mtype, serr
      integer(kind=kint) :: isEp
      real(kind=kreal)   :: naturalCoord(3), rot(3,3), mat(6,6), EPSTH(6)
      real(kind=kreal)   :: totaldisp(3,nn), elem(3,nn), elem1(3,nn), coordsys(3,3), tm(6,6)
      real(kind=kreal)   :: dstrain(6),dstress(6),dumstress(3,3),dum(3,3)
  