xslope 0.1.11__py3-none-any.whl → 0.1.13__py3-none-any.whl

xslope/fileio.py

@@ -29,8 +29,8 @@ def build_ground_surface(profile_lines):
     which represents the true ground surface.
 
     Parameters:
-        profile_lines (list of list of tuple): A list of profile lines, each represented
-            as a list of (x, y) coordinate tuples.
+        profile_lines (list of dict): A list of profile lines, each represented
+            as a dict with 'coords' key containing a list of (x, y) coordinate tuples.
 
     Returns:
         shapely.geometry.LineString: A LineString of the top surface, or an empty LineString
@@ -40,9 +40,12 @@ def build_ground_surface(profile_lines):
     if not profile_lines:
         return LineString([])
     
+    # Extract coordinate lists from profile line dicts
+    coord_lists = [line['coords'] for line in profile_lines]
+    
     # Step 1: Gather all points from all profile lines
     all_points = []
-    for line in profile_lines:
+    for line in coord_lists:
         all_points.extend(line)
     
     # Step 2: Group points by x-coordinate and find the highest y for each x
@@ -61,7 +64,7 @@ def build_ground_surface(profile_lines):
         
         # Check intersections with all profile lines
         is_topmost = True
-        for profile_line in profile_lines:
+        for profile_line in coord_lists:
             line = LineString(profile_line)
             if line.length == 0:
                 continue
@@ -125,10 +128,10 @@ def load_slope_data(filepath):
 
     try:
         template_version = main_df.iloc[4, 3]  # Excel row 5, column D
-        gamma_water = float(main_df.iloc[18, 3])  # Excel row 19, column D
-        tcrack_depth = float(main_df.iloc[19, 3])  # Excel row 20, column D
-        tcrack_water = float(main_df.iloc[20, 3])  # Excel row 21, column D
-        k_seismic = float(main_df.iloc[21, 3])  # Excel row 22, column D
+        gamma_water = float(main_df.iloc[7, 3])  # Excel row 8, column D
+        tcrack_depth = float(main_df.iloc[8, 3])  # Excel row 9, column D
+        tcrack_water = float(main_df.iloc[9, 3])  # Excel row 10, column D
+        k_seismic = float(main_df.iloc[10, 3])  # Excel row 11, column D
     except Exception as e:
         raise ValueError(f"Error reading static global values from 'main' tab: {e}")
 
@@ -139,37 +142,76 @@ def load_slope_data(filepath):
     max_depth = float(profile_df.iloc[1, 1])  # Excel B2 = row 1, column 1
 
     profile_lines = []
-
-    profile_data_blocks = [
-        {"header_row": 4, "data_start": 5, "data_end": 20},
-        {"header_row": 22, "data_start": 23, "data_end": 38}
-    ]
-    profile_block_width = 3
-
-    for block in profile_data_blocks:
-        for col in range(0, profile_df.shape[1], profile_block_width):
-            x_col, y_col = col, col + 1
+    
+    # New format: single data block, profile lines arranged horizontally
+    # First profile line: columns A:B, second: D:E, third: G:H, etc.
+    # Header row is row 4 (index 3), mat_id is in B5 (row 4, column 1)
+    # XY coordinates start in row 7 (index 6)
+    header_row = 3  # Excel row 4 (0-indexed)
+    mat_id_row = 4  # Excel row 5 (0-indexed)
+    coords_start_row = 7  # Excel row 8 (0-indexed)
+    
+    col = 0  # Start with column A (index 0)
+    while col < profile_df.shape[1]:
+        x_col = col
+        y_col = col + 1
+        
+        # Check if header row is empty (stop reading if empty)
+        try:
+            header_val = str(profile_df.iloc[header_row, x_col]).strip()
+            if not header_val or header_val.lower() == 'nan':
+                break  # No more profile lines
+        except:
+            break  # No more profile lines
+        
+        # Read mat_id from B5 (row 4, column 1) for this profile line
+        # Convert from 1-based to 0-based for internal use
+        try:
+            mat_id_val = profile_df.iloc[mat_id_row, y_col]
+            if pd.isna(mat_id_val):
+                mat_id = None
+            else:
+                # Convert to integer and subtract 1 to make it 0-based
+                mat_id = int(float(mat_id_val)) - 1
+                if mat_id < 0:
+                    mat_id = None  # Invalid mat_id
+        except (ValueError, TypeError):
+            mat_id = None
+        
+        # Read XY coordinates starting from row 7, stop at first empty row
+        coords = []
+        row = coords_start_row
+        while row < profile_df.shape[0]:
             try:
-                x_header = str(profile_df.iloc[block["header_row"], x_col]).strip().lower()
-                y_header = str(profile_df.iloc[block["header_row"], y_col]).strip().lower()
+                x_val = profile_df.iloc[row, x_col]
+                y_val = profile_df.iloc[row, y_col]
+                
+                # Stop at first empty row (both x and y are empty)
+                if pd.isna(x_val) and pd.isna(y_val):
+                    break
+                
+                # If at least one coordinate is present, try to convert
+                if pd.notna(x_val) and pd.notna(y_val):
+                    coords.append((float(x_val), float(y_val)))
             except:
-                continue
-            if x_header != 'x' or y_header != 'y':
-                continue
-            data = profile_df.iloc[block["data_start"]:block["data_end"], [x_col, y_col]]
-            data = data.dropna(how='all')
-            if data.empty:
-                continue
-            if data.iloc[0].isna().any():
-                continue
-            coords = data.dropna().apply(lambda r: (float(r.iloc[0]), float(r.iloc[1])), axis=1).tolist()
-            if len(coords) == 1:
-                raise ValueError("Each profile line must contain at least two points.")
-            if coords:
-                profile_lines.append(coords)
+                break
+            row += 1
+        
+        # Validate that we have at least 2 points
+        if len(coords) == 1:
+            raise ValueError(f"Each profile line must contain at least two points. Profile line starting at column {chr(65 + col)} has only one point.")
+        
+        if len(coords) >= 2:
+            # Store as dict with coords and mat_id
+            profile_lines.append({
+                'coords': coords,
+                'mat_id': mat_id
+            })
+        
+        # Move to next profile line (skip 3 columns: A->D, D->G, etc.)
+        col += 3
 
     # === BUILD GROUND SURFACE FROM PROFILE LINES ===
-
     ground_surface = build_ground_surface(profile_lines)
 
     # === BUILD TENSILE CRACK LINE ===
@@ -179,42 +221,65 @@ def load_slope_data(filepath):
         tcrack_surface = LineString([(x, y - tcrack_depth) for (x, y) in ground_surface.coords])
 
     # === MATERIALS (Optimized Parsing) ===
-    mat_df = xls.parse('mat', header=2)  # header=2 because the header row is row 3 in Excel
+    mat_df = xls.parse('mat', header=7)  # header=7 because the header row is row 8 in Excel (0-indexed row 7)
     materials = []
 
-    # Only process rows 4-15 (Excel), which are 0-indexed 0-11 in pandas 
-    for _, row in mat_df.iloc[0:12].iterrows():
-        # Check if the row is blank (columns 2-22, which are indices 1-21)
-        if row.iloc[1:22].isna().all():
-            continue
+    def _num(x):
+        v = pd.to_numeric(x, errors="coerce")
+        return float(v) if pd.notna(v) else 0.0
+
+    # Read materials row by row until we encounter an empty material name (Column B)
+    # Data starts at Excel row 9 (0-indexed row 0 after header=7)
+    for i in range(len(mat_df)):
+        row = mat_df.iloc[i]
+        
+        # Check if material name (Column B) is empty - stop reading if empty
+        material_name = row.get('name', '')
+        if pd.isna(material_name) or str(material_name).strip() == '':
+            break  # Stop reading when we encounter an empty material name
+        
+        # For seep workflows, 'g' (unit weight) and shear strength properties are not required.
+        # A material row is considered "missing" only if Excel columns C:X are empty.
+        # (Excel A:B are number and name; C:X contain the actual property fields.)
+        start_col = 2  # C
+        end_col = min(mat_df.shape[1], 24)  # X is column 24 (1-based) -> index 23, so slice end is 24
+        c_to_x_empty = True if start_col >= end_col else row.iloc[start_col:end_col].isna().all()
+        if c_to_x_empty:
+            # Excel row number: header is on row 8, first data row is row 9
+            excel_row = i + 9
+            raise ValueError(
+                "CRITICAL ERROR: Material row has empty property fields. "
+                f"Material '{material_name}' (Excel row {excel_row}) is blank in columns C:X."
+            )
+
         materials.append({
-            "name": row.get('name', ''),
-            "gamma": float(row.get('g', 0) or 0),
+            "name": str(material_name).strip(),
+            "gamma": _num(row.get("g", 0)),
             "option": str(row.get('option', '')).strip().lower(),
-            "c": float(row.get('c', 0) or 0),
-            "phi": float(row.get('f', 0) or 0),
-            "cp": float(row.get('cp', 0) or 0),
-            "r_elev": float(row.get('r-elev', 0) or 0),
-            "d": float(row.get('d', 0)) if pd.notna(row.get('d')) else 0,
-            "psi": float(row.get('ψ', 0)) if pd.notna(row.get('ψ')) else 0,
+            "c": _num(row.get('c', 0)),
+            "phi": _num(row.get('f', 0)),
+            "cp": _num(row.get('cp', 0)),
+            "r_elev": _num(row.get('r-elev', 0)),
+            "d": _num(row.get('d', 0)) if pd.notna(row.get('d')) else 0,
+            "psi": _num(row.get('ψ', 0)) if pd.notna(row.get('ψ')) else 0,
             "u": str(row.get('u', 'none')).strip().lower(),
-            "sigma_gamma": float(row.get('s(g)', 0) or 0),
-            "sigma_c": float(row.get('s(c)', 0) or 0),
-            "sigma_phi": float(row.get('s(f)', 0) or 0),
-            "sigma_cp": float(row.get('s(cp)', 0) or 0),
-            "sigma_d": float(row.get('s(d)', 0) or 0),
-            "sigma_psi": float(row.get('s(ψ)', 0) or 0),
-            "k1": float(row.get('k1', 0) or 0),
-            "k2": float(row.get('k2', 0) or 0),
-            "alpha": float(row.get('alpha', 0) or 0),
-            "kr0" : float(row.get('kr0', 0) or 0),
-            "h0" : float(row.get('h0', 0) or 0),
-            "E": float(row.get('E', 0) or 0),
-            "nu": float(row.get('n', 0) or 0)
+            "sigma_gamma": _num(row.get('s(g)', 0)),
+            "sigma_c": _num(row.get('s(c)', 0)),
+            "sigma_phi": _num(row.get('s(f)', 0)),
+            "sigma_cp": _num(row.get('s(cp)', 0)),
+            "sigma_d": _num(row.get('s(d)', 0)),
+            "sigma_psi": _num(row.get('s(ψ)', 0)),
+            "k1": _num(row.get('k1', 0)),
+            "k2": _num(row.get('k2', 0)),
+            "alpha": _num(row.get('alpha', 0)),
+            "kr0" : _num(row.get('kr0', 0)),
+            "h0" : _num(row.get('h0', 0)),
+            "E": _num(row.get('E', 0)),
+            "nu": _num(row.get('n', 0))
         })
 
     # === SEEPAGE ANALYSIS FILES ===
-    # Check if any materials use seepage analysis for pore pressure
+    # Check if any materials use seep analysis for pore pressure
     has_seep_materials = any(material["u"] == "seep" for material in materials)
     
     seep_mesh = None
@@ -222,111 +287,213 @@ def load_slope_data(filepath):
     seep_u2 = None
     
     if has_seep_materials:
-        # Read seepage file names directly from Excel cells L19, L20, L21
         try:
-            # Read the 'mat' sheet directly without header parsing
-            mat_raw_df = xls.parse('mat', header=None)
-            
-            # L19 = row 18, column 11 (0-indexed)
-            mesh_filename = str(mat_raw_df.iloc[18, 11]).strip()  # L19
-            solution1_filename = str(mat_raw_df.iloc[19, 11]).strip()  # L20  
-            solution2_filename = str(mat_raw_df.iloc[20, 11]).strip()  # L21
-            
-            # Validate required files
-            if not mesh_filename or mesh_filename.lower() == 'nan':
-                raise ValueError("CRITICAL ERROR: Mesh filename is required when using 'seep' pore pressure option but is blank in cell L19.")
-            if not solution1_filename or solution1_filename.lower() == 'nan':
-                raise ValueError("CRITICAL ERROR: Solution1 filename is required when using 'seep' pore pressure option but is blank in cell L20.")
-            
-            # Load mesh file
+            base, _ = os.path.splitext(filepath)
+            mesh_filename = f"{base}_mesh.json"
+            solution1_filename = f"{base}_seep.csv"
+            solution2_filename = f"{base}_seep2.csv"
+
+            missing_required = []
             if not os.path.exists(mesh_filename):
-                raise ValueError(f"CRITICAL ERROR: Mesh file '{mesh_filename}' not found.")
-            seep_mesh = import_mesh_from_json(mesh_filename)
-            
-            # Load solution1 file
+                missing_required.append(mesh_filename)
             if not os.path.exists(solution1_filename):
-                raise ValueError(f"CRITICAL ERROR: Solution1 file '{solution1_filename}' not found.")
-            solution1_df = pd.read_csv(solution1_filename)
-            # Skip the last row which contains the total flowrate comment
-            solution1_df = solution1_df.iloc[:-1]
-            seep_u = solution1_df["u"].to_numpy()
-            
-            # Load solution2 file if provided
-            if solution2_filename and solution2_filename.lower() != 'nan':
-                if not os.path.exists(solution2_filename):
-                    raise ValueError(f"CRITICAL ERROR: Solution2 file '{solution2_filename}' not found.")
-                solution2_df = pd.read_csv(solution2_filename)
-                # Skip the last row which contains the total flowrate comment
-                solution2_df = solution2_df.iloc[:-1]
-                seep_u2 = solution2_df["u"].to_numpy()
-                
-        except Exception as e:
-            if "CRITICAL ERROR" in str(e):
-                raise e
+                missing_required.append(solution1_filename)
+
+            if missing_required:
+                missing_list = ", ".join(f"'{path}'" for path in missing_required)
+                print(
+                    "WARNING: Seep pore pressure option selected but required seep files "
+                    f"were not found: {missing_list}. Continuing without seep data."
+                )
             else:
-                raise ValueError(f"Error reading seepage files: {e}")
+                seep_mesh = import_mesh_from_json(mesh_filename)
+                solution1_df = pd.read_csv(solution1_filename)
+                solution1_df = solution1_df.iloc[:-1]
+                seep_u = solution1_df["u"].to_numpy()
+
+                if os.path.exists(solution2_filename):
+                    solution2_df = pd.read_csv(solution2_filename)
+                    solution2_df = solution2_df.iloc[:-1]
+                    seep_u2 = solution2_df["u"].to_numpy()
+
+        except Exception as e:
+            print(f"WARNING: Error reading seepage files: {e}. Continuing without seep data.")
 
     # === PIEZOMETRIC LINE ===
-    piezo_df = xls.parse('piezo')
+    piezo_df = xls.parse('piezo', header=None)
     piezo_line = []
     piezo_line2 = []
 
-    # Read all data once (rows 4-18)
-    piezo_data = piezo_df.iloc[2:18].dropna(how='all')
+    # Read first piezometric line (columns A:B, starting at row 4, Excel row 4 = index 3)
+    # Keep reading until we encounter an empty row
+    start_row = 3  # Excel row 4 (0-indexed row 3)
+    x_col = 0  # Column A
+    y_col = 1  # Column B
     
-    if len(piezo_data) >= 2:
-        # Extract first table (A4:B18) - columns 0 and 1
+    row = start_row
+    while row < piezo_df.shape[0]:
         try:
-            piezo_data1 = piezo_data.dropna(subset=[piezo_data.columns[0], piezo_data.columns[1]], how='all')
-            if len(piezo_data1) < 2:
-                raise ValueError("First piezometric line must contain at least two points.")
-            piezo_line = piezo_data1.apply(lambda row: (float(row.iloc[0]), float(row.iloc[1])), axis=1).tolist()
-        except Exception:
-            raise ValueError("Invalid first piezometric line format.")
-
-        # Extract second table (D4:E18) - columns 3 and 4
+            x_val = piezo_df.iloc[row, x_col]
+            y_val = piezo_df.iloc[row, y_col]
+            
+            # Stop at first empty row (both x and y are empty)
+            if pd.isna(x_val) and pd.isna(y_val):
+                break
+            
+            # If at least one coordinate is present, try to convert
+            if pd.notna(x_val) and pd.notna(y_val):
+                piezo_line.append((float(x_val), float(y_val)))
+        except:
+            break
+        row += 1
+    
+    # Validate first piezometric line
+    if len(piezo_line) == 1:
+        raise ValueError("First piezometric line must contain at least two points.")
+    
+    # Read second piezometric line (columns D:E, starting at row 4, Excel row 4 = index 3)
+    # Keep reading until we encounter an empty row
+    x_col2 = 3  # Column D
+    y_col2 = 4  # Column E
+    
+    row = start_row
+    while row < piezo_df.shape[0]:
         try:
-            piezo_data2 = piezo_data.dropna(subset=[piezo_data.columns[3], piezo_data.columns[4]], how='all')
-            if len(piezo_data2) < 2:
-                raise ValueError("Second piezometric line must contain at least two points.")
-            piezo_line2 = piezo_data2.apply(lambda row: (float(row.iloc[3]), float(row.iloc[4])), axis=1).tolist()
-        except Exception:
-            # If second table reading fails, just leave piezo_line2 as empty list
-            piezo_line2 = []
-    elif len(piezo_data) == 1:
-        raise ValueError("Piezometric line must contain at least two points.")
+            x_val = piezo_df.iloc[row, x_col2]
+            y_val = piezo_df.iloc[row, y_col2]
+            
+            # Stop at first empty row (both x and y are empty)
+            if pd.isna(x_val) and pd.isna(y_val):
+                break
+            
+            # If at least one coordinate is present, try to convert
+            if pd.notna(x_val) and pd.notna(y_val):
+                piezo_line2.append((float(x_val), float(y_val)))
+        except:
+            break
+        row += 1
+    
+    # Validate second piezometric line (only if it has data)
+    if len(piezo_line2) == 1:
+        raise ValueError("Second piezometric line must contain at least two points if provided.")
 
     # === DISTRIBUTED LOADS ===
+    # Read first set from "dloads" tab
     dload_df = xls.parse('dloads', header=None)
     dloads = []
+    
+    # Start reading from column B (index 1), each distributed load uses 3 columns (X, Y, Normal)
+    # Keep reading to the right until we encounter an empty distributed load
+    start_row = 3  # Excel row 4 (0-indexed row 3)
+    col = 1  # Start with column B (index 1)
+    
+    while col < dload_df.shape[1]:
+        x_col = col
+        y_col = col + 1
+        normal_col = col + 2
+        
+        # Check if dataframe has enough rows before accessing start_row
+        if dload_df.shape[0] <= start_row:
+            break  # Not enough rows, stop reading
+
+        # Check if this distributed load block is empty (check first row for X coordinate)
+        if pd.isna(dload_df.iloc[start_row, x_col]):
+            break  # Stop reading when we encounter an empty distributed load
+        
+        # Read points for this distributed load, keep reading down until empty row
+        block_points = []
+        row = start_row
+        while row < dload_df.shape[0]:
+            try:
+                x_val = dload_df.iloc[row, x_col]
+                y_val = dload_df.iloc[row, y_col]
+                normal_val = dload_df.iloc[row, normal_col]
+                
+                # Stop at first empty row (all three values are empty)
+                if pd.isna(x_val) and pd.isna(y_val) and pd.isna(normal_val):
+                    break
+                
+                # If at least coordinates are present, try to convert
+                if pd.notna(x_val) and pd.notna(y_val):
+                    normal = float(normal_val) if pd.notna(normal_val) else 0.0
+                    block_points.append({
+                        "X": float(x_val),
+                        "Y": float(y_val),
+                        "Normal": normal
+                    })
+            except:
+                break
+            row += 1
+        
+        # Validate that we have at least 2 points
+        if len(block_points) == 1:
+            raise ValueError(f"Each distributed load must contain at least two points. Distributed load starting at column {chr(65 + col)} has only one point.")
+        
+        if len(block_points) >= 2:
+            dloads.append(block_points)
+        
+        # Move to next distributed load (skip 4 columns: 3 for the dload + 1 empty column)
+        col += 4
+    
+    # Read second set from "dloads (2)" tab
     dloads2 = []
-    dload_data_blocks = [
-        {"start_row": 3, "end_row": 13},
-        {"start_row": 16, "end_row": 26}
-    ]
-    dload_block_starts = [1, 5, 9, 13]
-
-    for block_idx, block in enumerate(dload_data_blocks):
-        for col in dload_block_starts:
-            section = dload_df.iloc[block["start_row"]:block["end_row"], col:col + 3]
-            section = section.dropna(how='all')
-            section = section.dropna(subset=[col, col + 1], how='any')
-            if len(section) >= 2:
+    try:
+        dload_df2 = xls.parse('dloads (2)', header=None)
+        
+        # Start reading from column B (index 1), each distributed load uses 3 columns (X, Y, Normal)
+        # Keep reading to the right until we encounter an empty distributed load
+        col = 1  # Start with column B (index 1)
+        
+        while col < dload_df2.shape[1]:
+            x_col = col
+            y_col = col + 1
+            normal_col = col + 2
+            
+            # Check if dataframe has enough rows before accessing start_row
+            if dload_df2.shape[0] <= start_row:
+                break  # Not enough rows, stop reading
+            
+            # Check if this distributed load block is empty (check first row for X coordinate)
+            if pd.isna(dload_df2.iloc[start_row, x_col]):
+                break  # Stop reading when we encounter an empty distributed load
+            
+            # Read points for this distributed load, keep reading down until empty row
+            block_points = []
+            row = start_row
+            while row < dload_df2.shape[0]:
                 try:
-                    block_points = section.apply(
-                        lambda row: {
-                            "X": float(row.iloc[0]),
-                            "Y": float(row.iloc[1]),
-                            "Normal": float(row.iloc[2])
-                        }, axis=1).tolist()
-                    if block_idx == 0:
-                        dloads.append(block_points)
-                    else:
-                        dloads2.append(block_points)
+                    x_val = dload_df2.iloc[row, x_col]
+                    y_val = dload_df2.iloc[row, y_col]
+                    normal_val = dload_df2.iloc[row, normal_col]
+                    
+                    # Stop at first empty row (all three values are empty)
+                    if pd.isna(x_val) and pd.isna(y_val) and pd.isna(normal_val):
+                        break
+                    
+                    # If at least coordinates are present, try to convert
+                    if pd.notna(x_val) and pd.notna(y_val):
+                        normal = float(normal_val) if pd.notna(normal_val) else 0.0
+                        block_points.append({
+                            "X": float(x_val),
+                            "Y": float(y_val),
+                            "Normal": normal
+                        })
                 except:
-                    raise ValueError("Invalid data format in distributed load block.")
-            elif len(section) == 1:
-                raise ValueError("Each distributed load block must contain at least two points.")
+                    break
+                row += 1
+            
+            # Validate that we have at least 2 points
+            if len(block_points) == 1:
+                raise ValueError(f"Each distributed load must contain at least two points. Distributed load starting at column {chr(65 + col)} has only one point.")
+            
+            if len(block_points) >= 2:
+                dloads2.append(block_points)
+            
+            # Move to next distributed load (skip 4 columns: 3 for the dload + 1 empty column)
+            col += 4
+    except (ValueError, KeyError):
+        # If "dloads (2)" tab doesn't exist, just leave dloads2 as empty list
+        pass
 
     # === CIRCLES ===
 
@@ -376,11 +543,16 @@ def load_slope_data(filepath):
     reinforce_df = xls.parse('reinforce', header=1)  # Header in row 2 (0-indexed row 1)
     reinforce_lines = []
     
-    # Process rows 3-22 (Excel) which are 0-indexed rows 0-19 in pandas after header=1
-    for i, row in reinforce_df.iloc[0:20].iterrows():
-        # Check if the row has coordinate data (x1, y1, x2, y2)
-        if pd.isna(row.iloc[1]) or pd.isna(row.iloc[2]) or pd.isna(row.iloc[3]) or pd.isna(row.iloc[4]):
-            continue  # Skip empty rows
+    # Process rows starting from row 3 (Excel) which is 0-indexed row 0 in pandas after header=1
+    # Keep reading until we encounter an empty value in column B
+    for i, row in reinforce_df.iterrows():
+        # Check if column B (x1 coordinate) is empty - stop reading if empty
+        if pd.isna(row.iloc[1]):
+            break  # Stop reading when column B is empty
+        
+        # Check if other required coordinates are present
+        if pd.isna(row.iloc[2]) or pd.isna(row.iloc[3]) or pd.isna(row.iloc[4]):
+            continue  # Skip rows with incomplete coordinate data
             
         # If coordinates are present, check for required parameters (Tmax, Lp1, Lp2)
         if pd.isna(row.iloc[5]) or pd.isna(row.iloc[7]) or pd.isna(row.iloc[8]):
@@ -508,76 +680,160 @@ def load_slope_data(filepath):
 
 
     # === SEEPAGE ANALYSIS BOUNDARY CONDITIONS ===
+    # Read first set from "seep bc" sheet
     seep_df = xls.parse('seep bc', header=None)
     seepage_bc = {"specified_heads": [], "exit_face": []}
-
-    # Specified Head #1
-    head1 = seep_df.iloc[2, 2] if seep_df.shape[1] > 2 and seep_df.shape[0] > 2 else None
-    coords1 = []
-    for i in range(4, 12):  # rows 5-12 (0-indexed 4-11)
-        if i >= seep_df.shape[0]:
-            break
-        x = seep_df.iloc[i, 1] if seep_df.shape[1] > 1 else None
-        y = seep_df.iloc[i, 2] if seep_df.shape[1] > 2 else None
-        if pd.notna(x) and pd.notna(y):
-            coords1.append((float(x), float(y)))
-    if head1 is not None and coords1:
-        seepage_bc["specified_heads"].append({"head": float(head1), "coords": coords1})
-
-    # Specified Head #2
-    head2 = seep_df.iloc[2, 5] if seep_df.shape[1] > 5 and seep_df.shape[0] > 2 else None
-    coords2 = []
-    for i in range(4, 12):
-        if i >= seep_df.shape[0]:
-            break
-        x = seep_df.iloc[i, 4] if seep_df.shape[1] > 4 else None
-        y = seep_df.iloc[i, 5] if seep_df.shape[1] > 5 else None
-        if pd.notna(x) and pd.notna(y):
-            coords2.append((float(x), float(y)))
-    if head2 is not None and coords2:
-        seepage_bc["specified_heads"].append({"head": float(head2), "coords": coords2})
-
-    # Specified Head #3
-    head3 = seep_df.iloc[2, 8] if seep_df.shape[1] > 8 and seep_df.shape[0] > 2 else None
-    coords3 = []
-    for i in range(4, 12):
-        if i >= seep_df.shape[0]:
-            break
-        x = seep_df.iloc[i, 7] if seep_df.shape[1] > 7 else None
-        y = seep_df.iloc[i, 8] if seep_df.shape[1] > 8 else None
-        if pd.notna(x) and pd.notna(y):
-            coords3.append((float(x), float(y)))
-    if head3 is not None and coords3:
-        seepage_bc["specified_heads"].append({"head": float(head3), "coords": coords3})
-
-    # Exit Face
+    
+    # Exit Face BC: starts at B5 (row 4, columns 1 and 2), continues down until empty x value
     exit_coords = []
-    for i in range(15, 23):  # rows 16-23 (0-indexed 15-22)
-        if i >= seep_df.shape[0]:
+    exit_start_row = 4  # Excel row 5 (0-indexed row 4)
+    exit_x_col = 1  # Column B
+    exit_y_col = 2  # Column C
+    
+    row = exit_start_row
+    while row < seep_df.shape[0]:
+        try:
+            x_val = seep_df.iloc[row, exit_x_col]
+            y_val = seep_df.iloc[row, exit_y_col]
+            
+            # Stop at first empty x value
+            if pd.isna(x_val):
+                break
+            
+            # If x is present, try to convert (y can be empty but we'll still add the point)
+            if pd.notna(x_val) and pd.notna(y_val):
+                exit_coords.append((float(x_val), float(y_val)))
+        except:
             break
-        x = seep_df.iloc[i, 1] if seep_df.shape[1] > 1 else None
-        y = seep_df.iloc[i, 2] if seep_df.shape[1] > 2 else None
-        if pd.notna(x) and pd.notna(y):
-            exit_coords.append((float(x), float(y)))
+        row += 1
     seepage_bc["exit_face"] = exit_coords
+    
+    # Specified Head BCs: start at columns E:F, then H:I, etc.
+    # Head value is in row 3 (index 2), XY values start at row 5 (index 4)
+    # Keep reading to the right until head value in row 3 is empty
+    head_row = 2  # Excel row 3 (0-indexed row 2)
+    data_start_row = 4  # Excel row 5 (0-indexed row 4)
+    col = 4  # Start with column E (index 4)
+    
+    while col < seep_df.shape[1]:
+        x_col = col
+        y_col = col + 1
+        head_col = col + 1  # Head value is in the Y column (F, I, L, etc.)
+        
+        # Check if head value in row 3 is empty - stop reading if empty
+        if seep_df.shape[0] <= head_row:
+            break
+        head_val = seep_df.iloc[head_row, head_col]
+        if pd.isna(head_val):
+            break  # Stop reading when head value is empty
+        
+        # Read XY coordinates starting from row 5, continue down until empty
+        coords = []
+        row = data_start_row
+        while row < seep_df.shape[0]:
+            try:
+                x_val = seep_df.iloc[row, x_col]
+                y_val = seep_df.iloc[row, y_col]
+                
+                # Stop at first empty x value
+                if pd.isna(x_val):
+                    break
+                
+                # If x is present, try to convert
+                if pd.notna(x_val) and pd.notna(y_val):
+                    coords.append((float(x_val), float(y_val)))
+            except:
+                break
+            row += 1
+        
+        if coords:  # Only add if we have coordinates
+            seepage_bc["specified_heads"].append({"head": float(head_val), "coords": coords})
+        
+        # Move to next specified head BC (skip 3 columns: E->H, H->K, etc.)
+        col += 3
+    
+    # Read second set from "seep bc (2)" sheet
+    seepage_bc2 = {"specified_heads": [], "exit_face": []}
+    try:
+        seep_df2 = xls.parse('seep bc (2)', header=None)
+        
+        # Exit Face BC: starts at B5 (row 4, columns 1 and 2), continues down until empty x value
+        exit_coords2 = []
+        row = exit_start_row
+        while row < seep_df2.shape[0]:
+            try:
+                x_val = seep_df2.iloc[row, exit_x_col]
+                y_val = seep_df2.iloc[row, exit_y_col]
+                
+                # Stop at first empty x value
+                if pd.isna(x_val):
+                    break
+                
+                # If x is present, try to convert
+                if pd.notna(x_val) and pd.notna(y_val):
+                    exit_coords2.append((float(x_val), float(y_val)))
+            except:
+                break
+            row += 1
+        seepage_bc2["exit_face"] = exit_coords2
+        
+        # Specified Head BCs: same structure as first sheet
+        col = 4  # Start with column E (index 4)
+        while col < seep_df2.shape[1]:
+            x_col = col
+            y_col = col + 1
+            head_col = col + 1  # Head value is in the Y column
+            
+            # Check if head value in row 3 is empty - stop reading if empty
+            if seep_df2.shape[0] <= head_row:
+                break
+            head_val = seep_df2.iloc[head_row, head_col]
+            if pd.isna(head_val):
+                break  # Stop reading when head value is empty
+            
+            # Read XY coordinates starting from row 5, continue down until empty
+            coords = []
+            row = data_start_row
+            while row < seep_df2.shape[0]:
+                try:
+                    x_val = seep_df2.iloc[row, x_col]
+                    y_val = seep_df2.iloc[row, y_col]
+                    
+                    # Stop at first empty x value
+                    if pd.isna(x_val):
+                        break
+                    
+                    # If x is present, try to convert
+                    if pd.notna(x_val) and pd.notna(y_val):
+                        coords.append((float(x_val), float(y_val)))
+                except:
+                    break
+                row += 1
+            
+            if coords:  # Only add if we have coordinates
+                seepage_bc2["specified_heads"].append({"head": float(head_val), "coords": coords})
+            
+            # Move to next specified head BC (skip 3 columns: E->H, H->K, etc.)
+            col += 3
+    except (ValueError, KeyError):
+        # If "seep bc (2)" sheet doesn't exist, just leave seepage_bc2 as empty
+        pass
 
     # === VALIDATION ===
  
     circular = len(circles) > 0
-    # Check if this is a seepage-only analysis (has seepage BCs but no slope stability surfaces)
+    # Check if this is a seep-only analysis (has seep BCs but no slope stability surfaces)
     has_seepage_bc = (len(seepage_bc.get("specified_heads", [])) > 0 or 
                      len(seepage_bc.get("exit_face", [])) > 0)
     is_seepage_only = has_seepage_bc and not circular and len(non_circ) == 0
     
-    # Only require circular/non-circular data if this is NOT a seepage-only analysis
+    # Only require circular/non-circular data if this is NOT a seep-only analysis
     if not is_seepage_only and not circular and len(non_circ) == 0:
         raise ValueError("Input must include either circular or non-circular surface data.")
     if not profile_lines:
         raise ValueError("Profile lines sheet is empty or invalid.")
     if not materials:
         raise ValueError("Materials sheet is empty.")
-    if len(materials) != len(profile_lines):
-        raise ValueError("Each profile line must have a corresponding material. You have " + str(len(materials)) + " materials and " + str(len(profile_lines)) + " profile lines.")
         
 
     # Add everything to globals_data
@@ -600,8 +856,9 @@ def load_slope_data(filepath):
     globals_data["dloads2"] = dloads2
     globals_data["reinforce_lines"] = reinforce_lines
     globals_data["seepage_bc"] = seepage_bc
+    globals_data["seepage_bc2"] = seepage_bc2
     
-    # Add seepage data if available
+    # Add seep data if available
     if has_seep_materials:
         globals_data["seep_mesh"] = seep_mesh
         globals_data["seep_u"] = seep_u