Bibliothèque pour suivi et éphémérides du soleil

Bonjour à toutes et à tous,

Comme le titre le dit, j'ai écrit une petite bibliothèque permettant de calculer la position du soleil dans le ciel pour une date et une heure donnée à une latitude et une longitude donnée en prenant en compte sa zone horaire.

Elle permet aussi de calculer les heures du lever, du zénith et du coucher de cet astre.

Cordialement.

Pierre

hello bonne idée le lien vers ton site fonctionne mais sur ton site, pas possible de charger ta bibliothèque

Désolé, normalement, c'est réparé. Apparemment, j'avais utilisé un caractère non autorisé dans le nom du fichier.

Cordialement.

Pierre

Dans la version V2, j'ai ajouté la prise en compte de l'altitude pour les lever et coucher du soleil.

Voir dans l'exemple joint cette influence entre le pied et le haut de la Tour Eiffel.

Cordialement.

Pierre

ChPr: Dans la version V2, j'ai ajouté la prise en compte de l'altitude pour les lever et coucher du soleil.

Voir dans l'exemple joint cette influence entre le pied et le haut de la Tour Eiffel.

Cordialement.

Pierre

bonjour sympa reste plus qu'a integrer de la recup de data GNSS (Lat/lon/Elev), une gestion d'affichage et tu mets au point un afficheur jour/nuit aeronautique concernant le lieu où tu es .... :grin:

Artouste: ... un afficheur jour/nuit aeronautique ...

Mis à part avoir fait du parapente avec pour seul instrument un variomètre, je n'y connais pas grand chose en instrumentation aéronautique. Donc, quelles sont les informations que donne un afficheur jour/nuit aéronautique ?

Cordialement.

Pierre

ChPr: Mis à part avoir fait du parapente avec pour seul instrument un variomètre, je n'y connais pas grand chose en instrumentation aéronautique. Donc, quelles sont les informations que donne un afficheur jour/nuit aéronautique ?

Cordialement.

Pierre

Bonsoir Pierre, Nickel comme bibliothèque. Pour finir, il ne reste plus qu'à introduire la "Réfraction atmosphérique" ;) @+

Mais si, elle y est à la ligne 93 de "Ephem_Soleil.cpp". Elle prend la valeur de -0.833°

case 0: // crépuscule sur l'horizon
      acr = cos(PI_180 * (90.833 + 0.015 * sqrt(alt)));

valeur à laquelle j'ai ajouté la prise en compte de l'altitude "alt" selon le document que tu m'avais mis en lien.

Cordialement.

Pierre

Un chose m’interpelle, cet angle : 0.833 correspond à la somme de la réfraction atmosphérique 34' (soit 0.566°) utilisée dans le document cité avec le demi-diamètre du soleil : 16' (0.266°) : 0.833 = 0.566 + 0.266

Mais jouant avec ces différentes valeurs, je n'arrive toujours pas à avoir de concordance entre mes cinq sources équivalentes et celle de IMCCE.

Cordialement.

Pierre

Bonjour,

ChPr: Mais si, elle y est à la ligne 93 de "Ephem_Soleil.cpp". Elle prend la valeur de -0.833°

Les 34' correspondent au décalage entre le centre et le bord supérieur du soleil. La réfraction est un phénomène physique qui se rajoute pour tenir compte de la position vraie par rapport à la position visible. En d'autre terme, on peut voir le soleil alors qu'il est déjà entièrement sous l'horizon. Ce phénomène n'est pas simple à mettre en œuvre car il est fonction de la nature des couches de l'atmosphère, de la température, de la pression (et peut être de l'âge du capitaine ;)) mais on peut arriver avec des valeurs appropriées à des résultats acceptables et proche de la réalité. @+

icare: Bonjour,Les 34' correspondent au décalage entre le centre et le bord supérieur du soleil. ...

Ben non, je cite ton document :

R est la réfraction à l'horizon. Les tables publiées dans l'annuaire du Bureau des longitudes utilisent la théorie de la réfraction de Radau qui conduit à R = 36' 36" mais l'on pourra utiliser la valeur R = 34' adoptée dans les Ephémérides Nautiques publiées par le Bureau des Longitudes et dans d'autres publications étrangères. 1/2 d est le demi-diamètre apparent de l'astre. On l'introduit dans la formule quand on calcule le lever et le coucher du bord supérieur du Soleil et de la Lune et non pas le lever et le coucher du centre de l'astre. On prend, aussi bien pour le Soleil que pour la Lune, 1/2 d = 16'.

Cordialement.

Pierre

Re,
J’avais fait un programme qui traitait tous les aspects de la lune et du soleil sous Delphi au siècle précédent.
Le programme était fait avec les algorithmes de Jean Meus et les contrôles avec le site IMCCE.
Il faut que je retrouve les sources.
@+

icare: ... sous Delphi au siècle précédent. ...

Toi aussi tu as fait du Delphi à outrance ;)

Cordialement.

Pierre

Ce que je pensais, a priori :

ChPr: Un chose m’interpelle, cet angle : 0.833 correspond à la somme de la réfraction atmosphérique 34' (soit 0.566°) utilisée dans le document cité avec le demi-diamètre du soleil : 16' (0.266°) : 0.833 = 0.566 + 0.266 ...

est expliqué clairement dans ce document Wikipédia :

"Pour le calcul des heures de lever et de coucher du Soleil la réfraction atmosphérique est également prise en compte. L'heure calculée ne correspond pas au moment où le Soleil atteint l'altitude 0°, mais au moment où il atteint l'altitude -50' : elle-même étant la somme de deux données: -16' pour tenir compte du rayon angulaire du Soleil, cette valeur correspondant au demi-diamètre apparent de l'astre exprimé en minutes (les tables astronomiques donnant généralement sa position héliocentrique) et -34' pour la réfraction atmosphérique à l'horizon."

Cordialement.

Pierre

Bonjour

Je n'ai malheureusement pas la dispo pour me replonger dans mes travaux de l'époque. Par contre je certifie plusieurs points :

  • Le site de l'IMCCE a évolué, ce qui est bien dommage car l'ancienne version permettait de récupérer facilement, pour un point donné (et altitude 0m), les éphémérides sur une longue période. Avec un résultat importable dans excel.
  • Mon algo est basé sur le franchissement du milieu du disque solaire à l'horizon (et non le bord supérieur), incluant la réfraction atmosphérique
  • Sur des milliers de tests comparatifs (plusieurs points du globe sur une année entière), et pour une latitude inférieure à 60° (Nord ou Sud), l'écart dans les calculs d'horaires lever/coucher de soleil était d'au maximum 3 secondes par rapport au site IMCCE de l'époque. Tout ceci chargé et vérifié dans excel

@ Bricoleau :

J'ai repris vos exemples et les ai à nouveau comparés à ceux trouvés par les quatre sources équivalentes de ce post en réponse #2. Selon ce qui est dit dans la réponse précédente de ce post, ces quatre sources se basent sur le haut du soleil et; vous, vous vous basez sur le centre.

Dans mon application (comme je suis maître des paramètres ;) ), pour voir, j'ai enlevé le demi-diamètre du soleil (16') et je tombe sur vos résultats.

Tous ces résultats sont donc concordants, encore faut-il préciser, comme vous le dites, sur quel point de l'astre l'on se base.

Cordialement.

Pierre

Pour ma part, j’ai utilisé le code de Bricoleau et l’ai intégré à un ensemble plus important.

En pièce jointe, l’écran relatif aux heures de lever, zénith, coucher et durée pour la journée en cours (infos reçues par DCF77).

Soleil.jpg