move main function from run_simulation.py to solifluction.py

add permafrost flag and model for heat_transfer

Author: saeb-faraji-gargari

source/package/soli3d/__init__.py

@@ -7,7 +7,6 @@
     interpolation,
     io_data_process,
     layer,
-    mass_conservation,
     momentum,
     phase_detect,
     solifluction,
@@ -29,7 +28,6 @@
     "interpolation",
     "io_data_process",
     "layer",
-    "mass_conservation",
     "momentum",
     "phase_detect",
     "solifluction",

source/package/soli3d/io_data_process.py

@@ -200,14 +200,29 @@ def load_daily_temperatures(csv_path: str) -> tuple[list[float], list[float]]:
     return days, temps
 
 
+def clean_bool(val: str, name: str = "unknown") -> bool:
+    value = val.split("#")[0].strip().lower()
+    true_values = ("true", "1", "yes", "y")
+    false_values = ("false", "0", "no", "n")
+
+    if value in true_values:
+        return True
+    if value in false_values:
+        return False
+
+    raise ValueError(
+        f"Invalid boolean value for '{name}': '{val}'. "
+        f"Expected one of {true_values + false_values}."
+    )
+
+
 def read_config_file(param_path: Path) -> tuple[
     int,  # number_of_iterations
     float,  # dt_momentum
     int,  # momentum_iteration_threshold
     float,  # dt_global_model
     float,  # dt_heat_transfer
     float,  # dt_mass_conservation
-    float,  # time_end_simulation
     float,  # write_intervals_time
     int,  # partition_shape_size
     float,  # h_mesh_step_value
@@ -247,11 +262,13 @@ def clean_int(val: str) -> int:
         input_variables["simulation"]["momentum_iteration_threshold"]
     )
 
-    time_end_simulation = clean_float(
-        input_variables["simulation"]["time_end_simulation"]
-    )
+    permafrost = clean_bool(input_variables["simulation"]["permafrost"], "permafrost")
+
+    # time_end_simulation = clean_float(
+    #     input_variables["simulation"]["time_end_simulation"]
+    # )
 
-    write_intervals_time = clean_float(
+    write_intervals_time: float = clean_float(
         input_variables["simulation"]["write_intervals_time"]
     )
 
@@ -301,7 +318,6 @@ def clean_int(val: str) -> int:
         dt_global_model,
         dt_heat_transfer,
         dt_mass_conservation,
-        time_end_simulation,
         write_intervals_time,
         partition_shape_size,
         h_mesh_step_value,
@@ -316,6 +332,7 @@ def clean_int(val: str) -> int:
         temps_temperature_file,
         results_pathname,
         slope_radian,
+        permafrost,
     )