An efficient multigrid approach for solving a discretized elliptic equation
whose boundary values are determined in part by integral relations is deve
loped, analyzed, and tested. The algorithm is motivated by a problem that i
s solved during integration of the 3D quasigeostrophic (QG) equations, whic
h model large-scale rotating stratified flows, where the integral constrain
ts represent mass conservation.