fP_prop_pp = compute_fP_prop_pp(vinf_array)
fP_prop_ab = compute_fP_prop_ab(vinf_array)
Q_array = np.array([1.53E3, 1.88E3, 2.77E3, 4.09E3, 5.82E3,
7.96E3, 10.5E3, 13.5E3, 16.8E3, 20.7E3,
24.9E3])
tps_mf_array = 0.091*np.power(Q_array, 0.51575)/100.0
ess_mf = 0.23
fP_aero = 1.0 - ess_mf - tps_mf_array
fig = plt.figure()
fig.set_size_inches([6.25,6.25])
rcParams['font.family'] = 'sans-serif'
rcParams['font.sans-serif'] = ['DejaVu Sans']
params = {'mathtext.default': 'regular' }
plt.rcParams.update(params)
plt.plot(vinf_array/1.0E3, fP_prop_pp, linestyle='-', linewidth=1.00, marker='s',ms=12, markerfacecolor="None", markeredgecolor='red', markeredgewidth=0.75, color='red', label='Propulsive')
plt.plot(vinf_array/1.0E3, fP_prop_ab, linestyle='-', linewidth=1.00, marker='o',ms=12, markerfacecolor="None", markeredgecolor='magenta', markeredgewidth=0.75, color='magenta', label='Aerobraking')
plt.plot(vinf_array/1.0E3, fP_aero, linestyle='-', linewidth=1.00, marker='^',ms=12, markerfacecolor="None", markeredgecolor='blue', markeredgewidth=0.75, color='blue', label='Aerocapture')
plt.xlabel("Arrival "+r'$V_\infty$'+r', km/s' ,fontsize=16)
plt.ylabel('Payload mass fraction, '+r'$f_P$' ,fontsize=16)
ax = plt.gca()
ax.tick_params(direction='in')
ax.yaxis.set_ticks_position('both')
ax.xaxis.set_ticks_position('both')
plt.legend(loc='upper right', fontsize=16, frameon=False)
ax.set_ylim(0,1.03)
plt.xticks(fontsize=16)
plt.yticks(fontsize=16)
ax.set_ylim(0,1.03)
ax.annotate('',
xy=(vinf_array[6]/1.0E3, fP_prop_ab[6]) ,
xytext=(vinf_array[6]/1.0E3, fP_aero[6]),
arrowprops=dict(arrowstyle="<-"), va="center", ha="center", fontsize=16)
plt.text(vinf_array[6]/1.0E3, 0.5*(fP_prop_ab[6]+fP_aero[6]),
str(int((fP_aero[6] - fP_prop_ab[6])*100/fP_prop_ab[6]))+'%' ,
fontsize=16, ha='center', va='center',
bbox=dict(facecolor='white', edgecolor='white', pad=0.0))
ax.annotate('',
xy=(vinf_array[8]/1.0E3, fP_prop_ab[8]) ,
xytext=(vinf_array[8]/1.0E3, fP_aero[8]),
arrowprops=dict(arrowstyle="<-"), va="center", ha="center", fontsize=16)
plt.text(vinf_array[8]/1.0E3, 0.5*(fP_prop_ab[8]+fP_aero[8]),
str(int((fP_aero[8] - fP_prop_ab[8])*100/fP_prop_ab[8]))+'%' ,
fontsize=16, ha='center', va='center',
bbox=dict(facecolor='white', edgecolor='white', pad=0.0))
ax.annotate('',
xy=(vinf_array[10]/1.0E3, 0.0) ,
xytext=(vinf_array[10]/1.0E3, fP_aero[10]),
arrowprops=dict(arrowstyle="<-"), va="center", ha="center", fontsize=16)
plt.text(vinf_array[10]/1.0E3, 0.5*(0.0+fP_aero[10]),
r'$\infty$' ,
fontsize=16, ha='center', va='center',
bbox=dict(facecolor='white', edgecolor='white', pad=0.0))
plt.savefig('../data/jsr-paper/titan/titan-mb.png', dpi= 300,bbox_inches='tight')
plt.savefig('../data/jsr-paper/titan/titan-mb.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../data/jsr-paper/titan/titan-mb.eps', dpi=300,bbox_inches='tight')
plt.show()