start a density flux routine based on advection