Add csv functions

Author: 1-Bart-1

src/polars.jl

@@ -165,4 +165,77 @@ function create_polars(; dat_path, polar_path, wind_vel, area, width, crease_fra
     serialize(polar_path, (alpha_range, delta_range, cl_matrix, cd_matrix, cm_matrix))
 
     toc()
-end
+end
+
+
+"""
+    write_aero_matrix(filepath::String, matrix::Matrix{Float64}, 
+                      alpha_range::Vector{Float64}, delta_range::Vector{Float64};
+                      label::String="C_l")
+
+Write an aerodynamic coefficient matrix to CSV with angle labels.
+The first row contains flap deflection angles, first column contains angles of attack.
+
+# Arguments
+- `filepath`: Path to output CSV file
+- `matrix`: Matrix of aerodynamic coefficients
+- `alpha_range`: Vector of angle of attack values in radians
+- `delta_range`: Vector of flap deflection angles in radians
+- `label`: Coefficient label for the header (default: "C_l")
+"""
+function write_aero_matrix(filepath::String, matrix::Matrix{Float64}, 
+                         alpha_range::Vector{Float64}, delta_range::Vector{Float64};
+                         label::String="C_l")
+    open(filepath, "w") do io
+        # Write header with delta values
+        println(io, "$label/delta," * join(["δ=$(round(rad2deg(δ), digits=1))°" for δ in delta_range], ","))
+        
+        # Write data rows with alpha values and coefficients
+        for i in eachindex(alpha_range)
+            row = "α=$(round(rad2deg(alpha_range[i]), digits=1))°," * join(matrix[i,:], ",")
+            println(io, row)
+        end
+    end
+end
+
+"""
+    read_aero_matrix(filepath::String) -> (Matrix{Float64}, Vector{Float64}, Vector{Float64})
+
+Read an aerodynamic coefficient matrix from CSV with angle labels.
+Returns the coefficient matrix and corresponding angle ranges.
+
+# Returns
+- `matrix`: Matrix of aerodynamic coefficients
+- `alpha_range`: Vector of angle of attack values in radians
+- `delta_range`: Vector of flap deflection angles in radians
+"""
+function read_aero_matrix(filepath::String)
+    lines = readlines(filepath)
+    
+    # Parse header to get delta values
+    header = split(lines[1], ',')
+    delta_values = map(header[2:end]) do δ_str
+        # Extract number between "δ=" and "°"
+        m = match(r"δ=(-?\d+\.?\d*)°", δ_str)
+        deg2rad(parse(Float64, m.captures[1]))
+    end
+    
+    # Initialize matrix
+    n_rows = length(lines) - 1  # Subtract header
+    n_cols = length(delta_values)
+    matrix = zeros(n_rows, n_cols)
+    alpha_values = zeros(n_rows)
+    
+    # Parse data rows
+    for (i, line) in enumerate(lines[2:end])
+        entries = split(line, ',')
+        # Extract alpha value
+        m = match(r"α=(-?\d+\.?\d*)°", entries[1])
+        alpha_values[i] = deg2rad(parse(Float64, m.captures[1]))
+        # Parse coefficient values
+        matrix[i,:] .= parse.(Float64, entries[2:end])
+    end
+    
+    return matrix, alpha_values, delta_values
+end
+

src/ram_geometry.jl

@@ -473,38 +473,26 @@ function RamAirWing(
     # Load or create polars
     (!endswith(dat_path, ".dat")) && (dat_path *= ".dat")
     (!isfile(dat_path)) && error("DAT file not found: $dat_path")
-    polar_path = dat_path[1:end-4] * "_polar.bin"
+    polar_path = dat_path[1:end-4] * "_polar.csv"
 
     (!endswith(obj_path, ".obj")) && (obj_path *= ".obj")
     (!isfile(obj_path)) && error("OBJ file not found: $obj_path")
-    info_path = obj_path[1:end-4] * "_info.bin"
 
-    if !ispath(info_path)
-        ! prn || @info "Reading $obj_path"
-        vertices, faces = read_faces(obj_path)
-        T_cad_body = center_to_com!(vertices, faces)
-        inertia_tensor = calculate_inertia_tensor(vertices, faces, mass, zeros(3))
+    ! prn || @info "Reading $obj_path"
+    vertices, faces = read_faces(obj_path)
+    T_cad_body = center_to_com!(vertices, faces)
+    inertia_tensor = calculate_inertia_tensor(vertices, faces, mass, zeros(3))
 
-        if align_to_principal
-            inertia_tensor, R_cad_body = calc_inertia_y_rotation(inertia_tensor)
-        else
-            R_cad_body = I(3)
-        end
-        circle_center_z, radius, gamma_tip = find_circle_center_and_radius(vertices)
-        le_interp, te_interp, area_interp = create_interpolations(vertices, circle_center_z, radius, gamma_tip, R_cad_body; interp_steps)
-
-        ! prn || @info "Writing $info_path"
-        serialize(info_path, (inertia_tensor, T_cad_body, R_cad_body, radius, gamma_tip, 
-            le_interp, te_interp, area_interp))
+    if align_to_principal
+        inertia_tensor, R_cad_body = calc_inertia_y_rotation(inertia_tensor)
+    else
+        R_cad_body = I(3)
     end
+    circle_center_z, radius, gamma_tip = find_circle_center_and_radius(vertices)
+    le_interp, te_interp, area_interp = create_interpolations(vertices, circle_center_z, radius, gamma_tip, R_cad_body; interp_steps)
 
-    ! prn || @info "Loading kite info from $info_path and polars from $polar_path"
+    ! prn || @info "Loading 2d polars from $polar_path"
     try
-        (inertia_tensor::Matrix, T_cad_body::Vector, R_cad_body::Matrix,
-            radius::Real, gamma_tip::Real, le_interp, te_interp, area_interp) = deserialize(info_path)
-
-        ((R_cad_body == I(3)) == align_to_principal) && throw(ArgumentError("Delete $info_path and try again."))
-
         if !ispath(polar_path)
             width = 2gamma_tip * radius
             area = area_interp(gamma_tip)