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
220 INPUT "HASTA QUE MAGNITUD LE INTERESA LAS EFEMERIDES "; MMM9
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 NUCOM = 149
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) + 10 * LOG(R1) / LOG(10) + k8
3670 P7 = 1 + P7
3680 IF P7 > P9 THEN 470
3690 IF E6 > 15 AND M9 < MMM9 AND AMI < A6 AND AMA > A6 AND DMI < D7 AND DMA > D7 AND BAND = 0 THEN PRINT nom$
3700 IF E6 > 15 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 > 15 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 > 15 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 > 15 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 > 15 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 GOTO 3900
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; " G="; 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

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