17 lines
584 B
Text
17 lines
584 B
Text
10 CLS
|
|
20 DEF FNR(x) = x * PI/180
|
|
30 DEF FND(x) = x * 180/PI
|
|
40 INPUT " Enter latitude (degrees): ";latitude
|
|
50 INPUT " Enter longitude (degrees): ";longitude
|
|
60 INPUT "Enter legal meridian (degrees): ";meridian
|
|
70 PRINT
|
|
80 PRINT " Time Sun hour angle Dial hour line angle"
|
|
90 FOR h = 6 TO 18
|
|
100 hra = 15*h-longitude+meridian-180
|
|
110 hla = FN D(ARCTAN(SIN(FN R(latitude))*TAN(FN R(hra))))
|
|
120 IF ABS(hra) > 90 THEN hla = hla+180*SGN(hra*latitude)
|
|
130 PRINT USING "##.##";h;TAB (10)
|
|
140 PRINT USING "####.###";hra;TAB (27)
|
|
150 PRINT USING "####.###";hla
|
|
160 NEXT h
|
|
170 END
|