<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" lang="en" xml:lang="en">
<head>
<title>sunriseset.cpp Source File</title>
<style type="text/css">
.cp-doNotDisplay { display: none; }
@media aural, braille, handheld, tty { .cp-doNotDisplay { display: inline; speak: normal; }}
.cp-edit { text-align: right; }
@media print, embossed { .cp-edit { display: none; }}
</style>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<meta http-equiv="Content-Style-Type" content="text/css" />
<link rel="meta" href="http://www.kde.org/labels.rdf" type="application/rdf+xml" title="ICRA labels" />
<meta name="trademark" content="KDE e.V." />
<meta name="description" content="K Desktop Environment Homepage, KDE.org" />
<meta name="MSSmartTagsPreventParsing" content="true" />
<meta name="robots" content="all" />
<meta name="no-email-collection" content="http://www.unspam.com/noemailcollection" />
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link rel="shortcut icon" href="../../favicon.ico" /><link rel="icon" href="../../favicon.ico" />
<link rel="stylesheet" media="screen" type="text/css" title="KDE Colors" href="../../kde.css" />
<link rel="stylesheet" media="print, embossed" type="text/css" href="../../print.css" />
<link rel="stylesheet" media="screen, aural, handheld, tty, braille" type="text/css" title="Flat" href="../../flat.css" />
<link rel="stylesheet" type="text/css" href="../../doxygen.css" />
<link rel="stylesheet" type="text/css" href="../../tabs.css" />
</head>
<body>
<ul class="cp-doNotDisplay">
<li><a href="#cp-content" accesskey="2">Skip to content</a></li>
<li><a href="#cp-menu" accesskey="5">Skip to link menu</a></li>
</ul>
<div id="container">
<div id="header">
<div id="header_top"><div><div>
<img alt ="" src="../../top-kde.jpg"/>
</div></div></div>
<div id="header_bottom">
<div id="location">
<ul>
<li><a href="http://api.kde.org">KDE API Reference</a></li>
<li><a href="../../index.html">kdepimlibs-4.10.5 API Reference</a></li>
</ul>
</div>
<div id="menu">
<ul>
<li><a href="http://www.kde.org/">KDE Home</a></li>
<li><a href="http://kde.org/contact/">Contact Us</a></li>
</ul>
</div>
</div>
</div>
<!-- End page header -->
<div id="body_wrapper">
<div id="body">
<!-- begin main content -->
<div id="right">
<div class="content">
<div id="main">
<div class="clearer"> </div>
<h2><a name="content"></a>KHolidays Library</h2>
<!-- BC -->
<div id="top">
<!-- Generated by Doxygen 1.8.3.1 -->
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_7d1f8f4788b4980602abcfcd1bfe22cd.html">kholidays</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">sunriseset.cpp</div> </div>
</div><!--header-->
<div class="contents">
<div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"> This file is part of the kholidays library.</span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment"></span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment"> Adapted from the Javascript found at http://www.esrl.noaa.gov/gmd/grad/solcalc</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> which is in the public domain.</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"></span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment"> Copyright (c) 2012 Allen Winter <winter@kde.org></span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment"></span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment"> This library is free software; you can redistribute it and/or</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment"> modify it under the terms of the GNU Library General Public</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment"> License as published by the Free Software Foundation; either</span></div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment"> version 2 of the License, or (at your option) any later version.</span></div>
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment"></span></div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment"> This library is distributed in the hope that it will be useful,</span></div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment"> but WITHOUT ANY WARRANTY; without even the implied warranty of</span></div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment"> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU</span></div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment"> Library General Public License for more details.</span></div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment"></span></div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment"> You should have received a copy of the GNU Library General Public License</span></div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment"> along with this library; see the file COPYING.LIB. If not, write to</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment"> the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,</span></div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="comment"> Boston, MA 02110-1301, USA.</span></div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="comment">*/</span></div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span> </div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="preprocessor">#include "sunriseset.h"</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="preprocessor">#include <cmath></span></div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="preprocessor">#ifndef KDE_USE_FINAL</span></div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="preprocessor"></span><span class="keyword">static</span> <span class="keywordtype">double</span> PI = 3.14159265358979323846;</div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="preprocessor"></span><span class="keyword">static</span> <span class="keywordtype">double</span> MaxLat = 89.99;</div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="keyword">static</span> <span class="keywordtype">double</span> MaxLong = 179.99;</div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> </div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span> <span class="keyword">using namespace </span>KHolidays;</div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="keyword">using namespace </span>SunRiseSet;</div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> </div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="keyword">static</span> <span class="keywordtype">double</span> degToRad( <span class="keywordtype">double</span> degree )</div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span> {</div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  <span class="keywordflow">return</span> ( degree * PI ) / 180.00;</div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> }</div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span> </div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="keyword">static</span> <span class="keywordtype">double</span> radToDeg( <span class="keywordtype">double</span> rad )</div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> {</div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  <span class="keywordflow">return</span> ( rad * 180.0 ) / PI;</div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> }</div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span> </div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="comment">//convert Julian Day to centuries since J2000.0.</span></div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcTimeJulianCent( <span class="keywordtype">int</span> jday )</div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span> {</div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="keywordflow">return</span> ( <span class="keywordtype">double</span>( jday ) - 2451545.0 ) / 36525.0;</div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span> }</div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span> </div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span> <span class="comment">//calculate the mean obliquity of the ecliptic (in degrees)</span></div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcMeanObliquityOfEcliptic( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span> {</div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  <span class="keywordtype">double</span> seconds = 21.448 - t * ( 46.8150 + t * ( 0.00059 - t * 0.001813 ) );</div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  <span class="keywordtype">double</span> e0 = 23.0 + ( 26.0 + ( seconds / 60.0 ) ) / 60.0;</div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  <span class="keywordflow">return</span> e0; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span> }</div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span> </div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span> <span class="comment">//calculate the corrected obliquity of the ecliptic (in degrees)</span></div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcObliquityCorrection( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span> {</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  <span class="keywordtype">double</span> e0 = calcMeanObliquityOfEcliptic( t );</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  <span class="keywordtype">double</span> omega = 125.04 - 1934.136 * t;</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  <span class="keywordtype">double</span> e = e0 + 0.00256 * cos( degToRad( omega ) );</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  <span class="keywordflow">return</span> e; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span> }</div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span> </div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span> <span class="comment">//calculate the Geometric Mean Longitude of the Sun (in degrees)</span></div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcGeomMeanLongSun( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span> {</div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <span class="keywordtype">double</span> L0 = 280.46646 + t * ( 36000.76983 + 0.0003032 * t );</div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  <span class="keywordflow">while</span> ( L0 > 360.0 ) {</div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  L0 -= 360.0;</div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  }</div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  <span class="keywordflow">while</span> ( L0 < 0.0 ) {</div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  L0 += 360.0;</div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  }</div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  <span class="keywordflow">return</span> L0; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span> }</div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span> </div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span> <span class="comment">//calculate the Geometric Mean Anomaly of the Sun (in degrees)</span></div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcGeomMeanAnomalySun( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span> {</div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  <span class="keywordtype">double</span> M = 357.52911 + t * ( 35999.05029 - 0.0001537 * t );</div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <span class="keywordflow">return</span> M; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span> }</div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span> </div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span> <span class="comment">//calculate the equation of center for the sun (in degrees)</span></div>
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcSunEqOfCenter( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00091"></a><span class="lineno"> 91</span> {</div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  <span class="keywordtype">double</span> m = calcGeomMeanAnomalySun( t );</div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <span class="keywordtype">double</span> mrad = degToRad( m );</div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <span class="keywordtype">double</span> sinm = sin( mrad );</div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  <span class="keywordtype">double</span> sin2m = sin( mrad + mrad );</div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <span class="keywordtype">double</span> sin3m = sin( mrad + mrad + mrad );</div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  <span class="keywordtype">double</span> C = sinm * ( 1.914602 - t * ( 0.004817 + 0.000014 * t ) ) +</div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  sin2m * ( 0.019993 - 0.000101 * t ) + sin3m * 0.000289;</div>
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  <span class="keywordflow">return</span> C; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span> }</div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span> </div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span> <span class="comment">//calculate the true longitude of the sun (in degrees)</span></div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcSunTrueLong( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span> {</div>
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  <span class="keywordtype">double</span> l0 = calcGeomMeanLongSun( t );</div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  <span class="keywordtype">double</span> c = calcSunEqOfCenter( t );</div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  <span class="keywordtype">double</span> O = l0 + c;</div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  <span class="keywordflow">return</span> O; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span> }</div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span> </div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span> <span class="comment">//calculate the apparent longitude of the sun (in degrees)</span></div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcSunApparentLong( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span> {</div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  <span class="keywordtype">double</span> o = calcSunTrueLong( t );</div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  <span class="keywordtype">double</span> omega = 125.04 - 1934.136 * t;</div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  <span class="keywordtype">double</span> lambda = o - 0.00569 - 0.00478 * sin( degToRad( omega ) );</div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  <span class="keywordflow">return</span> lambda; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span> }</div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span> </div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span> <span class="comment">//calculate the declination of the sun (in degrees)</span></div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcSunDeclination( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span> {</div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  <span class="keywordtype">double</span> e = calcObliquityCorrection( t );</div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  <span class="keywordtype">double</span> lambda = calcSunApparentLong( t );</div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span> </div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="keywordtype">double</span> sint = sin( degToRad( e ) ) * sin( degToRad( lambda ) );</div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  <span class="keywordtype">double</span> theta = radToDeg( asin( sint ) );</div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  <span class="keywordflow">return</span> theta; <span class="comment">// in degrees</span></div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span> }</div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span> </div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span> <span class="comment">//calculate the eccentricity of earth's orbit (unitless)</span></div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcEccentricityEarthOrbit( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span> {</div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  <span class="keywordtype">double</span> e = 0.016708634 - t * ( 0.000042037 + 0.0000001267 * t );</div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <span class="keywordflow">return</span> e; <span class="comment">// unitless</span></div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span> }</div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span> </div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span> <span class="comment">//calculate the difference between true solar time and mean solar time</span></div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span> <span class="comment">//(output: equation of time in minutes of time)</span></div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcEquationOfTime( <span class="keywordtype">double</span> t )</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span> {</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  <span class="keywordtype">double</span> epsilon = calcObliquityCorrection( t );</div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  <span class="keywordtype">double</span> l0 = calcGeomMeanLongSun( t );</div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  <span class="keywordtype">double</span> e = calcEccentricityEarthOrbit( t );</div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  <span class="keywordtype">double</span> m = calcGeomMeanAnomalySun( t );</div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span> </div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  <span class="keywordtype">double</span> y = tan( degToRad( epsilon ) / 2.0 );</div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  y *= y;</div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span> </div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  <span class="keywordtype">double</span> sin2l0 = sin( 2.0 * degToRad( l0 ) );</div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  <span class="keywordtype">double</span> sinm = sin( degToRad( m ) );</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  <span class="keywordtype">double</span> cos2l0 = cos( 2.0 * degToRad( l0 ) );</div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  <span class="keywordtype">double</span> sin4l0 = sin( 4.0 * degToRad( l0 ) );</div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  <span class="keywordtype">double</span> sin2m = sin( 2.0 * degToRad( m ) );</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span> </div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  <span class="keywordtype">double</span> Etime = y * sin2l0 - 2.0 * e * sinm + 4.0 * e *</div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  y * sinm * cos2l0 - 0.5 *</div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  y * y * sin4l0 - 1.25 * e * e * sin2m;</div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  <span class="keywordflow">return</span> radToDeg( Etime ) * 4.0; <span class="comment">// in minutes of time</span></div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span> }</div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span> </div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span> <span class="comment">//calculate the hour angle of the sun at sunrise for the latitude (in radians)</span></div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span> <span class="keyword">static</span> <span class="keywordtype">double</span> calcHourAngleSunrise( <span class="keywordtype">double</span> latitude, <span class="keywordtype">double</span> solarDec )</div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span> {</div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="keywordtype">double</span> latRad = degToRad( latitude );</div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  <span class="keywordtype">double</span> sdRad = degToRad( solarDec );</div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  <span class="keywordtype">double</span> HAarg = ( cos( degToRad( 90.833 ) ) /</div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  ( cos( latRad ) * cos( sdRad ) ) - tan( latRad ) * tan( sdRad ) );</div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  <span class="keywordtype">double</span> HA = acos( HAarg );</div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>  <span class="keywordflow">return</span> HA; <span class="comment">// in radians (for sunset, use -HA)</span></div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span> }</div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span> </div>
<div class="line"><a name="l00173"></a><span class="lineno"><a class="code" href="namespaceKHolidays_1_1SunRiseSet.html#ab139ae505de3c3b05c77046f06989287"> 173</a></span> QTime <a class="code" href="namespaceKHolidays_1_1SunRiseSet.html#ab139ae505de3c3b05c77046f06989287" title="Compute the sunrise time (UTC) for a date and Earth location.">KHolidays::SunRiseSet::utcSunrise</a>( <span class="keyword">const</span> QDate &date, <span class="keywordtype">double</span> latitude, <span class="keywordtype">double</span> longitude )</div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span> {</div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  latitude = qMax( qMin( MaxLat, latitude ), -MaxLat );</div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  longitude = qMax( qMin( MaxLong, longitude ), -MaxLong );</div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span> </div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  <span class="keywordtype">double</span> t = calcTimeJulianCent( date.toJulianDay() );</div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="keywordtype">double</span> eqTime = calcEquationOfTime( t );</div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <span class="keywordtype">double</span> solarDec = calcSunDeclination( t );</div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  <span class="keywordtype">double</span> hourAngle = calcHourAngleSunrise( latitude, solarDec );</div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  <span class="keywordtype">double</span> delta = longitude + radToDeg( hourAngle );</div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  QTime timeUTC;</div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  timeUTC = timeUTC.addSecs( ( 720 - ( 4.0 * delta ) - eqTime ) * 60 );</div>
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span>  <span class="keywordflow">return</span> QTime( timeUTC.hour(),</div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  timeUTC.second() > 29 ? timeUTC.minute() + 1 : timeUTC.minute(),</div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  0 );</div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span> }</div>
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span> </div>
<div class="line"><a name="l00190"></a><span class="lineno"><a class="code" href="namespaceKHolidays_1_1SunRiseSet.html#a2e5a51eb6fb5138417b88b3188bd9e99"> 190</a></span> QTime <a class="code" href="namespaceKHolidays_1_1SunRiseSet.html#a2e5a51eb6fb5138417b88b3188bd9e99" title="Compute the sunset time (UTC) for a date and Earth location.">KHolidays::SunRiseSet::utcSunset</a>( <span class="keyword">const</span> QDate &date, <span class="keywordtype">double</span> latitude, <span class="keywordtype">double</span> longitude )</div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span> {</div>
<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  latitude = qMax( qMin( MaxLat, latitude ), -MaxLat );</div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  longitude = qMax( qMin( MaxLong, longitude ), -MaxLong );</div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span> </div>
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <span class="keywordtype">double</span> t = calcTimeJulianCent( date.toJulianDay() );</div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="keywordtype">double</span> eqTime = calcEquationOfTime( t );</div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  <span class="keywordtype">double</span> solarDec = calcSunDeclination( t );</div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  <span class="keywordtype">double</span> hourAngle = -calcHourAngleSunrise( latitude, solarDec );</div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  <span class="keywordtype">double</span> delta = longitude + radToDeg( hourAngle );</div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  QTime timeUTC;</div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  timeUTC = timeUTC.addSecs( ( 720 - ( 4.0 * delta ) - eqTime ) * 60 );</div>
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <span class="keywordflow">return</span> QTime( timeUTC.hour(),</div>
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  timeUTC.second() > 29 ? timeUTC.minute() + 1 : timeUTC.minute(),</div>
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  0 );</div>
<div class="line"><a name="l00205"></a><span class="lineno"> 205</span> }</div>
</div><!-- fragment --></div><!-- contents -->
<div class="copyrights">
This file is part of the KDE documentation.<br/>
Documentation copyright © 1996-2013 The KDE developers.<br/>
Generated on Fri Jul 12 2013 16:48:07 by
<a href="http://www.doxygen.org/index.html">doxygen</a> 1.8.3.1 written
by <a href="mailto:dimitri@stack.nl">Dimitri van Heesch</a>, © 1997-2006
<p>
KDE's <a href="http://techbase.kde.org/Policies/Library_Documentation_Policy">Doxygen guidelines</a> are available online.
</p>
</div>
</div>
</div>
</div>
<div id="left">
<div class="menu_box">
<a name="cp-menu" /><div class="menutitle"><div>
<h2 id="cp-menu-project">KHolidays Library</h2>
<a href="#cp-skip-project" class="cp-doNotDisplay">Skip menu "KHolidays Library"</a>
</div></div>
<ul><li><a href="index.html">Main Page</a></li><li><a href="namespaces.html">Namespace List</a></li><li><a href="namespacemembers.html">Namespace Members</a></li><li><a href="classes.html">Alphabetical List</a></li><li><a href="annotated.html">Class List</a></li><li><a href="hierarchy.html">Class Hierarchy</a></li><li><a href="functions.html">Class Members</a></li><li><a href="files.html">File List</a></li><li><a href="pages.html">Related Pages</a></li></ul>
<!--
<h2>Class Picker</h2>
<div style="text-align: center;">
<form name="guideform">
<select name="guidelinks" style="width:100%;" onChange="window.location=document.guideform.guidelinks.options[document.guideform.guidelinks.selectedIndex].value">
<option value="annotated.html">-- Choose --</option>
<option value="classKHolidays_1_1HolidayParserDriver.html">kholidays::holidayparserdriver</option>, <option value="classKHolidays_1_1HolidayParserDriverPlan.html">kholidays::holidayparserdriverplan</option>, <option value="classKHolidays_1_1HolidayParserDriverPlanOld.html">kholidays::holidayparserdriverplanold</option>, <option value="classKHolidays_1_1HolidayParserPlan.html">kholidays::holidayparserplan</option>, <option value="classKHolidays_1_1HolidayScannerPlan.html">kholidays::holidayscannerplan</option>, <option value="classKHolidays_1_1location.html">kholidays::location</option>, <option value="classKHolidays_1_1LunarPhase.html">kholidays::lunarphase</option>, <option value="classKHolidays_1_1position.html">kholidays::position</option>, <option value="classKHolidays_1_1slice.html">kholidays::slice</option>,
</select>
</form>
</div>
-->
<div class="menu_box">
<a name="cp-menu" /><div class="menutitle"><div>
<h2 id="cp-menu-project">kdepimlibs-4.10.5 API Reference</h2>
<a href="#cp-skip-project" class="cp-doNotDisplay">Skip menu "kdepimlibs-4.10.5 API Reference"</a>
</div></div>
<div class="nav_list">
<ul>
<li><a href="../../akonadi/html/index.html">akonadi</a></li><li> <a href="../../akonadi/contact/html/index.html">contact</a></li><li> <a href="../../akonadi/kmime/html/index.html">kmime</a></li><li> <a href="../../akonadi/socialutils/html/index.html">socialutils</a></li><li><a href="../../kabc/html/index.html">kabc</a></li><li><a href="../../kalarmcal/html/index.html">kalarmcal</a></li><li><a href="../../kblog/html/index.html">kblog</a></li><li><a href="../../kcal/html/index.html">kcal</a></li><li><a href="../../kcalcore/html/index.html">kcalcore</a></li><li><a href="../../kcalutils/html/index.html">kcalutils</a></li><li><a href="../../kholidays/html/index.html">kholidays</a></li><li><a href="../../kimap/html/index.html">kimap</a></li><li><a href="../../kioslave/html/index.html">kioslave</a></li><li> <a href="../../kioslave/imap4/html/index.html">imap4</a></li><li> <a href="../../kioslave/mbox/html/index.html">mbox</a></li><li> <a href="../../kioslave/nntp/html/index.html">nntp</a></li><li><a href="../../kldap/html/index.html">kldap</a></li><li><a href="../../kmbox/html/index.html">kmbox</a></li><li><a href="../../kmime/html/index.html">kmime</a></li><li><a href="../../kontactinterface/html/index.html">kontactinterface</a></li><li><a href="../../kpimidentities/html/index.html">kpimidentities</a></li><li><a href="../../kpimtextedit/html/index.html">kpimtextedit</a></li><li><a href="../../kpimutils/html/index.html">kpimutils</a></li><li><a href="../../kresources/html/index.html">kresources</a></li><li><a href="../../ktnef/html/index.html">ktnef</a></li><li><a href="../../kxmlrpcclient/html/index.html">kxmlrpcclient</a></li><li><a href="../../mailtransport/html/index.html">mailtransport</a></li><li><a href="../../microblog/html/index.html">microblog</a></li><li><a href="../../qgpgme/html/index.html">qgpgme</a></li><li><a href="../../syndication/html/index.html">syndication</a></li><li> <a href="../../syndication/atom/html/index.html">atom</a></li><li> <a href="../../syndication/rdf/html/index.html">rdf</a></li><li> <a href="../../syndication/rss2/html/index.html">rss2</a></li>
</ul></div></div>
<!-- api_searchbox -->
</div>
</div>
<div class="clearer"></div>
</div>
<div class="clearer"></div>
</div>
<div id="end_body"></div>
<div id="footer"><div id="footer_text">
Report problems with this website to <a href="https://bugs.kde.org/enter_sysadmin_request.cgi?component=api.kde.org">our bug tracking system</a>.<br>
Contact the specific authors with questions and comments about the page contents.<p>
KDE<sup>®</sup> and <a href="/media/images/kde_gear_black.png">the K Desktop Environment<sup>®</sup> logo</a> are registered trademarks of <a href="http://ev.kde.org/" title="Homepage of the KDE non-profit Organization">KDE e.V.</a> |
<a href="http://www.kde.org/contact/impressum.php">Legal</a></p>
</div></div>
</div>
<!--
WARNING: DO NOT SEND MAIL TO THE FOLLOWING EMAIL ADDRESS! YOU WILL
BE BLOCKED INSTANTLY AND PERMANENTLY!
<a href="mailto:aaaatrap-45abe0e0c3bebc77@kde.org">Block me</a>
WARNING END
-->
</body>
</html>
