+(*)Set the default of USE_PSURF_IN_EOS to True

  Changed the default value of USE_PSURF_IN_EOS to True.  This could change
answers in some test cases with sea ice where USE_PSURF_IN_EOS is not set
explicitly, and there are changes to the MOM_parameter_doc files.

src/core/MOM.F90

@@ -1908,7 +1908,7 @@ subroutine initialize_MOM(Time, Time_init, param_file, dirs, CS, restart_CSp, &
                  default=3991.86795711963, scale=US%J_kg_to_Q)
     call get_param(param_file, "MOM", "USE_PSURF_IN_EOS", CS%use_p_surf_in_EOS, &
                  "If true, always include the surface pressure contributions "//&
-                 "in equation of state calculations.", default=.false.) !### Change the default.
+                 "in equation of state calculations.", default=.true.)
   endif
   if (use_EOS) call get_param(param_file, "MOM", "P_REF", CS%tv%P_Ref, &
                  "The pressure that is used for calculating the coordinate "//&