anchor.diff:      type(hecmwST_matrix)     :: hecMAT
       type(hecmwST_local_mesh) :: hecMESH
       integer(kind=kint) :: i, in, j, nodLOCAL(20), ip, inod
  -    integer(kind=kint) :: itype, iS, iE, ic_type, icel, isect, IMAT, in0, nn, NDOF
  +    integer(kind=kint) :: itype, iS, iE, ic_type, icel, isect, IMAT, in0, nn
       real(kind=kreal) :: temp(20), lumped(120), mass(20*6, 20*6), ecoord(3,20)
       real(kind=kreal) :: delta_time, surf, THICK, ALFA, BETA
   
  -    NDOF = hecMESH%n_dof
       beta = fstrHEAT%beta
   
       do itype = 1, hecMESH%n_elem_type
  
positive.code:         nn = hecmw_get_max_node(ic_type)
          do i = 1, nn
            nodLOCAL(i) = hecMESH%elem_node_item(in0+i)
            do j = 1, NDOF
              ecoord(j,i) = hecMESH%node(3*(nodLOCAL(i)-1)+j)
            enddo
            temp(i) = fstrHEAT%TEMP0(nodLOCAL(i))
  
negative.code: 
        !$omp parallel default(none), &
          !$omp&  private(icel,isect,IMAT,nn,temp,in0,i,j,nodLOCAL,ecoord,in,thick,surf,inod,lumped,mass), &
          !$omp&  shared(iS,iE,hecMESH,ic_type,hecMAT,fstrHEAT,delta_time)
        !$omp do
        do icel = iS, iE
          isect = hecMESH%section_ID(icel)
  