anchor.diff:          integer(kind=kint), allocatable :: id_spc(:)
           real(kind=kreal) :: a
           integer(kind=kint) :: cdsys_ID
  +
  +        tt = 0.d0
  +        tt0 = 0.d0
   !C
   !C set temperature
   !C
  
positive.code: 
  ! element loop
  !$omp parallel default(none), &
  !$omp&  private(icel,iiS,j,nodLOCAL,i,ecoord,tt0,ttn,tt,ddu,du,total_disp, &
  !$omp&          cdsys_ID,coords,thick,qf,isect,ihead,tmp), &
  !$omp&  shared(iS,iE,hecMESH,nn,fstrSOLID,ndof,hecMAT,ic_type,fstrPR, &
  !$omp&         strainEnergy,iter,tincr)
  !$omp do
            do icel = iS, iE
  
  
negative.code:         ndof = hecMAT%NDOF
          fstrSOLID%QFORCE=0.0d0
  
          tt0 = 0.d0
          ttn = 0.d0
          tt = 0.d0
  
  ! --------------------------------------------------------------------
  !      updated
  !         1. stress and strain  : ep^(k) = ep^(k-1)+dep^(k)
  