RCAIDE-GUI 1.0.0__py3-none-any.whl

tabs/solve/plots/mission/plot_flight_trajectory.py

@@ -0,0 +1,172 @@
+# RCAIDE_GUI/tabs/solve/plots/mission/plot_flight_trajectory.py
+
+# RCAIDE imports
+from RCAIDE.Framework.Core import Units
+from RCAIDE.Library.Plots.Common import set_axes, plot_style
+
+# Python imports
+import matplotlib.pyplot as plt
+import matplotlib.cm as cm
+import numpy as np
+
+# ------------------------------------------------------------------------------
+# Plot Flight Trajectory
+# ------------------------------------------------------------------------------
+def plot_flight_trajectory(results,
+                           line_color = 'bo-',
+                           line_color2 = 'rs--',
+                           save_figure = False,
+                           show_legend = True,
+                           save_filename = "Flight_Trajectory",
+                           file_type = ".png",
+                           width = 11, height = 7):
+    """
+    Creates a multi-panel visualization of aircraft trajectory and mission profile.
+
+    Parameters
+    ----------
+    results : Results
+        RCAIDE results data structure containing:
+        
+        - segments[i].conditions.frames.inertial.time[:,0]
+            Time history for each segment
+        - segments[i].conditions.frames.inertial.aircraft_range[:,0]
+            Range history for each segment
+        - segments[i].conditions.frames.inertial.position_vector[:,0:3]
+            3D position vector containing:
+                - [:,0]: x-position
+                - [:,1]: y-position
+                - [:,2]: z-position
+        - segments[i].tag
+            Name/identifier of each segment
+            
+    line_color : str, optional
+        Primary line style specification (default: 'bo-')
+        
+    line_color2 : str, optional
+        Secondary line style specification (default: 'rs--')
+        
+    save_figure : bool, optional
+        Flag for saving the figure (default: False)
+        
+    show_legend : bool, optional
+        Flag to display segment legend (default: True)
+        
+    save_filename : str, optional
+        Name of file for saved figure (default: "Flight_Trajectory")
+        
+    file_type : str, optional
+        File extension for saved figure (default: ".png")
+        
+    width : float, optional
+        Figure width in inches (default: 11)
+        
+    height : float, optional
+        Figure height in inches (default: 7)
+
+    Returns
+    -------
+    fig : matplotlib.figure.Figure
+        Handle to the generated figure containing four subplots:
+            - Range vs time
+            - Top view (x-y plane)
+            - Altitude vs time
+            - 3D trajectory
+
+    Notes
+    -----
+    Creates a four-panel plot showing:
+        1. Mission range profile
+        2. Top-down view of flight path
+        3. Altitude profile
+        4. Complete 3D trajectory
+    
+    **Major Assumptions**
+    
+    * Time is in minutes
+    * Range is in nautical miles
+    * Position coordinates are in meters
+    * Z-axis points downward in inertial frame
+    
+    **Definitions**
+    
+    'Range'
+        Total ground distance covered
+    'Ground Track'
+        Projection of flight path onto x-y plane
+    'Altitude Profile'
+        Variation of height with time
+    
+    See Also
+    --------
+    RCAIDE.Library.Plots.Mission.plot_flight_conditions : Detailed flight condition analysis
+    RCAIDE.Library.Plots.Mission.plot_aircraft_velocities : Aircraft speed profiles
+    """	     
+
+    # get plotting style 
+    ps      = plot_style()  
+
+    parameters = {'axes.labelsize': ps.axis_font_size,
+                  'xtick.labelsize': ps.axis_font_size,
+                  'ytick.labelsize': ps.axis_font_size,
+                  'axes.titlesize': ps.title_font_size}
+    plt.rcParams.update(parameters) 
+        
+    fig = plt.figure(save_filename)
+    fig.set_size_inches(width,height) 
+     
+    # get line colors for plots 
+    line_colors   = cm.inferno(np.linspace(0,0.9,len(results.segments)))    
+     
+    for i in range(len(results.segments)): 
+        time     = results.segments[i].conditions.frames.inertial.time[:,0] / Units.min
+        Range    = results.segments[i].conditions.frames.inertial.aircraft_range[:,0]/Units.nmi
+        x        = results.segments[i].conditions.frames.inertial.position_vector[:,0]  
+        y        = results.segments[i].conditions.frames.inertial.position_vector[:,1] 
+        z        = -results.segments[i].conditions.frames.inertial.position_vector[:,2] 
+
+        segment_tag  =  results.segments[i].tag
+        segment_name = segment_tag.replace('_', ' ')
+        
+        axes = plt.subplot(2,2,1)
+        axes.plot( time , Range, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width , label = segment_name)
+        axes.set_ylabel('Distance (nmi)')
+        axes.set_xlabel('Time (min)')
+        set_axes(axes)            
+
+        axes = plt.subplot(2,2,2)
+        axes.plot(x, y , line_color, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width )
+        axes.set_xlabel('x (m)')
+        axes.set_ylabel('y (m)')
+        set_axes(axes)
+
+        axes = plt.subplot(2,2,3)
+        axes.plot( time , z, line_color , color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width )
+        axes.set_ylabel('z (m)')
+        axes.set_xlabel('Time (min)')
+        set_axes(axes)   
+        
+        axes = plt.subplot(2,2,4, projection='3d') 
+        axes.scatter(x, y, z, marker='o',color =  line_colors[i])
+        axes.set_xlabel('x')
+        axes.set_ylabel('y')
+        axes.set_zlabel('z')
+        axes.set_box_aspect([1,1,1])
+        set_axes(axes)         
+        
+    if show_legend:        
+        leg =  fig.legend(bbox_to_anchor=(0.5, 0.95), loc='upper center', ncol = 4) 
+        leg.set_title('Flight Segment', prop={'size': ps.legend_font_size, 'weight': 'heavy'})    
+    
+    # Adjusting the sub-plots for legend 
+    fig.tight_layout()
+    fig.subplots_adjust(top=0.8)
+    
+    # set title of plot 
+    title_text    = 'Flight Trajectory'      
+    fig.suptitle(title_text)
+    
+    if save_figure:
+        plt.savefig(save_filename + file_type)   
+             
+    return fig         

tabs/solve/plots/stability/__init__.py

@@ -0,0 +1,5 @@
+# RCAIDE_GUI/tabs/solve/plots/stability/__init__.py
+
+from .plot_flight_forces_and_moments import plot_flight_forces_and_moments
+from .plot_longitudinal_stability    import plot_longitudinal_stability
+from .plot_lateral_stability         import plot_lateral_stability

tabs/solve/plots/stability/plot_flight_forces_and_moments.py

@@ -0,0 +1,163 @@
+# RCAIDE_GUI/tabs/solve/plots/stability/plot_flight_forces_and_moments.py
+
+# RCAIDE imports
+from RCAIDE.Framework.Core import Units
+from RCAIDE.Library.Plots.Common import set_axes, plot_style
+
+# Python imports
+import matplotlib.pyplot as plt
+import matplotlib.cm as cm
+import numpy as np
+
+# ------------------------------------------------------------------------------
+# Plot Flight Forces and Moments
+# ------------------------------------------------------------------------------
+def plot_flight_forces_and_moments(results,
+                             save_figure = False,
+                             show_legend = True,
+                             save_filename = "Aerodynamic_Wind_Frame_Forces_and_Moments",
+                             file_type = ".png",
+                             width = 12, height = 8):
+    """
+    Creates a multi-panel visualization of aircraft forces and moments in the inertial frame.
+
+    Parameters
+    ----------
+    results : Results
+        RCAIDE results data structure containing:
+            - segments[i].conditions.frames.inertial.time[:,0]
+                Time history for each segment
+            - segments[i].conditions.frames.inertial.total_force_vector[:,0:3]
+                Force components [X, Y, Z] in inertial frame
+            - segments[i].conditions.frames.inertial.total_moment_vector[:,0:3]
+                Moment components [L, M, N] in inertial frame
+            - segments[i].tag
+                Name/identifier of each segment
+            
+    save_figure : bool, optional
+        Flag for saving the figure (default: False)
+        
+    show_legend : bool, optional
+        Flag to display segment legend (default: True)
+        
+    save_filename : str, optional
+        Name of file for saved figure (default: "Inertial_Forces_and_Moments")
+        
+    file_type : str, optional
+        File extension for saved figure (default: ".png")
+        
+    width : float, optional
+        Figure width in inches (default: 12)
+        
+    height : float, optional
+        Figure height in inches (default: 8)
+
+    Returns
+    -------
+    fig : matplotlib.figure.Figure
+
+    Notes
+    -----
+    Creates visualization showing:
+        * All force components in inertial frame
+        * All moment components in inertial frame
+        * Time history for each segment
+        * Color-coded flight segments
+    
+    **Major Assumptions**
+        * Forces are in Newtons
+        * Moments are in Newton-meters
+        * Time is in minutes
+        * Forces/moments are total (aerodynamic + propulsive + gravity)
+
+    **Definitions**
+    
+    'Inertial Forces'
+        Forces in earth-fixed reference frame
+    'Inertial Moments'
+        Moments in earth-fixed reference frame
+    'Flight Segment'
+        Distinct portion of mission profile
+    
+    See Also
+    --------
+    RCAIDE.Library.Plots.Stability.plot_longitudinal_stability : Longitudinal dynamics analysis
+    RCAIDE.Library.Plots.Stability.plot_lateral_stability : Lateral-directional dynamics analysis
+    """
+
+    # get plotting style 
+    ps      = plot_style()  
+
+    parameters = {'axes.labelsize': ps.axis_font_size,
+                  'xtick.labelsize': ps.axis_font_size,
+                  'ytick.labelsize': ps.axis_font_size,
+                  'axes.titlesize': ps.title_font_size}
+    plt.rcParams.update(parameters)
+     
+    # get line colors for plots 
+    line_colors   = cm.inferno(np.linspace(0,0.9,len(results.segments)))     
+  
+    
+    fig   = plt.figure(save_filename)
+    fig.set_size_inches(width,height) 
+
+    axis_1 = plt.subplot(3,2,1)
+    axis_2 = plt.subplot(3,2,2)
+    axis_3 = plt.subplot(3,2,3)
+    axis_4 = plt.subplot(3,2,4)
+    axis_5 = plt.subplot(3,2,5)  
+    axis_6 = plt.subplot(3,2,6)    
+    
+    for i in range(len(results.segments)): 
+        time   = results.segments[i].conditions.frames.inertial.time[:,0] / Units.min 
+        X = results.segments[i].conditions.frames.wind.total_force_vector[:,0] # Changed to wind form inertial 
+        Y = results.segments[i].conditions.frames.wind.total_force_vector[:,1]
+        Z = results.segments[i].conditions.frames.wind.total_force_vector[:,2]
+        L = results.segments[i].conditions.frames.wind.total_moment_vector[:,0]
+        M = results.segments[i].conditions.frames.wind.total_moment_vector[:,1]
+        N = results.segments[i].conditions.frames.wind.total_moment_vector[:,2]
+        
+                       
+        segment_tag  =  results.segments[i].tag
+        segment_name = segment_tag.replace('_', ' ') 
+        axis_1.plot(time,X, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width, label = segment_name) 
+        axis_1.set_ylabel(r'X (N)') 
+        set_axes(axis_1)                
+        
+        axis_2.set_ylabel(r'L (Nm)')
+        axis_2.plot(time,L, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)  
+        set_axes(axis_2)     
+        
+        axis_3.plot(time,Y, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width) 
+        axis_3.set_ylabel(r'Y (N)') 
+        set_axes(axis_3)  
+        
+        axis_4.plot(time,M, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width) 
+        axis_4.set_ylabel(r'M (Nm)') 
+        set_axes(axis_4) 
+
+        axis_5.plot(time, Z, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_5.set_xlabel('Time (mins)')
+        axis_5.set_ylabel(r'Z (N)') 
+        set_axes(axis_5)    
+
+        axis_6.plot(time, N, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_6.set_xlabel('Time (mins)')
+        axis_6.set_ylabel(r'N (Nm)') 
+        set_axes(axis_6)   
+ 
+    if show_legend:    
+        leg =  fig.legend(bbox_to_anchor=(0.5, 0.95), loc='upper center', ncol = 4) 
+        leg.set_title('Flight Segment', prop={'size': ps.legend_font_size, 'weight': 'heavy'})    
+    
+    # Adjusting the sub-plots for legend 
+    fig.tight_layout()
+    fig.subplots_adjust(top=0.8) 
+    
+    # set title of plot 
+    title_text    = 'Intertial Forces and Moments'      
+    fig.suptitle(title_text)
+    
+    if save_figure:
+        plt.savefig(save_filename + file_type)   
+    return fig 

tabs/solve/plots/stability/plot_lateral_stability.py

@@ -0,0 +1,151 @@
+# RCAIDE/Library/Plots/Performance/Stability/plot_lateral_stability.py
+# 
+# 
+# Created:  Jul 2023, M. Clarke 
+
+# ----------------------------------------------------------------------------------------------------------------------
+#  IMPORT
+# ----------------------------------------------------------------------------------------------------------------------  
+from RCAIDE.Framework.Core import Units
+from RCAIDE.Library.Plots.Common import set_axes, plot_style
+import matplotlib.pyplot as plt
+import matplotlib.cm as cm
+import numpy as np  
+
+# ----------------------------------------------------------------------------------------------------------------------
+#  PLOTS
+# ----------------------------------------------------------------------------------------------------------------------   
+def plot_lateral_stability(results,
+                             save_figure = False,
+                             show_legend = True,
+                             save_filename = "Lateral_Stability",
+                             file_type = ".png",
+                             width = 11, height = 7):
+    """
+    Creates a multi-panel visualization of lateral-directional stability characteristics.
+
+    Parameters
+    ----------
+    results : Results
+        RCAIDE results data structure containing:
+            - segments[i].conditions.frames.inertial.time[:,0]
+                Time history for each segment
+            - segments[i].conditions.aerodynamics.angles.phi[:,0]
+                Bank angle history
+            - segments[i].conditions.control_surfaces.aileron.deflection[:,0]
+                Aileron deflection history
+            - segments[i].conditions.control_surfaces.rudder.deflection[:,0]
+                Rudder deflection history
+            - segments[i].tag
+                Name/identifier of each segment
+            
+    save_figure : bool, optional
+        Flag for saving the figure (default: False)
+        
+    show_legend : bool, optional
+        Flag to display segment legend (default: True)
+        
+    save_filename : str, optional
+        Name of file for saved figure (default: "Lateral_Stability")
+        
+    file_type : str, optional
+        File extension for saved figure (default: ".png")
+        
+    width : float, optional
+        Figure width in inches (default: 11)
+        
+    height : float, optional
+        Figure height in inches (default: 7)
+
+    Returns
+    -------
+    fig : matplotlib.figure.Figure
+        Handle to the generated figure containing three subplots:
+            * Bank angle vs time
+            * Aileron deflection vs time
+            * Rudder deflection vs time
+
+    Notes
+    -----
+    Creates visualization showing:
+        * Roll attitude response
+        * Lateral control inputs
+        * Directional control inputs
+        * Time history for each segment
+    
+    **Major Assumptions**
+        * Angles are in degrees
+        * Time is in minutes
+        * Positive deflections follow right-hand rule
+    
+    **Definitions**
+    
+    'Bank Angle'
+        Roll attitude relative to horizon
+    'Aileron Deflection'
+        Roll control surface position
+    'Rudder Deflection'
+        Yaw control surface position
+    
+    See Also
+    --------
+    RCAIDE.Library.Plots.Stability.plot_longitudinal_stability : Longitudinal stability analysis
+    RCAIDE.Library.Plots.Stability.plot_flight_forces_and_moments : Force/moment visualization
+    """  
+
+    # get plotting style 
+    ps      = plot_style()  
+
+    parameters = {'axes.labelsize': ps.axis_font_size,
+                  'xtick.labelsize': ps.axis_font_size,
+                  'ytick.labelsize': ps.axis_font_size,
+                  'axes.titlesize': ps.title_font_size}
+    plt.rcParams.update(parameters)
+     
+    # get line colors for plots 
+    line_colors   = cm.inferno(np.linspace(0,0.9,len(results.segments)))     
+     
+    fig   = plt.figure(save_filename)
+    fig.set_size_inches(width,height) 
+    axis_1 = plt.subplot(2,2,1) 
+    axis_2 = plt.subplot(2,2,2)  
+    axis_3 = plt.subplot(2,2,3)    
+    
+    for i in range(len(results.segments)): 
+        time     = results.segments[i].conditions.frames.inertial.time[:,0] / Units.min  
+        phi      = -results.segments[i].conditions.aerodynamics.angles.phi[:,0] / Units.deg          
+        delta_a  = results.segments[i].conditions.control_surfaces.aileron.deflection[:,0] / Units.deg  
+        delta_r  = results.segments[i].conditions.control_surfaces.rudder.deflection[:,0] / Units.deg   
+          
+        segment_tag  =  results.segments[i].tag
+        segment_name = segment_tag.replace('_', ' ')
+        
+        axis_1.plot(time, phi, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width, label = segment_name)
+        axis_1.set_ylabel(r'$Bank Angle (deg)$') 
+        set_axes(axis_1)     
+
+        axis_2.plot(time,delta_a , color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_2.set_xlabel('Time (mins)')
+        axis_2.set_ylabel(r'Aileron Defl. (deg)')
+        set_axes(axis_2)  
+
+        axis_3.plot(time,delta_r , color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_3.set_xlabel('Time (mins)')
+        axis_3.set_ylabel(r'Rudder Defl. (deg)')
+        set_axes(axis_3)         
+         
+    if show_legend:
+        leg =  fig.legend(bbox_to_anchor=(0.5, 0.95), loc='upper center', ncol = 4) 
+        leg.set_title('Flight Segment', prop={'size': ps.legend_font_size, 'weight': 'heavy'})    
+    
+    # Adjusting the sub-plots for legend
+    fig.tight_layout()
+    fig.subplots_adjust(top=0.8) 
+    
+    # set title of plot 
+    title_text    = 'Stability Coeffiicents'      
+    fig.suptitle(title_text)
+ 
+    if save_figure:
+        plt.savefig(save_filename + file_type)   
+    return fig 

tabs/solve/plots/stability/plot_longitudinal_stability.py

@@ -0,0 +1,173 @@
+# RCAIDE/Library/Plots/Performance/Stability/plot_longitudinal_stability.py
+# 
+# 
+# Created:  Jul 2023, M. Clarke 
+
+# ----------------------------------------------------------------------------------------------------------------------
+#  IMPORT
+# ----------------------------------------------------------------------------------------------------------------------  
+from RCAIDE.Framework.Core import Units
+from RCAIDE.Library.Plots.Common import set_axes, plot_style
+import matplotlib.pyplot as plt
+import matplotlib.cm as cm
+import numpy as np  
+
+# ----------------------------------------------------------------------------------------------------------------------
+#  PLOTS
+# ----------------------------------------------------------------------------------------------------------------------   
+def plot_longitudinal_stability(results,
+                             save_figure = False,
+                             show_legend = True,
+                             save_filename = "Longitudinal_Stability",
+                             file_type = ".png",
+                             width = 11, height = 7):
+    """
+    Creates a multi-panel visualization of longitudinal stability characteristics.
+
+    Parameters
+    ----------
+    results : Results
+        RCAIDE results data structure containing:
+            - segments[i].conditions.frames.inertial.time[:,0]
+                Time history for each segment
+            - segments[i].conditions.aerodynamics.angles.theta[:,0]
+                Pitch angle history
+            - segments[i].conditions.aerodynamics.angles.alpha[:,0]
+                Angle of attack history
+            - segments[i].conditions.control_surfaces.elevator.deflection[:,0]
+                Elevator deflection history
+            - segments[i].tag
+                Name/identifier of each segment
+            
+    save_figure : bool, optional
+        Flag for saving the figure (default: False)
+        
+    show_legend : bool, optional
+        Flag to display segment legend (default: True)
+        
+    save_filename : str, optional
+        Name of file for saved figure (default: "Longitudinal_Stability")
+        
+    file_type : str, optional
+        File extension for saved figure (default: ".png")
+        
+    width : float, optional
+        Figure width in inches (default: 11)
+        
+    height : float, optional
+        Figure height in inches (default: 7)
+
+    Returns
+    -------
+    fig : matplotlib.figure.Figure
+
+    Notes
+    -----
+    Creates visualization showing:
+        * Pitch attitude response
+        * Aerodynamic angle evolution
+        * Longitudinal control inputs
+        * Trajectory angle history
+        * Time history for each segment
+    
+    **Major Assumptions**
+        * Angles are in degrees
+        * Time is in minutes
+        * Positive deflections follow right-hand rule
+   
+    **Definitions**
+    
+    'Pitch Angle'
+        Nose-up/down attitude relative to horizon
+    'Angle of Attack'
+        Angle between velocity vector and body reference line
+    'Flight Path Angle'
+        Angle between velocity vector and horizon
+    'Elevator Deflection'
+        Pitch control surface position
+    
+    See Also
+    --------
+    RCAIDE.Library.Plots.Stability.plot_lateral_stability : Lateral-directional stability analysis
+    RCAIDE.Library.Plots.Stability.plot_flight_forces_and_moments : Force/moment visualization
+    """  
+
+    # get plotting style 
+    ps      = plot_style()  
+
+    parameters = {'axes.labelsize': ps.axis_font_size,
+                  'xtick.labelsize': ps.axis_font_size,
+                  'ytick.labelsize': ps.axis_font_size,
+                  'axes.titlesize': ps.title_font_size}
+    plt.rcParams.update(parameters)
+     
+    # get line colors for plots 
+    line_colors   = cm.inferno(np.linspace(0,0.9,len(results.segments)))     
+     
+    fig   = plt.figure(save_filename)
+    fig.set_size_inches(width,height)
+    
+
+    axis_1 = plt.subplot(3,2,1)    
+    axis_2 = plt.subplot(3,2,2) 
+    axis_3 = plt.subplot(3,2,3)  
+    axis_4 = plt.subplot(3,2,4)  
+    axis_5 = plt.subplot(3,2,5)
+    axis_6 = plt.subplot(3,2,6) 
+    
+    for i in range(len(results.segments)): 
+        time       = results.segments[i].conditions.frames.inertial.time[:,0] / Units.min
+        c_m        = results.segments[i].conditions.static_stability.coefficients.M[:,0]   
+        SM         = results.segments[i].conditions.static_stability.static_margin[:,0]  
+        delta_e    = results.segments[i].conditions.control_surfaces.elevator.deflection[:,0] / Units.deg
+        CM_delta_e = results.segments[i].conditions.static_stability.derivatives.CM_delta_e[:,0]
+        Cm_alpha   = results.segments[i].conditions.static_stability.derivatives.CM_alpha[:,0]
+        CL_alpha   = results.segments[i].conditions.static_stability.derivatives.Clift_alpha[:,0] 
+          
+        segment_tag  =  results.segments[i].tag
+        segment_name = segment_tag.replace('_', ' ')  
+        
+        axis_1.plot(time, c_m, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width, label = segment_name) 
+        axis_1.set_ylabel(r'$C_M$') 
+        set_axes(axis_1) 
+
+        axis_2.plot(time, Cm_alpha, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width) 
+        axis_2.set_ylabel(r'$C_M\alpha$') 
+        set_axes(axis_2) 
+        
+        axis_3.plot(time,SM , color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_3.set_xlabel('Time (mins)')
+        axis_3.set_ylabel(r'Static Margin (%)')
+        set_axes(axis_3)  
+
+        axis_4.plot(time,delta_e , color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_4.set_xlabel('Time (mins)')
+        axis_4.set_ylabel(r'Elevator Defl.n')  
+        set_axes(axis_4) 
+        
+        axis_5.plot(time,CM_delta_e , color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_5.set_xlabel('Time (mins)')
+        axis_5.set_ylabel(r'$C_M\delta_e$') 
+        set_axes(axis_5)
+        
+        axis_6.plot(time,CL_alpha, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_6.set_xlabel('Time (mins)')
+        axis_6.set_ylabel(r'$C_L\alpha$') 
+        set_axes(axis_6)    
+        
+    if show_legend:
+        leg =  fig.legend(bbox_to_anchor=(0.5, 0.95), loc='upper center', ncol = 4) 
+        leg.set_title('Flight Segment', prop={'size': ps.legend_font_size, 'weight': 'heavy'})     
+          
+    # Adjusting the sub-plots for legend
+    fig.tight_layout()
+    fig.subplots_adjust(top=0.8)
+    
+    # set title of plot 
+    title_text    = 'Stability Coefficients'      
+    fig.suptitle(title_text)
+ 
+    if save_figure:
+        plt.savefig(save_filename + file_type)   
+    return fig 
+ 

tabs/solve/plots/thermal_management/__init__.py

@@ -0,0 +1,7 @@
+# RCAIDE_GUI/tabs/solve/plots/thermal_management/__init__.py
+
+from .plot_thermal_management_performance             import plot_thermal_management_performance
+from .plot_wavy_channel_conditions                   import plot_wavy_channel_conditions
+from .plot_cross_flow_heat_exchanger_conditions      import plot_cross_flow_heat_exchanger_conditions
+from .plot_reservoir_conditions                      import plot_reservoir_conditions
+from .plot_air_cooled_conditions                     import plot_air_cooled_conditions