anchor.diff:      end if
       flow(:)=0d0
       !$omp parallel do default(shared)&
  -    !$omp& private(m,n,up,dn,dV,r,vf)
  +    !$omp& private(m,n,sameDirection,up,dn,dV,r,vf)
       do i=1,grid%nP
         m=grid%iEP(1,i)
         n=grid%iEP(2,i)
  
positive.code:       sensP(:,:)=0d0
      end if
      !$omp parallel do default(shared)&
      !$omp& private(m,n,up,dn,rhoUp,rhoDn,dRho,r,rhof,vMN,vMP,vPN,eps,flux)
      do i=1,grid%nP
        m=grid%iEP(1,i)
        n=grid%iEP(2,i)
  
negative.code:     end if
      flow(:,:)=0d0
      !$omp parallel do default(shared)&
      !$omp& private(m,n,sameDirection,up,dn,dV,r,vf)
      do i=1,grid%nP
        m=grid%iEP(1,i)
        n=grid%iEP(2,i)
  