periodic BC
```

```
w[nx+2,:] = w[2,:]
w[:,ny+2] = w[:,2]