Magnetopause location
The mp provides magnetopause location models
- Models:
Shue et al., (1998)
|
Calculates the magnetopause location following Shue et al. [1998]. |
Functions
- rbamlib.models.mp.S1998(theta, Bz, Pdyn)
Calculates the magnetopause location following Shue et al. [1998].
- Parameters:
theta (ndarray) – Angle \(\theta\) between the Earth-Sun line in radians (1D array).
Bz (ndarray) – IMF \(B_z\) in nT (1D array), must have the same size as \(P_{dyn}\).
Pdyn (ndarray) – Dynamic pressure \(P_{dyn}\) in nPa (1D array), must have the same size as \(B_z\).
- Returns:
r – Distance to the magnetopause from the Earth’s center in \(R_E\). The result will be a 2D array of shape (len(Pdyn), len(theta)).
- Return type:
ndarray
Notes
The magnetopause location is calculated using the empirical model proposed by Shue et al. (1998):
\[\begin{split}r_0(B_z) = \begin{cases} (11.4 + 0.013 \cdot B_z) \cdot P_{dyn}^{-1/6.6}, & \text{if} B_z \geq 0 \\ (11.4 + 0.14 \cdot B_z) \cdot P_{dyn}^{-1/6.6}, & \text{if} B_z < 0 \end{cases}\end{split}\]\[r = r_0 \cdot \left( \frac{2}{1 + \cos(\theta)} \right)^{\alpha}\]where:
\[\alpha = (0.58 - 0.010 \cdot B_z) \cdot (1 + 0.010 \cdot P_{dyn})\]