10 CLEAR 20 CLS 30 PRINT TAB(21); "************ EFEMERIDES ***********" 40 PRINT 50 PRINT " ___________________________________" 60 PRINT 70 PRINT " 2019 INFORTECNICA " 80 PRINT " LEPANT 422 " 90 PRINT " BARCELONA 08025 " 95 PRINT " Tel.93 433 42 43. " 100 PRINT " www.infortecnica.com" 110 PRINT " ___________________________________ " 115 PRINT 116 PRINT " SOFTWARE SERIE METEORITOS " 117 PRINT " USO PUBLICO PARA UN MEJOR CONOCIMIENTO DE LOS RADIANTES" 120 PRINT 130 PRINT "ESTE PROGRAMA DETERMINA EL RADIANTE DE UN METEOROIDE " 132 PRINT "A PARTIR DE LAS VELOCIDADES DE IMPACTO VX, VY VZ " 150 PRINT 160 TMOPT = 1 165 REM GOSUB 20000 170 INPUT " HORA DEL CALCULO EN T.U. DE O A 24.00-HORA.MIN) "; TOD 180 INPUT " DIA "; DAYB 190 INPUT " MES "; MONB 200 INPUT " AÑO "; YRB 210 INPUT " INTRODUCIR LATITUD DEL LUGAR (XX.XX N=+,S=-) "; LAT 220 INPUT " INTRODUCIR LONGITUD EN GRADOS Y FRACCION (XXX.XX O=+,E=-) "; LONGG 230 CLS 240 PRINT "Efemérides astronómicas calculadas para las "; TOD; " horas del día "; DAYB; " de "; MONB 250 PRINT "Los datos topocéntricos son los relativos a las coordenadas:" 260 PRINT "Longitud ="; LONG.; " grados. Latitud ="; LAT; " grados." 270 CP = 1 280 REM *** CALCULO DEL TIEMPO *** 285 TOD = TOD + RETH 290 TMNM$(1) = "GMT" 300 DAYL = DAYB'DIA LOCAL 310 IF TMOPT = 1 THEN DAYB = DAYB + INT(TOD) / 24 + (TOD - INT(TOD)) * 100 / 1440: GMT = TOD 320 IF GMT > 24 THEN GMT = GMT - 24 330 REM ******* CALCULO DEL DIA JULIANO ********* 340 IF MONB > 2 THEN YR = YRB: MON = MONB 350 IF MONB <= 2 THEN YR = YRB - 1: MON = MONB + 12 360 AA# = INT(YR / 100) 370 BB# = 2 - AA# + INT(AA# / 4) 380 JD# = INT(365.25 * YR) + INT(30.6001 * (MON + 1)) + DAYB + 1720994.5# + BB# 390 JDS# = INT(365.25 * YR) + INT(30.6001 * (MON + 1)) + INT(DAYB) + 1720994.5# + BB# 400 REM ******* CALCULO DE T EN J C HASTA 0.5 TE ENE 1900 ***** 410 T# = (JD# - 2415020#) / 36525! 420 TS# = (JDS# - 2415020#) / 36525! 430 REM *** CALCULO DEL TIEMPO SIDEREO *** 440 SIDT# = .276919398# + 100.0021359# * TS# + 1.075E-06 * TS# * TS# 450 SIDT# = (SIDT# - INT(SIDT#)) * 24 460 SIDT# = SIDT# + (DAYB - INT(DAYB)) * 24 * 1.002737908# 470 IF SIDT# > 24 THEN SIDT# = SIDT# - 24 480 LCST = SIDT# + (-LONGG / 15) 490 IF LCST > 24 THEN LCST = LCST - 24 500 IF LCST < 0 THEN LCST = LCST + 24 510 REM 520 REM 530 REM *** CALCULOS PARA EL SOL Y LA TIERRA *** 540 REM ******* CALCULO DE LA LONGITUD DEL SOL ********* 550 LSUN# = 279.69668# + 36000.76892# * T# + .0003025 * T# * T# 560 LSUN# = (LSUN# / 360 - INT(LSUN# / 360)) * 360 570 MSUN# = 358.47583# + 35999.04975# * T# - .00015 * T# * T# - .0000033 * T# * T# * T# 580 MSUN# = (MSUN# / 360 - INT(MSUN# / 360)) * 360 590 RCN = 57.2957795# ' DE RADIANES A GRADOS 600 YE = (TAN(23.442136# / 2 / RCN)) ^ 2 610 ES = 1.671834E-02 620 EQT = (YE * SIN(2 * LSUN# / RCN) - 2 * ES * SIN(MSUN# / RCN) + 4 * ES * YE * SIN(MSUN# / RCN) * COS(2 * LSUN# / RCN) - .5 * YE * YE * SIN(4 * LSUN# / RCN) - 1.25 * ES * ES * SIN(2 * MSUN#)) * RCN / 15 * 60 630 ECEARTH# = 1.675104E-02 - .0000418 * T# - 1.26E-07 * T# * T# 640 CSUN# = (1.91946 - .004789 * T# - .000014 * T# * T#) * SIN(MSUN# / RCN) + (.020094 - .0001 * T#) * SIN(2 * MSUN# / RCN) + .000293 * SIN(3 * MSUN# / RCN) 650 LSUNT# = LSUN# + CSUN# 660 REM ******* CALCULO DE LA LONGITUD DE LA TIERRA ******** 670 IF LSUNT# > 180 THEN LEARTH# = LSUNT# - 180 ELSE LEARTH# = LSUNT# + 180 680 NUSUN# = MSUN# + CSUN#'SUN'S TRUE ANOMALY 690 REM *** CALCULATE R.A. AND DEC. FOR SUN *** 700 TOP = SIN(LSUNT# / RCN) * COS(23.442136# / RCN) 710 BOT = COS(LSUNT# / RCN) 720 RASUN = ATN(TOP / BOT) * RCN 730 IF TOP > 0 AND BOT > 0 THEN RASUN = RASUN 740 IF TOP > 0 AND BOT < 0 THEN RASUN = RASUN + 180 750 IF TOP < 0 AND BOT > 0 THEN RASUN = RASUN + 360 760 IF TOP < 0 AND BOT < 0 THEN RASUN = RASUN + 180 770 RASUN = RASUN / 15 780 DX = SIN(23.442136# / RCN) * SIN(LSUNT# / RCN) 790 DECSUN = ATN(DX / SQR(1 - DX * DX)) * RCN 800 MINSUN = (RASUN - INT(RASUN)) * 60 810 DMINSUN = (DECSUN - FIX(DECSUN)) * 60 820 REM *** CALCULO DEL AZIMUT Y ALTURA DEL SOL *** 830 HSUN = (LCST - RASUN) * 15 840 AF = SIN(DECSUN / RCN) * SIN(LAT / RCN) + COS(DECSUN / RCN) * COS(LAT / RCN) * COS(HSUN / RCN) 850 ALTSUN = ATN(AF / SQR(1 - AF * AF)) * RCN 860 AZF = (SIN(DECSUN / RCN) - SIN(LAT / RCN) * SIN(ALTSUN / RCN)) / COS(LAT / RCN) / COS(ALTSUN / RCN) 870 AZSUN = (1.570796 - ATN(AZF / SQR(1 - AZF * AZF))) * RCN 880 IF SIN(HSUN / RCN) > 0 THEN AZSUN = 360 - AZSUN 890 TXSUN = .533128 * ((1 + ECEARTH# * COS(NUSUN# / RCN)) / (1 - ECEARTH# * ECEARTH#)) 900 REM *** CALCULO DE LA DISTANCIA TIERRA/SOL *** 910 REARTH# = 1.0000002# * (1 - ECEARTH# * ECEARTH#) / (1 + ECEARTH# * COS(NUSUN# / RCN)) 920 REM *** CALCULO DE LA SALIDA Y PUESTA DEL SOL *** 930 LET QT = -TAN(LAT / RCN) * TAN(DECSUN / RCN) 940 TRS = 1 / 15 * (1.570796 - ATN(QT / SQR(1 - QT * QT))) * RCN 950 GSTR = 24 + RASUN - TRS 960 IF GSTR > 24 THEN GSTR = GSTR - 24 970 GSTS = RASUN + TRS 980 IF GSTS > 24 THEN GSTS = GSTS - 24 990 GMTF = INT(GMT) + (GMT - INT(GMT)) / .6 - SIDT# 1000 GMTR = GSTR + GMTF 1010 IF GMTR < 0 THEN GMTR = GMTR + 24 1020 IF GMTR > 24 THEN GMTR = GMTR - 24 1030 GMTS = GSTS + GMTF 1040 IF GMTS < 0 THEN GMTS = GMTS + 24 1050 IF GMTS > 24 THEN GMTS = GMTS - 24 1060 IF TMOPT = 1 THEN RT = GMTR + LONGG / 15 1070 IF TMOPT = 1 THEN ST = GMTS + LONGG / 15 1080 IF TMOPT = 2 THEN RT = GMTR + (LONGG - 75) / 15 1090 IF TMOPT = 2 THEN ST = GMTS + (LONGG - 75) / 15 1100 IF TMOPT = 3 THEN RT = GMTR + (LONGG - 75) / 15 + 1 1110 IF TMOPT = 3 THEN ST = GMTS + (LONGG - 75) / 15 + 1 1120 IF TMOPT = 4 THEN RT = GMTR + (LONGG - 90) / 15 1130 IF TMOPT = 4 THEN ST = GMTS + (LONGG - 90) / 15 1140 IF TMOPT = 5 THEN RT = GMTR + (LONGG - 90) / 15 + 1 1150 IF TMOPT = 5 THEN ST = GMTS + (LONGG - 90) / 15 + 1 1160 IF TMOPT = 6 THEN RT = GMTR + (LONGG - 105) / 15 1170 IF TMOPT = 6 THEN ST = GMTS + (LONGG - 105) / 15 1180 IF TMOPT = 7 THEN RT = GMTR + (LONGG - 105) / 15 + 1 1190 IF TMOPT = 7 THEN ST = GMTS + (LONGG - 105) / 15 + 1 1200 IF TMOPT = 8 THEN RT = GMTR + (LONGG - 120) / 15 1210 IF TMOPT = 8 THEN ST = GMTS + (LONGG - 120) / 15 1220 IF TMOPT = 9 THEN RT = GMTR + (LONGG - 120) / 15 + 1 1230 IF TMOPT = 9 THEN ST = GMTS + (LONGG - 120) / 15 + 1 1240 REM *** CALCULO DE LA DURACION SOLAR *** 1250 TX = (COS(108 / RCN) - SIN(LAT / RCN) * SIN(DECSUN / RCN)) / COS(LAT / RCN) / COS(DECSUN / RCN): IF ABS(TX) > 1 THEN 1310 1260 TX1 = (1.570796 - ATN(TX / SQR(1 - TX * TX))) * RCN 1270 TY = -TAN(LAT / RCN) * TAN(DECSUN / RCN): IF TY > 1 THEN 1310 1280 TY1 = (1.570796 - ATN(TY / SQR(1 - TY * TY))) * RCN 1290 TL = ABS((TY1 - TX1) * 4 / 60) 1300 GOTO 1330 1310 TL = 24 - ST + RT 1320 REM 1330 R1 = RT: S1 = ST: IF CP > 2 THEN 1620 1340 PRINT 1350 PRINT 1360 PRINT "DATOS ASTRONOMICOS PARA"; DAYL; "/"; MONB; "/"; YRB; TOD; TMNM$(TMOPT); " LAT= "; LAT; " LONGG= "; LONGG 1370 PRINT 1380 PRINT " TIEMPO SIDEREO LOCAL = "; INT(LCST); " H "; CINT((LCST - INT(LCST)) * 60); " MIN" 1390 PRINT " ECUACION DEL TIEMPO MINUTOS = "; EQT; " MIN" 1400 PRINT " ASCENSION RECTA DEL SOL = "; INT(RASUN); " H "; CINT(MINSUN); " MIN" 1410 PRINT " DECLINACION DEL SOL = "; FIX(DECSUN); " GRD "; ABS(CINT(DMINSUN)); " MIN" 1420 PRINT " ALTURA DEL SOL = "; ALTSUN 1430 AZSUN = -(180 - AZSUN) 1440 PRINT " ACIMUT DEL SOL = "; AZSUN 1450 PRINT " SALIDA DEL SOL (TIEMPO UNIVERSAL) = "; INT(RT); " H "; INT((RT - INT(RT)) * 60); " MIN "; TMNM$(TMOPT) 1460 PRINT " PUESTA DEL SOL (TIEMPO UNIVERSAL) = "; INT(ST); " H "; INT((ST - INT(ST)) * 60); " MIN "; TMNM$(TMOPT) 1470 REARTH# = REARTH# * 1.496E+08 1480 PRINT 1490 XX = MONB + DAYB / 30 1500 VARA = 229.9794 + 29.6773 * XX - 25.5121 * XX * XX + 3.61772 * XX * XX * XX - .1408888 * XX * XX * XX * XX 1510 VARA = (VARA / 1000 + 1) * 1000 1520 VARB = 144.4 - 19.697 * XX + 13.18 * XX * XX - 1.78086 * XX * XX * XX + .067955 * XX * XX * XX * XX 1530 VARB = VARB / 1000 1540 VARC = 60.7957 - 25.6924 * XX + 16.5257 * XX * XX - 2.2392 * XX * XX * XX + .086189 * XX * XX * XX * XX 1550 VARC = VARC / 1000 1560 IDN = VARA / (2.7182818# ^ (VARB / SIN(ALTSUN / RCN))) 1570 PRINT 1620 REM INPUT "SI NUEVO CALCULO 1 "; NC IF NC = 1 THEN 170 3010 PRINT "SADEYA/INFORTECNICA" PRINT "RADIANTE DE UN METEORO A PARTIR DE VX,VY,VZ" PRINT TSG = LCST * 15 PRINT "TIEMPO SIDEREO GRADOS"; TSG INPUT "TIEMPO SIDEREO LONG 0 GRADOS"; TSIG PRINT "INTROD. VEL MET." 3110 INPUT "VX,VY,VZ"; vx, vy, vz vx = vx * -1 vy = vy * -1 vz = vz * -1 IF vx > 0 AND vy > 0 THEN c = 1 IF vx < 0 AND vy > 0 THEN c = 2 IF vx < 0 AND vy < 0 THEN c = 3 IF vx > 0 AND vy < 0 THEN c = 4 AXY = ATN(vx / vy) VXY = SQR(vx * vx + vy * vy) V = SQR(VXY * VXY + vz * vz) VI = SQR(V * V - 11.2 * 11.2) DEC = ATN(vz / VXY) RD = 180 / 3.141592 ar = AXY * RD IF c = 1 THEN ar = 90 - ar IF c = 2 THEN ar = 180 - (ar + 90) IF c = 3 THEN ar = 270 - ar IF c = 4 THEN ar = 360 - (ar + 90) DEC = DEC * RD ASR = ar + TSIG IF ASR > 360 THEN ASR = ASR - 360 PRINT "ANG.H0R,DEC,VIMP,VINF="; ar, DEC, V, VI PRINT "ASCENSION RECTA="; ASR PRINT "CUADRANTE=", c PRINT INPUT "PARA NUEVO CALCULO 1"; NC IF NC = 1 GOTO 3010 END