Demande J.A. Sauvaud
--------------------
== Objet
Calcul du parametre L pour Mercure en utilisant le modèle de Korth.
== Composant
1. Récupérer le modèle de Korth pour Mercure
PC>Mes documents>Bepi>Sauvaud>grl*.pdf
2. Récupérer algo Runge-Kutta pour intégrer les lignes de force du champ magnétique
3. Récupérer orbite Bepi (dans repère à déterminer)
4. 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))
