# Defines a differential equation describing changes in concentration with 
# respect to time due to advection.
  
macro advect(Δ, bc)
    return :( 
        function advect!(du,u,p,t)
            du = -p[1]*(Δ*bc*u)
        end
    )
  end
