<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=iso-8859-2"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><title>Re: ot: analemma</title><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.StylE-mailovZprvy18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=CS link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Děkuji všem. Pro mcu raději zvolím tabulku a nějakou tu aproximaci. Tabulka je <a href="http://mail.astrohk.cz/~mira/shc/advanced/analema/">http://mail.astrohk.cz/~mira/shc/advanced/analema/</a> viz </span><a href="http://mail.astrohk.cz/~mira/shc/advanced/analema/casova_rce.dat">casova_rce.dat</a><o:p></o:p></p><p class=MsoNormal>Míra<o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> hw-list-bounces@list.hw.cz [mailto:hw-list-bounces@list.hw.cz] <b>On Behalf Of </b>Ales Prochaska, Divesoft<br><b>Sent:</b> Monday, August 26, 2013 9:12 PM<br><b>To:</b> HW-news<br><b>Subject:</b> Re: ot: analemma<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p><span style='font-size:10.0pt'>{ Vypocet polohy Slunce }<br><br>uses<br> datim,graf0,m2;<br><br>type<br> real = extended;<br><br>const<br> lat = 50.0875;<br> lon = +14.4714;<br><br>{----------------------------------------------------------------------------}<br>function daynumber (Y,M,D : longint; UT : real) : real;<br>{ Plati pro roky 1901 - 2099 }<br>{----------------------------------------------------------------------------}<br>begin<br> daynumber:=367*Y-int(7*(Y+int((M+9)/12))/4)+int(275*M/9)+D-730530+UT/24;<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>function arcsin (x : real) : real;<br>{ ve stupnich }<br>{----------------------------------------------------------------------------}<br>begin<br> arcsin:=arctan(x/sqrt(1-sqr(x)))/pi*180;<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>function sin (x : real) : real;<br>{ ve stupnich }<br>{----------------------------------------------------------------------------}<br>begin<br> sin:=system.sin(x/180*pi);<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>function cos (x : real) : real;<br>{ ve stupnich }<br>{----------------------------------------------------------------------------}<br>begin<br> cos:=system.cos(x/180*pi);<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>function tan (x : real) : real;<br>{ ve stupnich }<br>{----------------------------------------------------------------------------}<br>begin<br> tan:=sin(x)/cos(x)<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>function arctan (x : real) : real;<br>{ ve stupnich }<br>{----------------------------------------------------------------------------}<br>begin<br> arctan:=system.arctan(x)/pi*180;<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>function atan2 (y,x : real) : real;<br>{ ve stupnich }<br>{----------------------------------------------------------------------------}<br>begin<br> if x>0 then atan2:=arctan(y/x)<br> else if x<0 then atan2:=arctan(y/x)+180<br> else if y>0 then atan2:=90<br> else if y<0 then atan2:=-90<br> else atan2:=0;<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>function norm (x : real) : real;<br>{ normalizuje uhel do intervalu (-180, 180> }<br>{----------------------------------------------------------------------------}<br>begin<br> { jak od blbyho (co byste taky cekali) }<br> while x>180 do x:=x-360;<br> while x<=-180 do x:=x+360;<br> norm:=x;<br> end;<br><br><br><br>{----------------------------------------------------------------------------}<br>procedure poloha ( lat : real; { sirka }<br> lon : real; { delka }<br> year : integer; { rok }<br> month : integer; { mesic }<br> day : integer; { den }<br> UT : real; { cas }<br> var az : real; { azimut }<br> var h : real); { vyska }<br>{ Vypocita azimut a vysku Slunce pro zadane misto v zadanem okamziku. }<br>{ Uhly jsou ve stupnich a desetinach, cas v hodinach a desetinach. }<br>{----------------------------------------------------------------------------}<br>var<br> d : real; { day number }<br> lst : real; { local sidereal time (ve stupnich) }<br> oblecl : real; { obliquity of the ecliptic }<br> w : real; { argument of perihelion }<br> M : real; { mean anomaly (ve stupnich) }<br> e : real; { eccentricity }<br> EA : real; { eccentric anomaly (ve stupnich) }<br> A,B : real; { pracovni promenne }<br> r : real; { distance }<br> v : real; { true anomaly }<br> slon : real; { Sun's true longitude }<br> sRA : real; { Sun's right ascension }<br> sDec : real; { Sun's declination }<br> t_sRA : real; { tan(sRA) }<br> s_sDec : real; { sin(sDec) }<br> HA : real; { hour angle }<br><br>begin {poloha}<br> d:=daynumber(year,month,day,UT);<br> lst:=norm(98.9818+0.985647352*d+UT*15+lon);<br> oblecl:=23.4393-3.563e-7*d;<br> w:=282.9494+4.70935e-5*d;<br> M:=356.0470+0.98566002585*d;<br> e:=0.016709-1.151e-9*d;<br> EA:=M+e*(180/pi)*sin(M)*(1+e*cos(M));<br> A:=cos(EA)-e;<br> B:=sqrt(1-e*e)*sin(EA);<br> r:=sqrt(A*A+B*B);<br> v:=atan2(B,A);<br> slon:=v+w;<br> t_sRA:=(sin(slon)*cos(oblecl))/cos(slon);<br> sRA:=atan2(sin(slon)*cos(oblecl),cos(slon));<br> s_sDec:=sin(oblecl)*sin(slon);<br> sDec:=arcsin(s_sDec);<br> HA:=lst-sRA;<br> az:=atan2(sin(HA),cos(HA)*sin(lat)-tan(sDec)*cos(lat))+180;<br> h:=arcsin(sin(lat)*sin(sDec)+cos(lat)*cos(sDec)*cos(HA));<br> end {poloha};<br><br><br><br>_______________________________________________<br>HW-list mailing list - sponsored by <a href="http://www.HW.cz">www.HW.cz</a><br><a href="mailto:Hw-list@list.hw.cz">Hw-list@list.hw.cz</a><br><a href="http://list.hw.cz/mailman/listinfo/hw-list">http://list.hw.cz/mailman/listinfo/hw-list</a></span> <o:p></o:p></p><p class=MsoNormal><br><br>__________ Informace od ESET NOD32 Antivirus, verze databaze 8730 (20130826) __________<br><br>Tuto zpravu proveril ESET NOD32 Antivirus.<br><br><a href="http://www.eset.cz">http://www.eset.cz</a><o:p></o:p></p></div> <BR>
<BR>
__________ Informace od ESET NOD32 Antivirus, verze databaze 8731 (20130826) __________<BR>
<BR>
Tuto zpravu proveril ESET NOD32 Antivirus.<BR>
<BR>
<A HREF="http://www.eset.cz">http://www.eset.cz</A><BR>
</body></html>