Couple using PISM's Neumann Boundary Condition

[citibeth]

Currently the coupling assumes a Dirichlet boundary condition for PISM (which is the default). Re-do this to use PISM's Neumann Boundary Condition, which should be preferable in a coupled case where fluxes are passed between the two models.