| |
I3XTY > TUTTI 19.05.09 20:27l 571 Lines 22674 Bytes #999 (0) @ ITA
BID : 40019_I3XTY
Read: GUEST IV3JER IZ6WQP
Subj: MSG n.4 di I3ZJV
Path: IZ3LSV<IK3GET<I3LUG<I3XTY<I3XTY
Sent: 090519/2103 @:I3XTY.TV.IVEN.ITA.EU $:40019_I3XTY Sally 4.1.232
Date: mar, 19 mag 2009 21:03:21 LT
From: I3XTY (Luigi)
To: TUTTI@ITA
Subject: MSG n.4 di I3ZJV
Inviato da: I3XTY@I3XTY.TV.IVEN.ITA.EU
From: I3ZJV@I3XTY.TV.IVEN.ITA.EU
To : ASTRA@ITA
PROGRAMMA DI CALCOLO DELLE EFFEMERIDI DEL SOLE , DELLA LUNA, E DEI PIANETI
PER UN GIORNO E UN'ORA DATI, E PER UN LUOGO DATO
(era in Applesoft tradotto in IBM PC)
(i numeri di sinistra devono essere uno di seguito all'aitro senza
nessun salto !! : Ricordarsi di questo particolare !!)
(Il basic si comanda con : alt F)
Per farlo scorrere, il programma domanda l'introduzione delle coordinate del
luogo di osservazione, la latitudine e la longitudine in gradi e decimali,
piu' il giorno sotto forma JJ,MM,AAAA, dato in cifre separate con delle vir-
gole, esempio : 3,5,1973, infine l'ora in tempo universale, ore e minuti se-
parati da una virgola.(La longitudine ad est di Greenwich e' negativa, ci si
deve mettere il meno davanti)
IL PROGRAMMA CALCOLA E RIPORTA :
- N il numero dei giorni trascorsi dopo il zero gennaio 1901 e il giorno della
settimana, questo per il controllo se e' vero,
- La longitudine, l'ascensione retta, la declinazione del sole, la sua ora UTC
di levata e tramonto per il luogo di osservazione. Le coordinate cartesiane
sono ugualmente calcolate ma non esposte, cioe' non si vedono
Il programma domanda l'introduzione delle due prime lettere del pianeta desi-
derato ( oppure LU per la luna, cioe' LU ME VE MA JU SA UR NE PL.
Con ET si calcola qualsiasi astro immettendo le sue coordnate eliocentriche
per avere quelle geocentriche immediatamente, cioe' per sapere dove puntare
il telescopio.Inoltre su SA si possono vedere dove si trovano Thetys, Dione,
Rhea e Titano, battendo prima OUI , poi O.
SU JU dovrebbero vedersi dove si trovano IO Europa Ganimede Callisto, ma do-
vrei sistemare ancora una riga. Ho controllato i dati col sole, vanno OK.
Con la luna vanno OK, va bene con tutti i satelliti)
PER CIASCUN PIANETA DA' :
- longitudine e latitudine
- ascensione retta e declinazione
- azimut e altezza
- diametro apparente e magnitudine
- ora di levata e del tramonto
Per Giove e Saturno, fornisce in piu',su richiesta il diagramma dei quattro
principali satelliti, battendo prima OUI indi O
PER LA LUNA FORNISCE .
- longitudine e latitudine
- ascensione retta e declinazione
- fase
- ora di levata e tramonto in ore UTC
L'entrata di ST in luogo delle due prime lettere provoca un END. Il programma
deve essere riinizializzato con un RUN per un nuovo calcolo, cioe' si deve
mettere da capo latitudine, longitudine negativa (per noi),data, ora UTC
Dopo aver mostrato le informazioni concernenti un pianeta oppure la luna, bat-
tendo un tasto qualsiasi provoca il ritorno per il calcolo di un nuovo pianeta
PER LO STESSO GIORNO E LA STESSA ORA.
NOTE TECNICHE
1)-i calcolatori tascabili ignorano,in generale i gradi,tutte le costanti fi-
gurano nel programma e, in particolare i dati delle linee da 510 a 586 sono
in radianti.
2)-Per Giove, Saturno, Urano e Nettuno, le costanti sono corrette per tener
conto delle perturbazioni a lungo periodo
3)-Le costanti di tutti i pianeti sono caricate da principio, nella parte co-
mune del programma.
4)-Il calcolo di N segue esattamente l'algoritmo indicato a pag. 34 ; non ac-
cetta che i dati gregoriani.(Questo dato e' esatto col mio controllo)
5)-il sotto-programma del calcolo delle coordinate eclittiche di un pianeta
segue rigorosamente l'algoritmo indicato nel capitolo 4 e 5 :
- Calcolo di L piu' M piu' u per la risoluzione dell'equazione di Keplero
- Calcolo di v con la formula 12b
- Calcolo di Omega piu b
- Calcolo di l con la formula 20
- Calcolo delle coordinate cartesiane eliocentriche
- Trasformazione in coordinate cartesiane geocentriche
- Calcolo di l e di b con le formule 24 e 25
(coloro che desiderano i dettagli dei capitoli 4 e 5 e le formule 12b
le formule 20 24 e 25 le chiedano, e la pag. 34: chiedi a I3ZJV)
6)-Il sotto-programma di risoluzione dell'equazione di Keplero si trova in
tre linee, il coefficiente KE che determina il numero delle iterazioni e'
inizializzato prima dell'entrata nel sotto-programma.
7)-La trasformazione delle coordinate eclittiche in coordinate equatoriali
utilizza le formule 1 e 3. La determinazione del quadrante corretto si fa
con la prova del segno del coseno di l.
8)-Un sotto-programma trasforma i radianti in gradi compresi fra zero e 360.
Un altro trasforma i radianti in ore, minuti e secondi.
9)-A causa del movimento rapido della luna, la sua posizione e' ricalcolata
per un'ora prossima della levata o tramonto, prima di calcolare questa in
modo definitivo. Tuttavia, la precisione non e' migliore che fra 5 e 10
minuti.(Io i3zjv direi di meno, cioe' c'e' maggior precisione)
10)-Gli utenti esperti con la programmazione BASIC non si priveranno di adat-
tare il programma presentato per il loro caso :per l'impiego in uno stes-
so luogo, la longitudine e latitudine del luogo saranno introdotte una so-
la volta; si potranno sopprimere i sotto-programmi corrispondenti alle
funzioni inutili, etc..
Abbiamo scelto il BASIC APPLESOFT per i motivi di facilita' delle prove; si
tratta d'altro canto di un linguaggio conosciuto. Gli utenti di COMMODORE
oppure di TRS non avranno alcuna difficolta' per fare gli adattamenti neces-
sari. Ci scusiamo di non poter soddisfare tutti i modi di BASIC attualmente
in uso.
Nota di I3ZJV :HO TRADOTTO IN IBM PC (ho fatto del mio meglio)
(coloro che desiderano le tabelle comparative fra IBM PC-M20-SIRIUS-RAINBOW
100-APPLE II- ZX 81-ZX SPECTRUM-TI 99/4A-VIC 20-CBM 64 me le chiedano)
Alla riga 750 la parentesi chiusa del 37 mostra una lineetta inclinata a de-
sta attaccata alla parentesi (lunga un millimetro, non so cosa sia).
(Il programma termina alla riga 3590).
(Se vedi una cifra con uno spazio vuoto, cos
e' sul libro).
Le righe 750 760 770 780 erano cosi con Applesoft :
750 CV = PEEK(37), CV
770 POKE 37,CV
770 POKE 36, 0
780 CALL -958
Io le ho modificate per il mio IBM PC cosi:
750 CV = PEEK(32530) - 1
760 POKE 32530, CV
770 POKE 32529, 0
780 GOSUB 950
(penso che la 780 non vadi bene, anche se il programma gira, infatti sui sa-
telliti di Saturno va bene, ma su quelli di JU Giove non va).
Alla riga 440 ho tolto NEXT K: perche' mi risultano in piu', o no ?
Alla riga 530 il numero 5.8338 e' distanziato cosi : 5 . 8338
Alla riga 540 il numero .048376 e' distanziato cosi : .04 8376
Alla riga 3520 mi vien fuori : 3520 "DIAMETRE APPARENT="; DA; "";""
ALLA riga 3560 il numero 80 va sostituito con B0
Alle righe 1770 1780 1790 1800 sostituire VTAB con LOCATE non e' esatto !
(ma che cosa ci si deve mettere ?)
I numeri da 10 fino a 3590 devono essere su ogni riga, non distanziati !
La riga 1762 mi risulta difettosa sui satelliti di JU
Alle tre righe 750 760 770 io ho fatto un programmino per trovarle, eccolo:
(per trovare il numero 32528)
Conclusione di I3zjv: il programma corre, a me mi serve per vedere dov'era
la luna nell'attimo dei terremoti, perche' questa bellina e'concausa di essi.
ma mi serve anche per gli astri !
Per te che vuoi sapere a che ora c'e' l'alta e bassa marea, potrai fare un
calcolo cos
fatto: Devi mettere le coordimate del luogo che calcoli per la
bassa o alta marea. Devi tenere conto che in circa un giorno ci sono due pe-
riodi di flusso e due di riflusso. Il periodo di questo movimento e' di mezza
giornata lunare, pari a 12 ore e 25 minuti. Nel novilunio e nel plenilunio
si hanno le maree piu' ampie:il presente programma ti da' il novilunio e
plenilunio: nel primo e ultimo quarto le maree piu' strette.Tenendo conto
del periodo 6 ore 12 minuti 30 secondi, per una sola volta osservi in quel
dato qth la bassa marea, quando capita, segnandoti quando cade o va posto
questo periodo di 6 ore 12 minuti 30 secondi, rispetto alle ore del primo,
o ultimo quarto per le basse maree, e invece per le alte osservi il ple-
nilunio o novilunio.Trovato come si "installano" queste 6 ore 12 minuti e 30
secondi, cioe' se sono a meta' un po' meno o di piu', allora tutte le altre
maree le potrai calcolare con esattezza. Questo e' un metodo empirico, ma con
un po' di perseveranza si potra' pervenire al calcolo esatto per quella data
localita'.(segnarsi i dati della luna-sole e osservare i ripetitivi)
Nel libro "cielo 2001" edizioni Drioli Editore, piazza Concordia 7, i-22030
Civiglio (Como), lire 36000+sped., telef. e fax 031-364049-
e-mail info@drioli.it
web www.drioli.it
troverai per ogni giorno dell'anno:
esempio :primo quarto ore 11 minuti 8 del tal giorno, avuti questi dati, che
vanno bene per 15 gradi di longitudine e 42 gradi di latitudine, tu ci metti
la latitudine e longitudine del qth che desideri, cos
avrai con esattezza,
l'inizio del primo quarto per quel dato qth.(sul libro "cielo 2001 ler ore
sono ore locali italiane, invece sul programma le ore sono in utc. La diffe-
renza e' di un'ora, cioe' all'ora sul libro diminuisci di un'ora per avere
l'ora utc).Il libro "cielo 2001 e' corredato anche con un cartoncino "ca-
lendario lunare 2001, dove hai i giorni, ma non le ore. Per queste le ore te
le devi trovare come sopradetto.
Per avere il librettino delle maree (Adriatico settentrionale) telefono :
041-5701438
(il programma termina alla riga 3590 RETURN)
(fatti una copia delle istruzioni, poi le levi dal programma, vedi quindi
se il programma gira)-Buon lavoro, 73 Carlo.
========================================================================
==
100 B"BASIC"
110 P=VARPTR(B
120 PRINT P : REM DESCRITTORE DI B
130 L=PEEK(P) : REM LUNGHEZZA B
140 PRINT L
150 BB=PEEK(P+1) : REM BYTE BASSO
160 PRINT BB
170 BA=PEEK(P+2) : REM BYTE ALTO
180 PRINT BA
190 :
200 ID=BA*256 + BB : REM INDIRIZZO B
210 PRINT ID
220 :
230 FOR K=IND TO IND+L-1
240 PRINT PEEK(K), CHRPEEK)(K))
240 NEXT K
RUN
Mi e' venuto fuori :
32528
5
0
180
224
116
255
255
255
Ho utilizzato 32528 aumentandolo di due 32530 32529 a similitudine di 37 36
dove c'era 35, ma credo che sia sbagliato. Ognuno ha un computer dove varia
la cifra.
========================================================================
====
========================================================================
====
PROGRAMMA PER IBM PC MODIFICATO DA I3ZJV: (ERA IN BASIC APPLESOFT)
(le parole sono in lingua francese)
10 CLS : LOCATE 13: PRINT "EPHEMERIDES"
20 LOCATE 13: PRINT "==========="
30 PRINT
40 PRINT "COORDONNEES DU LIEU D'OBSERVATION"
50 PRINT "EN DEGRES et DECIMALES"
60 PRINT "LONGITUDE NEGATIVE A L'EST DE GREENWICH": PRINT
70 PI = 3.14159265
80 INPUT "LATITUDE? ";LA:LA = LA * PI / 180
90 INPUT "LONGITUDE? ";LO
100 PRINT : INPUT "DATE JJ,MM,AAAA ";J,M,A
110 PRINT : INPUT "HEURE T.U.(HH,MM)";HE,MI:HS = (HE + MI / 60) / 24: PRINT
120 J = J + HS
130 N = A * 365 + 31 * (M - 1) + J
140 IF M > 2 GOTO 160
150 A = A - 1
160 N = N + INT (A / 4) - INT (A / 100) + INT (A / 400)
170 IF M < = 2 GOTO 190
180 N = N - INT ((M - 1) * .4 + 2.7)
190 N = N - 694325
200 PRINT "N=";N;
205 DIM L0(15): DIM LP(15): DIM A(15): DIM P0(7): DIM PP(7): DIM O0(7):
DIM OP(7): DIM E(7): DIM I(7)
210 DIM Z6):Z5) = "SAMEDI":Z6) = "DIMANCHE"
220 Z0) = "LUNDI":Z1) = MARDI":Z2) = "MERCREDI"
230 Z3) = "JEUDI":Z4) = "VENDREDI"
240 I = INT ((N / 7 - INT (N / 7)) * 7 + .005)
250 REM
260 PRINT SPC( 4)ZI): PRINT
270 REM*CALCUL DE POSITION*SOLEIL
280 DATA 4.8689,1.72027914E-2,4.9085,8.1856E-7,.01675104,1.00000023,3
290 READ L0,LP,P0,PP,E,A,KE
300 P = P0 + PP * N:LL - L0 + LP * N:M = LL = P:MS = M
310 GOSUB 420
320 V = 2 * ATN ( TAN (U / 2) * SQR ((1 + E) / (1 - E)))
330 R = A * (1 - E * cos (U)):L = V + P:LS = L
340 XS = R * COS (L):YS = R * SIN (L)
350 GOSUB 470
360 PRINT "LONGITUDE SOLEIL=";LD: PRINT
370 GOSUB 1170: GOSUB 1270
380 GOSUB 1080
390 AH = TS - AR
400 ZX = 90.85 * PI / 180: GOSUB 1510: GOSUB 1550
410 GOTO 510
420 U = M: REM EQUATION DE KEPLER
430 FOR K = 0 TO KE
440 U = M + E * SIN (U): RETURN
450 NEXT K
460 RETURN
470 LD = L * 180 / PI: REM TRANSFORMATION RADIAN-DEGRE
480 LD = (LD / 360 - INT (LD / 360)) * 360
490 IF LD < 0 THEN LD = LD + 360
500 LD = INT (LD * 10 + .5) / 10: RETURN
510 DATA 4.0117,7.14254534E-2,1.3249,7.4229E-7,.82304,5.6618E-7,.205615,
.1222,.387098
520 DATA 3.6086,2.79631195E-2,2.2716,6.5572E-7,1.3229,4.3668E-7,.006816,
.05923,.72333
530 DATA 2.1776,9.14676584E-3,5 . 8338,8.793E-7,.8516,3.712E-7,.093309,
3.2294E-2,1.52367 8
540 DATA 4.6879,1.4509868E-3,.2 289,857E-9,1.7358,483E-9,.04 8376,.02284,
5.202799
550 DATA 4.8567,5.8484028E-4,1.5 974,412E-9,1.9686,417E-9,.05 4311,435E-4,
9.552098
560 DATA 4.3224,205424E-9,2.9523,762E-9,1.2825,2.3824E-7,.047319,1.3482E-2,
19.21694
570 DATA 1.5223,105061E-9,.7637,393E-9,2.281,525E-9,.008262,3.1054E-2,
30.112912
580 DATA 1.6406,701214E-10,3.8978,6.672E-7,1.9034,6.672E-7,.250236,.29968,
39.438712
585 DATA 2.694,3.5515529,5.9059,5.771,1.76932338,9.3969,2.597,.87820859,
14.989,3.091,.37648689,26.3641
586 DATA .434,3.32830711,4.877,4.472,2.29571836,6.2464,6.081,1.39085442,
8.723,3.857,.39404325,20.213
590 FOR J = 0 TO 7
600 READ L0(J),LP(J),P0(J),PP(J),O0(J),OP(J),E(J),I(J),A(J): NEXT J
605 FOR I = 8 TO 15: READ L0(I),LP(I),A(I): NEXT
610 INPUT "INTRODUIRE LES DEUX PREMIERES LETTRES DE L'ASTRE CHERCHE ";P
620 IF P= "ME" GOTO 800
630 IF P= "VE" GOTO 810
640 IF P= "MA" GOTO 820
650 IF P= "JU" GOTO 830
660 IF P= "SA" GOTO 840
670 IF P= "UR" GOTO 850
680 IF P= "NE" GOTO 860
690 IF P= "PL" GOTO 870
700 IF P= "LU" GOTO 2000
720 IF P= "ST" THEN END
730 IF P= "ET" GOTO 2500
740 INPUT "INCONNU RECOMMENCEZ";P
750 CV = PEEK (32530) - 1
760 POKE 32530,CV
770 POKE 32529,0
780 GOSUB 950
790 GOTO 620
800 T= "MERCURE=":J = 0:KE = 5:D0 = 6.7:M0 = - .21:A0 = 3.8:B0 = - 3.4:C0 =
2: GOTO 880
810 T= "VENUS=":J = 1:KE = 3:D0 = 16.7:M0 = - 4.1:A0 = .1:B0 = 2.4:C0 =
- .65: GOTO 880
820 T= "MARS"=":J = 2:KE = 5:D0 = 9.4:M0 = - 1.36:A0 = 1.5: GOTO 880
830 T= "JUPITER=":J = 3:KE = 4:FG = 1:D0 = 196.9:M0 = - 9:A0 = 1.48:
GOTO 880
840 T= "SATURNE=":J = 4:KE = 4:FG = 2:D0 = 165.5:M0 = - 8.7:A0 = 1.7:
GOTO 880
850 T= "URANUS=":J = 5:KE = 4:D0 = 71.4:M0 = - 7: GOTO 880
860 T= "NEPTUNE=":J = 6:KE = 3:D0 = 68.3:M0 = - 7: GOTO 880
870 T= "PLUTON=":J = 7:KE = 7:D0 = 8.3:M0 = - 1.5
880 P = P0(J) + PP(J) * N:M = L0(J) + LP(J) * N - P
890 E = E(J): GOSUB 420
900 V = 2 * ATN ( TAN (U / 2) * SQR ((1 + E) / (1 - E)))
910 O = O0(J) + OP(J) * N:C = V + P - O
920 IF COS (C) = 0 THEN D = C: GOTO 950
930 D = ATN ( TAN (C) * COS (I(J)))
940 IF COS (C) < 0 THEN D = D + PI
950 LH = D + O
960 BS = ATN ( SIN (D) * TAN (I(J)))
970 RS = A(J) * (1 - E * COS (U))
980 XP = RS * COS (BS) * COS (LH) + XS:YP = RS * COS (BS) * SIN (LH) + YS:ZP
= RS * SIN (BS)
990 R = SQR (XP * XP + YP * YP):B = ATN (ZP / R):L = ATN (YP / XP)
1000 IF XP < 0 THEN L = L + PI
1010 GOSUB 470
1020 HOME : PRINT TAB( 17)T
1030 PRINT : PRINT "POSITION A ";HE;"H";MI;"TU"
1040 PRINT : PRINT "LONGITUDE ";LD;
1050 PRINT SPC( 2)"LATITUDE "; INT (B * 18000 / PI * .5) / 100
1060 PRINT
1070 GOSUB 1080: GOTO 1130
1080 REM *CALCUL *DU TEMPS SIDERAL
1090 RD = 1.7273 + 1.72027914E - 2 * N + HS * 2 * PI - LO * PI / 180
1100 GOSUB 1300
1110 T= RD
1120 RETURN
1130 GOSUB 1170: GOSUB 1270
1140 GOSUB 1360: GOSUB 3510
1150 GOSUB 1510: GOSUB 1550
1160 GET S IF S= "F" THEN END
1165 GOTO 610
1170 REM*CALCUL DE AR*ET DECLINAISON
1180 EP = .40927971
1190 SD = COS (EP) * SIN (B) + SIN (EP) * COS (B) * SIN (L)
1200 DE = ATN (SD) / SQR (1 - SD * SD))
1210 SR = - SIN (EP) * SIN (B) + COS (EP) * COS (B) * SIN (L)
1220 RD = ATN (SR / COS (B) / COS (L))
1230 IF COS (B) * COS (L) < 0 THEN RD = RD + PI
1240 IF RD < 0 THEN RD = RD + 2 * PI
1250 AR = RD: GOSUB 1300
1260 GOSUB 1320: RETURN
1270 PRINT "ASCENSION DROITE=";H;"H";M;"M";S;"S";
1280 PRINT SPC( 4)"DECLINAISON="; INT (DE * 18000 / PI + .5) / 100: PRINT
1290 RETURN
1300 RD = (RD / 2 / PI - INT ( RD / 2 / PI)) * 2 * PI
1305 IF RD < 0 THEN RD = RD + 2 * PI
1310 RETURN
1320 H = INT (RD / PI * 12)
1330 M = INT ((RD - H * PI / 12) * 720 / PI)
1340 S = INT ((RD - H * PI / 12 - M * PI / 720) * 43200 / PI)
1350 RETURN
1360 REM*CALCUL DE HAUTEUR*ET AZIMUT
1370 AH = TS - AR
1380 ZC = SIN (LA) * SIN (DE) + COS (LA) * COS (DE) * COS (AH)
1390 HT = ATN (ZC / SQR ( - ZC * ZC + 1))
1400 AS = COS (DE) * SIN (AH) / COS (HT)
1410 AC = ( - COS (LA) * SIN (DE) + SIN (LA) * COS (DE) * COS (AH)) / COS (HT)
1420 AZ = ATN (AS / AC)
1430 IF AC < 0 THEN AZ = AZ + PI
1440 LJ = L:L = AZ: GOSUB 470
1450 PRINT : PRINT "AZIMUT=";LD;
1460 HD = HT * 180 / PI:HD = INT (HD * 10 + .5) / 10
1470 PRINT SPC( 3)"HAUTEUR=";HD
1480 IF HD < 0 THEN PRINT "INVISIBLE A ";HE;"H";MI;"(TU)": GOTO 1500
1490 IF HD < 10 THEN PRINT "DIFFICILEMENT VISIBLE A L'HORIZON"
1500 RETURN
1510 REM*CALCUL*DU LEVER & COUCHER
1520 ZX = 90.58 * PI / 180
1530 HC = ( COS (ZX) - SIN (LA) * SIN (DE)) / COS (LA) / COS (DE)
1535 IF HC > = 1 THEN PRINT "CALCUL LEVER-COUCHER IMPOSSIBLE"
1540 HL = - ATN (HC / SQR ( - HC * HC + 1)) + PI / 2: RETURN
1550 HZ = 2 * PI * HS - AH:RD = HZ - HL
1560 GOSUB 1300: GOSUB 1320
1570 PRINT : PRINT "LEVER=";H;"H";M;"M";
1580 RD = HZ + HL
1590 GOSUB 1300: GOSUB 1320
1600 PRINT SPC( 4)"COUCHER=";H;"H";M;"M"
1610 IF FG = 1 OR FG = 2 GOTO 1650
1620 PRINT : RETURN
1650 PRINT : INPUT "SATELLITES (O OU N)? ";S
1660 IF S= "N" THEN FG = 0: PRINT: RETURN
1665 IF S< > "O" GOTO 1650
1690 IF FG = 2 GOTO 1720
1700 FOR I = 8 TO 11
1705 SS(I - 8) = L0(I) + LP(I) * (N - R / 173) - LJ
1710 EL(I - 8) = A(I) * SIN (SS(I - 8))
1715 NEXT
1716 GOTO 1750
1720 FOR I = 12 TO 15
1725 SS(I - 8) = L0(I) * LP(I) * (N - R / 173) - LJ
1730 EL(I - 8) = A(I) * SIN (SS(I - 8))
1735 NEXT
1740 CLS : LOCATE 15: GOTO 1821
1750 CLS : LOCATE 15
1761 FOR I = 0 TO 3
1762 LOCATE EL(I) * .72 + 20
1763 PRINT ".";
1764 NEXT
1769 LOCATE 20: PRINT "O"
1770 LOCATE 21: PRINT TAB( 20 + EL(0) * .8); "IO"
1780 LOCATE 22: PRINT TAB( 19 + EL(1) * .8); "EUROPE"
1790 LOCATE 23: PRINT TAB( 18 + EL(2) * .8); "GANYMEDE"
1800 LOCATE 24: PRINT TAB( 20 + EL(3) * .7); "CALLISTO"
1810 FG = 0: INPUT S IF S= "F" THEN END
1820 TEXT : GOTO 610
1821 PRINT TAB( 20)"O";
1822 FOR I = 4 TO 7
1823 LOCATE EL(I) * .72 + 20
1824 PRINT ".";
1825 NEXT
1826 LOCATE 20: PRINT "O"
1830 LOCATE 21: PRINT TAB( 18 + EL(4) * .8); "TETHYS"
1840 LOCATE 22: PRINT TAB( 18 + EL(5) * .8); "DIONE"
1850 LOCATE 23: PRINT TAB( 19 + EL(6) * .8); "RHEA"
1860 LOCATE 24: PRINT TAB( 18 + EL(7) * .8); "TITAN"
1870 FG = 0: GET S IF S= "F" THEN END
1880 SCREEN : GOTO 610
2000 REM* CALCUL LUNE *
2010 GOSUB 3000
2020 CLS : PRINT TAB( 19)"LUNE": PRINT TAB( 19)"----": PRINT
2030 PRINT "POSITION A ";HE;"H";MI;"(TU)": PRINT
2040 GOSUB 470
2050 PRINT "LONGITUDE=";LD;
2060 PRINT SPC( 4)"LATITUDE="; INT (B * 18000 / PI + .5) / 100: PRINT
2070 GOSUB 1180: GOSUB 1270
2080 AH = TS - AR
2090 REM*CALCUL DE LA PHASE*
2100 PH = L - LS:RD = PH: GOSUB 1300:PH = RD
2110 IF PH < PI / 24 GOTO 2190
2120 IF PH < 11 * PI / 24 THEN PRINT "ENTRE NL & PREMIER QUARTIER": GOTO 2200
2130 IF PH < 13 * PI / 24 THEN PRINT "PREMIER QUARTIER": GOTO 2200
2140 IF PH < 23 * PI / 24 THEN PRINT "ENTRE PQ & PLEINE LUNE": GOTO 2200
2150 IF PH < 25 * PI / 24 THEN PRINT "PLEINE LUNE": GOTO 2200
2160 IF PH < 35 * PI / 24 THEN PRINT "ENTRE PL ET DERNIER QUARTIER": GOTO 2200
2170 IF PH < 37 * PI / 24 THEN PRINT "DERNIER QUARTIER": GOTO 2200
2180 IF PH < 47 * PI / 24 THEN PRINT "ENTRE DERNIER QUARTIER ET NOUVELLE LUNE":
GOTO 2200
2190 PRINT "NOUVELLE LUNE"
2200 PRINT
2210 GOSUB 1520
2220 REM*FRACTION DE JOUR * AU LEVER & COUCHER
2230 HZ = 2 * PI * HS - AH - HL
2240 FJ = HZ / 2 / PI - INT (HZ / 2 / PI)
2245 HZ = 2 * PI * HS - AH + HL
2250 FK = HZ / 2 / PI - INT (HZ / 2 / PI)
2260 NI = N:HU = HS: REM MISE EN RESERVE
2270 N = INT (N) + FJ:HS = FJ
2280 GOSUB 3000: GOSUB 1090: GOSUB 1180
2290 GOSUB 2440
2300 ZX = 89.9 * PI / 180
2310 GOSUB 1530
2320 HZ = 2 * PI * FJ - AH - HL
2325 IF HZ < 0 or HZ > 2 * PI GOTO 2345
2330 RD = HZ: GOSUB 1300: GOSUB 1320
2340 PRINT : PRINT "LEVER=";H;"H";M;"M": GOTO 2350
2345 PRINT : PRINT "PAS DE LEVER"
2350 N = INT (NI) + FK:HS = FK
2360 GOSUB 3000: GOSUB 1090: GOSUB 1180
2370 GOSUB 2440
2380 GOSUB 1530
2390 HZ = 2 * PI * FK - AH + HL
2395 IF HZ < 0 OR HZ > 2 * PI GOTO 2415
2400 RD = HZ: GOSUB 1300: GOSUB 1320
2410 PRINT "COUCHER=";H;"H";M;"M": GOTO 2420
2415 PRINT : PRINT "PAS DE COUCHER"
2420 N = NI:HS = HU: REM RESTAURATION
2430 GOTO 1160
2440 AH = TS - AR
2445 IF AH > PI THEN AH = AH - 2 * PI
2450 IF AH < - PI THEN AH = AH + 2 * PI
2455 RETURN
2500 HOME : INPUT "ASC DROITE? ";AD
2510 PRINT : INPUT "DECLINAISON?";DC
2520 AR = AD * PI / 180
2530 DE = DC * PI / 180
2540 GOSUB 1080: GOSUB 1360
2550 GOTO 1160
3000 REM*LONGITUDE LUNE*
3010 L = .57999 + .229971503 * N
3020 O = 4.1867 - 9.2422E - 4 * N
3030 M = .3193 + .228027135 * N
3040 D = L - LL:F = L - O
3050 L = L + .1098 * SIN (M) + 3.728E - 3 * SIN (2 * M) + 1.149E - 2 *
SIN (2 * D)
3060 L = L - 3.329E - 3 * SIN (MS) + 2.224E - 2 * SIN (2 * D - M) -
1.945E - 3 * SIN ( 2 * F)
3070 L = L + 1.026E - 3 * SIN (2 * D - 2 * M) + 9.983E - 4 *
SIN (2 * D - M - MS)
3080 L = L + 9.903E - 4 * SIN (2 * D + M) + 8.011E - 4 * SIN (2 * D - MS)
3090 L = L + 7.156E - 4 * SIN (M - MS) - 5.323E - 4 * SIN (M + MS) -
6.074E - 4 * SIN (D)
3200 REM*LATITUDE LUNE*
3210 B = .08995 * SIN (F) + 49E- 4 * SIN (M + F) + 485E - 5 * SIN (M - F) +
.003 * SIN (2 * D - F)
3220 B = B + 97E - 5 * SIN (2 * D - M + F) + 8E - 4 * SIN (2 * D - M - F) +
57E - 5 * SIN (2 * D + F)
3230 RETURN
3500 REM* CALCUL DE DIAMETRE*ET MAGNITUDE
3510 DA = D0 / R:DA = INT (DA * 10) / 10
3520 PRINT "DIAMETRE APPARENT=";DA;"";""
3530 L = LH: GOSUB 470:LH = LD
3540 L = LJ: GOSUB 470:LK = LD
3550 PH = ABS (LK - LH): IF PH > 180 THEN PH = 360 - PH
3560 MG = M0 + 5 * LOG (RS * R) * .4343 + A0 * PH / 100 + B0 * PH *
PH / 10000 + C0 * PH * PH * PH / 1000000
3570 MG = INT (MG * 10 + .5) / 10:A0 = 0:B0 = 0:C0 = 0
3580 PRINT "MAGNITUDE=";MG
3590 RETURN
73 a tutti de Luigi.
li, mar 19 maggio 2009 21:01 LT (+2.00 UTC)
I3XTY Paese, Veneto
BBS I3XTY@I3XTY.TV.IVEN.ITA.EU
e_mail
http://
Messaggio inviato con Sally V4.1.232 software packet per Windows 32bit.
Message sent with Sally V4.1. Windows packet software.
Read previous mail | Read next mail
| |
