fixed jumps and normalization constant, compute visc only at each time step, not stage