Add kwarg to sort sections

Author: 1-Bart-1

src/body_aerodynamics.jl

@@ -114,7 +114,7 @@ function Base.setproperty!(obj::BodyAerodynamics, sym::Symbol, val)
 end
 
 """
-    reinit!(body_aero::BodyAerodynamics; init_aero, va, omega, refine_mesh, recompute_mapping)
+    reinit!(body_aero::BodyAerodynamics; init_aero, va, omega, refine_mesh, recompute_mapping, sort_sections)
 
 Initialize a BodyAerodynamics struct in-place by setting up panels and coefficients.
 
@@ -129,6 +129,8 @@ Initialize a BodyAerodynamics struct in-place by setting up panels and coefficie
   `deform!()` to preserve deformed geometry.
 - `recompute_mapping=true`: Whether to recompute the refined panel mapping.
   Set to `false` to skip mapping computation when it hasn't changed.
+- `sort_sections=true`: Whether to sort sections by spanwise position.
+  Set to `false` for REFINE wings where section order is determined by structural connectivity.
 
 # Returns
 nothing
@@ -138,12 +140,13 @@ function reinit!(body_aero::BodyAerodynamics;
     va=[15.0, 0.0, 0.0],
     omega=zeros(MVec3),
     refine_mesh=true,
-    recompute_mapping=true
+    recompute_mapping=true,
+    sort_sections=true
 )
     idx = 1
     vec = @MVector zeros(3)
     for wing in body_aero.wings
-        reinit!(wing; refine_mesh, recompute_mapping)
+        reinit!(wing; refine_mesh, recompute_mapping, sort_sections)
         panel_props = wing.panel_props
 
         # Create panels

src/wing_geometry.jl

@@ -302,7 +302,7 @@ function Wing(n_panels::Int;
 end
 
 """
-    reinit!(wing::AbstractWing; refine_mesh=true, recompute_mapping=true)
+    reinit!(wing::AbstractWing; refine_mesh=true, recompute_mapping=true, sort_sections=true)
 
 Reinitialize wing geometry and panel properties.
 
@@ -311,11 +311,13 @@ Reinitialize wing geometry and panel properties.
   `deform!()` to preserve deformed geometry while updating panel properties.
 - `recompute_mapping::Bool=true`: Whether to recompute the refined panel mapping.
   Set to `false` to skip mapping computation when it hasn't changed.
+- `sort_sections::Bool=true`: Whether to sort sections by spanwise position.
+  Set to `false` for REFINE wings where section order is determined by structural connectivity.
 """
-function reinit!(wing::AbstractWing; refine_mesh=true, recompute_mapping=true)
+function reinit!(wing::AbstractWing; refine_mesh=true, recompute_mapping=true, sort_sections=true)
     # Refine mesh unless explicitly disabled (e.g., to preserve deformation)
     if refine_mesh
-        refine_aerodynamic_mesh!(wing; recompute_mapping)
+        refine_aerodynamic_mesh!(wing; recompute_mapping, sort_sections)
     end
 
     # Calculate panel properties
@@ -589,19 +591,22 @@ function update_non_deformed_sections!(wing::AbstractWing)
 end
 
 """
-    refine_aerodynamic_mesh!(wing::AbstractWing; recompute_mapping=true)
+    refine_aerodynamic_mesh!(wing::AbstractWing; recompute_mapping=true, sort_sections=true)
 
 Refine the aerodynamic mesh of the wing based on spanwise panel distribution.
 
 # Keyword Arguments
 - `recompute_mapping::Bool=true`: Whether to recompute the refined panel mapping.
   Set to `false` to skip mapping computation when it hasn't changed.
+- `sort_sections::Bool=true`: Whether to sort sections by spanwise position (y-coordinate).
+  Set to `false` for REFINE wings where section order is determined by structural connectivity.
 
 Returns:
     Vector{Section}: List of refined sections
 """
-function refine_aerodynamic_mesh!(wing::AbstractWing; recompute_mapping=true)
-    sort!(wing.sections, by=s -> s.LE_point[2], rev=true)
+function refine_aerodynamic_mesh!(wing::AbstractWing; recompute_mapping=true, sort_sections=true)
+    # Only sort sections if requested (skip for REFINE wings with fixed structural order)
+    sort_sections && sort!(wing.sections, by=s -> s.LE_point[2], rev=true)
     n_sections = wing.n_panels + 1
     if length(wing.refined_sections) == 0
         if wing.spanwise_distribution == UNCHANGED || length(wing.sections) == n_sections