Impose homogeneous Dirichlet condition on the top and right surfaces, instead of homogeneous Neumann. It's a bit more natural, amd it matches the small example in the getdp documentation.