1 CLS
2 PRINT
3 PRINT "          ********** CALENDARIO Y TEMPERATURAS MEDIAS DE CATALU¥A ********** "
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
100 INPUT "NOMBRE DEL FICHERO PARA VOLCAR (LPT1: PARA IMPRESORA)"; FICH$
102 IF FICH$ = "" THEN PRINT "ERROR": GOTO 110
104 OPEN FICH$ FOR OUTPUT AS #1
106 PRINT "ALMANAQUE, PREVISION DE TEMPERATURAS DE CATALU¥A Y RADIACION ULTRAVIOLETA "
110 INPUT "A¥O, MES:"; A5, M1
120 INPUT "DIA INICIAL, DIA FINAL:"; D1, DI1
130 INPUT "INTERVALO EN DIAS"; INVAL
140 INPUT "LONGITUD GEOCENTRICA LOCAL (GRADOS)"; L2
150 INPUT "LATITUD GEOCENTRICA LOCAL (GRADOS)"; F
160 INPUT "INCREMENTO HORA OFICIAL="; INCHO
170 PRINT #1,
180 PRINT #1,
190 PRINT #1,
200 PRINT #1,
210 PRINT #1,
220 PRINT #1,
230 PRINT #1,
240 PRINT #1, "-----------------------------------------------------------"
250 PRINT #1,
260 PRINT #1, "CALENDARIO"
270 PRINT #1, "----------"
280 PRINT #1,
290 PRINT "CALENDARIO"
300 PRINT "----------"
310 MI5 = M1
320 DI5 = D1
330 AI5 = A5
340 M5 = 1: D5 = 1
350 GOSUB 550
360 JI0 = J0
370 J0 = 0: J2 = 0: M5 = MI5: D5 = DI5: A5 = AI5
380 PI = 3.1415926#
390 RAD = 180 / PI
400 J9 = 9
410 N7 = 7
420 D0 = 0
430 Q7 = 7
440 R9 = 180 / PI
450 S9 = 1
460 N9 = 2
470 X0 = 0
480 GOTO 1910
490 W7 = A9 / 360 - INT(A9 / 360)
500 B2 = 9 + W7
510 B3 = B2 - INT(B2)
520 A1 = B3 * 2 * PI
530 RETURN
540 Q7 = 18
550 IF M5 > 2 THEN 600
560 M6 = 12 + M5
570 A6 = A5 - 1
580 A5 = A6
590 M5 = M6
600 A = INT(A5 / 100)
610 B = 2 - A + INT(A / 4)
620 J0 = INT(365.25 * A5) + INT(30.6001 * (M5 + 1)) + D5 + 1720994.5#
630 IF J0 > 2299170.50001# THEN 660
640 IF Q7 = 18 THEN 680
650 RETURN
660 J0 = J0 + B
670 GOTO 640
680 PRINT "JD="; J0
690 END
700 PRINT "FECHA JULIANA=";
710 INPUT J0
720 J6 = J0 + .5
730 J5 = INT(J6)
740 A = INT((J5 - 1867216.25#) / (36524.25))
750 A1 = J5 + 1 + A - INT(A / 4)
760 IF J5 < 2299161.001# THEN 890
770 B = A1 + 1524
780 C = INT((B - 122.1) / (365.25))
790 D = INT(365.25 * C)
800 E = INT((B - D) / (30.6001))
810 D5 = B - D - INT(30.6001 * E) + (J6 - J5)
820 IF E > 13.5 THEN 850
830 M5 = E - 1
840 GOTO 860
850 M5 = E - 13
860 IF M5 > 2.5 THEN 920
870 A5 = C - 4715
880 GOTO 930
890 A1 = J5
900 GOTO 770
910 GOTO 930
920 A5 = C - 4716
930 PRINT #1, "A¤o"; A5; ", mes"; M5; "dia"; D5
940 RETURN
950 V7 = R
960 A5 = B
970 GOSUB 1030
980 Z = Y
990 V7 = X
1000 A5 = A
1010 GOSUB 1030
1020 RETURN
1030 R9 = 180 / PI
1040 A8 = (A5 / 360 - INT(A5 / 360)) * 360
1050 C1 = INT(A8 / 90) + 1
1060 A2 = (A8 / 90 - INT(A8 / 90)) * 90
1070 X1 = COS(A2 / R9) * V7
1080 Y1 = SIN(A2 / R9) * V7
1090 IF C1 > 3.99999 THEN 1500
1100 IF C1 > 2.99999 THEN 1530
1110 IF C1 > 1.99999 THEN 1560
1120 GOTO 1590
1130 X1 = ABS(X)
1140 Y1 = ABS(Y)
1150 V7 = SQR(X1 * X1 + Y1 * Y1)
1160 A = ATN(Y / X) * 180 / PI
1170 IF X > -1E-08 THEN 1210
1180 IF X < 0 THEN 1190
1190 A7 = 180 + A
1200 RETURN
1210 A7 = A
1220 RETURN
1230 PRINT "X=";
1240 INPUT X9
1250 PRINT "Y=";
1260 INPUT Y9
1270 PRINT "Z=";
1280 INPUT Z9
1290 GOSUB 1320
1300 PRINT "A"; A3; "B"; B7; "R"; V7
1310 END
1320 X = X9
1330 Y = Y9
1340 GOSUB 1130
1350 A3 = A7
1360 X = V7
1370 Y = Z9
1380 GOSUB 1130
1390 B7 = A7
1400 RETURN
1410 PRINT "DE POLARES A CARTESIANAS:A=";
1420 INPUT A
1430 PRINT "B=";
1440 INPUT B
1450 PRINT "R=";
1460 INPUT R
1470 GOSUB 950
1480 PRINT "X="; X; "Y="; Y; "Z="; Z
1490 END
1500 X = Y1
1510 Y = X1 * -1
1520 RETURN
1530 X = X1 * -1
1540 Y = Y1 * -1
1550 RETURN
1560 X = Y1 * -1
1570 Y = X1
1580 RETURN
1590 X = X1
1600 Y = Y1
1610 RETURN
1620 PRINT "DE GRADOS A GMS: GRD=";
1630 INPUT G7
1640 IF G7 < 0 THEN 1740
1650 G = INT(G7)
1660 M6 = (G7 - G) * 60
1670 M = INT(M6)
1680 S = (M6 - M) * 60
1690 G9 = G + M / 100 + S / 10000
1700 RETURN
1710 GOSUB 1620
1720 PRINT "GRD="; G7; "G="; G; "M="; M; "S="; S; "GMS="; G9
1730 END
1740 G7 = G7 + 360
1750 GOSUB 1650
1760 G9 = -1 * G9
1770 RETURN
1780 PRINT "G=";
1790 INPUT G
1800 PRINT "M=";
1810 INPUT M
1820 PRINT "S=";
1830 INPUT S
1840 G7 = G + M / 60 + S / 3600
1850 G9 = G + M / 100 + S / 10000
1860 RETURN
1870 PRINT "DE GMS A GRADOS";
1880 GOSUB 1780
1890 PRINT "GMS="; G9; "GRD="; G7
1900 END
1910 REM 7
1920 J1 = 1
1930 J2 = 0
1940 J8 = INCHO
1950 J4 = 1
1960 GOSUB 550
1970 IF J4 > J1 THEN 5520
1980 J4 = J4 + 1
1990 J7 = 0
2000 GOTO 3630
2010 R9 = 180 / PI
2020 L0 = L1 * R9
2030 GOTO 2810
2040 GOSUB 2060
2050 RETURN
2060 T = (J0 - 2415020!) / 36525!
2070 T3 = T
2080 L7 = 270.4344 + 481267.8831# * T - .0001133 * T * T
2090 L6 = 279.69668# + 36000.76892# * T + .0003025 * T * T
2100 A9 = L6
2110 GOSUB 490
2120 L6 = A1
2130 M6 = 358.47583# + 35999.0475# * T - .00015 * T * T
2140 A9 = M6
2150 GOSUB 490
2160 M6 = A1
2170 M7 = 296.1046 + 477198.8491# * T + .009192 * T * T
2180 C5 = 1.91946 - .004789 * T - .000014 * T * T
2190 C6 = PI / 180 * (C5 * SIN(M6) + (.020094 - .0001 * T) * SIN(2 * M6) + .000293 * SIN(3 * M6))
2200 N0 = 259.1833 - 1934.142 * T + .002078 * T * T
2210 L1 = C6 + L6
2220 V6 = M6 + C6
2230 E7 = 1.675104E-02 - .0000418 * T
2240 R6 = ((1.0000002# * (1 - E7 ^ 2)) / (1 + E7 * COS(V6)))
2250 N6 = 259.13 - 1934.142 * T
2260 A9 = N6
2270 GOSUB 490
2280 N6 = A1
2290 W = L1 - PI / 180 * (.00569 - .00479 * SIN(N6))
2300 N1 = -(17.2327 + .01737 * T) * SIN(N0 / R9) - 1.273 * SIN(2 * L6 / R9)
2310 N2 = 9.21 * COS(N0 / R9) + .552 * COS(2 * L6 / R9)
2320 O1 = 23.452294# - .0130125 * T + .00256 * COS(N6)
2330 O5 = 23.452294# - .0130125 * T
2340 L5 = L1 - PI / 180 * 1.396 * (T - .5)
2350 X1 = R6 * COS(W)
2360 X5 = R6 * COS(L5)
2370 Y7 = R6 * SIN(W) * COS(O1 / (180 / PI))
2380 Y5 = R6 * SIN(L5) * COS(O5 / (180 / PI))
2390 Z3 = R6 * SIN(W) * SIN(O1 / (180 / PI))
2400 Z5 = R6 * SIN(L5) * SIN(O5 / (180 / PI))
2410 R9 = 180 / PI
2420 V7 = R6
2430 A5 = L5 * R9
2440 GOSUB 1030
2450 X2 = X
2460 Y2 = Y
2470 V7 = R6
2480 A5 = W * R9
2490 GOSUB 1030
2500 X3 = X
2510 Y3 = Y
2520 X9 = X5
2530 Y9 = Y5
2540 Z9 = Z5
2550 GOSUB 1320
2560 W5 = A3
2570 D5 = B7
2580 X9 = X3
2590 Y9 = Y7
2600 Z9 = Z3
2610 GOSUB 1320
2620 W6 = A3
2630 D9 = B7
2640 DS9 = D9
2650 P1 = 8.793999 / R6
2660 S0 = .267 / R6
2670 G7 = W5 / 15
2680 GOSUB 1640
2690 W0 = G9
2700 G7 = W6 / 15
2710 GOSUB 1640
2720 W8 = G9
2730 G7 = D5
2740 GOSUB 1640
2750 D6 = G9
2760 G7 = D9
2770 GOSUB 1640
2780 D0 = G9
2790 RETURN
2800 END
2810 REM 1
2820 T1 = 99.691 + 36000.769# * T3
2830 T2 = (T1 - A3 - L2) / 360
2840 T2 = 24 - 24 * (T2 - INT(T2))
2850 TI2 = T2
2860 S = (Z + F + D9) / 2
2870 T9 = SIN((S - F) / R9) * SIN((S - D9) / R9)
2880 T8 = COS((S - Z) / R9) * COS(S / R9)
2890 T9 = SQR(ABS(T9 / T8)) * (T9 / T8) / ABS(T9 / T8)
2900 T9 = R9 * 2 * ATN(T9) / 15
2910 T7 = T2 - T9
2920 T8 = T2 + T9
2930 RETURN
2940 T7 = INT(T7) + (T7 - INT(T7)) * .6
2950 T8 = INT(T8) + (T8 - INT(T8)) * .6
2960 T2 = INT(T2) + (T2 - INT(T2)) * .6
2970 REM 2
2980 T6 = T + J7 * 1.0027 / 36525! / 24
2990 L9 = (270.434164# + 481267.8831# * T6) / RAD
3000 M8 = (358.475833# + 35999.0498# * T6) / RAD
3010 M9 = (296.10468# + 477198.8491# * T6) / RAD
3020 D8 = (350.737486# + 445267.1142# * T6) / RAD
3030 F8 = (11.250889# + 483202.0251# * T6) / RAD
3040 O8 = (259.183275# - 1934.142 * T6) / RAD
3050 E8 = 1 - .002495 * T6
3060 L7 = L9 + 6.28875 / RAD * SIN(M9) + 1.274018 / RAD * SIN(2 * D8 - M9) + .658309 / RAD * SIN(2 * D8)
3070 L7 = L7 + .213616 / RAD * SIN(2 * M9) - .185596 / RAD * SIN(M8) - .114336 / RAD * SIN(2 * F8)
3080 L7 = L7 + .05879 / RAD * SIN(2 * D8 - 2 * M9) + .0572 / RAD * SIN(2 * D8 - M8 - M9)
3090 L7 = L7 + .05332 / RAD * SIN(2 * D8 + M9) + .045874 / RAD * SIN(2 * D8 - M8)
3100 L7 = L7 + .041024 / RAD * SIN(M9 - M8) - .03472 / RAD * SIN(D8) - .030465 / RAD * SIN(M8 + M9)
3110 L7 = L7 / (360 / RAD)
3120 L7 = 360 / RAD * (L7 - INT(L7))
3130 B7 = 5.128189 / RAD * SIN(F8) + .280606 / RAD * SIN(M9 + F8) + .277693 / RAD * SIN(M9 - F8)
3140 B7 = B7 + .17324 / RAD * SIN(2 * D8 - F8) + .055413 / RAD * SIN(2 * D8 + F8 - M9) + .04627 / RAD * SIN(2 * D8 - F8 - M9)
3150 B7 = B7 + .03257 / RAD * SIN(2 * D8 + F8) + .0172 / RAD * SIN(2 * M9 + F8)
3160 C9 = COS(L7 - L1) * COS(B7)
3170 C8 = ATN(SQR(1 - C9 * C9) / C9)
3180 IF C9 > 0 THEN 3200
3190 C8 = 180 / RAD - C8
3200 I = 180 - C8 * RAD - .1468 * (1 - .0549 * SIN(M8) / (1 - .0167 * SIN(M7)) * SIN(C8))
3210 K = (1 + COS(I / RAD)) / 2
3220 P7 = .95072 / RAD + .05182 / RAD * COS(M9) + 9.530001E-03 / RAD * COS(2 * D8 - M9) + .00784 / RAD * COS(2 * D8)
3230 P7 = P7 + .002824 / RAD * COS(2 * M9)
3240 O3 = O1
3250 T7 = 99.691 + 36000.769# * T - L2 + J7 * 1.0027 * 15
3260 N = COS(L7) * COS(B7) - SIN(P7) * COS(F / RAD) * COS(T7 / RAD)
3270 L4 = (SIN(L7) * COS(B7) - SIN(P7) * (SIN(F / RAD) * SIN(O3 / RAD) + COS(F / RAD) * COS(O3 / RAD) * SIN(T7 / RAD))) / N
3280 L4 = ATN(L4) * RAD
3290 IF ABS(L4 - (L7 * RAD)) < 90 THEN 3320
3300 IF ABS(L4 - (L7 * RAD)) > 270 THEN 3320
3310 L4 = 180 + L4
3320 B4 = (COS(L4 / RAD) * (SIN(B7) - SIN(P7) * (SIN(F / RAD) * COS(O3 / RAD) - COS(F / RAD) * SIN(O3 / RAD) * SIN(T7 / RAD)))) / N
3330 B4 = ATN(B4) * RAD
3340 IF L4 > 0 THEN 3360
3350 L4 = 360 + L4
3360 A4 = (SIN(L4 / RAD) * COS(O3 / RAD) - TAN(B4 / RAD) * SIN(O3 / RAD)) / COS(L4 / RAD)
3370 A4 = ATN(A4) * RAD
3380 IF ABS(A4 - L4) < 90 THEN 3410
3390 IF ABS(A4 - L4) > 270 THEN 3410
3400 A4 = A4 + 180
3410 D3 = SIN(B4 / RAD) * COS(O3 / RAD) + COS(B4 / RAD) * SIN(O3 / RAD) * SIN(L4 / RAD)
3420 D4 = ATN(D3 / SQR(1 - D3 * D3)) * RAD
3430 IF A4 > 0 THEN 3450
3440 A4 = 360 + A4
3450 A3 = A4
3460 D9 = D4
3470 A7 = (SIN(L7) * COS(O3) - TAN(B7) * SIN(O3)) / COS(L7)
3480 A7 = ATN(A7) * RAD
3490 IF ABS(A7 - L7) < 90 THEN 3520
3500 IF ABS(A7 - L7) > 270 THEN 3520
3510 A7 = A7 + 180
3520 D3 = SIN(B7) * COS(O3) + COS(B7) * SIN(O3) * SIN(L7)
3530 D7 = ATN(D3 / SQR(1 - D3 * D3)) * RAD
3540 IF A7 > 0 THEN 3560
3550 A7 = 360 + A7
3560 A7 = A7 / 15
3570 A7 = INT(A7) + (A7 - INT(A7)) * .6
3580 A4 = A4 / 15
3590 A4 = INT(A4) + (A4 - INT(A4)) * .6
3600 RETURN
3610 PRINT A7; A4; D7; D4; "AR. Y DEC. GEO. Y TOPC."
3620 END
3630 GOSUB 2060
3640 Z = 108
3650 GOSUB 2810
3660 U1 = T7
3670 U2 = T8
3680 Z = 90.85
3690 GOSUB 2810
3700 U3 = T7
3710 U4 = T2
3720 U5 = T8
3730 GOTO 3860
3740 T = T3 + U3 / 24 / 36525!
3750 GOSUB 2080
3760 GOSUB 2810
3770 U3 = T7
3780 T = T3 + U4 / 24 / 36525!
3790 GOSUB 2080
3800 GOSUB 2810
3810 U4 = T2
3820 T = T3 + U5 / 24 / 36525!
3830 GOSUB 2080
3840 GOSUB 2810
3850 U5 = T8
3860 T = T3
3870 J7 = 0
3880 U9 = 0
3890 GOSUB 2970
3900 K9 = K
3910 GOSUB 2810
3920 U6 = T7
3930 U7 = T2
3940 U8 = T8
3950 J7 = U6
3960 GOSUB 2970
3970 GOSUB 2810
3980 U6 = T7
3990 J7 = U7
4000 GOSUB 2970
4010 GOSUB 2810
4020 U7 = T2
4030 J7 = T8
4040 GOSUB 2970
4050 GOSUB 2810
4060 U8 = T8
4070 U9 = U9 + 1
4080 IF U9 > 1 THEN 4100
4090 GOTO 3950
4100 U15 = U5: U13 = U3: U11 = U1
4110 U18 = U8: U16 = U6
4120 UC1 = U3 - U1
4130 U1 = INT(U1) + (U1 - INT(U1)) * .6 + J8
4140 U2 = INT(U2) + (U2 - INT(U2)) * .6 + J8
4150 U3 = INT(U3) + (U3 - INT(U3)) * .6 + J8
4160 U4 = INT(U4) + (U4 - INT(U4)) * .6 + J8
4170 U5 = INT(U5) + (U5 - INT(U5)) * .6 + J8
4180 U6 = U6 + J8
4190 IF U6 > 0 THEN 4230
4200 U6 = U6 + 24
4210 U6 = -1 * (INT(U6) + (U6 - INT(U6)) * .6)
4220 GOTO 4240
4230 U6 = INT(U6) + (U6 - INT(U6)) * .6
4240 U7 = INT(U7) + (U7 - INT(U7)) * .6 + J8
4250 U8 = INT(U8) + (U8 - INT(U8)) * .6 + J8
4260 GOSUB 720
4270 J9 = A5 + M5 / 100 + D5 / 10000
4280 J0 = J0 + J2
4290 DS = INT((J0 / 7 - INT(J0 / 7)) * 7)
4300 PRINT #1,
4310 GOSUB 5530
4320 PRINT #1,
4330 DB = J0 - JI0 + 1
4340 IF DS < 7 THEN 4370
4350 IF DS = 7 THEN DS1 = 1
4360 GOTO 4380
4370 DS1 = DS + 1
4380 NS = INT((DB - DS1) / 7 + 1.5)
4390 PRINT "Semana del a¤o n£mero "; NS
4400 PRINT #1, "Semana del a¤o n£mero "; NS
4410 PRINT #1,
4420 PRINT "Dias transcurridos desde el inicio del a¤o"; DB
4430 PRINT #1, "D¡as transcurridos desde el inicio del a¤o:"; DB
4440 PRINT #1,
4450 DA = 365
4460 IF (A5 / 4 - INT(A5 / 4)) = 0 THEN DA = 366
4470 PRINT "Dias que faltan hasta fin de a¤o"; DA - DB
4480 PRINT #1, "D¡as que faltan hasta fin de a¤o:"; DA - DB
4490 PRINT #1,
4590 PRINT "D¡a juliano a las 12 horas:"; J0
4600 PRINT #1, "D¡a juliano a las 12 horas:"; J0
4610 PRINT #1,
4620 TS1 = T1 / 360
4630 TS1 = 24 * (TS1 - INT(TS1))
4640 PRINT #1, "Tiempo sid‚reo de Greenwich (0h):"; INT(TS1); " horas "; INT(60 * (TS1 - (INT(TS1)))); " minutos"
4650 PRINT "Tiempo sid‚reo de Greenwich (0h):"; INT(TS1); " horas "; INT(60 * (TS1 - (INT(TS1)))); " minutos"
4660 PRINT #1,
4670 TS2 = (T1 + L2) / 360
4680 TS2 = 24 * (TS2 - INT(TS2))
4690 PRINT #1, "Tiempo sid‚reo local (0h):"; INT(TS2); " horas "; INT(60 * (TS2 - (INT(TS2)))); " minutos"
4700 PRINT "Tiempo sid‚reo local (0h):"; INT(TS2); " horas "; INT(60 * (TS2 - (INT(TS2)))); " minutos"
4710 PRINT #1,
4720 HS = U3 - U11 - J8
4730 PRINT "SOL Y LUNA"
4740 PRINT #1, "SOL Y LUNA"
4750 PRINT #1, "----------"
4760 EC1 = INT(UC1)
4770 FC1 = INT((UC1 - EC1) * 60)
4780 PRINT "La duraci¢n del crep£sculo es de "; EC1; " horas,"; FC1; " minutos."
4790 PRINT #1, "La duraci¢n del crep£sculo es de "; EC1; " horas,"; FC1; " minutos."
4800 PRINT #1,
4810 PRINT #1, "El crep£sculo matutino comienza a las "; INT(U1); " horas,"; INT((U1 - INT(U1)) * 100); " minutos."
4820 PRINT "El crep£sculo matutino comienza a las "; INT(U1); " horas,"; INT((U1 - INT(U1)) * 100); " minutos."
4830 PRINT #1,
4840 PRINT #1, "El orto solar comienza a las "; INT(U3); " horas,"; INT((U3 - INT(U3)) * 100); " minutos."
4850 PRINT "El orto solar comienza a las "; INT(U3); " horas,"; INT((U3 - INT(U3)) * 100); " minutos."
4860 PRINT #1,
4870 PRINT #1, "El ocaso solar finaliza a las "; INT(U5); " horas,"; INT((U5 - INT(U5)) * 100); " minutos."
4880 PRINT "El ocaso solar finaliza a las "; INT(U5); " horas,"; INT((U5 - INT(U5)) * 100); " minutos."
4890 PRINT #1,
4900 PRINT #1, "El crep£sculo vespertino finaliza a las "; INT(U2); " horas,"; INT((U2 - INT(U2)) * 100); " minutos."
4910 PRINT "El crep£sculo vespertino finaliza a las "; INT(U2); " horas,"; INT((U2 - INT(U2)) * 100); " minutos."
4920 PRINT #1,
4930 HS = U15 - U13
4940 PRINT "La duracion del dia solar es de "; INT(HS); "horas"; INT((HS - INT(HS)) * 60); "minutos"
4950 PRINT #1, "La duraci¢n del d¡a solar es de "; INT(HS); "horas"; INT((HS - INT(HS)) * 60); "minutos"
4960 PRINT #1,
4970 UR4 = INT(U4): UR5 = INT((U4 - INT(U4)) * 100)
4980 PRINT "El mediodia solar tendra lugar a las"; UR4; "horas"; UR5; "minutos"
4990 PRINT #1, "El mediod¡a solar tendr  lugar a las"; UR4; "horas"; UR5; "minutos"
5000 PRINT #1,
5010 AB = INT(((90 - 41.23) + DS9) + .5)
5020 PRINT "La m xima altura que alcanzara hoy el Sol sobre el horizonte sera de"; AB; "grados"
5030 PRINT #1, "La m xima altura que alcanzar  hoy el Sol sobre el horizonte ser  de"; AB; "grados"
5040 PRINT #1,
5050 GOTO 5370
5060 GOSUB 5930
5070 ME = M5 + D5 / 31 - 6
5080 IF ME > 6 THEN ME = ME - 12
5090 TL1 = ME * 3.01201152# + ME * ME * -.86100884# + ME * ME * ME * -.151525 + ME * ME * ME * ME * .010679761415# + ME * ME * ME * ME * ME * .0018947964# + 22.6
5100 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
5110 TL3 = ME * 3.0792 + ME * ME * -.78842 + ME * ME * ME * -.1635751# + ME * ME * ME * ME * .01112967944# + ME * ME * ME * ME * ME * .002171003# + 15.5579
5120 TB1 = ME * 2.9624 + ME * ME * -.5998 + ME * ME * ME * -.1586 + ME * ME * ME * ME * .00786 + ME * ME * ME * ME * ME * .0021 + 21.65
5130 TB2 = ME * 2.9137 + ME * ME * -.62315 + ME * ME * ME * -.1583 + ME * ME * ME * ME * .00787 + ME * ME * ME * ME * ME * .0021 + 25.21
5140 TB3 = ME * 2.9902 + ME * ME * -.57012 + ME * ME * ME * -.1576 + ME * ME * ME * ME * .00768 + ME * ME * ME * ME * ME * .00208 + 18.02
5150 TT1 = ME * 2.9076 + ME * ME * -.5222701 + ME * ME * ME * -.1479 + ME * ME * ME * ME * .00654 + ME * ME * ME * ME * ME * .0018 + 20.46
5160 TT2 = ME * 2.7891 + ME * ME * -.46717 + ME * ME * ME * -.1379 + ME * ME * ME * ME * .00545 + ME * ME * ME * ME * ME * .0016 + 23.59
5170 TT3 = ME * 3.0469 + ME * ME * -.57091 + ME * ME * ME * -.1607 + ME * ME * ME * ME * .00742 + ME * ME * ME * ME * ME * .0021 + 17.27
5180 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
5190 TG2 = ME * 3.0689 + ME * ME * -.68772 + ME * ME * ME * -.1711 + ME * ME * ME * ME * .00837 + ME * ME * ME * ME * ME * .0024 + 26.83
5200 TG3 = ME * 2.9784 + ME * ME * -.64505 + ME * ME * ME * -.1501 + ME * ME * ME * ME * .00854 + ME * ME * ME * ME * ME * .0018 + 14.93
5210 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."
5220 PRINT #1, "La temperatura media que corresponde a la fecha de hoy en Barcelona es de "; INT(TB1)
5230 PRINT #1,
5240 PRINT #1, "grados, en Tarragona es de "; INT(TT1); ",en Lleida "; INT(TL1); ", en Girona"; INT(TG1)
5250 PRINT #1,
5260 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"
5270 PRINT #1, "La m xima media para hoy en Barcelona es de"; INT(TB2); ",en Tarragona es de"; INT(TT2)
5280 PRINT #1,
5290 PRINT #1, ",en Lleida es de "; INT(TL2); ", en Girona es de "; INT(TG2); "grados"
5300 PRINT #1,
5310 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"
5320 PRINT #1, "En cuanto a la m¡nima media en Barcelona es de"; INT(TB3); ",en Tarragona es de"; INT(TT3)
5330 PRINT #1,
5340 PRINT #1, ",en Lleida es de "; INT(TL3); ", en Girona es de "; INT(TG3); "grados"
5350 PRINT #1,
5360 GOTO 5900
5370 LET AS0 = U7
5380 GOSUB 5690
5390 UR7 = INT(U7): UR8 = INT((U7 - INT(U7)) * 100)
5400 HL = U18 - U16
5410 PRINT "La duracion del dia lunar es de "; INT(HL); "horas"; INT((HL - INT(HL)) * 60); "minutos"
5420 PRINT #1, "La duraci¢n del d¡a lunar es de "; INT(HL); "horas"; INT((HL - INT(HL)) * 60); "minutos"
5430 PRINT #1,
5440 PRINT "El mediodia lunar a las "; UR7; "horas"; UR8; "minutos en la constelacion de "; C$
5450 PRINT #1, "El mediod¡a lunar a las "; UR7; "horas"; UR8; "minutos en la constelaci¢n de "; C$
5460 PRINT #1,
5470 GOSUB 5830
5480 PRINT "La Luna en "; R$; "con un disco iluminado del "; INT(K9 * 100); "por ciento"
5490 PRINT #1, "La Luna en "; R$; "con un disco iluminado del "; INT(K9 * 100); "por ciento"
5500 PRINT #1,
5510 GOTO 5060
5520 END
5530 DS = INT(DS + 1)
5540 IF DS = 7 THEN PRINT "Lunes"
5550 IF DS = 7 THEN PRINT #1, "Lunes"
5560 IF DS = 1 THEN PRINT "Martes"
5570 IF DS = 1 THEN PRINT #1, "Martes"
5580 IF DS = 2 THEN PRINT "Mi‚rcoles"
5590 IF DS = 2 THEN PRINT #1, "Mi‚rcoles"
5600 IF DS = 3 THEN PRINT "Jueves"
5610 IF DS = 3 THEN PRINT #1, "Jueves"
5620 IF DS = 4 THEN PRINT "Viernes"
5630 IF DS = 4 THEN PRINT #1, "Viernes"
5640 IF DS = 5 THEN PRINT "S bado"
5650 IF DS = 5 THEN PRINT #1, "S bado"
5660 IF DS = 6 THEN PRINT "Domingo"
5670 IF DS = 6 THEN PRINT #1, "Domingo"
5680 RETURN
5690 IF A4 > 1.4 AND A4 < 3.2 THEN C$ = "Aries"
5700 IF A4 > 3.2 AND A4 < 5.5 THEN C$ = "Taurus"
5710 IF A4 > 5.5 AND A4 < 8 THEN C$ = "Geminis"
5720 IF A4 > 8 AND A4 < 9.2 THEN C$ = "Cancer"
5730 IF A4 > 9.2 AND A4 < 11.35 THEN C$ = "Leo"
5740 IF A4 > 11.35 AND A4 < 14.2 THEN C$ = "Virgo"
5750 IF A4 > 14.2 AND A4 < 15.4 THEN C$ = "Libra"
5760 IF A4 > 15.4 AND A4 < 16.4 THEN C$ = "Escorpio"
5770 IF A4 > 16.4 AND A4 < 17.4 THEN C$ = "Ophiuco"
5780 IF A4 > 17.4 AND A4 < 20! THEN C$ = "Sagitario"
5790 IF A4 > 20! AND A4 < 22! THEN C$ = "Capricornio"
5800 IF A4 > 22! AND A4 < 24! THEN C$ = "Acuario"
5810 IF A4 > 0 AND A4 < 1.4 THEN C$ = "Piscis"
5820 RETURN
5830 IF U7 > U4 AND (U7 - U4) < 12 THEN R$ = "creciente "
5840 IF U7 < U4 THEN R$ = "menguante "
5850 IF K9 < .06 THEN R$ = "novilunio "
5860 IF K9 > .94 THEN R$ = "plenilunio "
5870 RETURN
5880 GOSUB 5930
5890 D1 = D1 + 1
5900 D1 = D1 + INVAL
5910 IF D1 > DI1 THEN END
5920 GOTO 170
5930 PRINT #1,
5940 PRINT #1, "IRRADIACION SOLAR"
5950 W = 20 / (COS(AB / (180 / 3.1416)))
5960 PRINT #1, "Irradiancia ultravioleta previsible en W/m*m ="; W
5970 PRINT #1,
5980 RETURN
5990 PRINT #1, "Tiempos m ximos en minutos recomendados de exposici¢n al Sol"
6000 PRINT #1, "------------------------------------------------------------"
6010 PRINT #1, "Pieles muy blancas: Primeras exposiciones"
6020 PRINT #1, "Ma¤anas:"; INT(325 / W * 2)
6030 PRINT #1, "Mediod¡a:"; INT(325 / W)
6040 PRINT #1, "Tardes:"; INT(325 / W * 2)
6050 PRINT #1,
6060 PRINT #1, "Pieles que pigmentan ligeramente: Primeras exposiciones"
6070 PRINT #1, "Ma¤anas:"; INT(325 / W * 4)
6080 PRINT #1, "Mediod¡a:"; INT(325 / W * 2)
6090 PRINT #1, "Tardes:"; INT(325 / W * 4)
6100 PRINT #1,
6110 PRINT #1, "Pigmentaci¢n normal: Primeras exposiciones"
6120 PRINT #1, "Ma¤anas:"; INT(325 / W * 6)
6130 PRINT #1, "Mediod¡a:"; INT(325 / W * 3)
6140 PRINT #1, "Tardes:"; INT(325 / W * 6)
6150 PRINT #1,
6160 PRINT #1, "Pieles de muy f cil pigmentaci¢n: Primeras exposiciones"
6170 PRINT #1, "Ma¤anas:"; INT(325 / W * 8)
6180 PRINT #1, "Mediod¡a:"; INT(325 / W * 4)
6190 PRINT #1, "Tardes:"; INT(325 / W * 8)
6200 PRINT #1,
6210 PRINT #1, "Pieles muy blancas: Exposiciones posteriores"
6220 PRINT #1, "Ma¤anas:"; INT(325 / W * 4)
6230 PRINT #1, "Mediod¡a:"; INT(325 / W * 2)
6240 PRINT #1, "Tardes:"; INT(325 / W * 4)
6250 PRINT #1,
6260 PRINT #1, "Pieles que pigmentan ligeramente: Exposiciones posteriores"
6270 PRINT #1, "Ma¤anas:"; INT(325 / W * 8)
6280 PRINT #1, "Mediod¡a:"; INT(325 / W * 4)
6290 PRINT #1, "Tardes:"; INT(325 / W * 8)
6300 PRINT #1,
6310 PRINT #1, "Pigmentaci¢n normal: Exposiciones posteriores"
6320 PRINT #1, "Ma¤anas:"; INT(325 / W * 12)
6330 PRINT #1, "Mediod¡a:"; INT(325 / W * 6)
6340 PRINT #1, "Tardes:"; INT(325 / W * 12)
6350 PRINT #1,
6360 PRINT #1, "Pieles de muy f cil pigmentaci¢n: Exposiciones posteriores"
6370 PRINT #1, "Ma¤anas:"; INT(325 / W * 16)
6380 PRINT #1, "Mediod¡a:"; INT(325 / W * 8)
6390 PRINT #1, "Tardes:"; INT(325 / W * 16)
6400 PRINT #1,
6410 PRINT #1, "Los tiempos m ximos se prolongan cuando se utilizan cremas protectoras"
6420 RETURN