lims = (vxlim, vzlim) = np.r_[0, 1.25], np.r_[0, -1.25]
tvec = np.linspace(0, 30, 351)
reload(plots)
from plots import phase_plotter as ppr
fig, ax = ppr(afdict,
pitch,
lims,
angle_rng,
tvec,
ngrid=201,
nseed=501,
nseed_skip=25,
quiver=False,
skip=25,
seed=False,
timer=True,
gamtest=gammas,
extrap=extrap,
traj=plots.ps_traj_dp5,
fig=None,
ax=None)
lab = 'sugar glider, ' + r'$\theta=$0' + u'\u00B0'
ax.text(.05, -1, lab, fontsize=16)
fig.savefig('Figures/figure5di_vpd0_sugar_glider.pdf', transparent=True)
# %%
# %% Velocity polar diagram, pitch = 0
afdict = dict(cli=cl_fun, cdi=cd_fun, clip=clprime_fun, cdip=cdprime_fun)
pitch = 0
arng = sn_angle_rng
extrap = (ale[0], ale[-1])
lims = (vxlim, vzlim) = np.r_[0, 1.25], np.r_[0, -1.25]
tvec = np.linspace(0, 30, 351)
reload(plots)
from plots import phase_plotter as ppr
fig, ax = ppr(afdict, pitch, lims, arng, tvec, ngrid=201,
nseed=501, nseed_skip=25, quiver=False, skip=10, seed=False,
timer=True, gamtest=gammas, extrap=extrap,
traj=plots.ps_traj_dp5, fig=None, ax=None)
lab = 'airfoil snake, ' + r'$\theta=$0' + u'\u00B0'
ax.text(.05, -1, lab, fontsize=16)
fig.savefig('Figures/figure5bi_vpd0_airfoil_snake.pdf', transparent=True)
# %% Velocity polar diagram, pitch = 0 with Z nullcline
afdict = dict(cli=cl_fun, cdi=cd_fun, clip=clprime_fun, cdip=cdprime_fun)
pitch = 0
arng = sn_angle_rng
extrap = (ale[0], ale[-1])