anchor.diff:     rot0_temp = 0.0D0
      rotu_temp = 0.0D0
      DO i=1,3
  -       rot0_temp(i) = Nuu0(i+3)
  +       rot0_temp(i) = Nuu0(i+6+3)
          rotu_temp(i) = uu0(i+3)
      ENDDO
      rot_temp = 0.0D0
  
positive.code:    rotu_temp = 0.0D0
     DO i=1,3
         E1(i) = E10(i) + uup(i)
         rotu_temp(i) = Nuuu(i+3)
     ENDDO
     rot_temp(:)  = 0.0D0
     rot0_temp(:) = 0.0D0
  
negative.code: 
     Nr = 0.0D0
     Nu_temp1 = 0.0D0
     Nu_temp1(1:3) = Nu(10:12)
     DO i=1,node_elem
         temp_id = (i - 1) * dof_node
         Nu_temp = 0.0D0
         DO k=1,3
  !           IF(i==1) Nu_temp1(k) = Nu(temp_id+k+3)
             Nu_temp(k) = Nu(temp_id+k+3)
         ENDDO
         Nr_temp = 0.0D0
  