anchor.diff:      allocate(f1c(DIMS,k))
       allocate(f2c(DIMS,DIMS,k))
       allocate(f3c(DIMS,k))
  +    !$omp workshare
       forall(i=1:k)
         u(:,i)=rhou(:,i)/rho(i)
         H(i)=(rhoE(i)+p(i))/rho(i)
  
positive.code:           dRhou(:,m)=dRhou(:,m)-flow(2:4)
            dRhoE(m)=dRhoE(m)-flow(5)
          end if
        end if
      end do
      deallocate(u,H)
      deallocate(f1c,f2c,f3c)
    end subroutine
  
negative.code:         &         dot_product(0.5d0*(f2c(:,3,m)+f2c(:,3,n)),grid%normP(:,i)),&
          &         dot_product(0.5d0*(f3c(:,m)+f3c(:,n)),grid%normP(:,i))]&
          &        -0.5d0*(dFEntropy(:)+dFAcoustic1(:)+dFAcoustic2(:)))
          !$omp critical
          if(n<=grid%nC)then
            dRho(m)=dRho(m)-flow(1)
            dRho(n)=dRho(n)+flow(1)
  