ghost points
s[1,:] = s[nx+1,:]
s[:,1] = s[:,ny+1]