1 CLS
2 PRINT
3 PRINT "          ********** CALENDARIO DE POLENES ALERGENOS ********** "
4 PRINT
5 PRINT "                      COPYRIGHT (C) 1995 INFORTECNICA      "
6 PRINT "                            PLAZA MARAGALL 4    "
7 PRINT "                            BARCELONA  08041    "
8 PRINT "                      VOZ 93 351 59 25, FAX 93 408 43 46 "
9 PRINT "                      www.infortecnica.com       "
10 PRINT
11 INPUT "NOMBRE DEL FICHERO PARA VOLCAR (LPT1: PARA IMPRESORA)"; FICH$
12 IF FICH$ = "" THEN PRINT "ERROR": GOTO 110
13 OPEN FICH$ FOR OUTPUT AS #1
19 PRINT "CALENDARIO DE PRESENCIA DE POLENES PRODUCTORES DE ALERGIAS  "
20 INPUT "A¥O, MES:"; A5, M1
30 INPUT "DIA INICIAL, DIA FINAL:"; D1, DI1
35 INPUT "INTERVALO EN DIAS"; INVAL
38 INVAL = INVAL - 1
60 PRINT #1,
100 PRINT #1,
110 PRINT #1,
120 PRINT #1,
130 PRINT #1, "-----------------------------------------------------------"
140 PRINT #1,
150 PRINT #1, "CALENDARIO"
160 PRINT #1, "----------"
170 PRINT #1,
180 PRINT "CALENDARIO"
190 PRINT "----------"
200 MI5 = M1
210 DI5 = D1
220 AI5 = A5
230 M5 = 1: D5 = 1
240 GOSUB 440
250 JI0 = J0
260 J0 = 0: J2 = 0: M5 = MI5: D5 = DI5: A5 = AI5
270 PI = 3.1415926#
280 RAD = 180 / PI
290 J9 = 9
300 N7 = 7
310 D0 = 0
320 Q7 = 7
330 R9 = 180 / PI
340 S9 = 1
350 N9 = 2
360 X0 = 0
370 GOTO 1800
380 W7 = A9 / 360 - INT(A9 / 360)
390 B2 = 9 + W7
400 B3 = B2 - INT(B2)
410 A1 = B3 * 2 * PI
420 RETURN
430 Q7 = 18
440 IF M5 > 2 THEN 490
450 M6 = 12 + M5
460 A6 = A5 - 1
470 A5 = A6
480 M5 = M6
490 A = INT(A5 / 100)
500 B = 2 - A + INT(A / 4)
510 J0 = INT(365.25 * A5) + INT(30.6001 * (M5 + 1)) + D5 + 1720994.5#
520 IF J0 > 2299170.50001# THEN 550
530 IF Q7 = 18 THEN 570
540 RETURN
550 J0 = J0 + B
560 GOTO 530
570 PRINT "JD="; J0
580 END
590 PRINT "FECHA JULIANA=";
600 INPUT J0
610 J6 = J0 + .5
620 J5 = INT(J6)
630 A = INT((J5 - 1867216.25#) / (36524.25))
640 A1 = J5 + 1 + A - INT(A / 4)
650 IF J5 < 2299161.001# THEN 780
660 B = A1 + 1524
670 C = INT((B - 122.1) / (365.25))
680 D = INT(365.25 * C)
690 E = INT((B - D) / (30.6001))
700 D5 = B - D - INT(30.6001 * E) + (J6 - J5)
710 IF E > 13.5 THEN 740
720 M5 = E - 1
730 GOTO 750
740 M5 = E - 13
750 IF M5 > 2.5 THEN 810
760 A5 = C - 4715
770 GOTO 820
780 A1 = J5
790 GOTO 660
800 GOTO 820
810 A5 = C - 4716
820 PRINT #1, "A¤o"; A5; ", mes"; M5; "dia"; D5
830 RETURN
840 V7 = R
850 A5 = B
860 GOSUB 920
870 Z = Y
880 V7 = X
890 A5 = A
900 GOSUB 920
910 RETURN
920 R9 = 180 / PI
930 A8 = (A5 / 360 - INT(A5 / 360)) * 360
940 C1 = INT(A8 / 90) + 1
950 A2 = (A8 / 90 - INT(A8 / 90)) * 90
960 X1 = COS(A2 / R9) * V7
970 Y1 = SIN(A2 / R9) * V7
980 IF C1 > 3.99999 THEN 1390
990 IF C1 > 2.99999 THEN 1420
1000 IF C1 > 1.99999 THEN 1450
1010 GOTO 1480
1020 X1 = ABS(X)
1030 Y1 = ABS(Y)
1040 V7 = SQR(X1 * X1 + Y1 * Y1)
1050 A = ATN(Y / X) * 180 / PI
1060 IF X > -1E-08 THEN 1100
1070 IF X < 0 THEN 1080
1080 A7 = 180 + A
1090 RETURN
1100 A7 = A
1110 RETURN
1120 PRINT "X=";
1130 INPUT X9
1140 PRINT "Y=";
1150 INPUT Y9
1160 PRINT "Z=";
1170 INPUT Z9
1180 GOSUB 1210
1190 PRINT "A"; A3; "B"; B7; "R"; V7
1200 END
1210 X = X9
1220 Y = Y9
1230 GOSUB 1020
1240 A3 = A7
1250 X = V7
1260 Y = Z9
1270 GOSUB 1020
1280 B7 = A7
1290 RETURN
1300 PRINT "DE POLARES A CARTESIANAS:A=";
1310 INPUT A
1320 PRINT "B=";
1330 INPUT B
1340 PRINT "R=";
1350 INPUT R
1360 GOSUB 840
1370 PRINT "X="; X; "Y="; Y; "Z="; Z
1380 END
1390 X = Y1
1400 Y = X1 * -1
1410 RETURN
1420 X = X1 * -1
1430 Y = Y1 * -1
1440 RETURN
1450 X = Y1 * -1
1460 Y = X1
1470 RETURN
1480 X = X1
1490 Y = Y1
1500 RETURN
1510 PRINT "DE GRADOS A GMS: GRD=";
1520 INPUT G7
1530 IF G7 < 0 THEN 1630
1540 G = INT(G7)
1550 M6 = (G7 - G) * 60
1560 M = INT(M6)
1570 S = (M6 - M) * 60
1580 G9 = G + M / 100 + S / 10000
1590 RETURN
1600 GOSUB 1510
1610 PRINT "GRD="; G7; "G="; G; "M="; M; "S="; S; "GMS="; G9
1620 END
1630 G7 = G7 + 360
1640 GOSUB 1540
1650 G9 = -1 * G9
1660 RETURN
1670 PRINT "G=";
1680 INPUT G
1690 PRINT "M=";
1700 INPUT M
1710 PRINT "S=";
1720 INPUT S
1730 G7 = G + M / 60 + S / 3600
1740 G9 = G + M / 100 + S / 10000
1750 RETURN
1760 PRINT "DE GMS A GRADOS";
1770 GOSUB 1670
1780 PRINT "GMS="; G9; "GRD="; G7
1790 END
1800 REM 7
1810 J1 = 1
1820 J2 = 0
1830 J8 = INCHO
1840 J4 = 1
1850 GOSUB 440
1860 IF J4 > J1 THEN 5430
1870 J4 = J4 + 1
1880 J7 = 0
1890 GOTO 3520
1900 R9 = 180 / PI
1910 L0 = L1 * R9
1920 GOTO 2700
1930 GOSUB 1950
1940 RETURN
1950 T = (J0 - 2415020!) / 36525!
1960 T3 = T
1970 L7 = 270.4344 + 481267.8831# * T - .0001133 * T * T
1980 L6 = 279.69668# + 36000.76892# * T + .0003025 * T * T
1990 A9 = L6
2000 GOSUB 380
2010 L6 = A1
2020 M6 = 358.47583# + 35999.0475# * T - .00015 * T * T
2030 A9 = M6
2040 GOSUB 380
2050 M6 = A1
2060 M7 = 296.1046 + 477198.8491# * T + .009192 * T * T
2070 C5 = 1.91946 - .004789 * T - .000014 * T * T
2080 C6 = PI / 180 * (C5 * SIN(M6) + (.020094 - .0001 * T) * SIN(2 * M6) + .000293 * SIN(3 * M6))
2090 N0 = 259.1833 - 1934.142 * T + .002078 * T * T
2100 L1 = C6 + L6
2110 V6 = M6 + C6
2120 E7 = 1.675104E-02 - .0000418 * T
2130 R6 = ((1.0000002# * (1 - E7 ^ 2)) / (1 + E7 * COS(V6)))
2140 N6 = 259.13 - 1934.142 * T
2150 A9 = N6
2160 GOSUB 380
2170 N6 = A1
2180 W = L1 - PI / 180 * (.00569 - .00479 * SIN(N6))
2190 N1 = -(17.2327 + .01737 * T) * SIN(N0 / R9) - 1.273 * SIN(2 * L6 / R9)
2200 N2 = 9.21 * COS(N0 / R9) + .552 * COS(2 * L6 / R9)
2210 O1 = 23.452294# - .0130125 * T + .00256 * COS(N6)
2220 O5 = 23.452294# - .0130125 * T
2230 L5 = L1 - PI / 180 * 1.396 * (T - .5)
2240 X1 = R6 * COS(W)
2250 X5 = R6 * COS(L5)
2260 Y7 = R6 * SIN(W) * COS(O1 / (180 / PI))
2270 Y5 = R6 * SIN(L5) * COS(O5 / (180 / PI))
2280 Z3 = R6 * SIN(W) * SIN(O1 / (180 / PI))
2290 Z5 = R6 * SIN(L5) * SIN(O5 / (180 / PI))
2300 R9 = 180 / PI
2310 V7 = R6
2320 A5 = L5 * R9
2330 GOSUB 920
2340 X2 = X
2350 Y2 = Y
2360 V7 = R6
2370 A5 = W * R9
2380 GOSUB 920
2390 X3 = X
2400 Y3 = Y
2410 X9 = X5
2420 Y9 = Y5
2430 Z9 = Z5
2440 GOSUB 1210
2450 W5 = A3
2460 D5 = B7
2470 X9 = X3
2480 Y9 = Y7
2490 Z9 = Z3
2500 GOSUB 1210
2510 W6 = A3
2520 D9 = B7
2530 DS9 = D9
2540 P1 = 8.793999 / R6
2550 S0 = .267 / R6
2560 G7 = W5 / 15
2570 GOSUB 1530
2580 W0 = G9
2590 G7 = W6 / 15
2600 GOSUB 1530
2610 W8 = G9
2620 G7 = D5
2630 GOSUB 1530
2640 D6 = G9
2650 G7 = D9
2660 GOSUB 1530
2670 D0 = G9
2680 RETURN
2690 END
2700 REM 1
2710 T1 = 99.691 + 36000.769# * T3
2720 T2 = (T1 - A3 - L2) / 360
2730 T2 = 24 - 24 * (T2 - INT(T2))
2740 TI2 = T2
2750 S = (Z + F + D9) / 2
2760 T9 = SIN((S - F) / R9) * SIN((S - D9) / R9)
2770 T8 = COS((S - Z) / R9) * COS(S / R9)
2780 T9 = SQR(ABS(T9 / T8)) * (T9 / T8) / ABS(T9 / T8)
2790 T9 = R9 * 2 * ATN(T9) / 15
2800 T7 = T2 - T9
2810 T8 = T2 + T9
2820 RETURN
2830 T7 = INT(T7) + (T7 - INT(T7)) * .6
2840 T8 = INT(T8) + (T8 - INT(T8)) * .6
2850 T2 = INT(T2) + (T2 - INT(T2)) * .6
2860 REM 2
2870 T6 = T + J7 * 1.0027 / 36525! / 24
2880 L9 = (270.434164# + 481267.8831# * T6) / RAD
2890 M8 = (358.475833# + 35999.0498# * T6) / RAD
2900 M9 = (296.10468# + 477198.8491# * T6) / RAD
2910 D8 = (350.737486# + 445267.1142# * T6) / RAD
2920 F8 = (11.250889# + 483202.0251# * T6) / RAD
2930 O8 = (259.183275# - 1934.142 * T6) / RAD
2940 E8 = 1 - .002495 * T6
2950 L7 = L9 + 6.28875 / RAD * SIN(M9) + 1.274018 / RAD * SIN(2 * D8 - M9) + .658309 / RAD * SIN(2 * D8)
2960 L7 = L7 + .213616 / RAD * SIN(2 * M9) - .185596 / RAD * SIN(M8) - .114336 / RAD * SIN(2 * F8)
2970 L7 = L7 + .05879 / RAD * SIN(2 * D8 - 2 * M9) + .0572 / RAD * SIN(2 * D8 - M8 - M9)
2980 L7 = L7 + .05332 / RAD * SIN(2 * D8 + M9) + .045874 / RAD * SIN(2 * D8 - M8)
2990 L7 = L7 + .041024 / RAD * SIN(M9 - M8) - .03472 / RAD * SIN(D8) - .030465 / RAD * SIN(M8 + M9)
3000 L7 = L7 / (360 / RAD)
3010 L7 = 360 / RAD * (L7 - INT(L7))
3020 B7 = 5.128189 / RAD * SIN(F8) + .280606 / RAD * SIN(M9 + F8) + .277693 / RAD * SIN(M9 - F8)
3030 B7 = B7 + .17324 / RAD * SIN(2 * D8 - F8) + .055413 / RAD * SIN(2 * D8 + F8 - M9) + .04627 / RAD * SIN(2 * D8 - F8 - M9)
3040 B7 = B7 + .03257 / RAD * SIN(2 * D8 + F8) + .0172 / RAD * SIN(2 * M9 + F8)
3050 C9 = COS(L7 - L1) * COS(B7)
3060 C8 = ATN(SQR(1 - C9 * C9) / C9)
3070 IF C9 > 0 THEN 3090
3080 C8 = 180 / RAD - C8
3090 I = 180 - C8 * RAD - .1468 * (1 - .0549 * SIN(M8) / (1 - .0167 * SIN(M7)) * SIN(C8))
3100 K = (1 + COS(I / RAD)) / 2
3110 P7 = .95072 / RAD + .05182 / RAD * COS(M9) + 9.530001E-03 / RAD * COS(2 * D8 - M9) + .00784 / RAD * COS(2 * D8)
3120 P7 = P7 + .002824 / RAD * COS(2 * M9)
3130 O3 = O1
3140 T7 = 99.691 + 36000.769# * T - L2 + J7 * 1.0027 * 15
3150 N = COS(L7) * COS(B7) - SIN(P7) * COS(F / RAD) * COS(T7 / RAD)
3160 L4 = (SIN(L7) * COS(B7) - SIN(P7) * (SIN(F / RAD) * SIN(O3 / RAD) + COS(F / RAD) * COS(O3 / RAD) * SIN(T7 / RAD))) / N
3170 L4 = ATN(L4) * RAD
3180 IF ABS(L4 - (L7 * RAD)) < 90 THEN 3210
3190 IF ABS(L4 - (L7 * RAD)) > 270 THEN 3210
3200 L4 = 180 + L4
3210 B4 = (COS(L4 / RAD) * (SIN(B7) - SIN(P7) * (SIN(F / RAD) * COS(O3 / RAD) - COS(F / RAD) * SIN(O3 / RAD) * SIN(T7 / RAD)))) / N
3220 B4 = ATN(B4) * RAD
3230 IF L4 > 0 THEN 3250
3240 L4 = 360 + L4
3250 A4 = (SIN(L4 / RAD) * COS(O3 / RAD) - TAN(B4 / RAD) * SIN(O3 / RAD)) / COS(L4 / RAD)
3260 A4 = ATN(A4) * RAD
3270 IF ABS(A4 - L4) < 90 THEN 3300
3280 IF ABS(A4 - L4) > 270 THEN 3300
3290 A4 = A4 + 180
3300 D3 = SIN(B4 / RAD) * COS(O3 / RAD) + COS(B4 / RAD) * SIN(O3 / RAD) * SIN(L4 / RAD)
3310 D4 = ATN(D3 / SQR(1 - D3 * D3)) * RAD
3320 IF A4 > 0 THEN 3340
3330 A4 = 360 + A4
3340 A3 = A4
3350 D9 = D4
3360 A7 = (SIN(L7) * COS(O3) - TAN(B7) * SIN(O3)) / COS(L7)
3370 A7 = ATN(A7) * RAD
3380 IF ABS(A7 - L7) < 90 THEN 3410
3390 IF ABS(A7 - L7) > 270 THEN 3410
3400 A7 = A7 + 180
3410 D3 = SIN(B7) * COS(O3) + COS(B7) * SIN(O3) * SIN(L7)
3420 D7 = ATN(D3 / SQR(1 - D3 * D3)) * RAD
3430 IF A7 > 0 THEN 3450
3440 A7 = 360 + A7
3450 A7 = A7 / 15
3460 A7 = INT(A7) + (A7 - INT(A7)) * .6
3470 A4 = A4 / 15
3480 A4 = INT(A4) + (A4 - INT(A4)) * .6
3490 RETURN
3500 PRINT A7; A4; D7; D4; "AR. Y DEC. GEO. Y TOPC."
3510 END
3520 GOSUB 1950
3530 Z = 108
3540 GOSUB 2700
3550 U1 = T7
3560 U2 = T8
3570 Z = 90.85
3580 GOSUB 2700
3590 U3 = T7
3600 U4 = T2
3610 U5 = T8
3620 GOTO 3750
3630 T = T3 + U3 / 24 / 36525!
3640 GOSUB 1970
3650 GOSUB 2700
3660 U3 = T7
3670 T = T3 + U4 / 24 / 36525!
3680 GOSUB 1970
3690 GOSUB 2700
3700 U4 = T2
3710 T = T3 + U5 / 24 / 36525!
3720 GOSUB 1970
3730 GOSUB 2700
3740 U5 = T8
3750 T = T3
3760 J7 = 0
3770 U9 = 0
3780 GOSUB 2860
3790 K9 = K
3800 GOSUB 2700
3810 U6 = T7
3820 U7 = T2
3830 U8 = T8
3840 J7 = U6
3850 GOSUB 2860
3860 GOSUB 2700
3870 U6 = T7
3880 J7 = U7
3890 GOSUB 2860
3900 GOSUB 2700
3910 U7 = T2
3920 J7 = T8
3930 GOSUB 2860
3940 GOSUB 2700
3950 U8 = T8
3960 U9 = U9 + 1
3970 IF U9 > 1 THEN 3990
3980 GOTO 3840
3990 U15 = U5: U13 = U3: U11 = U1
4000 U18 = U8: U16 = U6
4010 UC1 = U3 - U1
4020 U1 = INT(U1) + (U1 - INT(U1)) * .6 + J8
4030 U2 = INT(U2) + (U2 - INT(U2)) * .6 + J8
4040 U3 = INT(U3) + (U3 - INT(U3)) * .6 + J8
4050 U4 = INT(U4) + (U4 - INT(U4)) * .6 + J8
4060 U5 = INT(U5) + (U5 - INT(U5)) * .6 + J8
4070 U6 = U6 + J8
4080 IF U6 > 0 THEN 4120
4090 U6 = U6 + 24
4100 U6 = -1 * (INT(U6) + (U6 - INT(U6)) * .6)
4110 GOTO 4130
4120 U6 = INT(U6) + (U6 - INT(U6)) * .6
4130 U7 = INT(U7) + (U7 - INT(U7)) * .6 + J8
4140 U8 = INT(U8) + (U8 - INT(U8)) * .6 + J8
4150 GOSUB 610
4160 J9 = A5 + M5 / 100 + D5 / 10000
4170 J0 = J0 + J2
4180 DS = INT((J0 / 7 - INT(J0 / 7)) * 7)
4190 PRINT #1,
4200 GOSUB 5440
4210 PRINT #1,
4220 DB = J0 - JI0 + 1
4230 IF DS < 7 THEN 4260
4240 IF DS = 7 THEN DS1 = 1
4250 GOTO 4270
4260 DS1 = DS + 1
4270 NS = INT((DB - DS1) / 7 + 1.5)
4280 PRINT "Semana del a¤o n£mero "; NS
4290 PRINT #1, "Semana del a¤o n£mero "; NS
4300 PRINT #1,
4310 PRINT "Dias transcurridos desde el inicio del a¤o"; DB
4320 PRINT #1, "D¡as transcurridos desde el inicio del a¤o:"; DB
4330 PRINT #1,
4340 DA = 365
4350 IF (A5 / 4 - INT(A5 / 4)) = 0 THEN DA = 366
4360 PRINT "Dias que faltan hasta fin de a¤o"; DA - DB
4370 PRINT #1, "D¡as que faltan hasta fin de a¤o:"; DA - DB
4380 PRINT #1,
4390 GOTO 30000
4400 PRINT #1, "D¡as transcurridos desde el inicio del siglo:"; J0 - 2415019.5# + 1
4410 PRINT #1,
4420 PRINT "Dias que faltan hasta fin de siglo"; 36525! - (J0 - 2415019.5#) - 1
4430 PRINT #1, "D¡as que faltan hasta fin de siglo:"; 36525! - (J0 - 2415019.5#) - 1
4440 GOTO 30000
4470 PRINT #1,
4480 PRINT "D¡a juliano a las 12 horas:"; J0
4490 PRINT #1, "D¡a juliano a las 12 horas:"; J0
4500 PRINT #1,
4510 TS1 = T1 / 360
4520 TS1 = 24 * (TS1 - INT(TS1))
4530 PRINT #1, "Tiempo sid‚reo de Greenwich (0h):"; INT(TS1); " horas "; INT(60 * (TS1 - (INT(TS1)))); " minutos"
4540 PRINT "Tiempo sid‚reo de Greenwich (0h):"; INT(TS1); " horas "; INT(60 * (TS1 - (INT(TS1)))); " minutos"
4550 PRINT #1,
4560 TS2 = (T1 + L2) / 360
4570 TS2 = 24 * (TS2 - INT(TS2))
4580 PRINT #1, "Tiempo sid‚reo local (0h):"; INT(TS2); " horas "; INT(60 * (TS2 - (INT(TS2)))); " minutos"
4590 PRINT "Tiempo sid‚reo local (0h):"; INT(TS2); " horas "; INT(60 * (TS2 - (INT(TS2)))); " minutos"
4600 PRINT #1,
4610 HS = U3 - U11 - J8
4620 PRINT "SOL Y LUNA"
4630 PRINT #1, "SOL Y LUNA"
4640 PRINT #1, "----------"
4650 EC1 = INT(UC1)
4660 FC1 = INT((UC1 - EC1) * 60)
4670 PRINT "La duraci¢n del crep£sculo es de "; EC1; " horas,"; FC1; " minutos."
4680 PRINT #1, "La duraci¢n del crep£sculo es de "; EC1; " horas,"; FC1; " minutos."
4690 PRINT #1,
4700 PRINT #1, "El crep£sculo matutino comienza a las "; INT(U1); " horas,"; INT((U1 - INT(U1)) * 100); " minutos."
4710 PRINT "El crep£sculo matutino comienza a las "; INT(U1); " horas,"; INT((U1 - INT(U1)) * 100); " minutos."
4720 PRINT #1,
4730 PRINT #1, "El orto solar comienza a las "; INT(U3); " horas,"; INT((U3 - INT(U3)) * 100); " minutos."
4740 PRINT "El orto solar comienza a las "; INT(U3); " horas,"; INT((U3 - INT(U3)) * 100); " minutos."
4750 PRINT #1,
4760 PRINT #1, "El ocaso solar finaliza a las "; INT(U5); " horas,"; INT((U5 - INT(U5)) * 100); " minutos."
4770 PRINT "El ocaso solar finaliza a las "; INT(U5); " horas,"; INT((U5 - INT(U5)) * 100); " minutos."
4780 PRINT #1,
4790 PRINT #1, "El crep£sculo vespertino finaliza a las "; INT(U2); " horas,"; INT((U2 - INT(U2)) * 100); " minutos."
4800 PRINT "El crep£sculo vespertino finaliza a las "; INT(U2); " horas,"; INT((U2 - INT(U2)) * 100); " minutos."
4810 PRINT #1,
4820 HS = U15 - U13
4830 PRINT "La duracion del dia solar es de "; INT(HS); "horas"; INT((HS - INT(HS)) * 60); "minutos"
4840 PRINT #1, "La duraci¢n del d¡a solar es de "; INT(HS); "horas"; INT((HS - INT(HS)) * 60); "minutos"
4850 PRINT #1,
4860 UR4 = INT(U4): UR5 = INT((U4 - INT(U4)) * 100)
4870 PRINT "El mediodia solar tendra lugar a las"; UR4; "horas"; UR5; "minutos"
4880 PRINT #1, "El mediod¡a solar tendr  lugar a las"; UR4; "horas"; UR5; "minutos"
4890 PRINT #1,
4900 AB = INT(((90 - 41.23) + DS9) + .5)
4910 PRINT "La m xima altura que alcanzara hoy el Sol sobre el horizonte sera de"; AB; "grados"
4920 PRINT #1, "La m xima altura que alcanzar  hoy el Sol sobre el horizonte ser  de"; AB; "grados"
4930 PRINT #1,
4940 GOTO 5280
4950 GOSUB 10000
4980 ME = M5 + D5 / 31 - 6
4990 IF ME > 6 THEN ME = ME - 12
5000 TL1 = ME * 3.01201152# + ME * ME * -.86100884# + ME * ME * ME * -.151525 + ME * ME * ME * ME * .010679761415# + ME * ME * ME * ME * ME * .0018947964# + 22.6
5010 TL2 = ME * 2.989125 + ME * ME * -.926251 + ME * ME * ME * -.1453997# + ME * ME * ME * ME * 9.941039352999999D-03 + ME * ME * ME * ME * ME * .00174962292# + 29.65
5020 TL3 = ME * 3.0792 + ME * ME * -.78842 + ME * ME * ME * -.1635751# + ME * ME * ME * ME * .01112967944# + ME * ME * ME * ME * ME * .002171003# + 15.5579
5030 TB1 = ME * 2.9624 + ME * ME * -.5998 + ME * ME * ME * -.1586 + ME * ME * ME * ME * .00786 + ME * ME * ME * ME * ME * .0021 + 21.65
5040 TB2 = ME * 2.9137 + ME * ME * -.62315 + ME * ME * ME * -.1583 + ME * ME * ME * ME * .00787 + ME * ME * ME * ME * ME * .0021 + 25.21
5050 TB3 = ME * 2.9902 + ME * ME * -.57012 + ME * ME * ME * -.1576 + ME * ME * ME * ME * .00768 + ME * ME * ME * ME * ME * .00208 + 18.02
5060 TT1 = ME * 2.9076 + ME * ME * -.5222701 + ME * ME * ME * -.1479 + ME * ME * ME * ME * .00654 + ME * ME * ME * ME * ME * .0018 + 20.46
5070 TT2 = ME * 2.7891 + ME * ME * -.46717 + ME * ME * ME * -.1379 + ME * ME * ME * ME * .00545 + ME * ME * ME * ME * ME * .0016 + 23.59
5080 TT3 = ME * 3.0469 + ME * ME * -.57091 + ME * ME * ME * -.1607 + ME * ME * ME * ME * .00742 + ME * ME * ME * ME * ME * .0021 + 17.27
5090 TG1 = ME * 3.029 + ME * ME * -.66711 + ME * ME * ME * -.1608 + ME * ME * ME * ME * 8.500001E-03 + ME * ME * ME * ME * ME * .0021 + 20.89
5100 TG2 = ME * 3.0689 + ME * ME * -.68772 + ME * ME * ME * -.1711 + ME * ME * ME * ME * .00837 + ME * ME * ME * ME * ME * .0024 + 26.83
5110 TG3 = ME * 2.9784 + ME * ME * -.64505 + ME * ME * ME * -.1501 + ME * ME * ME * ME * .00854 + ME * ME * ME * ME * ME * .0018 + 14.93
5115 GOTO 5260
5120 PRINT "La temperatura media que corresponde a la fecha de hoy en Barcelona es de "; INT(TB1); "grados, en Tarragona es de "; INT(TT1); ",en Lleida es de "; INT(TL1); ", en Girona es de "; INT(TG1); "grados centigrados."
5130 PRINT #1, "La temperatura media que corresponde a la fecha de hoy en Barcelona es de "; INT(TB1)
5140 PRINT #1,
5150 PRINT #1, "grados, en Tarragona es de "; INT(TT1); ",en Lleida "; INT(TL1); ", en Girona"; INT(TG1)
5160 PRINT #1,
5170 PRINT "La maxima media para hoy en Barcelona es de "; INT(TB2); ",en Tarragona es de "; INT(TT2); ",en Lleida es de "; INT(TL2); ", en Girona es de "; INT(TG2); "grados"
5180 PRINT #1, "La m xima media para hoy en Barcelona es de"; INT(TB2); ",en Tarragona es de"; INT(TT2)
5190 PRINT #1,
5200 PRINT #1, ",en Lleida es de "; INT(TL2); ", en Girona es de "; INT(TG2); "grados"
5210 PRINT #1,
5220 PRINT "En cuanto a la minima media en Barcelona es de "; INT(TB3); ",en Tarragona es de "; INT(TT3); ",en Lleida es de "; INT(TL3); ", en Girona es de "; INT(TG3); "grados"
5230 PRINT #1, "En cuanto a la m¡nima media en Barcelona es de"; INT(TB3); ",en Tarragona es de"; INT(TT3)
5240 PRINT #1,
5250 PRINT #1, ",en Lleida es de "; INT(TL3); ", en Girona es de "; INT(TG3); "grados"
5260 PRINT #1,
5270 GOTO 5800
5280 LET AS0 = U7
5290 GOSUB 5600
5300 UR7 = INT(U7): UR8 = INT((U7 - INT(U7)) * 100)
5310 HL = U18 - U16
5320 PRINT "La duracion del dia lunar es de "; INT(HL); "horas"; INT((HL - INT(HL)) * 60); "minutos"
5330 PRINT #1, "La duraci¢n del d¡a lunar es de "; INT(HL); "horas"; INT((HL - INT(HL)) * 60); "minutos"
5340 PRINT #1,
5350 PRINT "El mediodia lunar a las "; UR7; "horas"; UR8; "minutos en la constelacion de "; C$
5360 PRINT #1, "El mediod¡a lunar a las "; UR7; "horas"; UR8; "minutos en la constelaci¢n de "; C$
5370 PRINT #1,
5380 GOSUB 5740
5390 PRINT "La Luna en "; R$; "con un disco iluminado del "; INT(K9 * 100); "por ciento"
5400 PRINT #1, "La Luna en "; R$; "con un disco iluminado del "; INT(K9 * 100); "por ciento"
5410 PRINT #1,
5420 GOTO 4950
5430 END
5440 DS = INT(DS + 1)
5450 IF DS = 7 THEN PRINT "Lunes"
5460 IF DS = 7 THEN PRINT #1, "Lunes"
5470 IF DS = 1 THEN PRINT "Martes"
5480 IF DS = 1 THEN PRINT #1, "Martes"
5490 IF DS = 2 THEN PRINT "Mi‚rcoles"
5500 IF DS = 2 THEN PRINT #1, "Mi‚rcoles"
5510 IF DS = 3 THEN PRINT "Jueves"
5520 IF DS = 3 THEN PRINT #1, "Jueves"
5530 IF DS = 4 THEN PRINT "Viernes"
5540 IF DS = 4 THEN PRINT #1, "Viernes"
5550 IF DS = 5 THEN PRINT "S bado"
5560 IF DS = 5 THEN PRINT #1, "S bado"
5570 IF DS = 6 THEN PRINT "Domingo"
5580 IF DS = 6 THEN PRINT #1, "Domingo"
5590 RETURN
5600 IF A4 > 1.4 AND A4 < 3.2 THEN C$ = "Aries"
5610 IF A4 > 3.2 AND A4 < 5.5 THEN C$ = "Taurus"
5620 IF A4 > 5.5 AND A4 < 8 THEN C$ = "Geminis"
5630 IF A4 > 8 AND A4 < 9.2 THEN C$ = "Cancer"
5640 IF A4 > 9.2 AND A4 < 11.35 THEN C$ = "Leo"
5650 IF A4 > 11.35 AND A4 < 14.2 THEN C$ = "Virgo"
5660 IF A4 > 14.2 AND A4 < 15.4 THEN C$ = "Libra"
5670 IF A4 > 15.4 AND A4 < 16.4 THEN C$ = "Escorpio"
5680 IF A4 > 16.4 AND A4 < 17.4 THEN C$ = "Ophiuco"
5690 IF A4 > 17.4 AND A4 < 20! THEN C$ = "Sagitario"
5700 IF A4 > 20! AND A4 < 22! THEN C$ = "Capricornio"
5710 IF A4 > 22! AND A4 < 24! THEN C$ = "Acuario"
5720 IF A4 > 0 AND A4 < 1.4 THEN C$ = "Piscis"
5730 RETURN
5740 IF U7 > U4 AND (U7 - U4) < 12 THEN R$ = "creciente "
5750 IF U7 < U4 THEN R$ = "menguante "
5760 IF K9 < .06 THEN R$ = "novilunio "
5770 IF K9 > .94 THEN R$ = "plenilunio "
5780 RETURN
5785 GOSUB 10000
5790 D1 = D1 + 1
5800 D1 = D1 + INVAL
5810 IF D1 > DI1 THEN END
5820 GOTO 60
10000 PRINT #1,
10010 PRINT #1, "IRRADIACION SOLAR"
10020 W = 20 / (COS(AB / (180 / 3.1416)))
10030 PRINT #1, "Irradiancia ultravioleta previsible en W/m*m ="; W
10035 PRINT #1,
10040 PRINT #1, "Tiempos m ximos en minutos recomendados de exposici¢n al Sol"
10050 PRINT #1, "------------------------------------------------------------"
10060 PRINT #1, "Pieles muy blancas: Primeras exposiciones"
10070 PRINT #1, "Ma¤anas:"; INT(325 / W * 2)
10080 PRINT #1, "Mediod¡a:"; INT(325 / W)
10090 PRINT #1, "Tardes:"; INT(325 / W * 2)
10095 PRINT #1,
10100 PRINT #1, "Pieles que pigmentan ligeramente: Primeras exposiciones"
10110 PRINT #1, "Ma¤anas:"; INT(325 / W * 4)
10120 PRINT #1, "Mediod¡a:"; INT(325 / W * 2)
10130 PRINT #1, "Tardes:"; INT(325 / W * 4)
10135 PRINT #1,
10140 PRINT #1, "Pigmentaci¢n normal: Primeras exposiciones"
10150 PRINT #1, "Ma¤anas:"; INT(325 / W * 6)
10160 PRINT #1, "Mediod¡a:"; INT(325 / W * 3)
10170 PRINT #1, "Tardes:"; INT(325 / W * 6)
10175 PRINT #1,
10180 PRINT #1, "Pieles de muy f cil pigmentaci¢n: Primeras exposiciones"
10190 PRINT #1, "Ma¤anas:"; INT(325 / W * 8)
10200 PRINT #1, "Mediod¡a:"; INT(325 / W * 4)
10210 PRINT #1, "Tardes:"; INT(325 / W * 8)
10215 PRINT #1,
10260 PRINT #1, "Pieles muy blancas: Exposiciones posteriores"
10270 PRINT #1, "Ma¤anas:"; INT(325 / W * 4)
10280 PRINT #1, "Mediod¡a:"; INT(325 / W * 2)
10290 PRINT #1, "Tardes:"; INT(325 / W * 4)
10295 PRINT #1,
10300 PRINT #1, "Pieles que pigmentan ligeramente: Exposiciones posteriores"
10310 PRINT #1, "Ma¤anas:"; INT(325 / W * 8)
10320 PRINT #1, "Mediod¡a:"; INT(325 / W * 4)
10330 PRINT #1, "Tardes:"; INT(325 / W * 8)
10335 PRINT #1,
10340 PRINT #1, "Pigmentaci¢n normal: Exposiciones posteriores"
10350 PRINT #1, "Ma¤anas:"; INT(325 / W * 12)
10360 PRINT #1, "Mediod¡a:"; INT(325 / W * 6)
10370 PRINT #1, "Tardes:"; INT(325 / W * 12)
10375 PRINT #1,
10380 PRINT #1, "Pieles de muy f cil pigmentaci¢n: Exposiciones posteriores"
10390 PRINT #1, "Ma¤anas:"; INT(325 / W * 16)
10400 PRINT #1, "Mediod¡a:"; INT(325 / W * 8)
10410 PRINT #1, "Tardes:"; INT(325 / W * 16)
10415 PRINT #1,
10420 PRINT #1, "Los tiempos m ximos se prolongan cuando se utilizan cremas protectoras"
10500 RETURN
30000 PRINT
30010 TD = DB
30020 PRINT #1,
30030 PRINT #1, "ANDALUCIA"
30040 IF TD > 92 AND TD < 183 THEN PRINT #1, "olivo, ";
30050 IF TD > 92 AND TD < 213 THEN PRINT #1, "gram¡neas, ";
30060 IF TD > 92 AND TD < 183 THEN PRINT #1, "encina, ";
30070 IF TD > 61 AND TD < 243 THEN PRINT #1, "parietaria, ";
30080 IF TD > 170 AND TD < 244 THEN PRINT #1, "artemisia, ";
30090 IF TD > 145 AND TD < 264 THEN PRINT #1, "chenopodium, ";
30100 IF TD > 244 AND TD < 274 THEN PRINT #1, "lili ceas, ";
30110 IF TD > 244 AND TD < 274 THEN PRINT #1, "fresno, ";
30120 IF TD > 244 AND TD < 274 THEN PRINT #1, "aliguste, ";
30130 IF TD > 122 AND TD < 183 THEN PRINT #1, "plantago, ";
30140 IF TD > 78 AND TD < 132 THEN PRINT #1, "pl tano de jard¡n, ";
30150 PRINT #1,
30160 PRINT #1, "ARAGON"
30170 IF TD > 122 AND TD < 183 THEN PRINT #1, "gram¡neas, ";
30180 IF TD > 122 AND TD < 167 THEN PRINT #1, "olivo, ";
30190 IF TD > 122 AND TD < 167 THEN PRINT #1, "encina, ";
30200 IF TD > 122 AND TD < 213 THEN PRINT #1, "plantago, ";
30210 IF TD > 117 AND TD < 155 THEN PRINT #1, "pl tano de jard¡n, ";
30220 IF TD > 60 AND TD < 91 THEN PRINT #1, "cipr‚s, ";
30230 PRINT #1,
30240 PRINT #1, "ASTURIAS"
30250 IF TD > 122 AND TD < 213 THEN PRINT #1, "gram¡neas, ";
30260 IF TD > 115 AND TD < 290 THEN PRINT #1, "plantago, ";
30270 PRINT #1,
30280 PRINT #1, "BALEARES"
30290 IF TD > 30 AND TD < 310 THEN PRINT #1, "parietaria, ";
30300 IF TD > 122 AND TD < 200 THEN PRINT #1, "gram¡neas, ";
30310 IF TD > 129 AND TD < 183 THEN PRINT #1, "olivo, ";
30320 IF TD > 183 AND TD < 244 THEN PRINT #1, "artemisia, ";
30330 IF TD > 100 AND TD < 132 THEN PRINT #1, "pl tano de jard¡n, ";
30340 PRINT #1,
30350 PRINT #1, "CANTABRIA"
30360 IF TD > 122 AND TD < 205 THEN PRINT #1, "gram¡neas, ";
30370 IF TD > 110 AND TD < 193 THEN PRINT #1, "roble, ";
30380 IF TD > 135 AND TD < 193 THEN PRINT #1, "encina, ";
30390 IF TD > 104 AND TD < 136 THEN PRINT #1, "pl tano de jard¡n, ";
30400 IF TD > 136 AND TD < 280 THEN PRINT #1, "ortiga-parietaria, ";
30410 IF TD > 112 AND TD < 245 THEN PRINT #1, "plantago, ";
30420 IF TD > 92 AND TD < 124 THEN PRINT #1, "abedul, ";
30430 PRINT #1,
30440 PRINT #1, "CASTILLA-LA MANCHA"
30450 IF TD > 126 AND TD < 183 THEN PRINT #1, "gram¡neas, ";
30460 IF TD > 126 AND TD < 183 THEN PRINT #1, "olivo, ";
30470 IF TD > 126 AND TD < 226 THEN PRINT #1, "chenopodium, ";
30480 IF TD > 126 AND TD < 226 THEN PRINT #1, "salsola, ";
30490 PRINT #1,
30500 PRINT #1, "CASTILLA-LEON"
30510 IF TD > 126 AND TD < 213 THEN PRINT #1, "gram¡neas, ";
30520 IF TD > 126 AND TD < 183 THEN PRINT #1, "olivo, ";
30530 IF TD > 104 AND TD < 183 THEN PRINT #1, "encina, ";
30540 IF TD > 135 AND TD < 183 THEN PRINT #1, "plantago, ";
30550 PRINT #1,
30560 PRINT #1, "CATALU¥A"
30570 IF TD > 60 AND TD < 300 THEN PRINT #1, "parietaria, ";
30580 IF TD > 106 AND TD < 183 THEN PRINT #1, "gram¡neas, ";
30590 IF TD > 122 AND TD < 173 THEN PRINT #1, "olivo, ";
30600 IF TD > 90 AND TD < 122 THEN PRINT #1, "pl tano de jard¡n, ";
30610 IF TD > 60 AND TD < 124 THEN PRINT #1, "abedul, ";
30620 IF TD > 60 AND TD < 140 THEN PRINT #1, "cipr‚s, ";
30630 IF TD > 80 AND TD < 153 THEN PRINT #1, "encina, ";
30640 PRINT #1,
30650 PRINT #1, "CANARIAS"
30660 IF TD > 120 AND TD < 183 THEN PRINT #1, "gram¡neas, ";
30670 IF TD > 160 AND TD < 200 THEN PRINT #1, "artemisia, ";
30680 PRINT #1,
30690 PRINT #1, "COMUNIDAD DE MADRID"
30700 IF TD > 122 AND TD < 183 THEN PRINT #1, "olivo, ";
30710 IF TD > 122 AND TD < 204 THEN PRINT #1, "plantago, ";
30720 IF TD > 60 AND TD < 91 THEN PRINT #1, "fresno, ";
30730 IF TD > 150 AND TD < 210 THEN PRINT #1, "gram¡neas, ";
30740 IF TD > 153 AND TD < 204 THEN PRINT #1, "chenopodium, ";
30750 IF TD > 122 AND TD < 153 THEN PRINT #1, "encina, ";
30760 PRINT #1,
30770 PRINT #1, "MURCIA"
30780 IF TD > 130 AND TD < 193 THEN PRINT #1, "gram¡neas, ";
30790 IF TD > 194 AND TD < 290 THEN PRINT #1, "artemisia, ";
30800 IF TD > 60 AND TD < 200 THEN PRINT #1, "parietaria, ";
30810 IF TD > 158 AND TD < 254 THEN PRINT #1, "chenopodium, ";
30820 IF TD > 158 AND TD < 254 THEN PRINT #1, "amaranio, ";
30830 IF TD > 91 AND TD < 122 THEN PRINT #1, "moral, ";
30840 IF TD > 128 AND TD < 188 THEN PRINT #1, "olivo, ";
30850 PRINT #1,
30860 PRINT #1, "COMUNIDAD DE NAVARRA"
30870 IF TD > 130 AND TD < 193 THEN PRINT #1, "gram¡neas, ";
30880 PRINT #1,
30890 PRINT #1, "COMUNIDAD VALENCIANA"
30900 IF TD > 60 AND TD < 270 THEN PRINT #1, "parietaria, ";
30910 IF TD > 91 AND TD < 244 THEN PRINT #1, "gram¡neas, ";
30920 IF TD > 122 AND TD < 183 THEN PRINT #1, "olivo, ";
30930 IF TD > 100 AND TD < 133 THEN PRINT #1, "pl tano de jard¡n, ";
30940 IF TD > 200 AND TD < 270 THEN PRINT #1, "artemisia, ";
30950 IF TD > 153 AND TD < 274 THEN PRINT #1, "amaranto, ";
30960 IF TD > 153 AND TD < 264 THEN PRINT #1, "chenopodium, ";
30970 PRINT #1,
30980 PRINT #1, "EXTREMADURA"
30990 IF TD > 122 AND TD < 183 THEN PRINT #1, "gram¡neas, ";
31000 IF TD > 122 AND TD < 183 THEN PRINT #1, "olivo, ";
31010 IF TD > 122 AND TD < 183 THEN PRINT #1, "encina, ";
31020 IF TD > 163 AND TD < 214 THEN PRINT #1, "plantago, ";
31030 IF TD > 110 AND TD < 234 THEN PRINT #1, "parietaria, ";
31040 PRINT #1,
31050 PRINT #1, "GALICIA"
31060 IF TD > 91 AND TD < 130 THEN PRINT #1, "abedul, ";
31070 IF TD > 153 AND TD < 214 THEN PRINT #1, "gram¡neas, ";
31080 IF TD > 183 AND TD < 274 THEN PRINT #1, "chenopodium, ";
31090 PRINT #1,
31100 PRINT #1, "LA RIOJA"
31110 IF TD > 163 AND TD < 214 THEN PRINT #1, "gram¡neas, ";
31120 PRINT #1,
31130 PRINT #1, "PAIS VASCO"
31140 IF TD > 168 AND TD < 204 THEN PRINT #1, "gram¡neas, ";
31150 IF TD > 153 AND TD < 244 THEN PRINT #1, "plantago, ";
31160 IF TD > 122 AND TD < 168 THEN PRINT #1, "fresno, ";
31170 IF TD > 62 AND TD < 98 THEN PRINT #1, "pino insignis, ";
32000 GOTO 5790