Jump to content

Horizontal sundial calculations: Difference between revisions

Horizontal sundial calculations (view source)

Revision as of 15:49, 19 January 2022

8,738 bytes added ,  2 years ago
m
→‎{{header|Phix}}: added syntax colouring, marked p2js compatible
(→‎{{header|AutoHotkey}}: marked incorrect (6AM,5AM..1AM,0PM))
m (→‎{{header|Phix}}: added syntax colouring, marked p2js compatible)
Line 2,169:
 
=={{header|Phix}}==
For better behaviour/proper input on pwa/p2js I suppose this should really be a pGUI app.
Copy of [[Horizontal_sundial_calculations#Euphoria|Euphoria]]
<!--<lang Phix>(phixonline)-->
<lang Phix>atom lat = prompt_number("Enter Latitude: ",{})
<span style="color: #008080;">with</span> <span style="color: #008080;">javascript_semantics</span>
atom lng = prompt_number("Enter Longitude: ",{})
<span style="color: #008080;">function</span> <span style="color: #000000;">prompt</span><span style="color: #0000FF;">(</span><span style="color: #004080;">string</span> <span style="color: #000000;">text</span><span style="color: #0000FF;">,</span> <span style="color: #004080;">atom</span> <span style="color: #000000;">v</span><span style="color: #0000FF;">)</span>
atom lm = prompt_number("Enter Legal Meridian: ",{})
<span style="color: #008080;">if</span> <span style="color: #7060A8;">platform</span><span style="color: #0000FF;">()!=</span><span style="color: #004600;">JS</span> <span style="color: #008080;">then</span>
puts(1,'\n')
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">"Enter %s (%g):"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">text</span><span style="color: #0000FF;">,</span><span style="color: #000000;">v</span><span style="color: #0000FF;">})</span>
<span style="color: #004080;">object</span> <span style="color: #000000;">answer</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">gets</span><span style="color: #0000FF;">(</span><span style="color: #000000;">0</span><span style="color: #0000FF;">)</span>
<span style="color: #7060A8;">puts</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span> <span style="color: #008000;">'\n'</span><span style="color: #0000FF;">)</span>
<span style="color: #008080;">if</span> <span style="color: #004080;">string</span><span style="color: #0000FF;">(</span><span style="color: #000000;">answer</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">then</span>
<span style="color: #000000;">v</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">to_number</span><span style="color: #0000FF;">(</span><span style="color: #7060A8;">trim</span><span style="color: #0000FF;">(</span><span style="color: #000000;">answer</span><span style="color: #0000FF;">),</span><span style="color: #000000;">v</span><span style="color: #0000FF;">)</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
<span style="color: #008080;">else</span>
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">"%s:%g\n"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">text</span><span style="color: #0000FF;">,</span><span style="color: #000000;">v</span><span style="color: #0000FF;">})</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
<span style="color: #008080;">return</span> <span style="color: #000000;">v</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">function</span>
<span style="color: #008080;">constant</span> <span style="color: #000000;">deg2rad</span> <span style="color: #0000FF;">=</span> <span style="color: #004600;">PI</span><span style="color: #0000FF;">/</span><span style="color: #000000;">180</span><span style="color: #0000FF;">,</span>
atom ha, hla
<span style="color: #000000;">rad2deg</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">180</span><span style="color: #0000FF;">/</span><span style="color: #004600;">PI</span>
<span style="color: #004080;">atom</span> <span style="color: #000000;">lat</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">prompt</span><span style="color: #0000FF;">(</span><span style="color: #008000;">"Latitude"</span><span style="color: #0000FF;">,-</span><span style="color: #000000;">4.95</span><span style="color: #0000FF;">),</span>
function Deg2Rad(atom degrees)
<span style="color: #000000;">lng</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">prompt</span><span style="color: #0000FF;">(</span><span style="color: #008000;">"Longitude"</span><span style="color: #0000FF;">,-</span><span style="color: #000000;">150.5</span><span style="color: #0000FF;">),</span>
return degrees * PI / 180
<span style="color: #000000;">mer</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">prompt</span><span style="color: #0000FF;">(</span><span style="color: #008000;">"Legal Meridian"</span><span style="color: #0000FF;">,-</span><span style="color: #000000;">150</span><span style="color: #0000FF;">),</span>
end function
<span style="color: #000000;">slat</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">sin</span><span style="color: #0000FF;">(</span><span style="color: #000000;">lat</span><span style="color: #0000FF;">*</span><span style="color: #000000;">deg2rad</span><span style="color: #0000FF;">)</span>
function Rad2Deg(atom radians)
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">"""
return radians * 180 / PI
end function
Sine of latitude: %g
Longitude offset: %g
function atan2(atom y, atom x)
return 2*arctan((sqrt(power(x,2)+power(y,2)) - x)/y)
Hour Sun hour angle Dial hour line angle
end function
-------- -------------- --------------------
"""</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">slat</span><span style="color: #0000FF;">,</span><span style="color: #000000;">lng</span><span style="color: #0000FF;">-</span><span style="color: #000000;">mer</span><span style="color: #0000FF;">})</span>
atom s_lat = sin(Deg2Rad(lat))
<span style="color: #008080;">for</span> <span style="color: #000000;">hour</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">-</span><span style="color: #000000;">6</span> <span style="color: #008080;">to</span> <span style="color: #000000;">6</span> <span style="color: #008080;">do</span>
puts(1,"Hour, Sun Hour Angle, Dial Hour Line Angle\n")
<span style="color: #004080;">atom</span> <span style="color: #000000;">ha</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">hour</span> <span style="color: #0000FF;">*</span> <span style="color: #000000;">15</span> <span style="color: #0000FF;">-</span> <span style="color: #000000;">lng</span> <span style="color: #0000FF;">+</span> <span style="color: #000000;">mer</span><span style="color: #0000FF;">,</span>
<span style="color: #000000;">s</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">sin</span><span style="color: #0000FF;">(</span><span style="color: #000000;">ha</span><span style="color: #0000FF;">*</span><span style="color: #000000;">deg2rad</span><span style="color: #0000FF;">),</span>
for hour = -6 to 6 do
<span style="color: #000000;">c</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">cos</span><span style="color: #0000FF;">(</span><span style="color: #000000;">ha</span><span style="color: #0000FF;">*</span><span style="color: #000000;">deg2rad</span><span style="color: #0000FF;">),</span>
ha = hour * 15 - lng + lm
<span style="color: #000000;">hla</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">atan2</span><span style="color: #0000FF;">(</span><span style="color: #000000;">slat</span><span style="color: #0000FF;">*</span><span style="color: #000000;">s</span><span style="color: #0000FF;">,</span><span style="color: #000000;">c</span><span style="color: #0000FF;">)*</span><span style="color: #000000;">rad2deg</span>
atom s = sin(Deg2Rad(ha))
<span style="color: #004080;">integer</span> <span style="color: #000000;">hour12</span> <span style="color: #0000FF;">=</span> <span style="color: #008080;">iff</span><span style="color: #0000FF;">(</span><span style="color: #000000;">hour</span><span style="color: #0000FF;"><=</span><span style="color: #000000;">0</span><span style="color: #0000FF;">?</span><span style="color: #000000;">hour</span><span style="color: #0000FF;">+</span><span style="color: #000000;">12</span><span style="color: #0000FF;">:</span><span style="color: #000000;">hour</span><span style="color: #0000FF;">)</span>
atom c = cos(Deg2Rad(ha))
<span style="color: #004080;">string</span> <span style="color: #000000;">am</span> <span style="color: #0000FF;">=</span> <span style="color: #008080;">iff</span><span style="color: #0000FF;">(</span><span style="color: #000000;">hour</span><span style="color: #0000FF;"><</span><span style="color: #000000;">0</span><span style="color: #0000FF;">?</span><span style="color: #008000;">"am"</span><span style="color: #0000FF;">:</span><span style="color: #008000;">"pm"</span><span style="color: #0000FF;">)</span>
hla = Rad2Deg(atan2(s_lat*s,c))
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">"%3d:00%s %12.3f %17.3f\n"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">hour12</span><span style="color: #0000FF;">,</span><span style="color: #000000;">am</span><span style="color: #0000FF;">,</span><span style="color: #000000;">ha</span><span style="color: #0000FF;">,</span><span style="color: #000000;">hla</span><span style="color: #0000FF;">})</span>
printf(1,"%3d %7.3f %7.3f\n",{hour+12,ha,hla})
<span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
end for
<!--</lang>-->
{} = wait_key()</lang>
{{out}}
<pre>
Enter Latitude: (-4.95):
Enter Longitude: (-150.5):
Enter Legal Meridian: (-150):
 
Sine of latitude: -0.086286
Longitude offset: -0.5
 
Hour, Sun Hourhour angle Angle, Dial Hourhour Lineline Angleangle
-------- -------------- --------------------
6 -89.500 84.225▼
7 6:00am -7489.500 17 84.283225
8 7:00am -5974.500 817.334283
9 8:00am -4459.500 4 8.847334
10 9:00am -2944.500 2 4.795847
11 10:00am -1429.500 1 2.278795
1211:00am -14.500 0.500 -01.043278
1312:00pm 150.500 -10.371043
14 1:00pm 30 15.500 -21.910371
15 2:00pm 45 30.500 -52.018910
16 3:00pm 60 45.500 -85.671018
17 4:00pm 75 60.500 -188.451671
18 5:00pm 90 75.500 -9518.775451
▲ 6:00pm -8990.500 84-95.225775
</pre>
 
Petelomax
7,820

edits

Retrieved from "https://rosettacode.org/wiki/Special:MobileDiff/150343"
Cookies help us deliver our services. By using our services, you agree to our use of cookies.
More information