anchor.diff:      end do
     end subroutine
     
  -  !> add Rhie-Chow flow of scalar s to r on polyFvGrid
  +  !> add Rhie-Chow advection of scalar s to r on polyFvGrid
     subroutine addRhieChowPolyScal(grid,s,p,gradP,rho,dt,r)
       use modPolyFvGrid
       class(polyFvGrid),intent(inout)::grid !< the grid
  
positive.code:   
  contains
    
    !> add Rhie-Chow flow of vector s to r on polyFvGrid
    subroutine addRhieChowPolyVect(grid,s,p,gradP,rho,dt,r)
      use modPolyFvGrid
      class(polyFvGrid),intent(inout)::grid !< the grid
  
negative.code:   ! constants
    integer,parameter::DIMS=3 !< dimensions
    
    !> generic adding Rhie-Chow flow/advection
    interface addRhieChow
      module procedure::addRhieChowMassFlowPoly
      module procedure::addRhieChowPolyVect
      module procedure::addRhieChowPolyScal
    end interface
  