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
resolution: