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 del Tempo Sidereo apparente

Calcolare il JD corrispondente alla data alle ore 0
T = (JD - 2415020.0) / 36525
θº = 0.276919398 + 100.0021359 * T + 0.000001075 * T²

Ts medio = Parte frazionaria di θº * 1.002737908 + orario UT (in ore decimali)

Si effetuano poi i calcoli per la nutazione e l'obliquità dell'ecclittica

Calcolare il JD corrispondente alla data ed alle ore UT

T = (JD - 2415020,0) / 36525

L = 279,6967° + 36000,7689° * T + 0,000303° * T²
L1 = 270,4342° + 481267,8831° * T - 0,001133° * T²
M = 358,4758° + 35999,0498° * T - 0,000150° * T²
M1 = 296,1046° + 477198,8491° * T + 0,009192° * T²
Ω = 259,1833° - 1934,1420° * T + 0,002078° * T²
Δψ = - (17,2327” + 0,01737” * T) * sen Ω - (1,2729” + 0,00013” * T) * sen (2L) + 0,2088” * sen (2Ω) - 0,2037” * sen (2L1) + (0,1261” - 0,00031” * T) * sen M + 0,0675” * sen M1 - (0,0497” - 0,00012” * T) * sen (2L + M) - 0,0342” * sen (2L1 - Ω) - 0,0261” * sen (2L1 + M1) + 0,0214” * sen (2L - M) - 0,0149” * sen (2L - 2L1 + M1) + 0,0124” * sen (2L - Ω) + 0,0114” * sen (2L1 - M1)
U = T / 100
ε = 23°26’21,448” - 0°00’4680,93” * (T / 100) - 0°00’01,55” * (T / 100)2 + 0°00’1999,25” * (T / 100)3 - 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

Correzione = Δψ * cos ε / 15

Ts apparente = Ts medio + Correzione