Calcolo del Giorno Giuliano

Unire la Data nel formato AAAA.MMGGdd
Esempio 21 Marzo 2012 alle ore 12:30 = 2012.032152
Se MM > 2 si prenda a = AAAA e m = MM
Se MM = 1 oppure 2 si prenda a = AAAA - 1 e m = MM + 12
Se il numero AAAA.MMGGdd è >= 1582.1015
Calcolare
A = Parte Intera di (a / 100)
B = 2 - A + Parte Intera di (A / 4)
Se invece AAAA.MMGGdd è < 1582.1015
allora A = B = 0

JD = Parte Intera di [365.25 * (a + 4716)] + Parte Intera di [30.6001 * (m + 1)] + GG.dd + B + 1524.5


Calcolo dell' Equazione del Tempo

Secoli Giuliani dall'epoca standard J2000 alla data inserita
T = (JDE - 2451545,0) ÷ 36525

Longitudine geometrica media del Sole
L = 280,46646° + 36000,76983° × (T) + 0,0003032° × (T)²

Anomalia media del Sole
M = 357,52911° + 35999,05029° × (T) - 0,0001537° × (T)²

Equazione del centro del Sole
C = [1,914602° - 0,004817° × (T) - 0,000014° × (T)²] × sen M° + [0,019993° - 0,000101° × (T)] × sen (2 × M°) + 0,000289° × sen (3 × M°)

Longitudine vera del Sole
LV = L + C

Longitudine apparente del Sole
LA = LV - 0,00569° - 0,00478° × sen (125,04° - 1934,136° × T)

Obliquità dell'ecclittica
ε = 23°26'21,448" - 0°00'4680,93" × (T ÷ 100) - 0°00'01,55" × (T ÷ 100)² + 0°00'1999,25" × (T ÷ 100)³ - 0°00'51,38" × (T ÷ 100)ˆ4 - 0°00'249,67" × (T ÷ 100)ˆ5 - 0°00'39,05" × (T ÷ 100)ˆ6 + 0°00'07,12" × (T ÷ 100)ˆ7 + 0°00'27,87" × (T ÷ 100)ˆ8 + 0°00'05,79" × (T ÷ 100)ˆ9 + 0°00'02,45" × (T ÷ 100)ˆ10

Eccentricità dell'orbita
Ec = 0,016708634 - 0,000042037 × (T) - 0,0000001267 × (T)²

Equazione del tempo
ET hms = {[tan (ε ÷ 2)]² × sen (2 × L) - 2 × Ec × sen M + (4 × Ec) × [tan (ε ÷ 2)]² × sen M × cos (2 × L) - (1 ÷ 2) × [tan (ε ÷ 2)]ˆ4 × sen (4 × L) - (5 ÷ 4) × (Ec)² × sen (2 × M)} × 180° ÷ 3,14159265359 ÷ 15