Calculateur des heures de prière.

rem ============================================================================
rem          Calculs_Astronomiques.bas Par Papydall
rem ============================================================================
'             Les formules utilisées dans ce programme sont issues du livre
'                        "Astronomical Algorithms" de Jean Meeus
rem ============================================================================

Init()

Gregorian_To_Hijri(jour,mois,annee)

message result$
terminate

end
rem ============================================================================
' déclaration des variables globales
SUB Init()
    dim heure,minute,seconde,jour,mois,annee, result$, yaum$
    Dim b, Y, m, a, jj, t, t2, Ct1, Ct, T2000, T2_2000, T3_2000, T4_2000, l, c, v
    dim Omega, longitude, d, e, h, longitude_soleil, excentricite, mercure
    Dim venus, longitude_moyenne_terre, anomalie_soleil
    Dim mars, jupiter, saturne, uranus, neptune
    Dim a1, a2, a3, Om, Longitude_noeud, W1, l_prime, Anomalie_lune, f, n
    dim petit_l, dzeta_lune, nu, nutation_en_longitude
    Dim nutation_en_obliquite, Obliquite, Correction, longitude_lune, l1, d1, d2
    dim latitude, latitude_lune, b1, p, distance_terre, p1, ascension, Bidouille
    dim declinaison, fraction_illuminee, Signe$, elongation_lune$
    Get_Date_heuree()
    Application_title "Calculs astronomiques par PAPYDALL"
    hide 0
END_SUB
rem ============================================================================
' Lire la date et l'heuree système
SUB Get_Date_heuree()
    heure = val(mid$(time$,1,2))
    minute = val(mid$(time$,4,2))
    seconde = val(mid$(time$,7,4))
    jour = val(mid$(date$,1,2))
    mois = val(mid$(date$,4,2))
    annee = val(mid$(date$,7,4))
END_SUB
rem ============================================================================
' Déterminer le Jour Julien à partir d'une date
SUB Jour_Julien(jour,mois,annee)
    dim_local a,b,d,m,y,JD
    d = jour : m = mois : y = annee
    if m < 3 then y = y - 1 : m = m + 12
    a = int(y/100) : b = 2 - a + int(a/4)
    JD = int(365.25*(y+4716)) + int(30.6001*(m+1)) + d + b - 1524.5
    result$ =  str$(JD)
END_SUB
rem ============================================================================
' Déterminer le nom du jour de la semaine en français dans result$
' et en Arabe dans yaum$
SUB Day_Of_Week(jour,mois,annee)
    dim_local j,m,a,jd, r
    j = jour : m = mois : a = annee
    Jour_Julien(j,m,a) : jd = val(result$) + 1.5
    r = mod(jd,7)
    select r
        case 0 : result$ = "DIMANCHE" : yaum$ = "AL-AHAD"
        case 1 : result$ = "LUNDI"    : yaum$ = "AL-ETHNINE
        case 2 : result$ = "MARDI"    : yaum$ = "AL-THOULATHA"
        case 3 : result$ = "MERCREDI" : yaum$ = "AL-ERBI'A"
        case 4 : result$ = "JEUDI"    : yaum$ = "AL-KHAMISS"
        case 5 : result$ = "VENDREDI" : yaum$ = "AL-JOUMOU'A"
        case 6 : result$ = "SAMEDI"   : yaum$ = "AL-SEBT"
    end_select
END_SUB
rem ============================================================================
' Conversion de date Grégorienne (ou julienne) vers date Héjirienne (musulmane)
SUB Gregorian_To_Hijri(jour,mois,annee)
    dim_local x,m,d,alfa,beta,a,b,c,e,c1,c2,dprime,q,r,s,k,o,j,h,jj,cl,dl,n,w
    dim day,month,year,m$,m1$,j$
' Si la date est donnée dans le calendrier grégorien, la convertir dans le
' calendrier julien
    x = annee : m = mois : d = jour
    if m < 3 then x = x - 1 : m = m + 12
    alfa = int(x/100) : beta = 2 - alfa + int(alfa/4)
    b = int(365.25*x) + int(30.6001*(m+1)) + d + 1722519 + beta
    c = int((b-122.1)/365.25) : d = int(365.25*c)
    e = int((b-d)/30.6001) : day = b-d-int(30.6001*e)
    if e < 14 then month = e-1
    if e > 13 then month = e-13
    if month > 2 then year = c-4716
    if month < 3 then year = c-4715
    if mod(x,4) = 0 then w = 1 : else : w = 2
    n = int(275*month/9) - w * int((month+9)/12) + day - 30
    a = x - 623 : b = int(a/4) : c = mod(a,4) : c1 = 365.2501*c : c2 = int(c1)
    if (c1 - c2) > 0.5 then c2 = c2 + 1
    dprime = 1461 * b + 170 + c2 : q = int(dprime/10631) : r = mod(dprime,10631)
    j = int(r/354) : k = mod(r,354) : o = int((11*j + 14)/30)
    h = 30*q + j + 1 : jj = k - o + n - 1
' jj est le numéro du jour dans l'année Héjirienne H
' On détermine si H est commune ou bissextile (KABISSA)
    cl = mod(h,30) : dl = mod((11*cl + 3),30)
    if dl < 19 then jj = jj - 354 : h = h + 1 : if jj < 0 then jj = jj + 354
    if dl > 18 then jj = jj - 355 : h = h + 1 : if jj < 0 then jj = jj + 355
    if jj = 0 then jj = 355 : h = h - 1
    s = int((jj-1)/29.5) : m = 1 + s : d = int(jj-29.5 * s)
    if jj = 355 then m = 12 : d = 30

    select m
        case 01 : m$ = "MOUHARRAM"
        case 02 : m$ = "SAFAR"
        case 03 : m$ = "RABI'AL-AWWAL"
        case 04 : m$ = "RABI'ATH-THANI"
        case 05 : m$ = "JOUMADA AL-OULA"
        case 06 : m$ = "JOUMADA AL-THANIA"
        case 07 : m$ = "RAJAB"
        case 08 : m$ = "CHAÂBAN"
        case 09 : m$ = "RAMADAN"
        case 10 : m$ = "CHAWWAL"
        case 11 : m$ = "DHOUL-QAÂDA"
        case 12 : m$ = "DHOUL-HIJJA"
    end_select
    
    select mois
       case 1 : m1$ = "JANVIER"
       case 2 : m1$ = "FEVRIER"
       case 3 : m1$ = "MARS"
       case 4 : m1$ = "AVRIL"
       case 5 : m1$ = "MAI"
       case 6 : m1$ = "JUIN"
       case 7 : m1$ = "JUILLET"
       case 8 : m1$ = "AOÛT"
       case 9 : m1$ = "SEPTEMBRE"
       case 10: m1$ = "OCTOBRE"
       case 11: m1$ = "NOVEMBRE"
       case 12: m1$ = "DECEMBRE"
    end_select
    
    Day_Of_Week(jour,mois,annee) : j$ = result$
    
    result$ = j$ + " " + str$(jour)+" "+ m1$ +" "+str$(annee)+ " correspond au : "
    result$ = chr$(13) + result$  + yaum$ + " " +  str$(d) +" "+ m$ + " " + str$(h) + " HIJRI"

END_SUB
rem ============================================================================