Code working, although +x boundary shows some oscillations. Will try periodic boundaries and Courand criteria