boundary condition for vorticity (Jensen) left and right
for j = 1:ny+1
w[1,j] = (-4.0*s[2,j]+0.5*s[3,j])/(dx*dx)
w[nx+1,j]= (-4.0*s[nx,j]+0.5*s[nx-1,j])/(dx*dx)
end