Hi Dabing Chen, thanks for the upate, but please note that ?advection.1D says that: "When solved dynamically, the euler method should be used, unless the first-order upstream method is used." Thomas P.