<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=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-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;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
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="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">I made an Ipython notebook to walk through Example 15a in “Astronomical Algorithms”, 2<sup>nd</sup> ed. [1] This was a debugging tool for me, and I hope it’s interesting to you too. I try to do things one step at a time and compare against
the “book answer” at each stage.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The process does expose some ambiguities in the Astronomical Algorithms text, at least to a non-astronomer like me. This mostly involves what units should be assigned to the arbitrary constants. Could I pick your brains on this topic?
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Assuming you said “yes”, there are a number of places where Meeus converts between fractional days and degrees. I’ve been trying to make his “fractional day” variables work out to “u.sday”, converting to degrees where instructed by multiplying
by 360 u.deg/u.sday. All is well and good until we hit “input block 21” on my notebook. He then wants to multiply his fractional day variable by 360.985647. I tried assigning the units u.deg/u.day, but then my sidereal time (theta0 in notebook) is off by a
degree.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I’m guessing that as long as theta0 is OK, all the subsequent calculations will be OK and the unit system will convert hourangle to real time for me. However, this seems an awful lot like subverting the unit system. Anyone have any ideas
what Meeus is doing, and maybe how I could express it better using astropy.units?
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks,<o:p></o:p></p>
<p class="MsoNormal">Bryce<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">[1] <a href="http://nbviewer.ipython.org/github/firelab/met_utils/blob/master/tests/ex15a.ipynb">
http://nbviewer.ipython.org/github/firelab/met_utils/blob/master/tests/ex15a.ipynb</a>
<o:p></o:p></p>
</div>
<br>
<br>
<br>
<br>
This electronic message contains information generated by the USDA solely for the intended recipients. Any unauthorized interception of this message or the use or disclosure of the information it contains may violate the law and subject the violator to civil
or criminal penalties. If you believe you have received this message in error, please notify the sender and delete the email immediately.
</body>
</html>