calculateur de marée pour HP48

Calculateur de Marée
Sur le graphe suivant, on représente en bleu le calcul simplifié avec le mode de la lune uniquement pour du1/1/2003 au 15/1/2003. En noir, il s'agit du calcul issu de l'ouvrage "Table des marées des grands ports du monde", éditions SHOM 1984 avec 20 modes (calculé avec Mathematica).


Les formules (issues de l'ouvrage "Table des marées des grands ports du monde", éditions SHOM 1984) sont codées pour la HP48 en ascii et en binaire (émulateur de HP48 "x48", "emu48" etc.).

Le directory est "SHOM"
remarquer l'élégance de la programmation: la phase est calculée par un produit scalaire "DOT":
par exemple, pour l'onde annuelle
ASA GSA [ 0 0 1 0 0 0 0 ] pa CTM EVAL
ASA GSA sont les coefficients, pa le vecteur des positions astronomiques et où CTM est le programme:
<< DOT - COS * 1000 / + DUP 4 DISP >>
fort non?
Les formules (issues de l'ouvrage "Table des marées des grands ports du monde", éditions SHOM 1984) sont codées (jv2003) sous Mathematica