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:                 nn = hecmw_get_max_node( ic_type )
  !C element loop
  !$omp parallel default(none), &
  !$omp&  private(icel,js,j,nodLOCAL,xx,yy,zz,tt,tt0,ecoord,edisp,isect,cdsys_ID,ihead,thick,mixflag,stiff,iflag,coords), &
  !$omp&  shared(iS,iE,hecMESH,temp,ref_temp,fstrSOLID,id_spc,force), &
  !$omp&  firstprivate(nn,ic_type)
  !$omp do
                  do icel = iS, iE
                          jS = hecMESH%elem_node_index(icel-1)
  
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)
  