Boundaries are updated by power law extapolations. Artificial viscosity is now set by the min(abs(v), cs) instead of just min(v, cs).