Fix advection multiple steps

[freifrauvonbleifrei]

fix the algorithm to do the right thing for multiple advection steps per run(); use std::memset to avoid overhead