debug: updating coeff primary

Author: RY4GIT

src/models/dCFE.py

@@ -66,6 +66,7 @@ def initialize(self):
         self.cfe_instance.refkdt = self.refkdt[:, 0]
         self.cfe_instance.satdk = self.satdk[:, 0]
         self.cfe_instance.reset_flux_and_states()
+        self.cfe_instance.reset_volume_tracking()
 
     def forward(self, x, t):  # -> (Tensor, Tensor):
         """

src/models/physics/bmi_cfe.py

@@ -411,7 +411,9 @@ def update_params(self, refkdt, satdk):
         self.Schaake_adjusted_magic_constant_by_soil_type = (
             self.refkdt * self.satdk / 2.0e-06
         )
-        self.soil_reservoir["coeff_primary"] = self.satdk
+        self.soil_reservoir["coeff_primary"] = (
+            self.satdk * self.soil_params["slop"] * self.time_step_size
+        )
         if self.verbose:
             print(
                 f"refkdt: {self.refkdt:.2f}; satdk: {self.satdk:.5f}; \