Demande J.A. Sauvaud
Objet
Calcul du parametre L pour Mercure en utilisant le modèle de Korth.
Composant
-
Récupérer le modèle de Korth pour Mercure
PC>Mes documents>Bepi>Sauvaud>grl*.pdf
-
Récupérer algo Runge-Kutta pour intégrer les lignes de force du champ magnétique
-
Récupérer orbite Bepi (dans repère à déterminer)
-
Construire ébauche algorithme en Fortran
Méthodologie
Pour chaque position de MIA (repère à determiner) :
-
Calculer B à partir du modèle de Korth (champ interne + externe)
-
Utilier Runge Kutta pour intégrer jusqu'à la surface de Mercure
-
Boucler avec sortie si R ⇐ 1
-
Prévoir de sortir en erreur si R > 12
-
Utiliser Runge Kutta avec le modèle interne uniquement pour caluler L
-
Sortie quand on trouve Bmin (dès que la dérivée de B change de sens)
⇒ on a alors L = R
-
On peut alors calculer latiture invariante = arcos (sqrt (1/L))
