fixing (?) courant array ranges for velocity interpolation in solver_pressure.hpp