Revisions of sub-ice pressure gradient fixes

  Refactored and revised the recently added code in MOM_PressureForce_FV.F90 to
reduce the number of calls to the equation of state routines, and corrected a
number of minor bugs in the original implementation.  Answers are bitwise
identical unless the new options to reset the pressure gradient calculations are
actively selected.