Created A and b for the continuity equation, all functions now use Nglob instead of Nz