Renamed the flux solvers in parameter output
Simplified computation of gravity rhs derivatives
Removed divergence correction, it is incorporated into the upwind solver now. Also simplified advection step