fix bug

birdpressure/Bird.py

@@ -45,7 +45,7 @@ def __init__(self, initial_density=950, speed_of_sound=1482.9, mass=1.81437, asp
         self.mass = mass
         self.k = k
         self.aspect_ratio = aspect_ratio
-        self.length = (16 *self.mass / (np.pi*self.initial_density)) ** 1 / 3
+        self.length = (16 *self.mass / (np.pi*self.initial_density)) ** (1 / 3)
         self.diameter = self.length / self.aspect_ratio
 
     def to_dict(self):

birdpressure/Timing.py

@@ -202,7 +202,7 @@ def plot_f_history(self,f_hist,t_hist):
 
         ax.plot(t_hist* 10 ** (3), f_hist * 10 ** (-3))
         ax.grid()
-        ax.set_xlabel(r'Time, [\mu s]', fontsize='xx-large')
+        ax.set_xlabel(r'Time, [$\mu s$]', fontsize='xx-large')
         ax.set_ylabel(r'Average Force [kN]', fontsize='xx-large')
         ax.set_xlim(xmin=0)
         ax.set_ylim(ymin=0)

tests/test_barber_p_dist.py

@@ -14,4 +14,4 @@ def test_barber_p_dist():
     barber = AnalBarber(bird_1,scenario, grid, eos)
 
     barber.plot_result(barber.P *10**(-6), 'Steady State Pressure', '[MPa]')
-    barber.plot_result(barber.c_p,'Steady State a Pressure Coefficient', '[-]')
+    barber.plot_result(barber.c_p,'Steady State a Pressure Coefficient', '[-]')

tests/test_bird_radii_impact_duration.py

@@ -39,8 +39,4 @@ def test_bird_radii_impact_duration():
     ax.set_ylabel(r'Impact Duaration, [$\mu$ s]',fontsize= 'xx-large')
     ax.set_xlim(xmin = 0)
     ax.set_ylim(ymin = 0)
-    plt.show()
-
-
-
-
+    plt.show()

tests/test_critical_aspect_ratio.py

@@ -45,4 +45,4 @@ def test_critical_aspect_ratio():
     ax.set_ylim(ymin = 0, ymax = critical_aspect[1])
     ax.set_xlim(xmin = 0)
     ax.set_ylim(ymin = 0)
-    plt.show()
+    plt.show()

tests/test_wilbeck_p_dist.py

@@ -28,5 +28,4 @@ def test_wilbeck_p_dist():
     ax.set_xlim(xmin = 0, xmax = 3.5)
     ax.set_ylim(ymin = 0)
     ax.legend()
-    plt.show()
-
+    plt.show()