Combine rebuild_global_mesh() and merge_meshes() functions

They were copy pasted from each other so might as well be one function for
consistency.
Add option to calculate predictor values even in non-temporally-adaptive Newton solve
Allow multiple meshes and clean up in copy_dof_pt_to_data_history_value(...)