Python ppr Examples

Example #1

File: script_sugar_glider.py Project: TheSochaLab/Global-dynamics-of-non-equilibrium-gliding-in-animals

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)

# %%

Example #2

File: script_airfoil_snake.py Project: TheSochaLab/Global-dynamics-of-non-equilibrium-gliding-in-animals

# %% 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])