RCAIDE-GUI 1.0.0__py3-none-any.whl

tabs/solve/plots/emissions/plot_emissions.py

@@ -0,0 +1,167 @@
+# RCAIDE_GUI/tabs/solve/plots/emissions/plot_emissions.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 Emissions
+# ------------------------------------------------------------------------------
+def plot_emissions(results,
+                    save_figure = False,
+                    show_legend = True,
+                    save_filename = "Emissions" ,
+                    file_type = ".png",
+                    width = 11, height = 5):
+    """
+    Generate plots showing CO2-equivalent emissions and emission indexes for various fuel species over mission segments.
+
+    Parameters
+    ----------
+    results : Data
+        Mission results data structure containing:
+        results.segments[i].conditions.emissions.total with fields:
+            - CO2 : array
+                Carbon dioxide emissions [kg]
+            - NOx : array
+                Nitrogen oxide emissions [kg]
+            - H2O : array
+                Water vapor emissions [kg]
+            - Contrails : array
+                Contrail formation impact [kg CO2e]
+            - Soot : array
+                Particulate emissions [kg]
+            - SO2 : array
+                Sulfur dioxide emissions [kg]
+
+    save_figure : bool, optional
+        Save figure to file if True, default False
+
+    show_legend : bool, optional
+        Display segment legend if True, default True
+
+    save_filename : str, optional
+        Name for saved figure file, default "CO2e_Emissions"
+
+    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
+        Figure showing stacked emissions contributions
+
+    Notes
+    -----
+    Creates a stacked area plot showing:
+        - Individual contributions from each emission type
+        - Cumulative total CO2-equivalent impact
+        - Breakdown by mission segment
+        - Time history of emissions
+
+    Different emission types are distinguished by fill colors
+    and segments use different shades from the inferno colormap.
+
+    **Definitions**
+
+    'CO2-equivalent (CO2e)'
+        Combined climate impact normalized to CO2
+    
+    'Global Warming Potential (GWP)'
+        Relative impact factor for different emissions
+    
+    'Contrail Impact'
+        Climate forcing from aviation-induced cloudiness
+    """
+ 
+    # 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)  
+
+    line_colors   = cm.inferno(np.linspace(0,0.9,len(results.segments)))
+    
+    cum_y0  = 0
+    cum_y1  = 0 
+    cum_y1_0 = 0  
+    
+    for i in range(len(results.segments)): 
+        time                = results.segments[i].conditions.frames.inertial.time[:, 0] / Units.min 
+        emissions_CO2       = results.segments[i].conditions.emissions.mass.CO2[:, 0] /1E3  
+        emissions_CO        = results.segments[i].conditions.emissions.mass.CO[:, 0]  /1E3  
+        emissions_NOx       = results.segments[i].conditions.emissions.mass.NOx[:, 0] /1E3  
+        emissions_H2O       = results.segments[i].conditions.emissions.mass.H2O[:, 0] /1E3   
+        EI_CO2              = results.segments[i].conditions.emissions.index.CO2[:, 0] 
+        EI_CO               = results.segments[i].conditions.emissions.index.CO[:, 0] 
+        EI_NOx              = results.segments[i].conditions.emissions.index.NOx[:, 0] 
+        EI_H2O              = results.segments[i].conditions.emissions.index.H2O[:, 0] 
+        EI_SO2              = results.segments[i].conditions.emissions.index.SO2[:, 0]        
+
+        cum_y0 = np.zeros_like(emissions_CO2)  
+        cum_y1 = cum_y1_0 + emissions_CO2 + emissions_NOx  + emissions_CO +  emissions_H2O  #+ emissions_Contrails +  emissions_Soot +  emissions_SO2    
+
+        segment_tag  =  results.segments[i].tag
+        segment_name = segment_tag.replace('_', ' ') 
+        axis_1 = plt.subplot(1,2,1)    
+        axis_1.fill_between(time, cum_y0, cum_y1, where=(cum_y0 < cum_y1), color= line_colors[i],  interpolate=True, label = segment_name)   
+        cum_y1_0 = cum_y1[-1]  
+        axis_1.set_ylabel(r'CO2e Emissions (Metric Tons)') 
+        axis_1.set_xlabel(r'Time (mins)') 
+        set_axes(axis_1)
+         
+        axis_2 = plt.subplot(1,2,2)
+        if i == 0: 
+            axis_2.plot(time, EI_CO2, color = line_colors[i], marker = ps.markers[0],markersize = ps.marker_size, linewidth = ps.line_width, label = r"$CO_2$") 
+            axis_2.plot(time, EI_CO , color = line_colors[i], marker = ps.markers[1],markersize = ps.marker_size, linewidth = ps.line_width, label = r"$CO$" ) 
+            axis_2.plot(time, EI_NOx, color = line_colors[i], marker = ps.markers[2],markersize = ps.marker_size, linewidth = ps.line_width, label = r"$NO_x$") 
+            axis_2.plot(time, EI_H2O, color = line_colors[i], marker = ps.markers[3],markersize = ps.marker_size, linewidth = ps.line_width, label = r"$H_2O$") 
+            axis_2.plot(time, EI_SO2, color = line_colors[i], marker = ps.markers[4],markersize = ps.marker_size, linewidth = ps.line_width, label = r"$SO_2$")
+        else:
+
+            axis_2.semilogy(time, EI_CO2, color = line_colors[i], marker = ps.markers[0],markersize = ps.marker_size, linewidth = ps.line_width) 
+            axis_2.semilogy(time, EI_CO , color = line_colors[i], marker = ps.markers[1],markersize = ps.marker_size, linewidth = ps.line_width)
+            axis_2.semilogy(time, EI_NOx, color = line_colors[i], marker = ps.markers[2],markersize = ps.marker_size, linewidth = ps.line_width) 
+            axis_2.semilogy(time, EI_H2O, color = line_colors[i], marker = ps.markers[3],markersize = ps.marker_size, linewidth = ps.line_width) 
+            axis_2.semilogy(time, EI_SO2, color = line_colors[i], marker = ps.markers[4],markersize = ps.marker_size, linewidth = ps.line_width)             
+    
+        axis_2.set_ylabel(r'Emission Index') 
+        axis_2.set_xlabel(r'Time (mins)') 
+        axis_2.minorticks_on()
+        axis_2.grid(which='major', linestyle='-', linewidth=0.5, color='grey')
+        axis_2.grid(which='minor', linestyle=':', linewidth=0.5, color='grey')
+        axis_2.grid(True)         
+                
+    if show_legend:
+        leg =  fig.legend(bbox_to_anchor=(0.5, 0.95), loc='upper center', ncol = 6) 
+        leg.set_title('Flight Segment', prop={'size': ps.legend_font_size, 'weight': 'heavy'})    
+    
+    # Adjusting the sub-plots for legend 
+    
+    # set title of plot 
+    title_text    = 'Emissions'     
+    fig.tight_layout() 
+    fig.suptitle(title_text)
+    
+    fig.subplots_adjust(top=0.7)
+    
+    if save_figure:
+        plt.savefig(save_filename + file_type)   
+    return fig 

tabs/solve/plots/energy/__init__.py

@@ -0,0 +1,12 @@
+# RCAIDE_GUI/tabs/solve/plots/energy/__init__.py
+
+from .plot_battery_module_conditions            import plot_battery_module_conditions
+from .plot_battery_cell_conditions              import plot_battery_cell_conditions
+from .plot_battery_degradation                  import plot_battery_degradation
+from .plot_battery_temperature                  import plot_battery_temperature
+from .plot_battery_module_C_rates               import plot_battery_module_C_rates
+from .plot_battery_pack_conditions              import plot_battery_pack_conditions
+from .plot_battery_ragone_diagram               import plot_battery_ragone_diagram
+from .plot_electric_propulsor_efficiencies      import plot_electric_propulsor_efficiencies
+from .plot_altitude_sfc_weight                  import plot_altitude_sfc_weight
+from .plot_propulsor_throttles                  import plot_propulsor_throttles

tabs/solve/plots/energy/plot_altitude_sfc_weight.py

@@ -0,0 +1,134 @@
+# RCAIDE_GUI/tabs/solve/plots/energy/plot_altitude_sfc_weight.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 Altitude SFC Weight
+# ------------------------------------------------------------------------------
+def plot_altitude_sfc_weight(results,
+                             save_figure = False,
+                             show_legend = True,
+                             save_filename = "Weight_Fuel_Consumption" ,
+                             file_type = ".png", 
+                             width = 11, height = 7):
+    """
+    Creates a four-panel plot showing throttle settings, vehicle weight, specific fuel consumption (SFC), 
+    and fuel consumption rate over time.
+
+    Parameters
+    ----------
+    results : Results
+        RCAIDE results structure containing segment data
+        
+    save_figure : bool, optional
+        Flag for saving the figure (default: False)
+        
+    show_legend : bool, optional
+        Flag for displaying plot legend (default: True)
+        
+    save_filename : str, optional
+        Name of file for saved figure (default: "Altitude_SFC_Weight")
+        
+    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
+        Figure handle containing the generated plots
+
+    Notes
+    -----
+    The function creates a 2x2 subplot containing:
+        2. Vehicle weight vs time
+        1. Fuel consumption vs time
+        3. Specific fuel consumption vs time
+        4. Fuel consumption rate vs time
+    
+    Each segment is plotted with a different color from the inferno colormap.
+  
+    **Definitions**
+    
+    'SFC'
+        Specific Fuel Consumption - measure of the fuel efficiency of an engine design 
+        with respect to thrust output
+    """
+ 
+    # 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)
+    axis_4 = plt.subplot(2,2,4)
+    
+    for i in range(len(results.segments)): 
+        time      = results.segments[i].conditions.frames.inertial.time[:, 0] / Units.min 
+        Weight    = (results.segments[i].conditions.weights.total_mass[:, 0] * 9.81)  / Units.lbf
+        mdot      = results.segments[i].conditions.weights.vehicle_mass_rate[:, 0]/ Units.lb
+        thrust    = abs(results.segments[i].conditions.frames.body.thrust_force_vector[:, 0])/ Units.lbf
+        fuel_mass = results.segments[i].conditions.energy.cumulative_fuel_consumption[:, 0]/ Units.lb
+        
+        axis_1.set_ylabel(r'Weight (lbf)')  
+        axis_2.set_ylabel(r'Fuel Consumption (lb)')
+        axis_3.set_ylabel(r'SFC (lb/lbf-hr)')
+        axis_4.set_ylabel(r'Fuel Rate (lb/s)')  
+        
+        sfc       = (mdot ) / (thrust ) * Units.hr 
+
+        segment_tag  =  results.segments[i].tag
+        segment_name = segment_tag.replace('_', ' ')       
+        axis_1.plot(time, Weight , color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width, label = segment_name )   
+        set_axes(axis_1) 
+
+        axis_2.plot(time, fuel_mass, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_2.set_xlabel('Time (mins)')
+        set_axes(axis_2) 
+
+        axis_3.plot(time, sfc, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_3.set_xlabel('Time (mins)')
+        set_axes(axis_3) 
+
+        axis_4.plot(time, mdot, color = line_colors[i], marker = ps.markers[0], linewidth = ps.line_width)
+        axis_4.set_xlabel('Time (mins)')
+        set_axes(axis_4)         
+        
+    
+    if show_legend:
+        leg =  fig.legend(bbox_to_anchor=(0.5, 0.95), loc='upper center', ncol = 4)    
+    
+    # Adjusting the sub-plots for legend 
+    fig.tight_layout()
+    fig.subplots_adjust(top=0.8)
+    
+    # set title of plot 
+    title_text    = 'Weight and Fuel Consumption'      
+    fig.suptitle(title_text)
+    
+    if save_figure:
+        plt.savefig(save_filename + file_type)   
+    return fig 

tabs/solve/plots/energy/plot_battery_cell_conditions.py

@@ -0,0 +1,154 @@
+# RCAIDE_GUI/tabs/solve/plots/energy/plot_battery_cell_conditions.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 Battery Cell Conditions
+# ------------------------------------------------------------------------------
+def plot_battery_cell_conditions(results,
+                                  save_figure = False,
+                                  show_legend = True,
+                                  save_filename = "Battery_Cell_Conditions_",
+                                  file_type = ".png",
+                                  width = 11, height = 7):
+    """
+    Creates a six-panel plot showing various battery cell-level conditions throughout flight.
+
+    Parameters
+    ----------
+    results : Results
+        RCAIDE results structure containing segment data and battery conditions
+        
+    save_figure : bool, optional
+        Flag for saving the figure (default: False)
+        
+    show_legend : bool, optional
+        Flag for displaying plot legend (default: True)
+        
+    save_filename : str, optional
+        Base name of file for saved figure (default: "Battery_Cell_Conditions_")
+        
+    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
+        Figure handle containing the generated plots
+
+    Notes
+    -----
+    The function creates a 3x2 subplot containing:
+        1. State of Charge (SOC) vs time
+        2. Cell energy vs time
+        3. Cell current vs time
+        4. Cell power vs time
+        5. Cell voltage vs time
+        6. Cell temperature vs time
+    
+    Each segment is plotted with a different color from the inferno colormap.
+    For multiple battery modules, each module is distinguished by different markers.
+    
+    **Definitions**
+    
+    'SOC'
+        State of Charge - the level of charge of a battery relative to its capacity
+    'Open Circuit Voltage'
+        The voltage across the battery terminals with no load applied
+    """ 
+    
+    # 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)))      
+    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 network in results.segments[0].analyses.vehicle.networks: 
+        busses  = network.busses 
+        for bus_i, bus in enumerate(busses):
+            for b_i, battery in enumerate(bus.battery_modules):
+                if b_i == 0 or bus.identical_battery_modules == False:
+                    for i in range(len(results.segments)):  
+                        time    = results.segments[i].conditions.frames.inertial.time[:,0] / Units.min    
+                        battery_conditions  = results.segments[i].conditions.energy.busses[bus.tag].battery_modules[battery.tag]
+                        cell_power          = battery_conditions.cell.power[:,0]
+                        cell_energy         = battery_conditions.cell.energy[:,0]
+                        cell_volts          = battery_conditions.cell.voltage_under_load[:,0]
+                        cell_volts_oc       = battery_conditions.cell.voltage_open_circuit[:,0]
+                        cell_current        = battery_conditions.cell.current[:,0]
+                        cell_SOC            = battery_conditions.cell.state_of_charge[:,0]   
+                        cell_temperature    = battery_conditions.cell.temperature[:,0]  
+                
+                        segment_tag  = results.segments[i].tag
+                        segment_name = segment_tag.replace('_', ' ') 
+    
+                        if b_i == 0 and i ==0:                    
+                            axis_1.plot(time, cell_SOC, color = line_colors[i], marker = ps.markers[bus_i], linewidth = ps.line_width, label = battery.tag)
+                        else:
+                            axis_1.plot(time, cell_SOC, color = line_colors[i], marker = ps.markers[bus_i], linewidth = ps.line_width)
+                        axis_1.set_ylabel(r'SOC')
+                        axis_1.set_ylim([0,1.1])
+                        set_axes(axis_1)     
+                         
+                        axis_2.plot(time, cell_energy/Units.Wh, color = line_colors[i], marker = ps.markers[bus_i], linewidth = ps.line_width)
+                        axis_2.set_ylabel(r'Energy (W-hr)')
+                        set_axes(axis_2) 
+                 
+                        axis_3.plot(time, cell_current, color = line_colors[i], marker = ps.markers[bus_i], linewidth = ps.line_width)
+                        axis_3.set_ylabel(r'Current (A)')
+                        set_axes(axis_3)  
+                 
+                        axis_4.plot(time, cell_power, color = line_colors[i], marker = ps.markers[bus_i], linewidth = ps.line_width)
+                        axis_4.set_ylabel(r'Power (W)')
+                        set_axes(axis_4)     
+                         
+                        axis_5.plot(time, cell_volts, color = line_colors[i], marker = ps.markers[bus_i], linewidth = ps.line_width) 
+                        axis_5.set_ylabel(r'Voltage (V)')
+                        set_axes(axis_5) 
+                 
+                        axis_6.plot(time, cell_temperature, color = line_colors[i], marker = ps.markers[bus_i], linewidth = ps.line_width)
+                        axis_6.set_ylabel(r'Temperature, $\degree$C')
+                        set_axes(axis_6)  
+                  
+    if show_legend:      
+        leg =  fig.legend(bbox_to_anchor=(0.5, 0.95), loc='upper center', ncol = 4)   
+    
+    # Adjusting the sub-plots for legend 
+    fig.tight_layout()
+    fig.subplots_adjust(top=0.8) 
+    
+    # set title of plot 
+    title_text   = 'Battery Cell Conditions'       
+    fig.suptitle(title_text) 
+    
+    if save_figure:
+        plt.savefig(save_filename + battery.tag + file_type)    
+    return fig  

tabs/solve/plots/energy/plot_battery_degradation.py

@@ -0,0 +1,157 @@
+# RCAIDE_GUI/tabs/solve/plots/energy/plot_battery_degradation.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 Battery Degradation
+# ------------------------------------------------------------------------------
+def plot_battery_degradation(results,
+                            save_figure = False,
+                            line_color = 'bo-',
+                            line_color2 = 'rs--',
+                            save_filename = "Battery_Degradation",
+                            file_type = ".png",
+                            width = 11, height = 7):
+    """
+    Creates a six-panel plot showing battery degradation metrics against various parameters.
+
+    Parameters
+    ----------
+    results : Results
+        RCAIDE results structure containing segment data and battery degradation metrics
+        
+    save_figure : bool, optional
+        Flag for saving the figure (default: False)
+        
+    line_color : str, optional
+        Matplotlib format string for first line style (default: 'bo-')
+        
+    line_color2 : str, optional
+        Matplotlib format string for second line style (default: 'rs--')
+        
+    save_filename : str, optional
+        Base name of file for saved figure (default: "Battery_Degradation")
+        
+    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
+        Figure handle containing the generated plots
+
+    Notes
+    -----
+    The function creates a 3x2 subplot showing:
+        Left column (Capacity fade vs):
+            1. Charge throughput (Ah)
+            2. Time (hours)
+            3. Time (days)
+
+        Right column (Resistance growth vs):
+            1. Charge throughput (Ah)
+            2. Time (hours)
+            3. Time (days)
+    
+    **Definitions**
+    
+    'Capacity Fade'
+        The loss of energy storage capacity over time/usage
+    'Resistance Growth'
+        The increase in internal resistance over time/usage
+    'Charge Throughput'
+        The cumulative amount of charge that has passed through the battery
+    """ 
+    # 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)
+    
+
+    for network in results.segments[0].analyses.vehicle.networks: 
+        busses  = network.busses
+        for bus in busses:
+            if bus.identical_battery_modules:
+                for i, battery in enumerate(bus.battery_modules):
+                    if i == 0:
+                        fig = plt.figure(save_filename + '_' + battery.tag)
+                        fig.set_size_inches(width,height)  
+                        num_segs          = len(results.segments)
+                        time_hrs          = np.zeros(num_segs)  
+                        capacity_fade     = np.zeros_like(time_hrs)
+                        resistance_growth = np.zeros_like(time_hrs)
+                        cycle_day         = np.zeros_like(time_hrs)
+                        charge_throughput = np.zeros_like(time_hrs)    
+                             
+                        for i in range(len(results.segments)): 
+                            time_hrs[i]    = results.segments[i].conditions.frames.inertial.time[-1,0]  / Units.hour   
+                            battery_conditions  = results.segments[i].conditions.energy.busses[bus.tag].battery_modules[battery.tag]
+                            cycle_day[i]          = battery_conditions.cell.cycle_in_day
+                            capacity_fade[i]      = battery_conditions.cell.capacity_fade_factor
+                            resistance_growth[i]  = battery_conditions.cell.resistance_growth_factor
+                            charge_throughput[i]  = battery_conditions.cell.charge_throughput[-1,0]  
+                 
+                        axis_1 = plt.subplot(3,2,1)
+                        axis_1.plot(charge_throughput, capacity_fade, color = ps.color , marker = ps.markers[0], linewidth = ps.line_width ) 
+                        axis_1.set_ylabel('$E/E_0$')
+                        axis_1.set_xlabel('Ah')
+                        set_axes(axis_1)      
+                    
+                        axis_2 = plt.subplot(3,2,3)
+                        axis_2.plot(time_hrs, capacity_fade, color = ps.color, marker = ps.markers[0], linewidth = ps.line_width ) 
+                        axis_2.set_ylabel('$E/E_0$')
+                        axis_2.set_xlabel('Time (hrs)')
+                        set_axes(axis_2)     
+                    
+                        axis_3 = plt.subplot(3,2,5)
+                        axis_3.plot(cycle_day, capacity_fade, color = ps.color, marker = ps.markers[0], linewidth = ps.line_width ) 
+                        axis_3.set_ylabel('$E/E_0$')
+                        axis_3.set_xlabel('Time (days)')
+                        set_axes(axis_3)     
+                    
+                        axis_4 = plt.subplot(3,2,2) 
+                        axis_4.plot(charge_throughput, resistance_growth, color = ps.color, marker = ps.markers[0], linewidth = ps.line_width )
+                        axis_4.set_ylabel('$R/R_0$')
+                        axis_4.set_xlabel('Ah')
+                        set_axes(axis_4)      
+                    
+                        axis_5 = plt.subplot(3,2,4) 
+                        axis_5.plot(time_hrs, resistance_growth, color = ps.color, marker = ps.markers[0], linewidth = ps.line_width )
+                        axis_5.set_ylabel('$R/R_0$')
+                        axis_5.set_xlabel('Time (hrs)')
+                        set_axes(axis_5)     
+                    
+                        axis_6 = plt.subplot(3,2,6) 
+                        axis_6.plot(cycle_day, resistance_growth, color = ps.color, marker = ps.markers[0], linewidth = ps.line_width )
+                        axis_6.set_ylabel('$R/R_0$')
+                        axis_6.set_xlabel('Time (days)')
+                        set_axes(axis_6)              
+    
+                        
+    # set title of plot 
+    title_text    = 'Battery Cell Degradation: ' + battery.tag   
+    fig.suptitle(title_text) 
+    
+    plt.tight_layout()    
+    if save_figure:    
+        fig.savefig(save_filename + '_'+ battery.tag + file_type) 
+
+    return fig 
+