// JavaScript Document
function jh(h){
if (h >= 24){h -= 24}
else {h = h}
return h;}
function abs_floor( x )
{
var ro;
if (x >= 0.0) ro = Math.floor(x);
else ro = Math.ceil(x);
return ro;
}
function mod2pi( x )
{
var bo = x / (2 * Math.PI);
var ao = (2 * Math.PI) * (bo - abs_floor(bo));
if (ao < 0) ao = (2 * Math.PI) + ao;
return ao;
}
function notn(x){
if (isNaN(x))x = " -- "
return x;
}
function aff(x){
if (x > 15)x = " -- ";
else x = x + "°";
return x;
}
function tronque(x) {
if (x>0.0) return(Math.floor(x));
else return Math.ceil(x);
}
//------------------------------------------------------------ DEFILEMENT pour des données en temps réel
//--------------- temps en millecondes pour la fonction "marche"
var time = 10000;
function marche()
{
//theTime = window.setTimeout("marche()", time);
date();
calculer();
}
function arret() {
//clearTimeout(theTime)
}
//------------------------------------------------------------ DATE
function date() {
var today = new Date();
jour = today.getUTCDate();
mois = today.getUTCMonth() +1;
annee = today.getUTCFullYear();
heure = today.getUTCHours();
minute = today.getUTCMinutes();
seconde = today.getUTCSeconds();
datas['day'] = jour;
datas['month'] = mois;
datas['year'] = annee;
datas['hour'] = heure;
datas['minut'] = minute;
datas['second'] = seconde;
} // date()
function joursem(jj){
var joursemain0 = (Math.round(jj)%7);
if (joursemain0 == 0)
{joursemain = "lundi"}
if (joursemain0 == 1)
{joursemain = "mardi"}
if (joursemain0 == 2)
{joursemain = "mercredi"}
if (joursemain0 == 3)
{joursemain = "jeudi"}
if (joursemain0 == 4)
{joursemain = "vendredi"}
if (joursemain0 == 5)
{joursemain = "samedi"}
if (joursemain0 == 6)
{joursemain = "dimanche"}
datas['joursem'] = joursemain;
}
//------------------------------------------------------------ CALCUL
function calculer()
{
var ct, T, T2, T3, R, A, B, C, D, F, E, G, H, M, N, I, J, W, X, Y, Z, a, b, c, d, e, f, l, om, n, m, l1, jj, jj2, r2d, y, m;
var jour = parseInt(datas['day']);
var mois = parseInt(datas['month']);
var annee = parseInt(datas['year']);
var heure = parseInt(datas['hour']);
var minute = parseInt(datas['minut']);
var seconde = parseInt(datas['second']);
if (mois>2) {
y = annee;
m = mois;
}
else {
y = annee - 1;
m = mois + 12;
}
with (Math)
{
a = floor(y / 100);
b = 2 - a + floor(a / 4);
//--------------- Jour Julien
var heure1 = heure + (minute / 60) +(seconde / 3600)
var ggg = 1
if( annee < 1582 ) ggg = 0;
if( annee < 1582 && mois < 10 ) ggg = 0;
if( annee <= 1582 && mois == 10 && jour < 5 ) ggg = 0;
var jd = - 1 * tronque(7 * (tronque((mois + 9) / 12) + annee) / 4);
var s = 1
if ((mois - 9) < 0) s = - 1;
a = abs(mois - 9);
var j1 = tronque(y + s * tronque(a / 7));
var j2 = - 1 * tronque(((tronque(j1 / 100) + 1) * 3 / 4));
var jd1 = jd + tronque((275 * mois / 9) + jour + (ggg * j2));
var jd2dat1 = jd1 + 1721027 + 2 * ggg + 367 * annee - 0.5 ;
var jj = jd2dat1 + (heure1 / 24);
datas['jourjulien'] = jj;
joursem(jj);
//--------------- Temps mesuré en siecle juliens de 36525 jours des éphémérides
// a partir du 0,5 janvier 1900 TE
T = (jj - 2415020.0) / 36525;
T2 = T * T;
T3 = T2 * T;
//--------------- Jour Julien
heure1 = heure + (minute / 60) +(seconde / 3600)
ggg = 1
if( annee < 1582 ) ggg = 0;
if( annee < 1582 && mois < 10 ) ggg = 0;
if( annee <= 1582 && mois == 10 && jour < 5 ) ggg = 0;
jd = - 1 * tronque(7 * (tronque((mois + 9) / 12) + annee) / 4);
s = 1
if ((mois - 9) < 0) s = - 1;
a = abs(mois - 9);
j1 = tronque(y + s * tronque(a / 7));
j2 = - 1 * tronque(((tronque(j1 / 100) + 1) * 3 / 4));
jd1 = jd + tronque((275 * mois / 9) + jour + (ggg * j2));
var jd2dat1 = jd1 + 1721027 + 2 * ggg + 367 * annee - 0.5 ;
jj = jd2dat1 + (heure1 / 24);
datas['jourjulien'] = jj;
joursem(jj);
//--------------- Temps mesuré en siecle juliens de 36525 jours des éphémérides
// a partir du 0,5 janvier 1900 TE
T = (jj - 2415020.0) / 36525;
T2 = T * T;
T3 = T2 * T;
//--------------- Jour Julien, DeltaT et Temps Terrestre
if (annee < 1620 | annee > 2010)
{
coeff = 86400;
T3=(annee - 2000) / 100
if (annee > -391 | annee <= + 946){dtAm = + 2177 + (497 * T3) + (44.1 * T3 * T3); dtA = dtAm}
if (annee > +948 | annee < + 1620){dtAm = + 102 + (102 * T3) + (25.3 * T3 * T3); dtA = dtAm}
}
else {
var dtA = eval('deltaT_' + annee);
coeff = 86400;
}
datas['dt'] = dtA;
jj2 = jj + (dtA / coeff);
datas['jourjulien2'] = jj2;
JD = datas['jourjulien2']-0;
Z = tronque(JD + 0.5);
F = JD + 0.5 - Z;
if (Z < 2299161) { A = Z; }
else
{
I = tronque((Z - 1867216.25) / 36524.25);
A = Z + 1 + I - tronque(I / 4);
}
B = A + 1524;
C = tronque((B - 122.1) / 365.25);
D = tronque(365.25 * C);
T = tronque((B - D) / 30.6001);
var RJ = B - D - tronque(30.6001 * T) + F;
var JJ = tronque(RJ);
var RH = (RJ - tronque(RJ)) * 24;
var Hre = tronque(RH);
var Mn = tronque((RH - Hre ) * 60);
var Sec = ((RH - Hre ) * 60 - Mn ) * 60;
if (T < 14)
{MM = T - 1}
else
{ if ((T == 14) || (T >= 15))
MM = T - 13 }
if (MM > 2) { AA = C - 4716 }
else
{ if ((MM == 1) || (MM == 2)) AA = C - 4715 }
datas['utday'] = JJ;
datas['utmonth'] = MM;
datas['utyear'] = AA;
datas['uth'] = Hre;
datas['utm'] = Mn;
datas['uts'] = floor(Sec);
//---------------
r2d = 180 / PI;
y = 2 * PI;
T = (jj2 - 2415020) / 36525;
T2 = T * T;
T3 = T2 * T;
//Temps mesuré en siecle juliens de 36525 jours des éphémérides
// a partir du 1,5 janvier 2000 TE
var T_2000 = (jj2 - 2451545) / 36525;
var T2_2000 = T_2000 * T_2000;
var T3_2000 = T2_2000 * T_2000;
var T4_2000 = T2_2000 * T2_2000;
//--------------- Longitude moyenne du Soleil, rapportée a l'équinoxe moyen de la date considérée
var Lo = 279.69668 + 36000.76892 * T + 0.0003025 * T2;
Lo = Lo / r2d;
//--------------- Anomalie moyenne du Soleil
M = (358.47583 + 35999.04975 * T) - 0.00015 * T2 - 0.0000033 * T3;
M = M / r2d;
anomalie_soleil = M;
//--------------- Excentricité de l'orbite terrestre
e = 0.01675104 - 0.0000418 * T - 0.000000126 * T2;
//--------------- Equation de centre C du Soleil
c = (1.919460 - 0.004789 * T - 0.000014 * T2) * sin(M) + (0.020094 - 0.000100 * T) * sin(2 * M) + 0.000293 * sin(3 * M);
c = c / r2d;
//--------------- Longitude vraie du Soleil rapportée a l'équinoxe moyen de la date considéré, Anomalie vraie
var Lov = Lo + c;
var v = M + c;
//--------------- Rayon vecteur
R = 1.0000002 *(1 - e * e) / (1 + e * cos(v));
//--------------- Longitude apparente du Soleil rapportée a l'équinoxe vrai de l'époque
var omega = 259.18 / r2d - 1934.142 / r2d * T;
var longapp = Lov - 0.00569 / r2d - 0.00479 / r2d * sin(omega);
var longapp0 = longapp - y * floor(longapp / y);
var longapp1 = longapp0 * r2d;
var d1 = ((longapp1 - floor(longapp1)) * 60);
var d2 = ((d1 - floor(d1)) * 60);
//---------------Corrections a appliquer a la longitude solaire et au rayon vecteur pour une meilleur précision :
//------- A et B sont les corrections dues a l'action de Vénus, C a celle de Jupiter,
//------- D a celle de la Lune, alors que E est une inégalité de longue période
A = 153.23 / r2d + 22518.7541 / r2d * T;
B = 216.57 / r2d + 45037.5082 / r2d * T;
C = 312.69 / r2d + 32964.3577 / r2d * T;
D = 350.74 / r2d + 445267.1142 / r2d * T - 0.00144 / r2d * T2;
E = 231.19 / r2d + 20.20 / r2d * T;
H = 353.40 / r2d + 65928.7155 / r2d * T;
var longitude = longapp + 0.00134 / r2d * cos(A) + 0.00154 / r2d * cos(B) + 0.002 / r2d * cos(C) + 0.00179 / r2d * sin(D) + 0.00178 / r2d * sin(E);
var longitude1 = longapp - y * floor(longapp /y);
var long1 = longitude1 * r2d;
var d1 = ((long1 - floor(long1)) * 60);
var d2 = ((d1- floor(d1)) * 60);
R = R + 0.00000543 * sin(A) + 0.00001575 * sin(B) + 0.00001627 * sin(C) + 0.00003076 * cos(D) + 0.00000927 * sin(H);
//--------------------------------------------- Lune
ct = 0.00000484814;
var longitude_soleil = longitude + 0.00134 / r2d * cos(a) + 0.00154 / r2d * cos(b) + 0.002 / r2d * cos(c) + 0.00179 / r2d * sin(d) + 0.00178 / r2d * sin(e);
longitude_soleil = longitude_soleil - y * floor(longitude_soleil / y);
e = 1 - 0.002516 * T_2000 - 0.0000074 * T2_2000
var excentricite = e * 0.01675104;
var mercure = (252.250906 + 149472.674636 * T_2000) / r2d;
mercure = mercure - y * floor(mercure /y);
var venus = (181.979801 + 58517.815676 * T_2000) / r2d;
venus = venus - y * floor(venus / y);
var longitude_moyenne_terre = (100.46645 + 35999.372854 * T_2000) / r2d;
longitude_moyenne_terre = longitude_moyenne_terre - y * floor(longitude_moyenne_terre /y);
var mars = (355.433275 + 19140.299331 * T_2000) / r2d;
mars = mars - y * floor(mars / y);
var jupiter = (34.351484 + 3034.905675 * T_2000) / r2d;
jupiter = jupiter - y * floor(jupiter / y);
var saturne = (50.077471 + 1222.113794 * T_2000) / r2d;
saturne = saturne - y * floor(saturne /y);
var uranus = ((314 + 3 / 60 + 18.01841 / 3600) + 1542481.19393 / 3600 * T_2000) / r2d;
uranus = uranus - y * floor(uranus / y);
var neptune = ((304 + 20 / 60 + 55.19575 / 3600) + 786550.32074 / 3600 * T_2000) / r2d;
neptune = neptune - y * floor(neptune / y);
a1=(119.75 + 131.849 * T_2000) / r2d;
a2=(53.09 + 479264.29 * T_2000) / r2d;
a3=(313.45 + 481266.484 * T_2000) / r2d;
om = (125.044555 - 1934.1361849 * T_2000 + 0.0020762 * T2_2000 + T3_2000 / 467410 - T4_2000 / 60616000) / r2d;
var longitude_noeud = om;
l = (218.3164591 + 481267.88134236 * T_2000 - 0.0013268 * T2_2000 + T3_2000 / 538841 - T4_2000 / 65194000) / r2d;
l = l - y * floor(l / y);
var w1 = l;
m = (357.5291092 + 35999.0502909 * T_2000 - 0.0001536 * T2_2000 + T3_2000 / 24490000) / r2d;
var l_prime = m;
var anomalie_lune = m;
f = (93.2720993 + 483202.0175273 * T_2000 - 0.0034029 * T2_2000 - T3_2000/3526000 + T4_2000 / 863310000) / r2d;
f = f - y * floor(f / y);
n = (134.9634114 + 477198.8676313 * T_2000 + 0.008997 * T2_2000 + T3_2000 / 69699 - T4_2000 / 14712000) /r2d;
n = n - y * floor(n / y);
petit_l = n;
d = (297.8502042 + 445267.1115168 * T_2000 - 0.00163 * T2_2000 + T3_2000 / 545868 - T4_2000 / 113065000) / r2d;
d = d - y * floor(d / y);
var dzeta_lune = w1 + 5029.0966 / 3600 / r2d * T_2000;
var nu = (-171996 - 174.2 * T_2000) * sin(om);
nu = nu + (-13187 - 1.6 * T_2000) * sin(-2 * d + 2 * f + 2 * om);
nu = nu + (-2274 - 0.2 * T_2000) * sin(2 * f + 2 * om);
nu = nu + (2062 + 0.2 * T_2000) * sin(2 * om);
nu = nu + (1426 - 3.4 * T_2000) * sin(m);
nu = nu + (712 + 0.1 * T_2000) * sin(n);
nu = nu + (-517 + 1.2 * T_2000) * sin(-2 * d + m + 2 * f + 2 * om);
nu = nu + (-386 - 0.4 * T_2000) * sin(2 * f + om);
nu = nu - 301 * sin(n + 2 * f + 2 * om);
nu = nu + (217 - 0.5 * T_2000) * sin(-2 * d - m + 2 * f + 2 * om);
nu = nu - 158 * sin(-2 * d + n);
nu = nu + (129 + 0.1 * T_2000) * sin (-2 * d + 2 * f + om);
nu = nu + 123 * sin (-n + 2 * f + 2 * om);
nu = nu + 63 * sin(2 * d);
nu = nu + (63 + 0.1 * T_2000) * sin (n + om);
nu = nu - 59 * sin(2 * d - n + 2 * f + 2 * om);
nu = nu - (58 - 0.1 * T_2000) * sin(-n + om);
nu = nu - 51 * sin(n + 2 * f + om);
nu = nu + 48 * sin(-2 * d + 2 * n);
nu = nu + 46 * sin(-2 * n + 2 * f + om);
nu = nu - 38 * sin(2 * (d + f + om));
nu = nu - 31 * sin(2 * (n + f + om));
nu = nu + 29 * sin(2 * n);
nu = nu + 29 * sin(-2 * d + n + 2 * f + 2 * om);
nu = nu + 26 * sin(2 * f);
nu = nu - 22 * sin(-2 * d + 2 * f);
nu = nu + 21 * sin(-n + 2 * f + om);
nu = nu + (17 - 0.1 * T_2000) * sin(2 * m);
nu = nu + 16 * sin(2 * d - n + om);
nu = nu - (16 + 0.1 * T_2000) * sin(-2 * d + 2 * m + 2 * f + 2 * om);
nu = nu - 15 * sin(m + om);
nu = nu - 13 * sin(-2 * d + n + om);
nu = nu - 12 * sin(-m - om);
nu = nu + 11 * sin(2 * n - 2 * f);
nu = nu - 10 * sin(2 * d - n + 2 * f + om);
nu = nu - 8 * sin(2 * d + n + 2 * f + 2 * om);
nu = nu + 7 * sin(m + 2 * f + 2 * om);
nu = nu - 7 * sin(-2 * d + m + n);
nu = nu - 7 * sin(-m + 2 * f + 2 * om);
nu = nu - 7 * sin(2 * d + 2 * f + om);
nu = nu + 6 * sin(2 * d + n);
nu = nu + 6 * sin(-2 * d + 2 * n + 2 * f + 2 * om);
nu = nu + 6 * sin(-2 * d + n + 2 * f + om);
nu = nu - 6 * sin(2 * d - 2 * n + om);
nu = nu - 6 * sin(2 * d + om);
nu = nu + 5 * sin(-m + n);
nu = nu - 5 * sin(-2 * d - m + 2 * f + om);
nu = nu - 5 * sin(-2 * d + om);
nu = nu - 5 * sin(2 * n + 2 * f + om);
nu = nu + 4 * sin(-2 * d + 2 * n + om);
nu = nu + 4 * sin(-2 * d + m + 2 * f + om);
nu = nu + 4 * sin(n - 2 * f);
nu = nu - 4 * sin(-d + n);
nu = nu - 4 * sin(-2 * d + m);
nu = nu - 4 * sin(d);
nu = nu + 3 * sin(n + 2 * f);
nu = nu - 3 * sin(-2 * n + 2 * f + 2 * om);
nu = nu - 3 * sin(-d - m - n);
nu = nu - 3 * sin(m + n);
nu = nu - 3 * sin(-m + n + 2 * f + 2 * om);
nu = nu - 3 * sin(2 * d - m - n + 2 * f + 2 * om);
nu = nu - 3 * sin(3 * n + 2 * f + 2 * om);
nu = nu - 3 * sin(2 * d - m + 2 * f + 2 * om);
var nutation_en_longitude = nu / 10000;
nu = (92025 + 8.9 * T_2000) * cos(om);
nu = nu + (5736 - 3.1 * T_2000) * cos(-2 * d + 2 * f + 2 * om);
nu = nu + (977 - 0.5 * T_2000) * cos(2 * f + 2 * om);
nu = nu + (-895 + 0.5 * T_2000) * cos(2 * om);
nu = nu + (54 - 0.1 * T_2000) * cos(m);
nu = nu - 7 * cos(n);
nu = nu + (224 - 0.6 * T_2000) * cos(-2 * d + m + 2 * f + 2 * om);
nu = nu + 200 * cos(2 * f + om);
nu = nu + (129 - 0.1 * T_2000) * cos(n + 2 * f + 2 * om);
nu = nu + (-95 + 0.3 * T_2000) * cos(-2 * d - m + 2 * f + 2 * om);
nu = nu - 70 * cos(-2 * d + 2 * f + om);
nu = nu - 53 * cos(-n + 2 * f + 2 * om);
nu = nu - 33 * cos(n + om);
nu = nu + 26 * cos(2 * d - n + 2 * f + 2 * om);
nu = nu + 32 * cos(-n + om);
nu = nu + 27 * cos(n + 2 * f + om);
nu = nu - 24 * cos(-2 * n + 2 * f + om);
nu = nu + 16 * cos(2 * (d + f + om));
nu = nu + 13 * cos(2 * (n + f + om));
nu = nu - 12 * cos(-2 * d + n + 2 * f + 2 * om);
nu = nu - 10 * cos(-n + 2 * f + om);
nu = nu - 8 * cos(2 * d - n + om);
nu = nu + 7 * cos(-2 * d + 2 * m + 2 * f + 2 * om);
nu = nu + 9 * cos(m + om);
nu = nu + 7 * cos(-2 * d + n + om);
nu = nu + 6 * cos(-m + om);
nu = nu + 5 * cos(2 * d - n + 2 * f + om);
nu = nu + 3 * cos(2 * d + n + 2 * f + 2 * om);
nu = nu - 3 * cos(m + 2 * f + 2 * om);
nu = nu + 3 * cos(-m + 2 * f + 2 * om);
nu = nu + 3 * cos(2 * d + 2 * f + om);
nu = nu - 3 * cos(-2 * d + 2 * n + 2 * f + 2 * om);
nu = nu - 3 * cos(-2 * d + n + 2 * f + om);
nu = nu + 3 * cos(2 * d - 2 * n + om);
nu = nu + 3 * cos(2 * d + om);
nu = nu + 3 * cos(-2 * d - m + 2 * f + om);
nu = nu + 3 * cos(-2 * d + om);
nu = nu + 3 * cos(2 * n + 2 * f + om);
var nutation_en_obliquite = nu / 10000;
var obliquite = (23 + 26 / 60 + 21.448 / 3600 - 46.815 / 3600 * T_2000 - 0.00059 / 3600 * T2_2000 + 0.001813 / 3600 * T3_2000) / r2d + nutation_en_obliquite * ct;
var correct = 22639.55 * sin (n);
correct = correct + 4586.43061 * sin (2 * d - n);
correct = correct + 2369.91227 * sin (2 * d);
correct = correct + 769.02326 * sin (2 * n);
correct = correct + 211.65487 * sin (2 * d - 2 * n);
correct = correct + 205.44315 * sin (2 * d - m - n);
correct = correct + 191.95575 * sin (2 * d + n);
correct = correct + 164.73458 * sin (2 * d - m);
correct = correct + 55.17801 * sin (2 * d - 2 * f);
correct = correct + 39.53393 * sin (n - 2 * f);
correct = correct + 38.42974 * sin (4 * d - n);
correct = correct + 36.12364 * sin (3 * n);
correct = correct + 30.77247 * sin (4 * d - 2 * n);
correct = correct + 17.95512 * sin (d + m);
correct = correct + 14.53078 * sin (2 * d - m + n);
correct = correct + 14.37964 * sin (2 * d + 2 * n);
correct = correct + 13.89903 * sin (4 * d);
correct = correct + 13.194 * sin (2 * d - 3 * n);
correct = correct + 8.60582 * sin (2 * d - m - 2 * n);
correct = correct + 8.05076 * sin (2 * d - 2 * m);
correct = correct + 7.37173 * sin (2 * d - 2 * m - n);
correct = correct + 4.37416 * sin (4 * d - m - n);
correct = correct + 2.73198 * sin (4 * d - m - 2 * n);
correct = correct + 2.48897 * sin (2 * d + m - 2 * n);
correct = correct + 2.14619 * sin (2 * d - m - 2 * f);
correct = correct + 1.977083 * sin (4 * d - m);
correct = correct + 1.26186 * sin (d + m + n);
correct = correct + 1.18682 * sin (4 * d - 3 * n);
correct = correct + 1.17704 * sin (2 * d - m + 2 * n);
correct = correct + 1.07773 * sin (d + m - n);
correct = correct + 1.05949 * sin (2 * d + 3 * n);
correct = correct + .94827 * sin (2 * d - 4 * n);
correct = correct + .75173 * sin (2 * d - 2 * m + n);
correct = correct + .57156 * sin (6 * d - 2 * n);
correct = correct + .47842 * sin (2 * d - m - 3 * n);
correct = correct + .42034 * sin (4 * f);
correct = correct + .41342 * sin (m + 2 * f);
correct = correct + .40423 * sin (3 * d);
correct = correct + .39451 * sin (6 * d - n);
correct = correct + .34966 * sin (d + m - 2 * n);
correct = correct + .33983 * sin (2 * d - 3 * m);
correct = correct + .30874 * sin (4 * d - 2 * m - n);
correct = correct + .30157 * sin (m - n - 2 * f);
correct = correct + .30086 * sin (4 * d - n - 2 * f);
correct = correct + .29422 * sin (2 * d - 2 * m - 2 * n);
correct = correct + .29255 * sin (6 * d - 3 * n);
correct = correct + .28251 * sin (4 * d - m + n);
correct = correct + .27377 * sin (3 * d + m - n);
correct = correct + .26338 * sin (m + n + 2 * f);
correct = correct + .25429 * sin (d + 2 * f);
correct = correct + .24697 * sin (2 * d - 3 * m - n);
correct = correct + .21853 * sin (4 * d + 2 * n);
correct = correct + .17903 * sin (2 * d - n - 2 * f);
correct = correct + .17624 * sin (2 * d + m - 3 * n);
correct = correct + .15781 * sin (4 * d - 2 * m - 2 * n);
correct = correct + .15227 * sin (4 * d - 2 * m);
correct = correct + .1499 * sin (3 * d + m);
correct = correct + .12616 * sin (6 * d);
correct = correct + .111 * sin (5 * n);
correct = correct + .09982 * sin (4 * d - m - 3 * n);
correct = correct + .0932 * sin (2 * d - m + 3 * n);
correct = correct + .09205 * sin (d + m + 2 * n);
correct = correct + .09092 * sin (n + 4 * f);
correct = correct + .09033 * sin (6 * d - m - 2 * n);
correct = correct + .08472 * sin (2 * d + m + n - 2 * f);
correct = correct + .077650 * sin (2 * d + 4 * n);
correct = correct + .07501 * sin (m - 2 * f);
correct = correct + .07142 * sin (6 * d - m - n);
correct = correct + .0685 * sin (2 * d - 5 * n);
correct = correct + .06742 * sin (2 * d + m - n + 2 * f);
correct = correct + .06541 * sin (2 * d + m + 2 * f);
correct = correct + .06507 * sin (3 * d - m);
correct = correct + .06439 * sin (2 * d - 2 * m + 2 * n);
correct = correct + .06314 * sin (2 * d - 2 * m - 2 * f);
correct = correct + .05165 * sin (m - 2 * n - 2 * f);
correct = correct + .0445 * sin (d + n + 2 * f);
correct = correct + .04338 * sin (m + 2 * n + 2 * f);
correct = correct + .04304 * sin (d - 2 * m);
correct = correct + .039 * sin (6* d - m - 3 * n);
correct = correct + .033 * sin (2 * d - 3 * m + n);
correct = correct + .03274 * sin (4 * d - m + 2 * n);
correct = correct + .02949 * sin (2 * d - m - 4 * n);
correct = correct + .02682 * sin (4 * d + m - 3 * n);
correct = correct + .02677 * sin (m + 2 * n - 2 * f);
correct = correct + .0251 * sin (6 * d - m);
correct = correct + .02429 * sin (m - 2 * n + 2 * f);
correct = correct + .02411 * sin (4 * d - 2 * m + n);
correct = correct + .02296 * sin (d + m - 3 * n);
correct = correct + .02289 * sin (4 * d - m - n - 2 * f);
correct = correct + .02285 * sin (6 * d + n);
correct = correct + .02244 * sin (3 * d + m + n);
correct = correct + .02149 * sin (4 * d + 3 * n);
correct = correct + .01993 * sin (2 * d - n + 4 * f);
correct = correct + .01819 * sin (2 * d + m - 4 * n);
correct = correct + .01741 * sin (4 * d - 3 * m - n);
correct = correct + .01605 * sin (2 * d + m + n + 2 * f);
correct = correct + .01598 * sin (d - n + 2 * f);
correct = correct + .01544 * sin (2 * d - 2 * m - 3 * n);
correct = correct + .01376 * sin (6 * d - 4 * n);
correct = correct + .01372 * sin (2 * d + 4 * f);
correct = correct + .01331 * sin (2 * d - 4 * m);
correct = correct + .01297 * sin (2 * n + 4 * f);
correct = correct + .01215 * sin (3 * d - n + 2 * f);
correct = correct + .00971 * sin (4 * d - 3 * m);
correct = correct + .00965 * sin (2 * d - 3 * m - 2 * n);
correct = correct + .00891 * sin (3 * d + m - 2 * f);
correct = correct + .00889 * sin (2 * d + m + 2 * n - 2 * f);
correct = correct + .00866 * sin (8 * d - 2 * n);
correct = correct + .0084 * sin (8 * d - 3 * n);
correct = correct + .00836 * sin (6 * d - 2 * m - 2 * n);
correct = correct + .00812 * sin (2 * d - 4 * m - n);
correct = correct + .00755 * sin (4 * d - 3 * m - 2 * n);
correct = correct + .00744 * sin (6 * d - 2 * m - n);
correct = correct + .0073 * sin (2 * d - m + 4 * n);
correct = correct + .00679 * sin (d + m + 3 * n);
correct = correct + .00666 * sin (4 * d - m - 2 * f);
correct = correct + .00665 * sin (6 * n);
correct = correct + .00662 * sin (4 * d - 2 * n - 2 * f);
correct = correct + .00623 * sin (m - 3 * n - 2 * f);
correct = correct + .00568 * sin (2 * d + 5 * n);
correct = correct + .0056 * sin (4 * d - 2 * m - 3 * n);
correct = correct + .0054 * sin (d + 2 * n + 2 * f);
correct = correct + .00538 * sin (2 * d - 2 * m + 3 * n);
correct = correct + .00526 * sin (m + 3 * n + 2 * f);
correct = correct + .00519 * sin (2 * m + 2 * f);
correct = correct + .00518 * sin (3 * d - 2 * m);
correct = correct + .00515 * sin (2 * d + 2 * m - n + 2 * f);
correct = correct + .00497 * sin (2 * d - 6 * n);
correct = correct + .00477 * sin (6 * d - m + n);
correct = correct + .00475 * sin (5 * d + m - n);
correct = correct + .00473 * sin (2 * m - n - 2 * f);
correct = correct + .00467 * sin (2 * d - 3 * n + 2 * f);
correct = correct + .00455 * sin (8 * d - n);
correct = correct + .00439 * sin (5 * d);
correct = correct + .00392 * sin (5 * d + m - 2 * n);
correct = correct + .00375 * sin (3 * d + 2 * f);
correct = correct + .00364 * sin (6 * d - 2 * n - 2 * f);
correct = correct + .00361 * sin (d + 2 * m - 2 * n);
correct = correct + .00353 * sin (4 * d + m - n + 2 * f);
correct = correct + .00344 * sin (2 * d + n + 4 * f);
correct = correct + .00336 * sin (4 * d - m + 3 * n);
correct = correct + .0033 * sin (3 * d - m + n);
correct = correct + .00324 * sin (8 * d - 4 * n);
correct = correct + .00318 * sin (6 * d + 2 * n);
correct = correct + .00312 * sin (6 * d - 2 * m - 3 * n);
correct = correct + .00298 * sin (3 * d - 2 * n + 2 * f);
correct = correct + .00295 * sin (2 * d - 3 * m + 2 * n);
correct = correct + .0029 * sin (4 * d - 2 * m + 2 * n);
correct = correct + .00289 * sin (d - 2 * n - 2 * f);
correct = correct + .00285 * sin (6 * d - 2 * m);
correct = correct + .00282 * sin (2 * d - 2 * n + 4 * f);
correct = correct + .0027 * sin (2 * m + n + 2 * f);
correct = correct + .00262 * sin (2 * d + m + 2 * n + 2 * f);
correct = correct + .00256 * sin (3 * d + m + 2 * n);
correct = correct + .00254 * sin (d - 3 * m);
correct = correct + .00229 * sin (d - 2 * m - n);
correct = correct + .0022 * sin (4 * d + m - 2 * n + 2 * f);
correct = correct + .00198 * sin (2 * d + m - 4 * f);
correct = correct + .00198 * sin (4 * d + 4 * n);
correct = correct + .00196 * sin (8 * d - m - 2 * n);
correct = correct + .00186 * sin (4 * d + m + 2 * f);
correct = correct + .00183 * sin (4 * d + m + n - 2 * f);
correct = correct + .00181 * sin (5 * d + m);
correct = correct + .00178 * sin (2 * d - m - 5 * n);
correct = correct + .00176 * sin (6 * d - m - 4 * n);
correct = correct + .00173 * sin (2 * d + m - 5 * n);
correct = correct + .0017 * sin (8 * d - m - 3 * n);
correct = correct + .00166 * sin (m + 3 * n - 2 * f);
correct = correct + .00163 * sin (2 * d - 3 * m - 2 * f);
correct = correct + .0016 * sin (4 * d - 3 * m + n);
correct = correct + .00155 * sin (d - m + 2 * f);
correct = correct + .00155 * sin (d + m - 4 * n);
correct = correct + .00153 * sin (3 * n + 4 * f);
correct = correct + .00139 * sin (8 * d);
correct = correct + .00133 * sin (2 * d - 4 * m + n);
correct = correct + .00123 * sin (d - 4 * f);
correct = correct + .00116 * sin (3 * d + m - n - 2 * f);
correct = correct + .00112 * sin (8 * d - m - n);
correct = correct + .00108 * sin (4 * d - 2 * m - n - 2 * f);
correct = correct + .00106 * sin (m - 3 * n + 2 * f);
correct = correct + .00102 * sin (5 * d - m);
correct = correct + .001 * sin (2 * m - 2 * n - 2 * f);
correct = correct + .00096 * sin (2 * d + 2 * m + 2 * f);
correct = correct - 666.44186 * sin (m);
correct = correct - 411.60287 * sin (2 * f);
correct = correct - 147.32654 * sin (m - n);
correct = correct - 124.98806 * sin (d);
correct = correct - 109.38419 * sin (m + n);
correct = correct - 45.10032 * sin (n + 2 * f);
correct = correct - 28.3981 * sin (2 * d + m - n);
correct = correct - 24.3591 * sin (2 * d + m);
correct = correct - 18.58467 * sin (d - n);
correct = correct - 9.67938 * sin (m - 2 * n);
correct = correct - 9.36601 * sin (2 * d - n + 2 * f);
correct = correct - 8.45308 * sin (d + n);
correct = correct - 7.63041 * sin (m + 2 * n);
correct = correct - 7.44804 * sin (2 * m);
correct = correct - 6.38325 * sin (2 * d + n - 2 * f);
correct = correct - 5.7417 * sin (2 * d + 2 * f);
correct = correct - 3.99767 * sin (2 * n + 2 * f);
correct = correct - 3.20968 * sin (3 * d - n);
correct = correct - 2.91464 * sin (2 * d + m + n);
correct = correct - 2.56813 * sin (2 * m - n);
correct = correct - 2.52138 * sin (2 * d + 2 * m - n);
correct = correct - 1.75296 * sin (d - 2 * n);
correct = correct - 1.43724 * sin (2 * d + m - 2 * f);
correct = correct - 1.37259 * sin (2 * n - 2 * f);
correct = correct - 1.22412 * sin (3 * d - 2 * n);
correct = correct - 1.16177 * sin (2 * m + n);
correct = correct - .99023 * sin (2 * d + n + 2 * f);
correct = correct - .6694 * sin (m - 3 * n);
correct = correct - .63523 * sin (4 * d + m - n);
correct = correct - .58399 * sin (d + 2 * n);
correct = correct - .58332 * sin (d - 2 * f);
correct = correct - .56065 * sin (2 * d - 2 * n - 2 * f);
correct = correct - .55694 * sin (d - m);
correct = correct - .54594 * sin (m + 3 * n);
correct = correct - .53572 * sin (2* d - 2 * n + 2 * f);
correct = correct - .4538 * sin (2 * d + 2 * n - 2 * f);
correct = correct - .42624 * sin (2 * d - m - n + 2 * f);
correct = correct - .38215 * sin (2 * d - m + 2 * f);
correct = correct - .37453 * sin (2 * d - m + n - 2 * f);
correct = correct - .35759 * sin (4 * d + m - 2 * n);
correct = correct - .32866 * sin (3 * n + 2 * f);
correct = correct - .29023 * sin (2 * d + m + 2 * n);
correct = correct - .28911 * sin (4 * d + m);
correct = correct - .25304 * sin (3 * d - 2 * f);
correct = correct - .2499 * sin (2 * d + 2 * m - 2 * n);
correct = correct - .23141 * sin (3 * d - m - n);
correct = correct - .20134 * sin (4 * d - n + 2 * f);
correct = correct - .19311 * sin (2 * m - 2 * n);
correct = correct - .18576 * sin (2 * d + 2 * m);
correct = correct - .16977 * sin (4 * d - 2 * n + 2 * f);
correct = correct - .13636 * sin (d - m - n);
correct = correct - .12812 * sin (d - 3 * n);
correct = correct - .12386 * sin (2 * d + 2 * n + 2 * f);
correct = correct - .12073 * sin (d - m + n);
correct = correct - .10136 * sin (3 * m);
correct = correct - .09154 * sin (2 * d - 3 * n - 2 * f);
correct = correct - .085 * sin (4 * d + 2 * f);
correct = correct - .08311 * sin (3 * d - m - 2 * n);
correct = correct - .08282 * sin (m + n - 2 * f);
correct = correct - .08049 * sin (m - n + 2 * f);
correct = correct - .08019 * sin (n - 4 * f);
correct = correct - .07518 * sin (2 * d - 4 * f);
correct = correct - .07373 * sin (2 * d - m + n + 2 * f);
correct = correct - .06601 * sin (4 * d + n - 2 * f);
correct = correct - .06513 * sin (2 * m + 2 * n);
correct = correct - .06103 * sin (2 * d - m - n - 2 * f);
correct = correct - .05725 * sin (5 * d - 2 * n);
correct = correct - .05684 * sin (3 * n - 2 * f);
correct = correct - .05142 * sin (3 * m - n);
correct = correct - .0507 * sin (4 * d + m + n);
correct = correct - .04702 * sin (m - 4 * n);
correct = correct - .04442 * sin (3 * d - 3 * n);
correct = correct - .04189 * sin (3 * d + m - 2 * n);
correct = correct - .04074 * sin (d + 3 * n);
correct = correct - .04012 * sin (d + n - 2 * f);
correct = correct - .03968 * sin (d + 2 * m);
correct = correct - .03947 * sin (m + 4 * n);
correct = correct - .03587 * sin (d + m + 2 * f);
correct = correct - .03514 * sin (4 * d + 2 * m - 2 * n);
correct = correct - .03336 * sin (2 * d + 3 * n - 2 * f);
correct = correct - .02979 * sin (3 * d - n - 2 * f);
correct = correct - .02887 * sin (2 * d - m + 2 * n - 2 * f);
correct = correct - .02804 * sin (2 * d - m - 2 * n - 2 * f);
correct = correct - .02676 * sin (2 * d + m + 3 * n);
correct = correct - .02602 * sin (4 * n + 2 * f);
correct = correct - .02391 * sin (4 * d - 2 * f);
correct = correct - .02379 * sin (d - n - 2 * f);
correct = correct - .02349 * sin (2 * d + 2 * m - 2 * f);
correct = correct - .02273 * sin (4 * d - m - n + 2 * f);
correct = correct - .02171 * sin (4 * d + 2 * m - n);
correct = correct - .02157 * sin (2 * d - m - 2 * n + 2 * f);
correct = correct - .01948 * sin (3 * d - m - 2 * f);
correct = correct - .01875 * sin (4 * d + m - n - 2 * f);
correct = correct - .01816 * sin (2 * d - 2 * m + 2 * f);
correct = correct - .01796 * sin (3 * m + n);
correct = correct - .01781 * sin (4 * d + n + 2 * f);
correct = correct - .01686 * sin (5 * d - 3 * n);
correct = correct - .01644 * sin (2 * d - 2 * m + n - 2 * f);
correct = correct - .01541 * sin (2 * d - 2 * m - n + 2 * f);
correct = correct - .01533 * sin (4 * d - m - 2 * n + 2 * f);
correct = correct - .01514 * sin (2 * m - 3 * n);
correct = correct - .01483 * sin (d - m + 2 * n);
correct = correct - .0135 * sin (5 * d - n);
correct = correct - .01343 * sin (2 * d + 2 * m + n);
correct = correct - .01332 * sin (2 * d + 3 * n + 2 * f);
correct = correct - .01282 * sin (6 * d + m - 2 * n);
correct = correct - .01281 * sin (d - m - 2 * f);
correct = correct - .01182 * sin (3 * d - 2 * m - n);
correct = correct - .01114 * sin (4 * d - m + 2 * f);
correct = correct - .01077 * sin (2 * d - 4 * n - 2 * f);
correct = correct - .01064 * sin (6 * d + m - n);
correct = correct - .01062 * sin (3 * d + n - 2 * f);
correct = correct - .01007 * sin (2 * d - m + 2 * n + 2 * f);
correct = correct - .0098 * sin (4 * d + 2 * n - 2 * f);
correct = correct - .00955 * sin (d - 4 * n);
correct = correct - .00944 * sin (2 * d + 2 * m - 3 * n);
correct = correct - .00934 * sin (4 * d - 3 * n + 2 * f);
correct = correct - .0085 * sin (2 * d - n - 4 * f);
correct = correct - .00849 * sin (d + 2 * m + n);
correct = correct - .00732 * sin (4 * d - m + n - 2 * f);
correct = correct - .00694 * sin (d - m - 2 * n);
correct = correct - .00693 * sin (5 * d - m - 2 * n);
correct = correct - .00668 * sin (4 * d + m + 2 * n);
correct = correct - .00659 * sin (d + m + n + 2 * f);
correct = correct - .00654 * sin (2 * d + 2 * m - n - 2 * f);
correct = correct - .00623 * sin (3 * d + m - 3 * n);
correct = correct - .00509 * sin (6 * d - 2 * n + 2 * f);
correct = correct - .00478 * sin (6 * d + m - 3 * n);
correct = correct - .00434 * sin (2 * d - 2 * m - n - 2 * f);
correct = correct - .00431 * sin (4 * d - 5 * n);
correct = correct - .00416 * sin (3 * m - 2 * n);
correct = correct - .00399 * sin (3 * d - 2 * m - 2 * n);
correct = correct - .00396 * sin (6 * d + m);
correct = correct - .00389 * sin (3 * d + 2 * n);
correct = correct - .00378 * sin (2 * d - 2 * m + n + 2 * f);
correct = correct - .00369 * sin (4 * d + 2 * m - 3 * n);
correct = correct - .00365 * sin (2 * d - m - 3 * n - 2 * f);
correct = correct - .00359 * sin (6 * d - n + 2 * f);
correct = correct - .00355 * sin (2 * m - 2 * f);
correct = correct - .00354 * sin (4 * n - 2 * f);
correct = correct - .00346 * sin (2 * d + m - 2 * n - 2 * f);
correct = correct - .00341 * sin (2 * m + 3 * n);
correct = correct - .00335 * sin (5 * d - n - 2 * f);
correct = correct - .00332 * sin (m - 5 * n);
correct = correct - .003 * sin ( d + 2 * m - n);
correct = correct - .00297 * sin (3 * d - m - 3 * n);
correct = correct - .00287 * sin (m + 5 * n);
correct = correct - .00287 * sin (6 * d - 3 * n + 2 * f);
correct = correct - .00286 * sin (2 * d - m - 4 * f);
correct = correct - .00285 * sin (d + 4 * n);
correct = correct - .00274 * sin (4 * d + 2 * n + 2 * f);
correct = correct - .00251 * sin (4 * d - m + n + 2 * f);
correct = correct - .00247 * sin (2 * d + 4 * n - 2 * f);
correct = correct - .00236 * sin (2 * d + m + 4 * n);
correct = correct - .00232 * sin (2 * d - m + 3 * n - 2 * f);
correct = correct - .00228 * sin (2 * d + m - n - 2 * f);
correct = correct - .00214 * sin (6 * d - 2 * f);
correct = correct - .00212 * sin (d - m + n - 2 * f);
correct = correct - .00208 * sin (4 * d + 2 * m);
correct = correct - .00201 * sin (5 * n + 2 * f);
correct = correct - .002 * sin (2 * d + 2 * m + n - 2 * f);
correct = correct - .00191 * sin (3 * d + 2 * m);
correct = correct - .00189 * sin (3 * d - m - n - 2 * f);
correct = correct - .00189 * sin (5 * d - m - 3 * n);
correct = correct - .00188 * sin (2 * d + 3 * m - n);
correct = correct - .00174 * sin (3 * d - 4 * n);
correct = correct - .0016 * sin (4 * d - 2 * m - n + 2 * f);
correct = correct - .00157 * sin (d + m + n - 2 * f);
correct = correct - .00154 * sin (5 * d - m - n);
correct = correct - .00149 * sin (d - m + 3 * n);
correct = correct - .00142 * sin (d - 2 * n + 2 * f);
correct = correct - .00138 * sin (3 * d + m - n + 2 * f);
correct = correct - .00137 * sin (5 * d - 2 * f);
correct = correct - .00133 * sin (2 * d - 2 * m + 2 * n - 2 * f);
correct = correct - .00132 * sin (6 * d + 2 * f);
correct = correct - .00131 * sin (2 * d + 4 * n + 2 * f);
correct = correct - .00128 * sin (4 * m);
correct = correct - .00127 * sin (3 * d + 2 * m - n);
correct = correct - .00121 * sin (4 * d - m + 2 * n - 2 * f);
correct = correct - .00119 * sin (2 * m - 4 * n);
correct = correct - .00117 * sin (2 * d - m + 3 * n + 2 * f);
correct = correct - .00116 * sin (2 * d + m - 3 * n - 2 * f);
correct = correct - .00111 * sin (2 * d - 2 * m - 2 * n - 2 * f);
correct = correct - .00111 * sin (2 * d - 5 * n - 2 * f);
correct = correct - .00109 * sin (4 * d + 3 * n - 2 * f);
correct = correct - .00108 * sin (4 * m - n);
correct = correct - .00102 * sin (d + 2 * m + 2 * n);
correct = correct - .00102 * sin (3 * d - 2 * m - 2 * f);
correct = correct - .001 * sin (d - m - n - 2 * f);
correct = correct - .00098 * sin (7 * d - 3 * n);
correct = correct + 14.2488 * sin (18 * venus - 16 * longitude_moyenne_terre - n + 26.54261 / r2d);
correct = correct + 1.1431 * sin (2 * longitude_moyenne_terre - 2 * jupiter + 2 * d - n + 180.11977 / r2d);
correct = correct + 0.9011 * sin (4 * longitude_moyenne_terre - 8 * mars + 3 * jupiter + 285.98707 / r2d);
correct = correct + 0.8216 * sin (venus - longitude_moyenne_terre + 180.00988 / r2d);
correct = correct + 0.7881 * sin (18 * venus - 16 * longitude_moyenne_terre - 2 * n + 26.54324 / r2d);
correct = correct + 0.7393 * sin (18 * venus - 16 * longitude_moyenne_terre + 26.54560 / r2d);
correct = correct + 0.6437 * sin (3 * venus - 3 * longitude_moyenne_terre + 2 * d - n + 179.98144 / r2d);
correct = correct + 0.6388 * sin (longitude_moyenne_terre - jupiter + 1.22890 / r2d);
correct = correct + 0.5634 * sin (10 * venus - 3 * longitude_moyenne_terre - n + 333.30551 / r2d);
correct = correct + 0.4453 * sin (2 * longitude_moyenne_terre - 3 * jupiter + 2 * d - n + 10.07001 / r2d);
correct = correct + 0.3436 * sin (2 * venus - 3 * longitude_moyenne_terre + 269.95393 / r2d);
correct = correct + 0.3246 * sin (longitude_moyenne_terre - 2 * mars + 318.13776 / r2d);
correct = correct + 0.3016 * sin (2 * venus - 2 * longitude_moyenne_terre + 0.20448 / r2d);
correct = correct + 7.06304 * sin (dzeta_lune - f + 0.00094 / r2d);
correct = correct + 0.49331 * sin (dzeta_lune + petit_l - f + 0.00127 / r2d);
correct = correct + 0.49141 * sin (dzeta_lune - petit_l - f + 0.00127 / r2d);
correct = correct + 0.36061 * sin (dzeta_lune + f + 0.00071 / r2d);
correct = correct + 0.09642 * sin (dzeta_lune + 2 * d - f + 0.0009 / r2d);
correct = correct + 0.06569 * sin (dzeta_lune - 2 * d - f + 0.001 / r2d);
correct = correct + 0.06456 * sin (dzeta_lune + 2 * d - petit_l - f + 0.00042 / r2d);
correct = correct + 0.05036 * sin (dzeta_lune - petit_l + f + 0.00051 / r2d);
correct = correct + 0.04962 * sin (dzeta_lune - 2*d + petit_l - f + 0.00029 / r2d);
correct = correct + 0.04746 * sin (dzeta_lune - 2 * d + f + 0.00076 / r2d);
correct = correct + 0.03838 * sin (dzeta_lune + petit_l + f + 0.0007 / r2d);
correct = correct + 0.03638 * sin (2 * dzeta_lune - 2 * f + 180 / r2d);
correct = correct + 0.03402 * sin (dzeta_lune + 2 * petit_l - f + 0.00126 / r2d);
correct = correct + 0.03279 * sin (dzeta_lune - 2 * petit_l - f + 0.00128 / r2d);
correct = correct + 0.02206 * sin (2 * d - petit_l);
correct = correct + 0.01492 * sin (dzeta_lune - 3 * f + 180.00086 / r2d);
correct = correct + 0.01234 * sin (dzeta_lune + 2 * d + petit_l - f + 0.00102 / r2d);
l = l + correct / 3600 / r2d + nutation_en_longitude * ct;
longitude_lune = l;
l1 = l * r2d;
d1 = ((l1 - floor(l1)) * 60);
d2 = ((d1 - floor(d1)) * 60);
var sg = strSigne(l1);
datas['longlune'] = "
" + (floor(l1) % 30) + "° " + sg + " " + zero(floor(d1)) + "'";
if (GetObject("longilune") !== null)
{
var node = GetObject("longilune");
node.innerHTML = datas['longlune'] + '
' +
datas['day'] +'/' + datas['month'] +'/' + datas['year'] + ' ' +
datas['hour'] + ':'+datas['minut']+ ' GMT
Longitude tropicale de la Lune'+'';
}
}
} // calculer