Fix rebase integration and parameter metadata

Author: Fei Yang

docs/parameters.yaml

@@ -22,7 +22,7 @@ parameters:
     type: Three Integers
     description: |
       Replicate the input STRU by Na, Nb, and Nc along its lattice vectors for distributed MDCell workflows. The default is 1 1 1, which preserves the input structure.
-    default_value: "1 1 1"
+    default_value: 1 1 1
     unit: ""
     availability: ""
   - name: calculation
@@ -1465,7 +1465,7 @@ parameters:
     category: Molecular dynamics
     type: Integer
     description: |
-      The output frequency of OUT.${suffix}/MD_dump in molecular dynamics calculations, which includes lattice and atomic information. Set to 0 to disable MD_dump output; scalar MD progress remains printed to the terminal and running_md.log every step.
+      The output frequency of OUT.${suffix}/MD_dump in molecular dynamics calculations, which includes lattice and atomic information. Set to 0 to disable MD_dump output.
     default_value: "1"
     unit: ""
     availability: ""

source/source_io/module_parameter/read_input_item_md.cpp

@@ -194,7 +194,7 @@ Note that md_tlast is only used in NVT/NPT simulations. If md_tlast is unset or
         item.description = "Whether to output the TOTAL-FORCE table in OUT.${suffix}/running_md.log for MDCell molecular dynamics. This does not affect force calculation or molecular dynamics integration.";
         item.default_value = "True";
         item.unit = "";
-        item.availability = "";
+        item.set_availability("");
         read_sync_bool(input.mdp.md_out_force);
         this->add_item(item);
     }

source/source_md/test/lj_pot_test.cpp

@@ -123,16 +123,16 @@ TEST_F(LJ_pot_test, mdcell_stress_includes_external_pressure)
 {
     ModuleESolver::ESolver_LJ p_esolver;
     MDCell mdcell = Setcell::setup_mdcell(ucell, param);
-    p_esolver.before_all_runners(mdcell, param.inp);
+    Input_para input = param.inp;
+    p_esolver.before_all_runners(mdcell, input);
     p_esolver.runner(mdcell, 0);
 
-    Input_para& global_input = const_cast<Input_para&>(PARAM.inp);
-    const double saved_press1 = global_input.press1;
-    const double saved_press2 = global_input.press2;
-    const double saved_press3 = global_input.press3;
-    global_input.press1 = 1.0;
-    global_input.press2 = 2.0;
-    global_input.press3 = 3.0;
+    const double saved_press1 = input.press1;
+    const double saved_press2 = input.press2;
+    const double saved_press3 = input.press3;
+    input.press1 = 1.0;
+    input.press2 = 2.0;
+    input.press3 = 3.0;
 
     p_esolver.cal_stress(mdcell, stress);
 
@@ -141,9 +141,9 @@ TEST_F(LJ_pot_test, mdcell_stress_includes_external_pressure)
     EXPECT_NEAR(stress(1, 1), p_esolver.lj_virial(1, 1) - 2.0 / unit_transform, doublethreshold);
     EXPECT_NEAR(stress(2, 2), p_esolver.lj_virial(2, 2) - 3.0 / unit_transform, doublethreshold);
 
-    global_input.press1 = saved_press1;
-    global_input.press2 = saved_press2;
-    global_input.press3 = saved_press3;
+    input.press1 = saved_press1;
+    input.press2 = saved_press2;
+    input.press3 = saved_press3;
 }
 
 TEST_F(LJ_pot_test, RcutSearchRadius)