Re: stima orari alba e tramonto a date coordinate geografiche (più data calendario)

From: Angelo Vannini <iz4dzd_at_gmail.com>
Date: Mon, 17 Aug 2020 01:16:27 -0700 (PDT)

Il giorno domenica 16 agosto 2020 alle 23:30:02 UTC+2 Soviet_Mario ha scritto:
> On 16/08/20 20:51, Elio Fabri wrote:
> > Soviet_Mario ha scritto:
> >> si, l'avevo capito, infatti sono alla ricerca di una nuova
> >> fonte ...
> >> sperando che qui qualcuno sia appassionato di astronomia
> >> anche fina
> >> a sé stessa
> > Va bene, tu cerca. Io a tempo perso vedo se riesco a dare
> > una formula
> > semplice.
> > E' un'occasiome per ristudiare (sui miei appunti) cose che ho
> > insegnato tanti anni fa :-)
> Elio, come dire : solo se la cosa ti diverte. Per quel che
> mi riguarda cercavo una formula come surrogato di una tavola
> dati pesante.
>
> Ma era una posizione sciocca in principio : visto che non mi
> serve precisione, Asmagorod mi ha già fatto notare che va
> benissimo "sotto-campionare" un calendario completo, ad es.
> con un dato a settimana, e poi interpolare (ecco boh magari
> non proprio linearmente, ma imho pure linearmente non è che
> si commettano errori tanto ingenti in una settimana).
>
> Ergo forse, utilitaristicamente parlando, mi sarebbe alla
> fine più utile una fonte reputabile per estrarre i dati,
> almeno alba e tramonto.
> L'altezza del sole l'avevo pescata solo per quel problema,
> ma senza molto sperare di poterla sfangare con poco.
> >
> >> non so se possa servire. Ho un antenna ancora montata, una
> >> parabola enorme (che ho montato da solo beccando il
> >> satellite con un
> >> cicalino ... ore e ore di tira e molla viti
> > Ho capito. Del GPS non sai proprio niente :-)
> no niente, salvo che esistono satelliti che "ti triangolano"
> (non so se sia giusto)
> >
> >> il KA-SAT penso fosse un geostazionario a 36000 km.
> > Appunto!
> >
> >> non so, ma ho come l'impressione che sia un compito al di
> >> sopra
> >> delle mie capacità ...
> > Premesso quanto ti ho già detto sul mistero delle tue
> > capacità, ora
> > provo a dirti qualcosa sul GPS.
> > Intanto non è *un* satellite ma una *costellazione*
> si questo lo sapevo (so che anche i Cinesi ora hanno la loro
> costellazione, e i russi il Glonass)
> > (attualmente una
> > trentina) su orbite a circa 20000 km dalla sup. terrestre,
> > con periodi
> > di 12 ore.
> > Quindi ben lontani dall'essere geostazionari.
> > Dai un'occhiata a "gps" in wiki (inglese): c'è una bella
> > animazione
> > con la Terra che gira e i satelliti che orbitano. Nella
> > versione
> > italiana la figura non c'è :-(
> ci guarderò, le figure le capisco !
> (beh, mica tutte, ma questa magari si)
> >
> > Dato che il concetto base è che ogni punto della Terra possa
> > vedere
> > simultaneamente almeno 4 satelliti (di più è meglio)
> > scordati il tuo
> > parabolone: le antenne per GPS sono *omnidirezionali*.
> > Non ti sto a dire come fa il sistema a ottenere la posizione
> > del
> > ricevitore, perché non è questa l'idea.
> > Come ho già detto, si tratta di ottenere invece da ciascun
> > satellite
> > dati sufficienti a calcolare la *sua* posizione nel cielo
> > dell'antenna
> > ricevente, in particolare al momento della sua apparizione o
> > scomparsa
> > dietro gli ostacoli (tipicamente montagne, ma anche altro:
> > per es. ho
> > sperimentato che basta un bosco di abeti per occultare i
> > satelliti).
> >
> > Questi dati ci sono mel messaggio che ogni satellite
> > trasmette, anche
> > se ora ho dimenticato come.
> > Qui entra in ballo il software del ricevitore.
> > Per inciso ho verificato l'ovvio: esistono antenne GPS che si
> > connettono a una porta USB di un PC. Quindi non occorre
> > smartphone.
> > Ed esistono vari software. Ne ho scaricato uno free per
> > Linux, anche
> > se non potrò usarlo a meno che non mi compri un'antenna (fra
> > i 10 e i
> > 30€).
> manco tanto !!! Cavoli, così la cosa mi ingustisce. Sempre
> che il software faccia lui i "ragionamenti" che mi spieghi
> > Vorrei vedere se mi permette di "ficcare il naso" nei dati,
> > come mi
> > permetteva quello di tanti anni fa, o se invece fa tutto lui
> > fino al
> > calcolo finale della posizione del ricevitore.
> che per me, come dire, sarebbe la soluzione perfetta.
> >
> >> In effetti potrei fare un pannellino solare
> > Su questo non dico niente, sia perché nomn ho esperienza,
> > sia perché
> > m'interessa meno :-)
> beh ma non c'è poi molto da dire : avere una quaterna anche
> solo di fotocelle agli angoli estremi del rettangolo
> pannellato è un modo empirico per stimare l'entità
> dell'illuminazione.
> Se non gli ho dato inizialmente molti punti, è perché non
> sapendo che "cervello" dovrò comprare, non so nemmeno quante
> porte logiche abbia e come gestire.
> Sull'elettronica ho cercato di volare bassissimo, il minimo
> indispendabile : sensori di T, e ... boh, relays ? per
> comandare pompa e elettrovalvole.
>
> Poi si è intromesso a forza bruta il display e/o i led, e
> pure l'orologio (che però penso si potrà connettere con meno
> problemi del resto, nascendo digitale)
> >
> >> guarda, devo un po' capire bene cosa possa fare quel che già
> >> ho in ambiti (GPS) che non ho mai indagato.
> > Mi sa che qui ti sei mangiato qualche parola...
> no mi sono espresso proprio male. Intendevo che ipotizzavo
> che magari possedessi già qualche hardware inutilizzato. Ma
> da come hai spiegato no : l'antenna GPS quella non l'ho mai
> avuta.
> >
> >> Non so se la griglia cellulare offra di suo servizi simili,
> >> o se possa semplicemente, appoggiandosi ai satelliti,
> >> indirettamente offrire questo servizio GPS.
> > Secondo me te lo puoi scordare.
> in effetti ... come potrebbe la rete cellulare interrogare
> il GPS, se essa stessa non ti può localizzare con una
> risoluzione particolarmente buona ? (specie quella
> verticale, penso), sull'orizzontale mi pare che la
> risoluzione sia un metro dove c'è buona copertura, e forse
> qualche metro qui dove è un disastro totale ... magari però
> mi sbaglio o ricordo male.
> >
>
>
> --
> 1) Resistere, resistere, resistere.
> 2) Se tutti pagano le tasse, le tasse le pagano tutti
> Soviet_Mario - (aka Gatto_Vizzato)
Source:
        Almanac for Computers, 1990
        published by Nautical Almanac Office
        United States Naval Observatory
        Washington, DC 20392

Inputs:
        day, month, year: date of sunrise/sunset
        latitude, longitude: location for sunrise/sunset
        zenith: Sun's zenith for sunrise/sunset
          offical = 90 degrees 50'
          civil = 96 degrees
          nautical = 102 degrees
          astronomical = 108 degrees
        
        NOTE: longitude is positive for East and negative for West
        NOTE: the algorithm assumes the use of a calculator with the
        trig functions in "degree" (rather than "radian") mode. Most
        programming languages assume radian arguments, requiring back
        and forth convertions. The factor is 180/pi. So, for instance,
        the equation RA = atan(0.91764 * tan(L)) would be coded as RA
        = (180/pi)*atan(0.91764 * tan((pi/180)*L)) to give a degree
        answer with a degree input for L.


1. first calculate the day of the year

        N1 = floor(275 * month / 9)
        N2 = floor((month + 9) / 12)
        N3 = (1 + floor((year - 4 * floor(year / 4) + 2) / 3))
        N = N1 - (N2 * N3) + day - 30

2. convert the longitude to hour value and calculate an approximate time

        lngHour = longitude / 15
        
        if rising time is desired:
          t = N + ((6 - lngHour) / 24)
        if setting time is desired:
          t = N + ((18 - lngHour) / 24)

3. calculate the Sun's mean anomaly
        
        M = (0.9856 * t) - 3.289

4. calculate the Sun's true longitude
        
        L = M + (1.916 * sin(M)) + (0.020 * sin(2 * M)) + 282.634
        NOTE: L potentially needs to be adjusted into the range [0,360) by adding/subtracting 360

5a. calculate the Sun's right ascension
        
        RA = atan(0.91764 * tan(L))
        NOTE: RA potentially needs to be adjusted into the range [0,360) by adding/subtracting 360

5b. right ascension value needs to be in the same quadrant as L

        Lquadrant = (floor( L/90)) * 90
        RAquadrant = (floor(RA/90)) * 90
        RA = RA + (Lquadrant - RAquadrant)

5c. right ascension value needs to be converted into hours

        RA = RA / 15

6. calculate the Sun's declination

        sinDec = 0.39782 * sin(L)
        cosDec = cos(asin(sinDec))

7a. calculate the Sun's local hour angle
        
        cosH = (cos(zenith) - (sinDec * sin(latitude))) / (cosDec * cos(latitude))
        
        if (cosH > 1)
          the sun never rises on this location (on the specified date)
        if (cosH < -1)
          the sun never sets on this location (on the specified date)

7b. finish calculating H and convert into hours
        
        if if rising time is desired:
          H = 360 - acos(cosH)
        if setting time is desired:
          H = acos(cosH)
        
        H = H / 15

8. calculate local mean time of rising/setting
        
        T = H + RA - (0.06571 * t) - 6.622

9. adjust back to UTC
        
        UT = T - lngHour
        NOTE: UT potentially needs to be adjusted into the range [0,24) by adding/subtracting 24

10. convert UT value to local time zone of latitude/longitude
        
        localT = UT + localOffset
Received on Mon Aug 17 2020 - 10:16:27 CEST