1 PRINT "PREPARACION DE BASE DE DATOS -ELEM4.TXT- HOME P. J.WALKER(1)"
2 PRINT "ADAPTACION DE -TP.TXT-(2), ADAPTACION -ELEM4.TXT- SOFT. BISQUE (3)"
3 INPUT "INICIO PROG.(4)"; TTPP
4 IF TTPP = 1 THEN 30000
5 IF TTPP = 2 THEN 40000
6 IF TTPP = 3 THEN 50000
10 KEY OFF
40 CLS
50 PRINT "INFORTECNICA 2002"
60 PRINT "Lepanto 422"
70 PRINT "Barcelona  08025"
80 PRINT "Tel. 934334243 / Fax. 934334239"
81 PRINT "Website - www.sadeya.org"
82 PRINT
83 PRINT "INFORTECNICA SERVICIO DE CALCULO Y SOFTWARE:"
84 PRINT "ESTUDIOS ESTADISTICOS, ASTRONOMIA, ROBOTICA,"
85 PRINT "ESTUDIOS DE MERCADO, INGENIERIA, MECANICA, TEC. MILITAR,"
86 PRINT "ARQUITECTURA, CALEFACCION, ANTENAS PARABOLICAS, MATEMATICA,"
87 PRINT "PROGRAMACION LINEAL, AERONAUTICA, ASTRONAUTICA, GEODESIA, TECNOLOGIA, "
88 PRINT "RESISTENCIA DE MATERIALES, FISICA, METEOROLOGIA, CARTOGRAFIA,"
89 PRINT "Solicite informaci¢n sin compromiso.   (pulse una tecla)"
100 PRINT
110 PRINT
130 OPEN "EFC2000.DAT" FOR APPEND AS 1
140 A$ = INPUT$(1): FOR A = 1 TO 8000: NEXT A
150 CLS
160 DEFDBL A-Z
170 DEFSNG I
180 PRINT
190 nuco = 0
200 PRINT " ***EFEMERIDES ECUATORIALES DE GRUPOS DE COMETAS***"
210 PRINT
212 INPUT "FIGURAR ELEMENTOS SI(1), NO(2) "; ELESI
214 INPUT "ELONGACION MINIMA PARA EFEMERIDES "; ELO15
216 INPUT "DISTANCIA A LA TIERRA MAXIMA PARA EFEMERIDES "; DIMI
220 INPUT "HASTA QUE MAGNITUD LE INTERESA LAS EFEMERIDES "; MMM9
225 INPUT "NUMERO DE COMETAS A ESTUDIAR "; NUCOM
230 PRINT "NUMERO DE EFEMERIDES PARA CADA COMETA";
240 INPUT PPP9
250 PRINT "INTERVALO EN DIAS=";
260 INPUT PPP4
270 P4 = PPP4
280 INPUT "PARA EL TOTAL DEL CIELO PULSE (1). PARA UN SECTOR DETERMINADO (2)"; RIF
290 AMA = 90: AMI = -90: DMA = 100: DMI = -100
300 IF RIF = 1 THEN 360
310 PRINT "PODRA ESTUDIAR SI EXISTEN COMETAS EN UN SECTOR DEL CIELO."
320 INPUT "MARQUE LOS LIMITES DEL SECTOR. ASCENSION RECTA MAXIMA (HH.MM)"; AMA
330 INPUT "ASCENSION RECTA MINIMA (HH.MM)"; AMI
340 INPUT "DECLINACION MAXIMA (GG.MM)"; DMA
350 INPUT "DECLINACION MINIMA (GG.MM)"; DMI
360 PRINT "TIEMPO DE LA EFEMERIDE INICIAL ";
370 CLS
380 GOSUB 620
390 TT9 = J0
400 JJ1 = J0
410 REM
420 REM NUCOM=150
430 PRINT "EFEMERIDES ECUATORIALES EQUIN. 2000,0"
440 PRINT #1, "EFEMERIDES ECUATORIALES EQUIN. 2000,0"
450 PRINT "Fecha     A.R.(2000)Dec.    Delta    r.    Elon.  F.A. Mag"
460 PRINT #1, "Fecha     A.R.(2000)Dec.    Delta    r.    Elon.  F.A. Mag"
470 nuco = nuco + 1: J0 = TT9: BAND = 0
480 PRINT #1,
490 P7 = 0
500 IF nuco > NUCOM THEN GOTO 3940
510 GOSUB 4801
520 PI = 3.141592654#
530 R9 = 180 / PI
540 P7 = 0
550 GOTO 2860
560 W7 = A9 / 360 - INT(A9 / 360)
570 B2 = 9 + W7
580 B3 = B2 - INT(B2)
590 A1 = B3 * 2 * PI
600 RETURN
610 Q7 = 18
620 PRINT "A¥O=";
630 INPUT a5
640 PRINT "MES=";
650 INPUT M5
660 PRINT "DIA=";
670 INPUT d5
680 IF M5 > 2 THEN 730
690 M6 = 12 + M5
700 A6 = a5 - 1
710 a5 = A6
720 M5 = M6
730 A0 = INT(a5 / 100)
740 B0 = 2 - A0 + INT(A0 / 4)
750 J0 = INT(365.25 * a5) + INT(30.6001 * (M5 + 1)) + d5 + 1720994.5#
760 IF J0 > 2299170.50001# THEN 780
770 RETURN
780 J0 = J0 + B0
790 GOTO 770
800 PRINT "JD="; J0
810 CLOSE 1: END
820 PRINT "FECHA JULIANA=";
830 INPUT J0
840 J6 = J0 + .5
850 J5 = INT(J6)
860 A0 = INT((J5 - 1867216.25#) / (36524.25#))
870 A1 = J5 + 1 + A0 - INT(A0 / 4)
880 IF J5 < 2299161.001# THEN 1010
890 B0 = A1 + 1524
900 C0 = INT((B0 - 122.1) / (365.25))
910 D0 = INT(365.25 * C0)
920 E0 = INT((B0 - D0) / (30.6001))
930 d5 = B0 - D0 - INT(30.6001 * E0) + (J6 - J5)
940 IF E0 > 13.5 THEN 970
950 M5 = E0 - 1
960 GOTO 980
970 M5 = E0 - 13
980 IF M5 > 2.5 THEN 1040
990 a5 = C0 - 4715
1000 GOTO 1050
1010 A1 = J5
1020 GOTO 890
1030 GOTO 1030
1040 a5 = C0 - 4716
1050 RETURN
1060 V7 = R
1070 K5 = B
1080 GOSUB 1140
1090 Z = Y
1100 V7 = X
1110 K5 = A
1120 GOSUB 1140
1130 RETURN
1140 R9 = 180 / PI
1150 A8 = (K5 / 360 - INT(K5 / 360)) * 360
1160 C1 = INT(A8 / 90) + 1
1170 A2 = (A8 / 90 - INT(A8 / 90)) * 90
1180 X1 = COS(A2 / R9) * V7
1190 Y1 = SIN(A2 / R9) * V7
1200 IF C1 > 3.99999 THEN 1610
1210 IF C1 > 2.99999 THEN 1640
1220 IF C1 > 1.99999 THEN 1670
1230 GOTO 1700
1240 X1 = ABS(X)
1250 Y1 = ABS(Y)
1260 V7 = SQR(X1 * X1 + Y1 * Y1)
1270 A = ATN(Y / X) * 180 / PI
1280 IF X > -1E-08 THEN 1320
1290 IF X < 0 THEN 1300
1300 K5 = 180 + A
1310 RETURN
1320 K5 = A
1330 RETURN
1340 PRINT "X=";
1350 INPUT X9
1360 PRINT "Y=";
1370 INPUT Y9
1380 PRINT "Z=";
1390 INPUT Z9
1400 GOSUB 1430
1410 PRINT "A"; A3; "B"; B7; "R"; V7
1420 CLOSE 1: END
1430 X = X9
1440 Y = Y9
1450 GOSUB 1240
1460 A3 = K5
1470 X = V7
1480 Y = Z9
1490 GOSUB 1240
1500 B7 = K5
1510 RETURN
1520 PRINT "DE POLARES A CARTESIANAS:A=";
1530 INPUT A
1540 PRINT "B=";
1550 INPUT B
1560 PRINT "R=";
1570 INPUT R
1580 GOSUB 1060
1590 PRINT "X="; X; "Y="; Y; "Z="; Z
1600 CLOSE 1: END
1610 X = Y1
1620 Y = X1 * -1
1630 RETURN
1640 X = X1 * -1
1650 Y = Y1 * -1
1660 RETURN
1670 X = Y1 * -1
1680 Y = X1
1690 RETURN
1700 X = X1
1710 Y = Y1
1720 RETURN
1730 PRINT "DE GRADOS A GMS: GRD=";
1740 INPUT G7
1750 IF G7 < 0 THEN 1850
1760 G = INT(G7)
1770 M6 = (G7 - G) * 60
1780 M = INT(M6)
1790 S = (M6 - M) * 60
1800 G9 = G + M / 100 + S / 10000
1810 RETURN
1820 GOSUB 1730
1830 PRINT "GRD="; G7; "G="; G; "M="; M; "S="; S; "GMS="; G9
1840 CLOSE 1: END
1850 G7 = G7 + 360
1860 GOSUB 1760
1870 G9 = -1 * G9
1880 RETURN
1890 PRINT "G=";
1900 INPUT G
1910 PRINT "M=";
1920 INPUT M
1930 PRINT "S=";
1940 INPUT S
1950 G7 = G + M / 60 + S / 3600
1960 G9 = G + M / 100 + S / 10000
1970 RETURN
1980 PRINT "DE GMS A GRADOS";
1990 GOSUB 1890
2000 PRINT "GMS="; G9; "GRD="; G7
2010 CLOSE 1: END
2020 PRINT "COORDENADAS DEL SOL";
2030 GOSUB 620
2040 GOSUB 2110
2050 R9 = 180 / PI
2060 L0 = L1 * R9
2070 PRINT "1.VER.="; L0; "R="; R6; "AR.APR.="; W8; "DEC.APR.="; D0
2080 CLOSE 1: END
2090 GOSUB 2110
2100 RETURN
2110 T = (J0 - 2415020#) / 36525!
2120 L7 = 270.4344# + 481267.8831# * T - .0001133 * T ^ 2
2130 L6 = 279.69668# + 36000.76892# * T + .0003025 * T ^ 2
2140 A9 = L6
2150 GOSUB 560
2160 L6 = A1
2170 M6 = 358.47583# + 35999.04975# * T - .00015 * T ^ 2
2180 A9 = M6
2190 GOSUB 560
2200 M6 = A1
2210 M7 = 296.1046# + 477198.8491# * T + .009192 * T ^ 2
2220 C5 = 1.91946 - .004789 * T - .000014 * T ^ 2
2230 C6 = PI / 180 * (C5 * SIN(M6) + (.020094 - .0001 * T) * SIN(2 * M6) + .000293 * SIN(3 * M6))
2240 N0 = 259.1833# - 1934.142# * T + .002078 * T ^ 2
2250 L1 = C6 + L6
2260 V6 = M6 + C6
2270 E7 = .01675104# - .0000418 * T
2280 R6 = ((1.0000002# * (1 - E7 ^ 2)) / (1 + E7 * COS(V6)))
2290 N6 = 259.13 - 1934.142# * T
2300 A9 = N6
2310 GOSUB 560
2320 N6 = A1
2330 W = L1 - PI / 180 * (.00569 - .00479 * SIN(N6))
2340 R9 = 180 / PI
2350 N1 = -(17.2327 + .01737 * T) * SIN(N0 / R9) - 1.273 * SIN(2 * L6 / R9)
2360 N2 = 9.21 * COS(N0 / R9) + .552 * COS(2 * L6 / R9)
2370 O1 = 23.452294# - .0130125 * T + .00256 * COS(N6)
2380 O5 = 23.452294# - .0130125 * 1
2390 L5 = L1 - PI / 180 * 1.396 * (T - 1)
2400 X1 = R6 * COS(W)
2410 X5 = R6 * COS(L5)
2420 Y7 = R6 * SIN(W) * COS(O1 / (180 / PI))
2430 Y5 = R6 * SIN(L5) * COS(O5 / (180 / PI))
2440 Z3 = R6 * SIN(W) * SIN(O1 / (180 / PI))
2450 Z5 = R6 * SIN(L5) * SIN(O5 / (180 / PI))
2460 R9 = 180 / PI
2470 V7 = R6
2480 K5 = L5 * R9
2490 GOSUB 1140
2500 X2 = X
2510 Y2 = Y
2520 V7 = R6
2530 K5 = W * R9
2540 GOSUB 1140
2550 X3 = X
2560 Y3 = Y
2570 X9 = X5
2580 Y9 = Y5
2590 Z9 = Z5
2600 GOSUB 1430
2610 W5 = A3
2620 d5 = B7
2630 X9 = X3
2640 Y9 = Y7
2650 Z9 = Z3
2660 GOSUB 1430
2670 W6 = A3
2680 D9 = B7
2690 P1 = 8.793999 / R6
2700 S0 = .267 / R6
2710 G7 = W5 / 15
2720 GOSUB 1750
2730 W0 = G9
2740 G7 = W6 / 15
2750 GOSUB 1750
2760 W8 = G9
2770 G7 = d5
2780 GOSUB 1750
2790 D6 = G9
2800 G7 = D9
2810 GOSUB 1750
2820 D0 = G9
2830 RETURN
2840 D7 = D7 * -1
2850 GOTO 3570
2860 GOSUB 5401
2870 GOSUB 5901
2880 GOSUB 6401
2890 GOSUB 6901
2900 GOSUB 7401
2910 IF E8 = 1 THEN 2930
2920 GOTO 2940
2930 E8 = .999999
2940 S8 = Q8 / (1 - E8)
2950 GOSUB 7901
2960 P9 = PPP9
2970 P4 = PPP4
2980 GOTO 3790
2990 H7 = 0
3000 T = (T9 - T8) - H7
3010 Q = Q8
3020 E0 = E8
3030 GOTO 3950
3040 V1 = V
3050 R1 = R
3060 K5 = V1 * R9 + P8
3070 V7 = R1
3080 GOSUB 1140
3090 G0 = X
3100 H0 = Y
3110 V7 = H0
3120 K5 = I8
3130 GOSUB 1140
3140 H0 = X
3150 I0 = Y
3160 X = G0
3170 Y = H0
3180 GOSUB 1240
3190 K5 = K5 + N8
3200 GOSUB 1140
3210 G0 = X
3220 H0 = Y
3230 G1 = G0 + X2
3240 H1 = H0 + Y2
3250 I1 = I0
3260 X = H1
3270 Y = I1
3280 GOSUB 1240
3290 K5 = K5 + O5
3300 GOSUB 1140
3310 H2 = X
3320 I2 = Y
3330 G2 = G1
3340 X9 = G2
3350 Y9 = H2
3360 Z9 = I2
3370 GOSUB 1430
3380 A4 = A3
3390 D4 = B7
3400 G7 = A4 / 15
3410 GOSUB 1750
3420 W0 = G9
3430 G7 = D4
3440 GOSUB 1750
3450 IF H7 > 0 THEN 3480
3460 H7 = .005775 * V7
3470 GOTO 3000
3480 D1 = G9 - .0001 * COS((A4 + 168.75) / (180 / PI)) * SIN(G9 / (180 / PI))
3490 D7 = INT(ABS(D4)) + .6 * (ABS(D4) - INT(ABS(D4)))
3500 U8 = V7
3510 A4 = A4 - .0001 * SIN((A4 + 168.75) / (180 / PI)) / COS(D1 / (180 / PI))
3520 W9 = A4 / 360 + 9
3530 W9 = 24 * (W9 - INT(W9))
3540 A6 = INT(W9) + .6 * (W9 - INT(W9))
3550 GOSUB 840
3560 IF D4 < 0 THEN 2840
3570 T5 = a5 + M5 / 100 + d5 / 10000
3580 E6 = (R6 * R6 + U8 * U8 - R1 * R1) / (2 * R6 * U8)
3590 E6 = R9 * ATN(SQR(1 - E6 * E6) / E6)
3600 IF E6 > 0 THEN 3620
3610 E6 = 180 + E6
3620 F6 = (R1 * R1 + U8 * U8 - R6 * R6) / (2 * R1 * U8)
3630 F6 = R9 * ATN(SQR(1 - F6 * F6) / F6)
3640 IF F6 > 0 THEN 3660
3650 F6 = 180 + F6
3660 M9 = 5 * LOG(U8) / LOG(10) + kk8 * LOG(R1) / LOG(10) + k8
3670 P7 = 1 + P7
3680 IF P7 > P9 THEN 470
3690 IF E6 > ELO15 AND M9 < MMM9 AND AMI < A6 AND AMA > A6 AND DMI < D7 AND DMA > D7 AND BAND = 0 THEN PRINT NOM$
3700 IF E6 > ELO15 AND M9 < MMM9 AND AMI < A6 AND AMA > A6 AND DMI < D7 AND DMA > D7 AND BAND = 0 THEN PRINT #1, NOM$
3710 IF E6 > ELO15 AND M9 < MMM9 AND AMI < A6 AND AMA > A6 AND DMI < D7 AND DMA > D7 AND BAND = 0 THEN PRINT #1, MOM$
3720 IF E6 > ELO15 AND M9 < MMM9 AND AMI < A6 AND AMA > A6 AND DMI < D7 AND DMA > D7 THEN BAND = 1
3730 MM9 = INT(M9 * 10) / 10
3735 M = 1000
3736 D = 10000
3740 IF E6 > ELO15 AND DIMI > U8 AND M9 < MMM9 AND AMI < A6 AND AMA > A6 AND DMI < D7 AND DMA > D7 THEN PRINT USING "####.#### ##.#### +##.### ###.### ###.### ###.# ###.# ##.#"; T5; INT(A6 * D) / D; INT(D7 * M) / M; INT(U8 * M) / M; INT(R1 * M) / M; _
 INT(E6 * 10) / 10; INT(F6 * 10) / 10; M9
3750 FF6 = INT(F6 * 10) / 10
3760 IF E6 > ELO15 AND DIMI > U8 AND M9 < MMM9 AND AMI < A6 AND AMA > A6 AND DMI < D7 AND DMA > D7 THEN PRINT #1, USING "####.#### ##.#### +##.### ###.### ###.### ###.# ###.# ##.#"; T5; INT(A6 * D) / D; INT(D7 * M) / M; INT(U8 * M) / M; INT(R1 * M)  _
/ M; INT(E6 * 10) / 10; F6; M9
3770 J0 = P4 + J0
3780 GOTO 3900
3790 GOSUB 8201
3800 GOSUB 5225
3810 T8 = JJ0
3820 T8 = JJ0
3830 IF ELESI = 2 THEN 3900
3835 PRINT
3836 PRINT #1,
3840 PRINT "EFEMERIDES ECUATORIALES DEL COMETA "; NOM$; " EQUIN. 2000,0"
3850 PRINT #1, "EFEMERIDES ECUATORIALES DEL COMETA "; NOM$; " EQUIN. 2000,0"
3860 PRINT "T="; T8; " PERI="; P8; " NODO="; N8; " INCL="; I8; " EXC="; E8; " Q="; Q8; " H="; k8
3870 PRINT #1, "T="; T8; " PERI="; P8; " NODO="; N8; " INCL="; I8; " EXC="; E8; " Q="; Q8; " G="; k8
3880 PRINT "Fecha     A.R.(2000)Dec.    Delta    r.    Elon.  F.A. Mag"
3890 PRINT #1, "Fecha     A.R.(2000)Dec.    Delta    r.    Elon.  F.A. Mag"
3900 J1 = J0
3910 T9 = J0
3920 GOSUB 2110
3930 GOTO 2990
3940 CLOSE 1: END
3950 K = .01720209895#
3960 IF E0 > .95 THEN 4250
3970 IF E0 > .9999999899999999# THEN 4700
3980 A1 = Q / (1 - E0)
3990 M = K * T * A1 ^ (-1.5)
4000 F = SGN(M)
4010 M = ABS(M) / (2 * PI)
4020 M = (M - INT(M)) * 2 * PI * F
4030 IF M < 0 THEN 4050
4040 GOTO 4060
4050 M = M + 2 * PI
4060 F = 1
4070 IF M > PI THEN 4090
4080 GOTO 4100
4090 F = -1
4100 IF M > PI THEN 4120
4110 GOTO 4130
4120 M = 2 * PI - M
4130 E = PI / 2
4140 D = PI / 4
4150 FOR I1 = 1 TO 23
4160 M1 = E - E0 * SIN(E)
4170 E = E + SGN(M - M1) * D
4180 D = D / 2
4190 NEXT I1
4200 V = SQR((1 + E0) / (1 - E0))
4210 E = E * F
4220 V = 2 * ATN(V * SIN(E / 2) / COS(E / 2))
4230 R = A1 * (1 - E0 * COS(E))
4240 GOTO 4660
4250 A = SQR((1 + 9 * E0) / 10)
4260 B = 5 * (1 - E0) / (1 + 9 * E0)
4270 C = SQR(5 * (1 + E0) / (1 + 9 * E0))
4280 B1 = 3 * A * K * T / SQR(2 * Q * Q * Q)
4290 B2 = 1
4300 W1 = B2 * B1
4310 B3 = ATN(2 / W1)
4320 T1 = SIN(B3 / 2) / COS(B3 / 2)
4330 S1 = SGN(T1)
4340 T1 = ABS(T1)
4350 T2 = T1 ^ (1 / 3) * S1
4360 G = ATN(T2)
4370 S = 2 * COS(2 * G) / SIN(2 * G)
4380 A2 = B * S * S
4390 B0 = B2
4400 B2 = 0
4410 IF ABS(A2) > .3 THEN 3970
4420 B2 = B2 + 1 * A2 ^ 0
4430 B2 = B2 + 0 * A2 ^ 1
4440 B2 = B2 - .017142857# * A2 * A2
4450 B2 = B2 - 3.809524E-03 * A2 ^ 3
4460 B2 = B2 - 1.104267E-03 * A2 ^ 4
4470 B2 = B2 - 3.67358E-04 * A2 ^ 5
4480 B2 = B2 - 1.31675E-04 * A2 ^ 6
4490 B2 = B2 - 4.9577E-05 * A2 ^ 7
4500 IF ABS(B2 - B0) > 1E-08 THEN 4300
4510 C1 = 0
4520 C1 = C1 + 1 * A2 ^ 0
4530 C1 = C1 - .8 * A2 ^ 1
4540 C1 = C1 + .04571429# * A2 ^ 2
4550 C1 = C1 + .0152381 * A2 ^ 3
4560 C1 = C1 + .0056282 * A2 ^ 4
4570 C1 = C1 + 2.18783E-03 * A2 ^ 5
4580 C1 = C1 + 8.7905E-04 * A2 ^ 6
4590 C1 = C1 + 3.6155E-04 * A2 ^ 7
4600 C1 = SQR(1 / C1)
4610 V1 = C * C1 * S
4620 D1 = 1 / (1 + A2 * C1 * C1)
4630 V = 2 * ATN(V1)
4640 R = Q * D1 * (1 + V1 * V1)
4650 GOTO 4660
4660 IF V < 0 THEN 4680
4670 GOTO 4690
4680 V = V + 2 * PI
4690 GOTO 3040
4700 PRINT "FUERA DE LIMITE"
4710 CLOSE 1: END
4720 GOTO 3920

REM campo nombre
4801 numerito = 1
4802 OPEN "nombre.txt" FOR INPUT AS #2
4803  DO WHILE NOT EOF(2)
4804     INPUT #2, nombre$
4805     IF nuco = numerito THEN NOM$ = nombre$
4806     numerito = numerito + 1
4807  LOOP
4808 CLOSE #2
REM fin campo nombre

5224 RETURN
5225 GOSUB 8201
5230 IF M5 > 2 THEN 5280
5240 M6 = 12 + M5
5250 A6 = a5 - 1
5260 a5 = A6
5270 M5 = M6
5280 A0 = INT(a5 / 100)
5290 B0 = 2 - A0 + INT(A0 / 4)
5300 JJ0 = INT(365.25 * a5) + INT(30.6001 * (M5 + 1)) + d5 + 1720994.5#
5310 IF JJ0 > 2299170.50001# THEN 5330
5320 RETURN
5330 JJ0 = JJ0 + B0
5340 GOTO 5320

REM campo perihelio
5401 numerito = 1
5402 OPEN "peri.txt" FOR INPUT AS #2
5403  DO WHILE NOT EOF(2)
5404     INPUT #2, valor
5405     IF nuco = numerito THEN P8 = valor
5406     numerito = numerito + 1
5407  LOOP
5408 CLOSE #2
REM fin campo perihelio

5849 RETURN

REM campo nodo
5901 numerito = 1
5902 OPEN "nodo.txt" FOR INPUT AS #2
5903  DO WHILE NOT EOF(2)
5904     INPUT #2, valor
5905     IF nuco = numerito THEN N8 = valor
5906     numerito = numerito + 1
5907  LOOP
5908 CLOSE #2
REM fin campo nodo

6329 RETURN

REM campo inclinacion
6401 numerito = 1
6402 OPEN "inclina.txt" FOR INPUT AS #2
6403  DO WHILE NOT EOF(2)
6404     INPUT #2, valor
6405     IF nuco = numerito THEN I8 = valor
6406     numerito = numerito + 1
6407  LOOP
6408 CLOSE #2
REM fin campo inclinacion

6809 RETURN

REM campo E
6901 numerito = 1
6902 OPEN "campo_e.txt" FOR INPUT AS #2
6903  DO WHILE NOT EOF(2)
6904     INPUT #2, valor
6905     IF nuco = numerito THEN E8 = valor
6906     numerito = numerito + 1
6907  LOOP
6908 CLOSE #2
REM fin campo E

7300 RETURN

REM campo Q
7401 numerito = 1
7402 OPEN "campo_q.txt" FOR INPUT AS #2
7403  DO WHILE NOT EOF(2)
7404     INPUT #2, valor
7405     IF nuco = numerito THEN Q8 = valor
7406     numerito = numerito + 1
7407  LOOP
7408 CLOSE #2
REM fin campo Q

7800 RETURN

REM campo magnitudes
7901 numerito = 1
7902 OPEN "magnitud.txt" FOR INPUT AS #2
7903  DO WHILE NOT EOF(2)
7904     INPUT #2, valor
7905     IF nuco = numerito THEN k8 = valor
7906     numerito = numerito + 1
7907  LOOP
7908 CLOSE #2
REM fin campo magnitudes

REM campo magnitudesb
7951 numerito = 1
7952 OPEN "magnitub.txt" FOR INPUT AS #2
7953  DO WHILE NOT EOF(2)
7954     INPUT #2, valor
7955     IF nuco = numerito THEN kk8 = valor
7956     numerito = numerito + 1
7957  LOOP
7958 CLOSE #2
REM fin campo magnitudes

8199 RETURN

REM campo a¤o
8201 numerito = 1
8202 OPEN "a¤o.txt" FOR INPUT AS #2
8203  DO WHILE NOT EOF(2)
8204     INPUT #2, valor
8205     IF nuco = numerito THEN a5 = valor
8206     numerito = numerito + 1
8207  LOOP
8208 CLOSE #2
REM fin campo a¤o

REM campo mes
8601 numerito = 1
8602 OPEN "mes.txt" FOR INPUT AS #2
8603  DO WHILE NOT EOF(2)
8604     INPUT #2, valor
8605     IF nuco = numerito THEN M5 = valor
8606     numerito = numerito + 1
8607  LOOP
8608 CLOSE #2
REM fin campo mes

REM campo dia
9101 numerito = 1
9102 OPEN "dia.txt" FOR INPUT AS #2
9103  DO WHILE NOT EOF(2)
9104     INPUT #2, valor
9105     IF nuco = numerito THEN d5 = valor
9106     numerito = numerito + 1
9107  LOOP
9108 CLOSE #2
REM fin campo dia

9600 RETURN

30000 OPEN "nombre.txt" FOR OUTPUT AS #1
30010 OPEN "elem4.txt" FOR INPUT AS #2
30020 OPEN "peri.txt" FOR OUTPUT AS #3
30030 OPEN "nodo.txt" FOR OUTPUT AS #4
30040 OPEN "inclina.txt" FOR OUTPUT AS #5
30050 OPEN "campo_q.txt" FOR OUTPUT AS #6
30060 OPEN "campo_e.txt" FOR OUTPUT AS #7
30070 OPEN "a¤o.txt" FOR OUTPUT AS #8
30080 OPEN "mes.txt" FOR OUTPUT AS #9
30090 OPEN "dia.txt" FOR OUTPUT AS #10
30100 OPEN "magnitud.txt" FOR OUTPUT AS #11
30105 OPEN "magnitub.txt" FOR OUTPUT AS #21
30110 OPEN "tp.txt" FOR OUTPUT AS #12
30120 IF EOF(2) THEN 30250
30130 INPUT #2, NOM$, TP$, Q8, E8, P8, N8, I8, A8, PERIOD$, MPC$
30140 PRINT #1, NOM$
30150 PRINT #3, P8
30160 PRINT #4, N8
30170 PRINT #5, I8
30180 PRINT #6, Q8
30190 PRINT #7, E8
30200 PRINT #12, TP$
30210 PRINT #9, M5
30220 PRINT #10, d5
30230 REM PRINT #11, g8
30240 GOTO 30120
30250 CLOSE #2: END



40000 OPEN "tp.txt" FOR INPUT AS #2
40008 OPEN "a¤o.txt" FOR OUTPUT AS #8
40009 OPEN "mes.txt" FOR OUTPUT AS #9
40010 OPEN "dia.txt" FOR OUTPUT AS #10
40020 IF EOF(2) THEN 50
40030 INPUT #2, a5, M5, d5
40037 PRINT #8, a5
40038 PRINT #9, M5
40039 PRINT #10, d5
40048 GOTO 40020
40050 CLOSE #2: END

50000 OPEN "nom.txt" FOR OUTPUT AS #1
50010 OPEN "elem4.txt" FOR INPUT AS #2
50020 OPEN "p8.txt" FOR OUTPUT AS #3
50030 OPEN "n8.txt" FOR OUTPUT AS #4
50040 OPEN "i8.txt" FOR OUTPUT AS #5
50050 OPEN "q8.txt" FOR OUTPUT AS #6
50060 OPEN "e8.txt" FOR OUTPUT AS #7
50070 OPEN "a5.txt" FOR OUTPUT AS #8
50080 OPEN "m5.txt" FOR OUTPUT AS #9
50090 OPEN "d5.txt" FOR OUTPUT AS #10
50100 OPEN "magnitud.txt" FOR OUTPUT AS #11
50105 OPEN "magnitub.txt" FOR OUTPUT AS #21
50110 OPEN "tp.txt" FOR OUTPUT AS #12
50120 OPEN "y5.txt" FOR OUTPUT AS #13
50130 OPEN "s8.txt" FOR OUTPUT AS #14
50140 REM OPEN "g8$.txt" FOR OUTPUT AS #15
50150 IF EOF(2) THEN 50280
50160 INPUT #2, NOM$, eq, TP$, Q8, E8, P8, N8, I8, g8, gg8, HH8$
50170 PRINT #1, NOM$
50180 PRINT #3, P8
50190 PRINT #4, N8
50200 PRINT #5, I8
50210 PRINT #6, Q8
50220 PRINT #7, E8
50230 PRINT #12, TP$
50240 PRINT #9, M5
50250 REM PRINT #15, g8$
50260 PRINT #11, g8
50265 PRINT #21, gg8
50270 GOTO 50150
50280 CLOSE #2: END