anchor.diff:          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)
  
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:         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
  