module advection
  integer, parameter :: NADV = 2, NADVS=1 ! add'l boundary points
end module advection