XL 2021 calcul de l'heure quand le soleil est au Zénith

Nicolas JACQUIN · 12 Août 2024

Bonjour à tous,

Avec chatGPT, j'ai fait une approche sur le calcul de l'heure quand le soleil est au Zénith, mais pas bon, mon format reçu n'est pas celui souhaité.

Le résultat est en cellule "B20"

Je ne sais pas si c'est juste une histoire de format ou de conversion, je devrais avoir 14h18m à une ou deux minutes prêt.

VB:

Function ZenithTime(ByVal Date_ As Date, ByVal Longitude As Double) As Double

            A = Year(Range("B1"))
            M = Month(Range("B1"))
            J = ActiveCell.Value
            vdate = DateSerial(A, M, J)

    ' Variables
    Dim n As Double ' Nombre de jours depuis le 1er janvier de l'année
    Dim B As Double ' Correction pour l'équation du temps
    Dim EoT As Double ' Équation du temps en minutes
    Dim SolarNoon As Double ' Midi solaire en heures

    ' Nombre de jours depuis le début de l'année
    n = DateSerial(Year(A), Month(M), Day(J)) - DateSerial(Year(A), 1, 1) + 1
  
    ' Calcul de l'équation du temps (EoT)
    B = (n - 81) * (360 / 365)
    EoT = 9.87 * Sin(2 * Application.WorksheetFunction.Radians(B)) - 7.53 * Cos(Application.WorksheetFunction.Radians(B)) - 1.5 * Sin(Application.WorksheetFunction.Radians(B))
  
    ' Calcul de l'heure du zénith solaire (en heures)
    SolarNoon = 12 + (4 * Longitude - EoT) / 60
  
    ' Retourner l'heure du zénith en heures décimales
    ZenithTime = SolarNoon
End Function

Sub testzenith()

Cells(20, 2) = ZenithTime(vdate, -3.36667)

End Sub

J'ai pris ma position sur Lorient.

Merci à tous
Nicolas

vgendron · 12 Août 2024

ca donnerait ceci:

VB:

'***********************************************
Function ZenithTime(ByVal vdate As Date, ByVal Longitude As Double) As Date
    ' Variables
    Dim n As Double ' Nombre de jours depuis le 1er janvier de l'année
    Dim B As Double ' Correction pour l'équation du temps
    Dim EoT As Double ' Équation du temps en minutes
    Dim SolarNoon As Double ' Midi solaire en heures
    Dim DateHe As Date 'date du passage à l'heure d'été
    Dim DateHh As Date 'date du passage à l'heure d'hiver
  
    DateHe = DateSerial(2024, 3, 31)
    DateHh = DateSerial(2024, 10, 28)
    If vdate >= DateHe And vdate < DateHh Then
        He = TimeSerial(2, 0, 0)
    Else
        He = TimeSerial(1, 0, 0)
    End If
  
    A = Year(vdate)...

vgendron · 12 Août 2024

Bonjour

déjà, dans la fonction Zenith, il y a une erreur pour calculer le nombre de jour depuis le début d'année

VB:

' Nombre de jours depuis le début de l'année
    n = vdate - DateSerial(A, 1, 1) + 1

ensuite, il faut peut etre mettre le format de ta cellule au format hh:mm:ss
mais la, ca donne un zenith le soir... comme quoi.. le soleil et la Bretagne...:-D
ou c'est juste que la formule n'est pas bonne non plus..?

Nicolas JACQUIN · 12 Août 2024

vgendron à dit:
comme quoi.. le soleil et la Bretagne...:-D

Bonjour, ça sais sûre. 🤣 🤣

Mais oui y a un truc mais 🙁🙁

Merci
Nico

Nicolas JACQUIN · 12 Août 2024

J'ai fait une correction, mais pas mieux

VB:

'***********************************************
Function ZenithTime(ByVal Date_ As Date, ByVal Longitude As Double) As Double

            A = Year(Range("B1"))
            M = Month(Range("B1"))
            J = ActiveCell.Value
            vdate = DateSerial(A, M, J)

    ' Variables
    Dim n As Double ' Nombre de jours depuis le 1er janvier de l'année
    Dim B As Double ' Correction pour l'équation du temps
    Dim EoT As Double ' Équation du temps en minutes
    Dim SolarNoon As Double ' Midi solaire en heures

    ' Nombre de jours depuis le début de l'année
    n = vdate - DateSerial(A, 1, 1) + 1
   
    ' Calcul de l'équation du temps (EoT)
    B = (n - 81) * (360 / 365)
    EoT = 9.87 * Sin(2 * Application.WorksheetFunction.Radians(B)) - 7.53 * Cos(Application.WorksheetFunction.Radians(B)) - 1.5 * Sin(Application.WorksheetFunction.Radians(B))
   
    ' Calcul de l'heure du zénith solaire (en heures)
    SolarNoon = 12 + (4 * Longitude - EoT) / 60
    Debug.Print SolarNoon
   
    ' Retourner l'heure du zénith en heures décimales
    ZenithTime = SolarNoon
End Function

Sub testzenith()

Cells(20, 2) = ZenithTime(vdate, -3.36667) 'Longitude Lorient = -3.36667

End Sub

vgendron · 12 Août 2024

Regarde ici

VB:

Function ZenithTime(ByVal Date_ As Date, ByVal Longitude As Double) As Date
            Pi = WorksheetFunction.Pi()
            A = Year(Range("B1"))
            M = Month(Range("B1"))
            J = ActiveCell.Value
            vdate = DateSerial(A, M, J)

    ' Variables
    Dim n As Double ' Nombre de jours depuis le 1er janvier de l'année
    Dim B As Double ' Correction pour l'équation du temps
    Dim EoT As Double ' Équation du temps en minutes
    Dim SolarNoon As Double ' Midi solaire en heures

    ' Nombre de jours depuis le début de l'année
    n = vdate - DateSerial(A, 1, 1) + 1
    
    ' Calcul de l'équation du temps (EoT)
    'Source du calcul: https://www.techno-science.net/definition/6567.html
    
    B = (2 * Pi * (n - 81) / 364)
    EoT = -9.87 * Sin(2 * B) + 7.53 * Cos(B) + 1.5 * Sin(B)
    
    ' Calcul de l'heure du zénith solaire (en heures) 'c'est ici que ca doit coincer
    SolarNoon = 12 + (4 * Longitude + EoT) * 60 '/60 ?
    
    ' Retourner l'heure du zénith en heures décimales
    ZenithTime = Format(SolarNoon, "h:mm:ss")
End Function

Sub testzenith()
    Cells(20, 2) = ZenithTime(vdate, -3.36667)
End Sub

la fonction EoT me semble correcte
maintenant, c'est le SolarNoon, et les unités à vérifier

Nicolas JACQUIN · 12 Août 2024

Bonjour Jacky, J'ai entrée dans le code la longitude de Lorient, se qui devrait me donner 14h18m environ, après je suis pas à 5 minutes prêt, le mieux c'est d'être au plus proche 😉

Nicolas JACQUIN · 12 Août 2024

re vgendron,

moi j'ai pris comme base horaire sur ce site,

Lever du Soleil et Coucher du Soleil | Lorient

Lorient, France - Un tableau et des graphiques du lever et du coucher du soleil aujourd'hui et pour toute l'année 2025.

meteogram.fr

on est déjà plus proche mais quand même pas ça.

edit, je suis à 16h34m au lieu de 14h18m

Merci

TooFatBoy · 12 Août 2024

Salut Nico 😉

J'ai supprimé mon message car le camarade vgendron semble avoir modifié quelque chose, donc je suppose que ça devrait mieux fonctionner. 👍

vgendron · 12 Août 2024

en calculant l'EOT comme suit, ??

VB:

SolarNoon = 12 + 4 * Longitude + EoT

je n'ai pas réussi à trouver de formule claire pour le solarNoon.. en as tu une? (autre que chatgpt)

vgendron · 12 Août 2024

TooFatBoy à dit:
Salut Nico 😉

J'ai supprimé mon message car le camarade vgendron semble avoir modifié quelque chose, donc je suppose que ça devrait mieux fonctionner. 👍

pas sur :-D

Nicolas JACQUIN · 12 Août 2024

Non, rien de plus pour aider

vgendron · 12 Août 2024

alors, voici ce que j'ai trouvé
à mon avis, il y a bien un pb d'unité
la correction donnée par l'EoT est en minutes
tout comme la correction due à la longitude

VB:

'***********************************************
Function ZenithTime(ByVal Date_ As Date, ByVal Longitude As Double) As Date
            Pi = WorksheetFunction.Pi()
            A = Year(Range("B1"))
            M = Month(Range("B1"))
            J = ActiveCell.Value
            vdate = DateSerial(A, M, J)

    ' Variables
    Dim n As Double ' Nombre de jours depuis le 1er janvier de l'année
    Dim B As Double ' Correction pour l'équation du temps
    Dim EoT As Double ' Équation du temps en minutes
    Dim SolarNoon As Double ' Midi solaire en heures

    ' Nombre de jours depuis le début de l'année
    n = vdate - DateSerial(A, 1, 1) + 1
    
    ' Calcul de l'équation du temps (EoT)
    'Source du calcul: https://www.techno-science.net/definition/6567.html
    
    B = (2 * Pi * (n - 81) / 364) 'pi est en radian ==> B est déjà en radian
    EoT = -9.87 * Sin(2 * B) + 7.53 * Cos(B) + 1.5 * Sin(B)
    
    ' Calcul de l'heure du zénith solaire (en heures) 'c'est ici que ca doit coincer
    SolarNoon = 12 + 4 * Longitude + EoT ' * 60 + EoT '* 60 '/60 ?
    
    'https://webetab.ac-bordeaux.fr/Pedagogie/Physique/TPE/midi.htm
    'HEURE LEGALE = HEURE SOLAIRE + CORRECTION LONGITUDE + 1 h (ou 2 l'été) + CORRECTION "EQUATION DU TEMPS"
    HS = TimeSerial(12, 0, 0)
    CL = TimeSerial(0, 4 * Longitude, 0)
    He = TimeSerial(2, 0, 0)
    
    HL = HS + CL + He + TimeSerial(0, EoT, 0)
    
    ' Retourner l'heure du zénith en heures décimales
    'ZenithTime = Format(SolarNoon, "h:mm:ss")
    ZenithTime = Format(HL, "hh:mm:ss")
End Function

à vérifier et confirmer sur plusieurs dates, et plusieurs positions..

Nicolas JACQUIN · 12 Août 2024

re, merci du retour et de l'intérêt,
j'ai une bonne 1/2 heure de décalage
avec le code comme ça:

VB:

Function ZenithTime(ByVal Date_ As Date, ByVal longitude As Double) As Date

            Pi = WorksheetFunction.Pi()
            A = Year(Range("B1"))
            M = Month(Range("B1"))
            J = ActiveCell.Value
            vdate = DateSerial(A, M, J)
            
longitude = -3.36667 'Lorient

    ' Variables
    Dim n As Double ' Nombre de jours depuis le 1er janvier de l'année
    Dim B As Double ' Correction pour l'équation du temps
    Dim EoT As Double ' Équation du temps en minutes
    Dim SolarNoon As Double ' Midi solaire en heures

    ' Nombre de jours depuis le début de l'année
    n = vdate - DateSerial(A, 1, 1) + 1
    
    ' Calcul de l'équation du temps (EoT)
    'Source du calcul: https://www.techno-science.net/definition/6567.html
    
    B = (2 * Pi * (n - 81) / 364) 'pi est en radian ==> B est déjà en radian
    EoT = -9.87 * Sin(2 * B) + 7.53 * Cos(B) + 1.5 * Sin(B)
    
    ' Calcul de l'heure du zénith solaire (en heures) 'c'est ici que ca doit coincer
    SolarNoon = 12 + 4 * longitude + EoT ' * 60 + EoT '* 60 '/60 ?
    
    'https://webetab.ac-bordeaux.fr/Pedagogie/Physique/TPE/midi.htm
    'HEURE LEGALE = HEURE SOLAIRE + CORRECTION LONGITUDE + 1 h (ou 2 l'été) + CORRECTION "EQUATION DU TEMPS"
    HS = TimeSerial(12, 0, 0)
    CL = TimeSerial(0, 4 * longitude, 0)
    He = TimeSerial(2, 0, 0)
    
    HL = HS + CL + He + TimeSerial(0, EoT, 0)
    
    ' Retourner l'heure du zénith en heures décimales
    'ZenithTime = Format(SolarNoon, "h:mm:ss")
    ZenithTime = Format(HL, "hh:mm:ss")
End Function

Sub testzenith()

longitude = -3.36667

    Cells(20, 2) = ZenithTime(vdate, longitude)
    
End Sub

Et ce sur plusieurs dates.

vgendron · 12 Août 2024

bon, on doit pas utiliser la meme chose

1) le code de la fonction purgée des commentaires ET des calculs de vdate qui sont inutiles puisque la fonction est appelée avec une date...

VB:

Function ZenithTime(ByVal vdate As Date, ByVal Longitude As Double) As Date
    ' Variables
    Dim n As Double ' Nombre de jours depuis le 1er janvier de l'année
    Dim B As Double ' Correction pour l'équation du temps
    Dim EoT As Double ' Équation du temps en minutes
    Dim SolarNoon As Double ' Midi solaire en heures

    A = Year(vdate)
    Pi = WorksheetFunction.Pi()
    ' Nombre de jours depuis le début de l'année
    n = vdate - DateSerial(A, 1, 1) + 1
   
    ' Calcul de l'équation du temps (EoT)
    'Source du calcul: https://www.techno-science.net/definition/6567.html
   
    B = (2 * Pi * (n - 81) / 364) 'pi est en radian ==> B est déjà en radian
    EoT = -9.87 * Sin(2 * B) + 7.53 * Cos(B) + 1.5 * Sin(B)
   
    ' Calcul de l'heure du zénith solaire (en heures) 'c'est ici que ca doit coincer
    SolarNoon = 12 + 4 * Longitude + EoT  ' * 60 + EoT '* 60 '/60 ?
   
    'https://webetab.ac-bordeaux.fr/Pedagogie/Physique/TPE/midi.htm
    'HEURE LEGALE = HEURE SOLAIRE + CORRECTION LONGITUDE + 1 h (ou 2 l'été) + CORRECTION "EQUATION DU TEMPS"
    HS = TimeSerial(12, 0, 0)
    CL = TimeSerial(0, 4 * Longitude, 0)
    He = TimeSerial(2, 0, 0)
   
    HL = HS + CL + He + TimeSerial(0, EoT, 0)
   
    ' Retourner l'heure du zénith en heures
    ZenithTime = Format(HL, "hh:mm:ss")
End Function

Sub testzenith()
    MsgBox ZenithTime(DateSerial(2024, 8, 31), 3.36667)
End Sub

avec ca, j'ai le meme résultat que sur ton site: https://meteogram.fr/soleil/france/lorient/
à UNE minute près

piga25 · 12 Août 2024

Bonjour,
Avec le calcul du jour Julien pour le 12 août 2024 j'ai 13:55

VB:

Function ZenithSolaire(Latitude As Double, Longitude As Double, DateChoisie As Date) As String
    Dim JourJulien As Double
    Dim T As Double
    Dim Déclinaison As Double
    Dim EquationTemps As Double
    Dim HeureLocale As Double
  
    ' Calcul du jour Julien
    JourJulien = DateChoisie + 2415018.5 + (TimeValue("12:00:00") / 24)
  
    ' Calcul du temps en siècles depuis J2000.0
    T = (JourJulien - 2451545#) / 36525
  
    ' Calcul de la déclinaison solaire (en radians)
    Déclinaison = 0.006918 - 0.399912 * Cos(T) + 0.070257 * Sin(T) - 0.006758 * Cos(2 * T) + 0.000907 * Sin(2 * T) - 0.002697 * Cos(3 * T) + 0.00148 * Sin(3 * T)
  
    ' Calcul de l'équation du temps (en minutes)
    EquationTemps = 229.18 * (0.000075 + 0.001868 * Cos(T) - 0.032077 * Sin(T) - 0.014615 * Cos(2 * T) - 0.040849 * Sin(2 * T))
  
    ' Calcul de l'heure locale du zénith solaire
    HeureLocale = 12 + Range("B3").Value + (Longitude / 15) - (EquationTemps / 60)
  
    ' Conversion de l'heure locale en heures et minutes
    Dim Heure As Integer
    Dim Minute As Integer
    Heure = Int(HeureLocale)
    Minute = Round((HeureLocale - Heure) * 60)
  
    ' Formatage de l'heure du zénith solaire en chaîne de caractères
    ZenithSolaire = Format(Heure, "00") & ":" & Format(Minute, "00")
End Function

Sub CalculerZenith()
    Dim Latitude As Double
    Dim Longitude As Double
    Dim DateChoisie As Date
    Dim HeureZenith As String
  
    ' Exemple de latitude et de longitude (Lorient, France)
    Latitude = 47.482524
    Longitude = -3.3702449
  
    ' Date choisie
    DateChoisie = Range("A2").Value 'DateValue("12/08/2024")
  
    ' Appel de la fonction
    HeureZenith = ZenithSolaire(Latitude, Longitude, DateChoisie)
  
    ' Affichage de l'heure du zénith
    Range("D2") = HeureZenith
    MsgBox "L'heure du zénith solaire le " & DateChoisie & " est à " & HeureZenith
End Sub

XL 2021 calcul de l'heure quand le soleil est au Zénith

Nicolas JACQUIN

XLDnaute Accro

Pièces jointes

vgendron

vgendron

XLDnaute Barbatruc

Nicolas JACQUIN

XLDnaute Accro

Nicolas JACQUIN

XLDnaute Accro

vgendron

XLDnaute Barbatruc

Nicolas JACQUIN

XLDnaute Accro

Nicolas JACQUIN

XLDnaute Accro

Lever du Soleil et Coucher du Soleil | Lorient

TooFatBoy

XLDnaute Barbatruc

vgendron

XLDnaute Barbatruc

vgendron

XLDnaute Barbatruc

Nicolas JACQUIN

XLDnaute Accro

vgendron

XLDnaute Barbatruc

Nicolas JACQUIN

XLDnaute Accro

vgendron

XLDnaute Barbatruc

piga25

XLDnaute Barbatruc

Pièces jointes

Discussions similaires

XL 2021 calcul de l'heure quand le soleil est au Zénith

XLDnaute Accro

Pièces jointes

XLDnaute Barbatruc

XLDnaute Accro

XLDnaute Accro

XLDnaute Barbatruc

XLDnaute Accro

XLDnaute Accro

XLDnaute Barbatruc

XLDnaute Barbatruc

XLDnaute Barbatruc

XLDnaute Accro

XLDnaute Barbatruc

XLDnaute Accro

XLDnaute Barbatruc

XLDnaute Barbatruc

Pièces jointes

Discussions similaires

Privacy & Transparency

Privacy & Transparency