irfpy.mars.bowshock

Bow shock model by Vignes et al.

• BowshockVignes provides default bow shock model

• BowshockVignesScaled provides a bow shock model of your own scaling

class irfpy.mars.bowshock.BowshockVignes

Bases: object

Vignes et al, (TBC) model for bow shock position.

>>> bsx, bsr = BowshockVignes.xr()

Conic curve (hyperbola) model with a formulation of

$$\begin{array}{rl}x& =R\cos \theta +x_0\\ y& =R\sin \theta \end{array}$$

with

$$R=\frac{L}{1+ϵ\cos \theta }$$

The parameters are

$$\begin{array}{rl}L& =2.04\\ ϵ& =1.03\\ x_0& =0.64\end{array}$$

in the unit of Martian radii.

L = 2.04

e = 1.03

x0 = 0.64

classmethod xr()

Return the bow shock position.

Returns:

The array with shape of (2, 180). Unit is in Rm.

>>> print(BowshockVignes.xr()[:, 50])
[1.49948598 0.77115753]

classmethod inside(x, y, z)

Return True if s/c is in bow shock.

>>> print(BowshockVignes.inside(1.2, 0, 0))
True
>>> print(BowshockVignes.inside(3, 0, 0))
False

irfpy.mars.bowshock.BowshockVignesScaled(scale=1)

Vignes et al., while scaled by changing the parameter L.

Parameters:

scale – Default 1.

Returns:

A Vignes model with scaled by the given scale.

10% bigger model is obtained as follows

>>> vig10 = BowshockVignesScaled(1.1)    # 10% bigger model
>>> x10, r10 = vig10.xr()
>>> plt.plot(x10, r10)      
>>> vig10.inside(0, 2.7, 0)
np.True_

>>> x00, r00 = BowshockVignes.xr()
>>> plt.plot(x00, r00)      
>>> BowshockVignes.inside(0, 2.7, 0)
np.False_