anchor.diff:        real(kind=kreal) :: G,Gg,K,Kg,Kth, exp_n,mu_0,mu_n,dq_n,dtau, theta
         real(kind=kreal) :: ina(1), outa(2), EE, PP, ddt
   
  -      real(kind=kreal) :: devstrain(6), en(6)
  +      real(kind=kreal) :: devstrain(6), en(6), D(6,6)
         type(tTable), pointer  :: dicval
         logical :: ierr
   
  
positive.code:           PP = outa(2)
          endif
          if( matl%mtype>VISCOELASTIC ) ddt=trsinc(temp,tempn, matl%mtype, matl%variables)*dt
        else
          EE = matl%variables(M_YOUNGS)
          PP = matl%variables(M_POISSON)
  
negative.code: !>  \author                date              version
  !>  X.Yuan(Advancesoft)    2010/10/06        original
  !>  X.Yuan                 2012/09/18        add trs function(WLF, ARRHENIUS)
  !>  X.Yuan                 2016/09/26        abnormal processing in TRS
  !
  !======================================================================!
  module mViscoElastic
  