irfpy.venus.bowshock

Collection of bow shock model

Reference:

Martinecz et al. (2008): “Location of the bow shock and ion composition boundaries at Venus–initial determinations from Venus Express ASPERA-4” PSS 56, 780-784.

Code author: Moa Persson

irfpy.venus.bowshock.xr_martinecz08(scale=1.0)

Return the (x,r) coordinates of Martinecz 08 model.

Returns:

A tuple, (x, r), where x is the array of x coordinates and r is the array of $\sqrt{y^2+z^2}$. The unit is Rv.

The parameter is $L=1.303,ϵ=1.056,x_0=0.788$. Then, $r=L/(1+ϵ\cos \theta )$, where $(r,\theta )$ are measured at $(x_0,0,0)$.

>>> from irfpy.venus import bowshock
>>> x, r = bowshock.xr_martinecz08()
>>> print(x.shape, r.shape)
(161,) (161,)

You can plot the boundary as

import matplotlib.pyplot as plt
plt.plot(x, r, 'r')

irfpy.venus.bowshock.inside_martinecz08(x, y, z, scale=1.0)

Return True if the given position is inside Martinecz shock model.

Parameters:

z (x, y,) – Position. The unit should be Rv.

Returns:

True if the given position is inside the Martinecz model.

>>> inside_martinecz08(1, 0, 0)
np.True_
>>> inside_martinecz08(2, 0, 0)
np.False_
>>> inside_martinecz08(1, 1, 0)
np.True_
>>> inside_martinecz08(1, 1, 1)
np.False_