+(*)Non-Boussinesq default for Z_INIT_REMAP_GENERAL

  Changed the default value of Z_INIT_REMAP_GENERAL to true for fully
non-Boussinesq configurations.  All existing fully non-Boussinesq cases that use
INIT_LAYERS_FROM_Z_FILE = True use this setting, and it is likely that such
cases will not work at all if this is false.  This change reduces the parameters
that need to be changed to go between equivalent Boussinesq and non-Boussinesq
configurations to just BOUSSINESQ and SEMI_BOUSSINESQ.  The previous default
(false) is being retained for any Boussinesq or semi-Bousssinesq cases so that
all answers and output are bitwise identical in those cases, but by default some
non-Boussinesq solutions change answers, and there are changes to the
MOM_parameter_doc files for those non-Boussinesq cases.

src/initialization/MOM_state_initialization.F90

@@ -2558,7 +2558,7 @@ subroutine MOM_temp_salt_initialize_from_Z(h, tv, depth_tot, G, GV, US, PF, just
   call get_param(PF, mdl, "Z_INIT_REMAP_GENERAL", remap_general, &
                  "If false, only initializes to z* coordinates. "//&
                  "If true, allows initialization directly to general coordinates.", &
-                 default=.false., do_not_log=just_read)
+                 default=.not.(GV%Boussinesq.or.GV%semi_Boussinesq) , do_not_log=just_read)
   call get_param(PF, mdl, "Z_INIT_REMAP_FULL_COLUMN", remap_full_column, &
                  "If false, only reconstructs profiles for valid data points. "//&
                  "If true, inserts vanished layers below the valid data.", &