Out of bound in simulation.cpp:velocityStep
The h+1
is past the end
real vp2 = std::pow(h < ni-1 // \todo recode loops for speed
? (svphi(i,h,j)+svphi(i,h+1,j)+svphi(i-1,h,j)+svphi(i-1,h+1,j))/4
: (svphi(i,h,j)+svphi(i-1,h,j)+svphi(i-1,h+1,j))/3,
2);