Barotropic Algorithm - Pop-alpha
Simultaneously solve for: free surface height          and both velocities,
Invert using iterative CG routine
smoothing within each iteration is too costly!
momentum forcing terms
momentum forcing terms
What if we eliminate just this one smoothing step?
resolution: