Is your model stable? How could you utilize the larger grid size and longer time step? Write down the coefficients of your algorithm? What would negative coefficient mean? What kind of initial and boundary conditions do you require to satisfy? How do you ensure that no contaminant will numerically leak out of upstream boundary?