diff --git a/scripts/2_disruption_analysis.py b/scripts/2_disruption_analysis.py
index 63e2134..a44e34a 100644
--- a/scripts/2_disruption_analysis.py
+++ b/scripts/2_disruption_analysis.py
@@ -22,6 +22,20 @@ def intersect_features_with_raster(
     features,
     flood_type,
 ):
+    """Intersects vector features with a raster dataset and computes flood depth for
+    each intersected feature.
+
+    Parameters:
+        raster_path (str): Path to the raster file containing flood data.
+        raster_key (str): Identifier for the raster dataset.
+        features (GeoDataFrame): Vector features (e.g., linestrings).
+        flood_type (str): Type of flood (e.g., "surface" or "river").
+
+    Returns:
+        GeoDataFrame: Intersected features with flood depth values,
+        reprojected to EPSG:27700.
+
+    """
     print(f"Intersecting features with raster {raster_key}...")
     # read the raster data: depth (meter)
     raster = io.read_raster_band_data(raster_path)
@@ -84,6 +98,26 @@ def compute_damage_level_on_flooded_roads(
     road_label,
     fldDepth,
 ):
+    """
+    Computes the damage level of roads based on flood type, road classification,
+    and flood depth.
+
+    Parameters:
+        fldType (str): Type of flood, either "surface" or "river".
+        road_classification (str): Motorway, A Road or B Road.
+        trunk_road (bool): Whether the road is a trunk road (True/False).
+        road_label (str): Label of the road, e.g., road, tunnel, bridge.
+        fldDepth (float): Depth of the flood in metres.
+
+    Returns:
+        str: Damage level as one of the following categories:
+            - "no": No damage
+            - "minor": Minor damage
+            - "moderate": Moderate damage
+            - "extensive": Extensive damage
+            - "severe": Severe damage
+    """
+
     depth = fldDepth * 100  # cm
     if fldType == "surface":
         if road_label == "tunnel" and (
@@ -325,11 +359,21 @@ def features_with_damage(
 
 
 def main(depth_thres):
-    """Inputs:
-    - base scenario: edge_flows_32p.gpq
-    - network: GB_road_links_with_bridges.gpq
-    - flood event map: Thames Lloyd's RDS (RASTER)
-    - clip mask: Thames Lloyd's RDS (VECTOR)
+    """Main function
+
+    Parameters
+    ----------
+    edge_flows_32p.gpq: base scenario output
+    GB_road_links_with_bridges.gpq: network element
+    Thames Lloyd's RDS (RASTER): JBA Flood Map
+    Thames Lloyd's RDS (Vector): Clip file
+
+    Returns
+    -------
+    intersections_x.pq: feature intersections with flood depth and damage level
+    road_links_x.gpq: with maxmimum flood depth and damage level by aggregating the
+        corresponding intersections.
+
     """
     # base scenario simulation results
     base_scenario_links = gpd.read_parquet(
diff --git a/scripts/3_damage_analysis.py b/scripts/3_damage_analysis.py
index 4496c04..bfc7852 100644
--- a/scripts/3_damage_analysis.py
+++ b/scripts/3_damage_analysis.py
@@ -138,10 +138,6 @@ def compute_damage_values(
         - max_cost (float): The maximum estimated damage cost.
         - mean_cost (float): The mean estimated damage cost.
 
-    Raises:
-    -------
-    AssertionError
-        If `road_label` is "bridge" but `bridge_width` is not provided.
     """
 
     def compute_bridge_damage(length, width, flood_type, damage_level):
@@ -357,11 +353,19 @@ def format_intersections(intersections, road_links):
 
 
 def main(depth_thres):
-    """Inputs:
-    - damage ratios
-    - damage values
-    - road links
-    - intersections in module 2
+    """Main function
+
+    Parameters
+    ----------
+    damage_ratio_road_flood.xlsx: damage curves
+    damage_cost_road_flood_uk.xlsx: asset damage values
+    GB_road_links_with_bridges.gpq: network links
+    intersections: module 2 output.
+
+    Returns
+    -------
+    intersections_with_damages (min, max, mean)
+
     """
     # damage curves
     """
diff --git a/scripts/4_rerouting_and_recovery.py b/scripts/4_rerouting_and_recovery.py
index 3fe7d2c..a900f88 100644
--- a/scripts/4_rerouting_and_recovery.py
+++ b/scripts/4_rerouting_and_recovery.py
@@ -92,13 +92,22 @@ def ordinary_road_recovery(
 
 
 def main(depth_thres, number_of_cpu):
-    """Inputs:
-    - model parameters
-    - road_links with added attributes:
+    """Main function:
+
+    Parameters
+    ----------
+    Model Parameters
+    odpfc_32p.pq: base scenario output
+    road_links_x.gpq: disruption analysis output
         [disruption analysis] road_label, flood_depth_max, damage_level_max,
         [base scenario analysis] current_capacity, current_speed,
         [config] free_flow_speed, min_flow_speeds, max_speed, initial_flow_speeds
-    - odpfc
+
+    Returns
+    -------
+    Daily edge flows during the recovery period (D-0 to D-110).
+    Isolated trips after daily recovery (D-0 to D-110).
+
     """
     # bridge recovery rates
     with open(base_path / "parameters" / "capt_minor.json", "r") as f: