Ysurac / FlightAirMap

newest.php

Braces +1 added lines, -1 removed lines

@@ -18,7 +18,7 @@
   $limit_start = 0;
   $limit_end = 25;
   $absolute_difference = 25;
-}  else {
+} else {
 	$limit_explode = explode(",", $_GET['limit']);
 	$limit_start = $limit_explode[0];
 	$limit_end = $limit_explode[1];

Spacing +5 added lines, -5 removed lines

@@ -6,7 +6,7 @@ 
 
 if (isset($_POST['category']))
 {
-	$category = filter_input(INPUT_POST,'category',FILTER_SANITIZE_STRING);
+	$category = filter_input(INPUT_POST, 'category', FILTER_SANITIZE_STRING);
 	header('Location: '.$globalURL.'/newest/'.$category);
 }
 
@@ -14,18 +14,18 @@ 
 require_once('header.php');
 
 //calculuation for the pagination
-if(!isset($_GET['limit']))
+if (!isset($_GET['limit']))
 {
   $limit_start = 0;
   $limit_end = 25;
   $absolute_difference = 25;
-}  else {
+} else {
 	$limit_explode = explode(",", $_GET['limit']);
 	$limit_start = $limit_explode[0];
 	$limit_end = $limit_explode[1];
 }
 
-$sort = filter_input(INPUT_GET,'sort',FILTER_SANITIZE_STRING);
+$sort = filter_input(INPUT_GET, 'sort', FILTER_SANITIZE_STRING);
 
 $absolute_difference = abs($limit_start - $limit_end);
 $limit_next = $limit_end + $absolute_difference;
@@ -36,7 +36,7 @@ 
 {
 	$category = "aircraft";
 } else {
-	$category = filter_input(INPUT_GET,'category',FILTER_SANITIZE_STRING);
+	$category = filter_input(INPUT_GET, 'category', FILTER_SANITIZE_STRING);
 }
 
 $page_url = $globalURL.'/newest/'.$category;

country-sub-menu.php

Spacing +6 added lines, -6 removed lines

@@ -3,9 +3,9 @@ 
 </span>
 <div class="sub-menu sub-menu-container">
 	<ul class="nav nav-pills">
-		<li><a href="<?php print $globalURL; ?>/country/<?php print $country; ?>" <?php if (strtolower($current_page) == "country-detailed"){ print 'class="active"'; } ?>><?php echo _("Detailed"); ?></a></li>
+		<li><a href="<?php print $globalURL; ?>/country/<?php print $country; ?>" <?php if (strtolower($current_page) == "country-detailed") { print 'class="active"'; } ?>><?php echo _("Detailed"); ?></a></li>
 		<li class="dropdown">
-		    <a class="dropdown-toggle <?php if(strtolower($current_page) == "country-statistics-aircraft" || strtolower($current_page) == "country-statistics-registration" || strtolower($current_page) == "country-statistics-manufacturer"){ print 'active'; } ?>" data-toggle="dropdown" href="#">
+		    <a class="dropdown-toggle <?php if (strtolower($current_page) == "country-statistics-aircraft" || strtolower($current_page) == "country-statistics-registration" || strtolower($current_page) == "country-statistics-manufacturer") { print 'active'; } ?>" data-toggle="dropdown" href="#">
 		      <?php echo _("Aircraft"); ?> <span class="caret"></span>
 		    </a>
 		    <ul class="dropdown-menu" role="menu">
@@ -15,7 +15,7 @@ 
 		    </ul>
 		</li>
 		<li class="dropdown">
-		    <a class="dropdown-toggle <?php if(strtolower($current_page) == "country-statistics-airline" || strtolower($current_page) == "country-statistics-airline-country"){ print 'active'; } ?>" data-toggle="dropdown" href="#">
+		    <a class="dropdown-toggle <?php if (strtolower($current_page) == "country-statistics-airline" || strtolower($current_page) == "country-statistics-airline-country") { print 'active'; } ?>" data-toggle="dropdown" href="#">
 		      <?php echo _("Airline"); ?> <span class="caret"></span>
 		    </a>
 		    <ul class="dropdown-menu" role="menu">
@@ -24,7 +24,7 @@ 
 		    </ul>
 		</li>
 		<li class="dropdown">
-		    <a class="dropdown-toggle <?php if(strtolower($current_page) == "country-statistics-departure-airport" || strtolower($current_page) == "country-statistics-departure-airport-country" || strtolower($current_page) == "country-statistics-arrival-airport" || strtolower($current_page) == "country-statistics-arrival-airport-country"){ print 'active'; } ?>" data-toggle="dropdown" href="#">
+		    <a class="dropdown-toggle <?php if (strtolower($current_page) == "country-statistics-departure-airport" || strtolower($current_page) == "country-statistics-departure-airport-country" || strtolower($current_page) == "country-statistics-arrival-airport" || strtolower($current_page) == "country-statistics-arrival-airport-country") { print 'active'; } ?>" data-toggle="dropdown" href="#">
 		      <?php echo _("Airport"); ?> <span class="caret"></span>
 		    </a>
 		    <ul class="dropdown-menu" role="menu">
@@ -34,7 +34,7 @@ 
 			  <li><a href="<?php print $globalURL; ?>/country/statistics/arrival-airport-country/<?php print $country; ?>"><?php echo _("Arrival Airport by Country"); ?></a></li>
 		    </ul>
 		</li>
-		<li><a href="<?php print $globalURL; ?>/country/statistics/route/<?php print $country; ?>" <?php if (strtolower($current_page) == "country-statistics-route"){ print 'class="active"'; } ?>><?php echo _("Route"); ?></a></li>
-		<li><a href="<?php print $globalURL; ?>/country/statistics/time/<?php print $country; ?>" <?php if (strtolower($current_page) == "country-statistics-time"){ print 'class="active"'; } ?>><?php echo _("Time"); ?></a></li>
+		<li><a href="<?php print $globalURL; ?>/country/statistics/route/<?php print $country; ?>" <?php if (strtolower($current_page) == "country-statistics-route") { print 'class="active"'; } ?>><?php echo _("Route"); ?></a></li>
+		<li><a href="<?php print $globalURL; ?>/country/statistics/time/<?php print $country; ?>" <?php if (strtolower($current_page) == "country-statistics-time") { print 'class="active"'; } ?>><?php echo _("Time"); ?></a></li>
 	</ul>
 </div>
\ No newline at end of file

registration-statistics-arrival-airport.php

Braces +3 added lines, -1 removed lines

@@ -8,7 +8,9 @@
 if ($registration != '') {
 	$spotter_array = $Spotter->getSpotterDataByRegistration($registration, "0,1", $sort);
 	$aircraft_array = $Spotter->getAircraftInfoByRegistration($registration);
-} else $spotter_array=array();
+} else {
+	$spotter_array=array();
+}
 
 if (!empty($spotter_array))
 {

Spacing +7 added lines, -7 removed lines

@@ -3,16 +3,16 @@ 
 require_once('require/class.Spotter.php');
 require_once('require/class.Language.php');
 $Spotter = new Spotter();
-$sort=filter_input(INPUT_GET,'sort',FILTER_SANITIZE_STRING);
-$registration = filter_input(INPUT_GET,'registration',FILTER_SANITIZE_STRING);
+$sort = filter_input(INPUT_GET, 'sort', FILTER_SANITIZE_STRING);
+$registration = filter_input(INPUT_GET, 'registration', FILTER_SANITIZE_STRING);
 if ($registration != '') {
 	$spotter_array = $Spotter->getSpotterDataByRegistration($registration, "0,1", $sort);
 	$aircraft_array = $Spotter->getAircraftInfoByRegistration($registration);
-} else $spotter_array=array();
+} else $spotter_array = array();
 
 if (!empty($spotter_array))
 {
-	$title = sprintf(_("Most Common Arrival Airports of aircraft with registration %s"),$registration);
+	$title = sprintf(_("Most Common Arrival Airports of aircraft with registration %s"), $registration);
 
 	require_once('header.php');
 		print '<div class="info column">';
@@ -25,7 +25,7 @@ 
 	include('registration-sub-menu.php');
 	print '<div class="column">';
 	print '<h2>'._("Most Common Arrival Airports").'</h2>';
-	print '<p>'.sprintf(_("The statistic below shows all arrival airports of flights with aircraft registration <strong>%s</strong>."),$registration).'</p>';
+	print '<p>'.sprintf(_("The statistic below shows all arrival airports of flights with aircraft registration <strong>%s</strong>."), $registration).'</p>';
 	$airport_airport_array = $Spotter->countAllArrivalAirportsByRegistration($registration);
 	print '<script type="text/javascript" src="'.$globalURL.'/js/d3.min.js"></script>';
 	print '<script type="text/javascript" src="'.$globalURL.'/js/topojson.v2.min.js"></script>';
@@ -34,7 +34,7 @@ 
 	print '<script>';
 	print 'var series = [';
 	$airport_data = '';
-	foreach($airport_airport_array as $airport_item)
+	foreach ($airport_airport_array as $airport_item)
 	{
 		$airport_data .= '[ "'.$airport_item['airport_arrival_icao_count'].'", "'.$airport_item['airport_arrival_name'].' ('.$airport_item['airport_arrival_icao'].')",'.$airport_item['airport_arrival_latitude'].','.$airport_item['airport_arrival_longitude'].'],';
 	}
@@ -88,7 +88,7 @@ 
 	print '</thead>';
 	print '<tbody>';
 	$i = 1;
-	foreach($airport_airport_array as $airport_item)
+	foreach ($airport_airport_array as $airport_item)
 	{
 		print '<tr>';
 		print '<td><strong>'.$i.'</strong></td>';

registration-statistics-arrival-airport-country.php

Braces +3 added lines, -1 removed lines

@@ -8,7 +8,9 @@
 if ($registration != '') {
 	$spotter_array = $Spotter->getSpotterDataByRegistration($registration, "0,1", $sort);
 	$aircraft_array = $Spotter->getAircraftInfoByRegistration($registration);
-} else $spotter_array=array();
+} else {
+	$spotter_array=array();
+}
 
 if (!empty($spotter_array))
 {

Spacing +7 added lines, -7 removed lines

@@ -3,16 +3,16 @@ 
 require_once('require/class.Spotter.php');
 require_once('require/class.Language.php');
 $Spotter = new Spotter();
-$sort = filter_input(INPUT_GET,'sort',FILTER_SANITIZE_STRING);
-$registration = filter_input(INPUT_GET,'registration',FILTER_SANITIZE_STRING);
+$sort = filter_input(INPUT_GET, 'sort', FILTER_SANITIZE_STRING);
+$registration = filter_input(INPUT_GET, 'registration', FILTER_SANITIZE_STRING);
 if ($registration != '') {
 	$spotter_array = $Spotter->getSpotterDataByRegistration($registration, "0,1", $sort);
 	$aircraft_array = $Spotter->getAircraftInfoByRegistration($registration);
-} else $spotter_array=array();
+} else $spotter_array = array();
 
 if (!empty($spotter_array))
 {
-	$title = sprintf(_("Most Common Arrival Airports by Country of aircraft with registration %s"),$registration);
+	$title = sprintf(_("Most Common Arrival Airports by Country of aircraft with registration %s"), $registration);
 	require_once('header.php');
 	print '<div class="info column">';
 	print '<h1>'.$registration.' - '.$aircraft_array[0]['aircraft_name'].' ('.$aircraft_array[0]['aircraft_icao'].')</h1>';
@@ -24,7 +24,7 @@ 
 	include('registration-sub-menu.php');
 	print '<div class="column">';
 	print '<h2>'._("Most Common Arrival Airports by Country").'</h2>';
-	print '<p>'.sprintf(_("The statistic below shows all arrival airports by Country of origin of flights with aircraft registration <strong>%s</strong>."),$registration).'</p>';
+	print '<p>'.sprintf(_("The statistic below shows all arrival airports by Country of origin of flights with aircraft registration <strong>%s</strong>."), $registration).'</p>';
 	$airport_country_array = $Spotter->countAllArrivalAirportCountriesByRegistration($registration);
 	print '<script type="text/javascript" src="'.$globalURL.'/js/d3.min.js"></script>';
 	print '<script type="text/javascript" src="'.$globalURL.'/js/topojson.v2.min.js"></script>';
@@ -32,7 +32,7 @@ 
 	print '<div id="chartCountry" class="chart" width="100%"></div><script>';
 	print 'var series = [';
 	$country_data = '';
-	foreach($airport_country_array as $airport_item)
+	foreach ($airport_country_array as $airport_item)
 	{
 		$country_data .= '[ "'.$airport_item['arrival_airport_country_iso3'].'",'.$airport_item['airport_arrival_country_count'].'],';
 	}
@@ -79,7 +79,7 @@ 
 		print '</thead>';
 		print '<tbody>';
 		$i = 1;
-		foreach($airport_country_array as $airport_item)
+		foreach ($airport_country_array as $airport_item)
 		{
 			print '<tr>';
 			print '<td><strong>'.$i.'</strong></td>';

polar-geojson.php

Braces +3 added lines, -1 removed lines

@@ -35,7 +35,9 @@
 		$first = '';
 		foreach($data as $value => $key) {
 			$final_coord = $Common->getCoordfromDistanceBearing($initial_latitude,$initial_longitude,$value*22.5,$key);
-			if ($first == '') $first = '['.round($final_coord['longitude'],5).','.round($final_coord['latitude'],5).']';
+			if ($first == '') {
+				$first = '['.round($final_coord['longitude'],5).','.round($final_coord['latitude'],5).']';
+			}
 			$output .= '['.$final_coord['longitude'].','.$final_coord['latitude'].'],';
 		}
 		$output .= $first;

Spacing +6 added lines, -6 removed lines

@@ -17,10 +17,10 @@ 
 header('Content-Type: text/javascript');
 
 
-$polar = $Stats->getStatsSource('polar',date('Y'),date('m'),date('d'));
+$polar = $Stats->getStatsSource('polar', date('Y'), date('m'), date('d'));
 $output = '{"type": "FeatureCollection","features": [';
 if (!empty($polar)) {
-	foreach($polar as $eachpolar) {
+	foreach ($polar as $eachpolar) {
 		$data = json_decode($eachpolar['source_data']);
 		$name = $eachpolar['source_name'];
 		$coord = $Location->getLocationInfobySourceName($name);
@@ -33,15 +33,15 @@ 
 			$initial_longitude = $globalCenterLongitude;
 		}
 		$first = '';
-		foreach($data as $value => $key) {
-			$final_coord = $Common->getCoordfromDistanceBearing($initial_latitude,$initial_longitude,$value*22.5,$key);
-			if ($first == '') $first = '['.round($final_coord['longitude'],5).','.round($final_coord['latitude'],5).']';
+		foreach ($data as $value => $key) {
+			$final_coord = $Common->getCoordfromDistanceBearing($initial_latitude, $initial_longitude, $value*22.5, $key);
+			if ($first == '') $first = '['.round($final_coord['longitude'], 5).','.round($final_coord['latitude'], 5).']';
 			$output .= '['.$final_coord['longitude'].','.$final_coord['latitude'].'],';
 		}
 		$output .= $first;
 		$output .= ']]}},';
 	}
-	$output  = substr($output, 0, -1);
+	$output = substr($output, 0, -1);
 }
 $output .= ']}';
 print $output;

require/libs/Predict/Predict.php

Doc Comments +3 added lines, -17 removed lines

@@ -146,20 +146,6 @@ 
 
     /** Predict first pass after a certain time.
      *
-     *  @param Predict_Sat $sat   The satellite data.
-     *  @param Predict_QTH $qth   The observer's location data.
-     *  @param float       $start Starting time.
-     *  @param int         $maxdt The maximum number of days to look ahead (0 for no limit).
-     *
-     *  @return Predict_Pass or NULL if there was an error.
-     *
-     * This function will find the first upcoming pass with AOS no earlier than
-     * t = start and no later than t = (start+maxdt).
-     *
-     *  note For no time limit use maxdt = 0.0
-     *
-     *  note the data in sat will be corrupt (future) and must be refreshed
-     *       by the caller, if the caller will need it later on
      */
     public function get_pass(Predict_Sat $sat_in, Predict_QTH $qth, $start, $maxdt)
     {
@@ -407,7 +393,7 @@ 
      *  @param Predict_QTH $qth   The observer's location (QTH) data.
      *  @param float       $start The julian date where calculation should start.
      *  @param int         $maxdt The upper time limit in days (0.0 = no limit)
-     *  @return The julain date of the next AOS or 0.0 if the satellite has no AOS.
+     *  @return double julain date of the next AOS or 0.0 if the satellite has no AOS.
      *
      * This function finds the time of AOS for the first coming pass taking place
      * no earlier that start.
@@ -563,7 +549,7 @@ 
      *  @param Predict_QTH $qth The QTH observer location data.
      *  @param float       $start The time where calculation should start. (Julian Date)
      *  @param int         $maxdt The upper time limit in days (0.0 = no limit)
-     *  @return The time (julian date) of the next LOS or 0.0 if the satellite has no LOS.
+     *  @return double time (julian date) of the next LOS or 0.0 if the satellite has no LOS.
      *
      * This function finds the time of LOS for the first coming pass taking place
      * no earlier that start.
@@ -648,7 +634,7 @@ 
      *  @param Predict_Sat $sat   The satellite to find AOS for.
      *  @param Predict_QTH $qth   The ground station.
      *  @param float       $start Start time, prefereably now.
-     *  @return The time of the previous AOS or 0.0 if the satellite has no AOS.
+     *  @return double time of the previous AOS or 0.0 if the satellite has no AOS.
      *
      * This function can be used to find the AOS time in the past of the
      * current pass.

Indentation +791 added lines, -791 removed lines

@@ -52,824 +52,824 @@
  */
 class Predict
 {
-    const de2ra    =  1.74532925E-2;   /* Degrees to Radians */
-    const pi       =  3.1415926535898; /* Pi */
-    const pio2     =  1.5707963267949; /* Pi/2 */
-    const x3pio2   =  4.71238898;      /* 3*Pi/2 */
-    const twopi    =  6.2831853071796; /* 2*Pi  */
-    const e6a      =  1.0E-6;
-    const tothrd   =  6.6666667E-1;    /* 2/3 */
-    const xj2      =  1.0826158E-3;    /* J2 Harmonic */
-    const xj3      = -2.53881E-6;      /* J3 Harmonic */
-    const xj4      = -1.65597E-6;      /* J4 Harmonic */
-    const xke      =  7.43669161E-2;
-    const xkmper   =  6.378135E3;      /* Earth radius km */
-    const xmnpda   =  1.44E3;          /* Minutes per day */
-    const km2mi    =  0.621371;        /* Kilometers per Mile */
-    const ae       =  1.0;
-    const ck2      =  5.413079E-4;
-    const ck4      =  6.209887E-7;
-    const __f      =  3.352779E-3;
-    const ge       =  3.986008E5;
-    const __s__    =  1.012229;
-    const qoms2t   =  1.880279E-09;
-    const secday   =  8.6400E4;        /* Seconds per day */
-    const omega_E  =  1.0027379;
-    const omega_ER =  6.3003879;
-    const zns      =  1.19459E-5;
-    const c1ss     =  2.9864797E-6;
-    const zes      =  1.675E-2;
-    const znl      =  1.5835218E-4;
-    const c1l      =  4.7968065E-7;
-    const zel      =  5.490E-2;
-    const zcosis   =  9.1744867E-1;
-    const zsinis   =  3.9785416E-1;
-    const zsings   = -9.8088458E-1;
-    const zcosgs   =  1.945905E-1;
-    const zcoshs   =  1;
-    const zsinhs   =  0;
-    const q22      =  1.7891679E-6;
-    const q31      =  2.1460748E-6;
-    const q33      =  2.2123015E-7;
-    const g22      =  5.7686396;
-    const g32      =  9.5240898E-1;
-    const g44      =  1.8014998;
-    const g52      =  1.0508330;
-    const g54      =  4.4108898;
-    const root22   =  1.7891679E-6;
-    const root32   =  3.7393792E-7;
-    const root44   =  7.3636953E-9;
-    const root52   =  1.1428639E-7;
-    const root54   =  2.1765803E-9;
-    const thdt     =  4.3752691E-3;
-    const rho      =  1.5696615E-1;
-    const mfactor  =  7.292115E-5;
-    const __sr__   =  6.96000E5;      /*Solar radius - kilometers (IAU 76)*/
-    const AU       =  1.49597870E8;   /*Astronomical unit - kilometers (IAU 76)*/
-
-    /* visibility constants */
-    const SAT_VIS_NONE     = 0;
-    const SAT_VIS_VISIBLE  = 1;
-    const SAT_VIS_DAYLIGHT = 2;
-    const SAT_VIS_ECLIPSED = 3;
-
-    /* preferences */
-    public $minEle     = 10; // Minimum elevation
-    public $timeRes    = 10; // Pass details: time resolution
-    public $numEntries = 20; // Pass details: number of entries
-    public $threshold  = -6; // Twilight threshold
-
-    /**
-     *  Predict the next pass.
-     *
-     * This function simply wraps the get_pass function using the current time
-     * as parameter.
-     *
-     * Note: the data in sat will be corrupt (future) and must be refreshed
-     *       by the caller, if the caller will need it later on (eg. if the caller
-     *       is GtkSatList).
-     *
-     * @param Predict_Sat $sat   The satellite data.
-     * @param Predict_QTH $qth   The observer data.
-     * @param int         $maxdt The maximum number of days to look ahead.
-     *
-     * @return Predict_Pass Pointer instance or NULL if no pass can be
-     *         found.
-     */
-    public function get_next_pass(Predict_Sat $sat, Predict_QTH $qth, $maxdt)
-    {
-        /* get the current time and call the get_pass function */
-        $now = Predict_Time::get_current_daynum();
-
-        return $this->get_pass($sat, $qth, $now, $maxdt);
-    }
-
-    /** Predict first pass after a certain time.
-     *
-     *  @param Predict_Sat $sat   The satellite data.
-     *  @param Predict_QTH $qth   The observer's location data.
-     *  @param float       $start Starting time.
-     *  @param int         $maxdt The maximum number of days to look ahead (0 for no limit).
-     *
-     *  @return Predict_Pass or NULL if there was an error.
-     *
-     * This function will find the first upcoming pass with AOS no earlier than
-     * t = start and no later than t = (start+maxdt).
-     *
-     *  note For no time limit use maxdt = 0.0
-     *
-     *  note the data in sat will be corrupt (future) and must be refreshed
-     *       by the caller, if the caller will need it later on
-     */
-    public function get_pass(Predict_Sat $sat_in, Predict_QTH $qth, $start, $maxdt)
-    {
-        $aos = 0.0;    /* time of AOS */
-        $tca = 0.0;    /* time of TCA */
-        $los = 0.0;    /* time of LOS */
-        $dt = 0.0;     /* time diff */
-        $step = 0.0;   /* time step */
-        $t0 = $start;
-        $tres = 0.0;   /* required time resolution */
-        $max_el = 0.0; /* maximum elevation */
-        $pass = null;
-        $detail = null;
-        $done = false;
-        $iter = 0;      /* number of iterations */
-        /* FIXME: watchdog */
-
-        /*copy sat_in to a working structure*/
-        $sat         = clone $sat_in;
-        $sat_working = clone $sat_in;
-
-        /* get time resolution; sat-cfg stores it in seconds */
-        $tres = $this->timeRes / 86400.0;
-
-        /* loop until we find a pass with elevation > SAT_CFG_INT_PRED_MIN_EL
+	const de2ra    =  1.74532925E-2;   /* Degrees to Radians */
+	const pi       =  3.1415926535898; /* Pi */
+	const pio2     =  1.5707963267949; /* Pi/2 */
+	const x3pio2   =  4.71238898;      /* 3*Pi/2 */
+	const twopi    =  6.2831853071796; /* 2*Pi  */
+	const e6a      =  1.0E-6;
+	const tothrd   =  6.6666667E-1;    /* 2/3 */
+	const xj2      =  1.0826158E-3;    /* J2 Harmonic */
+	const xj3      = -2.53881E-6;      /* J3 Harmonic */
+	const xj4      = -1.65597E-6;      /* J4 Harmonic */
+	const xke      =  7.43669161E-2;
+	const xkmper   =  6.378135E3;      /* Earth radius km */
+	const xmnpda   =  1.44E3;          /* Minutes per day */
+	const km2mi    =  0.621371;        /* Kilometers per Mile */
+	const ae       =  1.0;
+	const ck2      =  5.413079E-4;
+	const ck4      =  6.209887E-7;
+	const __f      =  3.352779E-3;
+	const ge       =  3.986008E5;
+	const __s__    =  1.012229;
+	const qoms2t   =  1.880279E-09;
+	const secday   =  8.6400E4;        /* Seconds per day */
+	const omega_E  =  1.0027379;
+	const omega_ER =  6.3003879;
+	const zns      =  1.19459E-5;
+	const c1ss     =  2.9864797E-6;
+	const zes      =  1.675E-2;
+	const znl      =  1.5835218E-4;
+	const c1l      =  4.7968065E-7;
+	const zel      =  5.490E-2;
+	const zcosis   =  9.1744867E-1;
+	const zsinis   =  3.9785416E-1;
+	const zsings   = -9.8088458E-1;
+	const zcosgs   =  1.945905E-1;
+	const zcoshs   =  1;
+	const zsinhs   =  0;
+	const q22      =  1.7891679E-6;
+	const q31      =  2.1460748E-6;
+	const q33      =  2.2123015E-7;
+	const g22      =  5.7686396;
+	const g32      =  9.5240898E-1;
+	const g44      =  1.8014998;
+	const g52      =  1.0508330;
+	const g54      =  4.4108898;
+	const root22   =  1.7891679E-6;
+	const root32   =  3.7393792E-7;
+	const root44   =  7.3636953E-9;
+	const root52   =  1.1428639E-7;
+	const root54   =  2.1765803E-9;
+	const thdt     =  4.3752691E-3;
+	const rho      =  1.5696615E-1;
+	const mfactor  =  7.292115E-5;
+	const __sr__   =  6.96000E5;      /*Solar radius - kilometers (IAU 76)*/
+	const AU       =  1.49597870E8;   /*Astronomical unit - kilometers (IAU 76)*/
+
+	/* visibility constants */
+	const SAT_VIS_NONE     = 0;
+	const SAT_VIS_VISIBLE  = 1;
+	const SAT_VIS_DAYLIGHT = 2;
+	const SAT_VIS_ECLIPSED = 3;
+
+	/* preferences */
+	public $minEle     = 10; // Minimum elevation
+	public $timeRes    = 10; // Pass details: time resolution
+	public $numEntries = 20; // Pass details: number of entries
+	public $threshold  = -6; // Twilight threshold
+
+	/**
+	 *  Predict the next pass.
+	 *
+	 * This function simply wraps the get_pass function using the current time
+	 * as parameter.
+	 *
+	 * Note: the data in sat will be corrupt (future) and must be refreshed
+	 *       by the caller, if the caller will need it later on (eg. if the caller
+	 *       is GtkSatList).
+	 *
+	 * @param Predict_Sat $sat   The satellite data.
+	 * @param Predict_QTH $qth   The observer data.
+	 * @param int         $maxdt The maximum number of days to look ahead.
+	 *
+	 * @return Predict_Pass Pointer instance or NULL if no pass can be
+	 *         found.
+	 */
+	public function get_next_pass(Predict_Sat $sat, Predict_QTH $qth, $maxdt)
+	{
+		/* get the current time and call the get_pass function */
+		$now = Predict_Time::get_current_daynum();
+
+		return $this->get_pass($sat, $qth, $now, $maxdt);
+	}
+
+	/** Predict first pass after a certain time.
+	 *
+	 *  @param Predict_Sat $sat   The satellite data.
+	 *  @param Predict_QTH $qth   The observer's location data.
+	 *  @param float       $start Starting time.
+	 *  @param int         $maxdt The maximum number of days to look ahead (0 for no limit).
+	 *
+	 *  @return Predict_Pass or NULL if there was an error.
+	 *
+	 * This function will find the first upcoming pass with AOS no earlier than
+	 * t = start and no later than t = (start+maxdt).
+	 *
+	 *  note For no time limit use maxdt = 0.0
+	 *
+	 *  note the data in sat will be corrupt (future) and must be refreshed
+	 *       by the caller, if the caller will need it later on
+	 */
+	public function get_pass(Predict_Sat $sat_in, Predict_QTH $qth, $start, $maxdt)
+	{
+		$aos = 0.0;    /* time of AOS */
+		$tca = 0.0;    /* time of TCA */
+		$los = 0.0;    /* time of LOS */
+		$dt = 0.0;     /* time diff */
+		$step = 0.0;   /* time step */
+		$t0 = $start;
+		$tres = 0.0;   /* required time resolution */
+		$max_el = 0.0; /* maximum elevation */
+		$pass = null;
+		$detail = null;
+		$done = false;
+		$iter = 0;      /* number of iterations */
+		/* FIXME: watchdog */
+
+		/*copy sat_in to a working structure*/
+		$sat         = clone $sat_in;
+		$sat_working = clone $sat_in;
+
+		/* get time resolution; sat-cfg stores it in seconds */
+		$tres = $this->timeRes / 86400.0;
+
+		/* loop until we find a pass with elevation > SAT_CFG_INT_PRED_MIN_EL
             or we run out of time
             FIXME: we should have a safety break
         */
-        while (!$done) {
-            /* Find los of next pass or of current pass */
-            $los = $this->find_los($sat, $qth, $t0, $maxdt); // See if a pass is ongoing
-            $aos = $this->find_aos($sat, $qth, $t0, $maxdt);
-            /* sat_log_log(SAT_LOG_LEVEL_MSG, "%s:%s:%d: found aos %f and los %f for t0=%f", */
-            /*          __FILE__,  */
-            /*          __FUNCTION__, */
-            /*          __LINE__, */
-            /*          aos, */
-            /*          los,  */
-            /*          t0); */
-            if ($aos > $los) {
-                // los is from an currently happening pass, find previous aos
-                $aos = $this->find_prev_aos($sat, $qth, $t0);
-            }
-
-            /* aos = 0.0 means no aos */
-            if ($aos == 0.0) {
-                $done = true;
-            } else if (($maxdt > 0.0) && ($aos > ($start + $maxdt)) ) {
-                /* check whether we are within time limits;
+		while (!$done) {
+			/* Find los of next pass or of current pass */
+			$los = $this->find_los($sat, $qth, $t0, $maxdt); // See if a pass is ongoing
+			$aos = $this->find_aos($sat, $qth, $t0, $maxdt);
+			/* sat_log_log(SAT_LOG_LEVEL_MSG, "%s:%s:%d: found aos %f and los %f for t0=%f", */
+			/*          __FILE__,  */
+			/*          __FUNCTION__, */
+			/*          __LINE__, */
+			/*          aos, */
+			/*          los,  */
+			/*          t0); */
+			if ($aos > $los) {
+				// los is from an currently happening pass, find previous aos
+				$aos = $this->find_prev_aos($sat, $qth, $t0);
+			}
+
+			/* aos = 0.0 means no aos */
+			if ($aos == 0.0) {
+				$done = true;
+			} else if (($maxdt > 0.0) && ($aos > ($start + $maxdt)) ) {
+				/* check whether we are within time limits;
                     maxdt = 0 mean no time limit.
                 */
-                $done = true;
-            } else {
-                //los = find_los (sat, qth, aos + 0.001, maxdt); // +1.5 min later
-                $dt = $los - $aos;
+				$done = true;
+			} else {
+				//los = find_los (sat, qth, aos + 0.001, maxdt); // +1.5 min later
+				$dt = $los - $aos;
 
-                /* get time step, which will give us the max number of entries */
-                $step = $dt / $this->numEntries;
+				/* get time step, which will give us the max number of entries */
+				$step = $dt / $this->numEntries;
 
-                /* but if this is smaller than the required resolution
+				/* but if this is smaller than the required resolution
                     we go with the resolution
                 */
-                if ($step < $tres) {
-                    $step = $tres;
-                }
-
-                /* create a pass_t entry; FIXME: g_try_new in 2.8 */
-                $pass = new Predict_Pass();
-
-                $pass->aos      = $aos;
-                $pass->los      = $los;
-                $pass->max_el   = 0.0;
-                $pass->aos_az   = 0.0;
-                $pass->los_az   = 0.0;
-                $pass->maxel_az = 0.0;
-                $pass->vis      = '---';
-                $pass->satname  = $sat->nickname;
-                $pass->details  = array();
-
-                /* iterate over each time step */
-                for ($t = $pass->aos; $t <= $pass->los; $t += $step) {
-
-                    /* calculate satellite data */
-                    $this->predict_calc($sat, $qth, $t);
-
-                    /* in the first iter we want to store
+				if ($step < $tres) {
+					$step = $tres;
+				}
+
+				/* create a pass_t entry; FIXME: g_try_new in 2.8 */
+				$pass = new Predict_Pass();
+
+				$pass->aos      = $aos;
+				$pass->los      = $los;
+				$pass->max_el   = 0.0;
+				$pass->aos_az   = 0.0;
+				$pass->los_az   = 0.0;
+				$pass->maxel_az = 0.0;
+				$pass->vis      = '---';
+				$pass->satname  = $sat->nickname;
+				$pass->details  = array();
+
+				/* iterate over each time step */
+				for ($t = $pass->aos; $t <= $pass->los; $t += $step) {
+
+					/* calculate satellite data */
+					$this->predict_calc($sat, $qth, $t);
+
+					/* in the first iter we want to store
                         pass->aos_az
                     */
-                    if ($t == $pass->aos) {
-                        $pass->aos_az = $sat->az;
-                        $pass->orbit  = $sat->orbit;
-                    }
-
-                    /* append details to sat->details */
-                    $detail             = new Predict_PassDetail();
-                    $detail->time       = $t;
-                    $detail->pos->x     = $sat->pos->x;
-                    $detail->pos->y     = $sat->pos->y;
-                    $detail->pos->z     = $sat->pos->z;
-                    $detail->pos->w     = $sat->pos->w;
-                    $detail->vel->x     = $sat->vel->x;
-                    $detail->vel->y     = $sat->vel->y;
-                    $detail->vel->z     = $sat->vel->z;
-                    $detail->vel->w     = $sat->vel->w;
-                    $detail->velo       = $sat->velo;
-                    $detail->az         = $sat->az;
-                    $detail->el         = $sat->el;
-                    $detail->range      = $sat->range;
-                    $detail->range_rate = $sat->range_rate;
-                    $detail->lat        = $sat->ssplat;
-                    $detail->lon        = $sat->ssplon;
-                    $detail->alt        = $sat->alt;
-                    $detail->ma         = $sat->ma;
-                    $detail->phase      = $sat->phase;
-                    $detail->footprint  = $sat->footprint;
-                    $detail->orbit      = $sat->orbit;
-                    $detail->vis        = $this->get_sat_vis($sat, $qth, $t);
-
-                    /* also store visibility "bit" */
-                    switch ($detail->vis) {
-                        case self::SAT_VIS_VISIBLE:
-                            $pass->vis[0] = 'V';
-                            break;
-                        case self::SAT_VIS_DAYLIGHT:
-                            $pass->vis[1] = 'D';
-                            break;
-                        case self::SAT_VIS_ECLIPSED:
-                            $pass->vis[2] = 'E';
-                            break;
-                        default:
-                            break;
-                    }
-
-                    // Using an array, no need to prepend and reverse the list
-                    // as gpredict does
-                    $pass->details[] = $detail;
-
-                    // Look up apparent magnitude if this is a visible pass
-                    if ($detail->vis === self::SAT_VIS_VISIBLE) {
-                        $apmag = $sat->calculateApparentMagnitude($t, $qth);
-                        if ($pass->max_apparent_magnitude === null || $apmag < $pass->max_apparent_magnitude) {
-                            $pass->max_apparent_magnitude = $apmag;
-                        }
-                    }
-
-                    /* store elevation if greater than the
+					if ($t == $pass->aos) {
+						$pass->aos_az = $sat->az;
+						$pass->orbit  = $sat->orbit;
+					}
+
+					/* append details to sat->details */
+					$detail             = new Predict_PassDetail();
+					$detail->time       = $t;
+					$detail->pos->x     = $sat->pos->x;
+					$detail->pos->y     = $sat->pos->y;
+					$detail->pos->z     = $sat->pos->z;
+					$detail->pos->w     = $sat->pos->w;
+					$detail->vel->x     = $sat->vel->x;
+					$detail->vel->y     = $sat->vel->y;
+					$detail->vel->z     = $sat->vel->z;
+					$detail->vel->w     = $sat->vel->w;
+					$detail->velo       = $sat->velo;
+					$detail->az         = $sat->az;
+					$detail->el         = $sat->el;
+					$detail->range      = $sat->range;
+					$detail->range_rate = $sat->range_rate;
+					$detail->lat        = $sat->ssplat;
+					$detail->lon        = $sat->ssplon;
+					$detail->alt        = $sat->alt;
+					$detail->ma         = $sat->ma;
+					$detail->phase      = $sat->phase;
+					$detail->footprint  = $sat->footprint;
+					$detail->orbit      = $sat->orbit;
+					$detail->vis        = $this->get_sat_vis($sat, $qth, $t);
+
+					/* also store visibility "bit" */
+					switch ($detail->vis) {
+						case self::SAT_VIS_VISIBLE:
+							$pass->vis[0] = 'V';
+							break;
+						case self::SAT_VIS_DAYLIGHT:
+							$pass->vis[1] = 'D';
+							break;
+						case self::SAT_VIS_ECLIPSED:
+							$pass->vis[2] = 'E';
+							break;
+						default:
+							break;
+					}
+
+					// Using an array, no need to prepend and reverse the list
+					// as gpredict does
+					$pass->details[] = $detail;
+
+					// Look up apparent magnitude if this is a visible pass
+					if ($detail->vis === self::SAT_VIS_VISIBLE) {
+						$apmag = $sat->calculateApparentMagnitude($t, $qth);
+						if ($pass->max_apparent_magnitude === null || $apmag < $pass->max_apparent_magnitude) {
+							$pass->max_apparent_magnitude = $apmag;
+						}
+					}
+
+					/* store elevation if greater than the
                         previously stored one
                     */
-                    if ($sat->el > $max_el) {
-                        $max_el         = $sat->el;
-                        $tca            = $t;
-                        $pass->maxel_az = $sat->az;
-                    }
-
-                    /*     g_print ("TIME: %f\tAZ: %f\tEL: %f (MAX: %f)\n", */
-                    /*           t, sat->az, sat->el, max_el); */
-                }
-
-                /* calculate satellite data */
-                $this->predict_calc($sat, $qth, $pass->los);
-                /* store los_az, max_el and tca */
-                $pass->los_az = $sat->az;
-                $pass->max_el = $max_el;
-                $pass->tca    = $tca;
-
-                /* check whether this pass is good */
-                if ($max_el >= $this->minEle) {
-                    $done = true;
-                } else {
-                    $done = false;
-                    $t0 = $los + 0.014; // +20 min
-                    $pass = null;
-                }
-
-                $iter++;
-            }
-        }
-
-        return $pass;
-    }
-
-    /**
-     * Calculate satellite visibility.
-     *
-     * @param Predict_Sat $sat     The satellite structure.
-     * @param Predict_QTH $qth     The QTH
-     * @param float       $jul_utc The time at which the visibility should be calculated.
-     *
-     * @return int The visiblity constant, 0, 1, 2, or 3 (see above)
-     */
-    public function get_sat_vis(Predict_Sat $sat, Predict_QTH $qth, $jul_utc)
-    {
-        /* gboolean sat_sun_status;
+					if ($sat->el > $max_el) {
+						$max_el         = $sat->el;
+						$tca            = $t;
+						$pass->maxel_az = $sat->az;
+					}
+
+					/*     g_print ("TIME: %f\tAZ: %f\tEL: %f (MAX: %f)\n", */
+					/*           t, sat->az, sat->el, max_el); */
+				}
+
+				/* calculate satellite data */
+				$this->predict_calc($sat, $qth, $pass->los);
+				/* store los_az, max_el and tca */
+				$pass->los_az = $sat->az;
+				$pass->max_el = $max_el;
+				$pass->tca    = $tca;
+
+				/* check whether this pass is good */
+				if ($max_el >= $this->minEle) {
+					$done = true;
+				} else {
+					$done = false;
+					$t0 = $los + 0.014; // +20 min
+					$pass = null;
+				}
+
+				$iter++;
+			}
+		}
+
+		return $pass;
+	}
+
+	/**
+	 * Calculate satellite visibility.
+	 *
+	 * @param Predict_Sat $sat     The satellite structure.
+	 * @param Predict_QTH $qth     The QTH
+	 * @param float       $jul_utc The time at which the visibility should be calculated.
+	 *
+	 * @return int The visiblity constant, 0, 1, 2, or 3 (see above)
+	 */
+	public function get_sat_vis(Predict_Sat $sat, Predict_QTH $qth, $jul_utc)
+	{
+		/* gboolean sat_sun_status;
         gdouble  sun_el;
         gdouble  threshold;
         gdouble  eclipse_depth;
         sat_vis_t vis = SAT_VIS_NONE; */
 
-        $eclipse_depth  = 0.0;
-        $zero_vector    = new Predict_Vector();
-        $obs_geodetic   = new Predict_Geodetic();
-
-        /* Solar ECI position vector  */
-        $solar_vector = new Predict_Vector();
-
-        /* Solar observed az and el vector  */
-        $solar_set = new Predict_ObsSet();
-
-        /* FIXME: could be passed as parameter */
-        $obs_geodetic->lon   = $qth->lon * self::de2ra;
-        $obs_geodetic->lat   = $qth->lat * self::de2ra;
-        $obs_geodetic->alt   = $qth->alt / 1000.0;
-        $obs_geodetic->theta = 0;
-
-        Predict_Solar::Calculate_Solar_Position($jul_utc, $solar_vector);
-        Predict_SGPObs::Calculate_Obs($jul_utc, $solar_vector, $zero_vector, $obs_geodetic, $solar_set);
-
-        if (Predict_Solar::Sat_Eclipsed($sat->pos, $solar_vector, $eclipse_depth)) {
-            /* satellite is eclipsed */
-            $sat_sun_status = false;
-        } else {
-            /* satellite in sunlight => may be visible */
-            $sat_sun_status = true;
-        }
-
-        if ($sat_sun_status) {
-            $sun_el = Predict_Math::Degrees($solar_set->el);
-
-            if ($sun_el <= $this->threshold && $sat->el >= 0.0) {
-                $vis = self::SAT_VIS_VISIBLE;
-            } else {
-                $vis = self::SAT_VIS_DAYLIGHT;
-            }
-        } else {
-            $vis = self::SAT_VIS_ECLIPSED;
-        }
-
-        return $vis;
-    }
-
-    /** Find the AOS time of the next pass.
-     *  @author Alexandru Csete, OZ9AEC
-     *  @author John A. Magliacane, KD2BD
-     *  @param Predict_Sat $sat   The satellite data.
-     *  @param Predict_QTH $qth   The observer's location (QTH) data.
-     *  @param float       $start The julian date where calculation should start.
-     *  @param int         $maxdt The upper time limit in days (0.0 = no limit)
-     *  @return The julain date of the next AOS or 0.0 if the satellite has no AOS.
-     *
-     * This function finds the time of AOS for the first coming pass taking place
-     * no earlier that start.
-     * If the satellite is currently within range, the function first calls
-     * find_los to get the next LOS time. Then the calculations are done using
-     * the new start time.
-     *
-     */
-    public function find_aos(Predict_Sat $sat, Predict_QTH $qth, $start, $maxdt)
-    {
-        $t = $start;
-        $aostime = 0.0;
-
-
-        /* make sure current sat values are
+		$eclipse_depth  = 0.0;
+		$zero_vector    = new Predict_Vector();
+		$obs_geodetic   = new Predict_Geodetic();
+
+		/* Solar ECI position vector  */
+		$solar_vector = new Predict_Vector();
+
+		/* Solar observed az and el vector  */
+		$solar_set = new Predict_ObsSet();
+
+		/* FIXME: could be passed as parameter */
+		$obs_geodetic->lon   = $qth->lon * self::de2ra;
+		$obs_geodetic->lat   = $qth->lat * self::de2ra;
+		$obs_geodetic->alt   = $qth->alt / 1000.0;
+		$obs_geodetic->theta = 0;
+
+		Predict_Solar::Calculate_Solar_Position($jul_utc, $solar_vector);
+		Predict_SGPObs::Calculate_Obs($jul_utc, $solar_vector, $zero_vector, $obs_geodetic, $solar_set);
+
+		if (Predict_Solar::Sat_Eclipsed($sat->pos, $solar_vector, $eclipse_depth)) {
+			/* satellite is eclipsed */
+			$sat_sun_status = false;
+		} else {
+			/* satellite in sunlight => may be visible */
+			$sat_sun_status = true;
+		}
+
+		if ($sat_sun_status) {
+			$sun_el = Predict_Math::Degrees($solar_set->el);
+
+			if ($sun_el <= $this->threshold && $sat->el >= 0.0) {
+				$vis = self::SAT_VIS_VISIBLE;
+			} else {
+				$vis = self::SAT_VIS_DAYLIGHT;
+			}
+		} else {
+			$vis = self::SAT_VIS_ECLIPSED;
+		}
+
+		return $vis;
+	}
+
+	/** Find the AOS time of the next pass.
+	 *  @author Alexandru Csete, OZ9AEC
+	 *  @author John A. Magliacane, KD2BD
+	 *  @param Predict_Sat $sat   The satellite data.
+	 *  @param Predict_QTH $qth   The observer's location (QTH) data.
+	 *  @param float       $start The julian date where calculation should start.
+	 *  @param int         $maxdt The upper time limit in days (0.0 = no limit)
+	 *  @return The julain date of the next AOS or 0.0 if the satellite has no AOS.
+	 *
+	 * This function finds the time of AOS for the first coming pass taking place
+	 * no earlier that start.
+	 * If the satellite is currently within range, the function first calls
+	 * find_los to get the next LOS time. Then the calculations are done using
+	 * the new start time.
+	 *
+	 */
+	public function find_aos(Predict_Sat $sat, Predict_QTH $qth, $start, $maxdt)
+	{
+		$t = $start;
+		$aostime = 0.0;
+
+
+		/* make sure current sat values are
             in sync with the time
         */
-        $this->predict_calc($sat, $qth, $start);
-
-        /* check whether satellite has aos */
-        if (($sat->otype == Predict_SGPSDP::ORBIT_TYPE_GEO) ||
-            ($sat->otype == Predict_SGPSDP::ORBIT_TYPE_DECAYED) ||
-            !$this->has_aos($sat, $qth)) {
-
-            return 0.0;
-        }
-
-        if ($sat->el > 0.0) {
-            $t = $this->find_los($sat, $qth, $start, $maxdt) + 0.014; // +20 min
-        }
-
-        /* invalid time (potentially returned by find_los) */
-        if ($t < 0.1) {
-            return 0.0;
-        }
-
-        /* update satellite data */
-        $this->predict_calc($sat, $qth, $t);
-
-        /* use upper time limit */
-        if ($maxdt > 0.0) {
-
-            /* coarse time steps */
-            while (($sat->el < -1.0) && ($t <= ($start + $maxdt))) {
-                $t -= 0.00035 * ($sat->el * (($sat->alt / 8400.0) + 0.46) - 2.0);
-                $this->predict_calc($sat, $qth, $t);
-            }
-
-            /* fine steps */
-            while (($aostime == 0.0) && ($t <= ($start + $maxdt))) {
-
-                if (abs($sat->el) < 0.005) {
-                    $aostime = $t;
-                } else {
-                    $t -= $sat->el * sqrt($sat->alt) / 530000.0;
-                    $this->predict_calc($sat, $qth, $t);
-                }
-            }
-        } else {
-            /* don't use upper time limit */
-
-            /* coarse time steps */
-            while ($sat->el < -1.0) {
-
-                $t -= 0.00035 * ($sat->el * (($sat->alt / 8400.0) + 0.46) - 2.0);
-                $this->predict_calc($sat, $qth, $t);
-            }
-
-            /* fine steps */
-            while ($aostime == 0.0) {
-
-                if (abs($sat->el) < 0.005) {
-                    $aostime = $t;
-                } else {
-                    $t -= $sat->el * sqrt($sat->alt) / 530000.0;
-                    $this->predict_calc($sat, $qth, $t);
-                }
-
-            }
-        }
-
-        return $aostime;
-    }
-
-    /** SGP4SDP4 driver for doing AOS/LOS calculations.
-     *  @param Predict_Sat $sat The satellite data.
-     *  @param Predict_QTH $qth The QTH observer location data.
-     *  @param float       $t   The time for calculation (Julian Date)
-     *
-     */
-    public function predict_calc(Predict_Sat $sat, Predict_QTH $qth, $t)
-    {
-        $obs_set      = new Predict_ObsSet();
-        $sat_geodetic = new Predict_Geodetic();
-        $obs_geodetic = new Predict_Geodetic();
-
-        $obs_geodetic->lon   = $qth->lon * self::de2ra;
-        $obs_geodetic->lat   = $qth->lat * self::de2ra;
-        $obs_geodetic->alt   = $qth->alt / 1000.0;
-        $obs_geodetic->theta = 0;
-
-        $sat->jul_utc = $t;
-        $sat->tsince = ($sat->jul_utc - $sat->jul_epoch) * self::xmnpda;
-
-        /* call the norad routines according to the deep-space flag */
-        $sgpsdp = Predict_SGPSDP::getInstance($sat);
-        if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
-            $sgpsdp->SDP4($sat, $sat->tsince);
-        } else {
-            $sgpsdp->SGP4($sat, $sat->tsince);
-        }
-
-        Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
-
-        /* get the velocity of the satellite */
-        $sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
-        $sat->velo = $sat->vel->w;
-        Predict_SGPObs::Calculate_Obs($sat->jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
-        Predict_SGPObs::Calculate_LatLonAlt($sat->jul_utc, $sat->pos, $sat_geodetic);
-
-        while ($sat_geodetic->lon < -self::pi) {
-            $sat_geodetic->lon += self::twopi;
-        }
-
-        while ($sat_geodetic->lon > (self::pi)) {
-            $sat_geodetic->lon -= self::twopi;
-        }
-
-        $sat->az = Predict_Math::Degrees($obs_set->az);
-        $sat->el = Predict_Math::Degrees($obs_set->el);
-        $sat->range = $obs_set->range;
-        $sat->range_rate = $obs_set->range_rate;
-        $sat->ssplat = Predict_Math::Degrees($sat_geodetic->lat);
-        $sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
-        $sat->alt = $sat_geodetic->alt;
-        $sat->ma = Predict_Math::Degrees($sat->phase);
-        $sat->ma *= 256.0 / 360.0;
-        $sat->phase = Predict_Math::Degrees($sat->phase);
-
-        /* same formulas, but the one from predict is nicer */
-        //sat->footprint = 2.0 * xkmper * acos (xkmper/sat->pos.w);
-        $sat->footprint = 12756.33 * acos(self::xkmper / (self::xkmper + $sat->alt));
-        $age = $sat->jul_utc - $sat->jul_epoch;
-        $sat->orbit = floor(($sat->tle->xno * self::xmnpda / self::twopi +
-                        $age * $sat->tle->bstar * self::ae) * $age +
-                        $sat->tle->xmo / self::twopi) + $sat->tle->revnum - 1;
-    }
-
-    /** Find the LOS time of the next pass.
-     *  @author Alexandru Csete, OZ9AEC
-     *  @author John A. Magliacane, KD2BD
-     *  @param Predict_Sat $sat The satellite data.
-     *  @param Predict_QTH $qth The QTH observer location data.
-     *  @param float       $start The time where calculation should start. (Julian Date)
-     *  @param int         $maxdt The upper time limit in days (0.0 = no limit)
-     *  @return The time (julian date) of the next LOS or 0.0 if the satellite has no LOS.
-     *
-     * This function finds the time of LOS for the first coming pass taking place
-     * no earlier that start.
-     * If the satellite is currently out of range, the function first calls
-     * find_aos to get the next AOS time. Then the calculations are done using
-     * the new start time.
-     * The function has a built-in watchdog to ensure that we don't end up in
-     * lengthy loops.
-     *
-     */
-    public function find_los(Predict_Sat $sat, Predict_QTH $qth, $start, $maxdt)
-    {
-        $t = $start;
-        $lostime = 0.0;
-
-
-        $this->predict_calc($sat, $qth, $start);
-
-        /* check whether satellite has aos */
-        if (($sat->otype == Predict_SGPSDP::ORBIT_TYPE_GEO) ||
-            ($sat->otype == Predict_SGPSDP::ORBIT_TYPE_DECAYED) ||
-            !$this->has_aos ($sat, $qth)) {
-
-            return 0.0;
-        }
-
-        if ($sat->el < 0.0) {
-            $t = $this->find_aos($sat, $qth, $start, $maxdt) + 0.001; // +1.5 min
-        }
-
-        /* invalid time (potentially returned by find_aos) */
-        if ($t < 0.01) {
-            return 0.0;
-        }
-
-        /* update satellite data */
-        $this->predict_calc($sat, $qth, $t);
-
-        /* use upper time limit */
-        if ($maxdt > 0.0) {
-
-            /* coarse steps */
-            while (($sat->el >= 1.0) && ($t <= ($start + $maxdt))) {
-                $t += cos(($sat->el - 1.0) * self::de2ra) * sqrt($sat->alt) / 25000.0;
-                $this->predict_calc($sat, $qth, $t);
-            }
-
-            /* fine steps */
-            while (($lostime == 0.0) && ($t <= ($start + $maxdt)))  {
-
-                $t += $sat->el * sqrt($sat->alt) / 502500.0;
-                $this->predict_calc($sat, $qth, $t);
-
-                if (abs($sat->el) < 0.005) {
-                    $lostime = $t;
-                }
-            }
-        } else {
-        /* don't use upper limit */
-
-            /* coarse steps */
-            while ($sat->el >= 1.0) {
-                $t += cos(($sat->el - 1.0) * self::de2ra) * sqrt($sat->alt) / 25000.0;
-                $this->predict_calc($sat, $qth, $t);
-            }
-
-            /* fine steps */
-            while ($lostime == 0.0) {
-
-                $t += $sat->el * sqrt($sat->alt) / 502500.0;
-                $this->predict_calc($sat, $qth, $t);
-
-                if (abs($sat->el) < 0.005)
-                    $lostime = $t;
-            }
-        }
-
-        return $lostime;
-    }
-
-    /** Find AOS time of current pass.
-     *  @param Predict_Sat $sat   The satellite to find AOS for.
-     *  @param Predict_QTH $qth   The ground station.
-     *  @param float       $start Start time, prefereably now.
-     *  @return The time of the previous AOS or 0.0 if the satellite has no AOS.
-     *
-     * This function can be used to find the AOS time in the past of the
-     * current pass.
-     */
-    public function find_prev_aos(Predict_Sat $sat, Predict_QTH $qth, $start)
-    {
-        $aostime = $start;
-
-        /* make sure current sat values are
+		$this->predict_calc($sat, $qth, $start);
+
+		/* check whether satellite has aos */
+		if (($sat->otype == Predict_SGPSDP::ORBIT_TYPE_GEO) ||
+			($sat->otype == Predict_SGPSDP::ORBIT_TYPE_DECAYED) ||
+			!$this->has_aos($sat, $qth)) {
+
+			return 0.0;
+		}
+
+		if ($sat->el > 0.0) {
+			$t = $this->find_los($sat, $qth, $start, $maxdt) + 0.014; // +20 min
+		}
+
+		/* invalid time (potentially returned by find_los) */
+		if ($t < 0.1) {
+			return 0.0;
+		}
+
+		/* update satellite data */
+		$this->predict_calc($sat, $qth, $t);
+
+		/* use upper time limit */
+		if ($maxdt > 0.0) {
+
+			/* coarse time steps */
+			while (($sat->el < -1.0) && ($t <= ($start + $maxdt))) {
+				$t -= 0.00035 * ($sat->el * (($sat->alt / 8400.0) + 0.46) - 2.0);
+				$this->predict_calc($sat, $qth, $t);
+			}
+
+			/* fine steps */
+			while (($aostime == 0.0) && ($t <= ($start + $maxdt))) {
+
+				if (abs($sat->el) < 0.005) {
+					$aostime = $t;
+				} else {
+					$t -= $sat->el * sqrt($sat->alt) / 530000.0;
+					$this->predict_calc($sat, $qth, $t);
+				}
+			}
+		} else {
+			/* don't use upper time limit */
+
+			/* coarse time steps */
+			while ($sat->el < -1.0) {
+
+				$t -= 0.00035 * ($sat->el * (($sat->alt / 8400.0) + 0.46) - 2.0);
+				$this->predict_calc($sat, $qth, $t);
+			}
+
+			/* fine steps */
+			while ($aostime == 0.0) {
+
+				if (abs($sat->el) < 0.005) {
+					$aostime = $t;
+				} else {
+					$t -= $sat->el * sqrt($sat->alt) / 530000.0;
+					$this->predict_calc($sat, $qth, $t);
+				}
+
+			}
+		}
+
+		return $aostime;
+	}
+
+	/** SGP4SDP4 driver for doing AOS/LOS calculations.
+	 *  @param Predict_Sat $sat The satellite data.
+	 *  @param Predict_QTH $qth The QTH observer location data.
+	 *  @param float       $t   The time for calculation (Julian Date)
+	 *
+	 */
+	public function predict_calc(Predict_Sat $sat, Predict_QTH $qth, $t)
+	{
+		$obs_set      = new Predict_ObsSet();
+		$sat_geodetic = new Predict_Geodetic();
+		$obs_geodetic = new Predict_Geodetic();
+
+		$obs_geodetic->lon   = $qth->lon * self::de2ra;
+		$obs_geodetic->lat   = $qth->lat * self::de2ra;
+		$obs_geodetic->alt   = $qth->alt / 1000.0;
+		$obs_geodetic->theta = 0;
+
+		$sat->jul_utc = $t;
+		$sat->tsince = ($sat->jul_utc - $sat->jul_epoch) * self::xmnpda;
+
+		/* call the norad routines according to the deep-space flag */
+		$sgpsdp = Predict_SGPSDP::getInstance($sat);
+		if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
+			$sgpsdp->SDP4($sat, $sat->tsince);
+		} else {
+			$sgpsdp->SGP4($sat, $sat->tsince);
+		}
+
+		Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
+
+		/* get the velocity of the satellite */
+		$sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
+		$sat->velo = $sat->vel->w;
+		Predict_SGPObs::Calculate_Obs($sat->jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
+		Predict_SGPObs::Calculate_LatLonAlt($sat->jul_utc, $sat->pos, $sat_geodetic);
+
+		while ($sat_geodetic->lon < -self::pi) {
+			$sat_geodetic->lon += self::twopi;
+		}
+
+		while ($sat_geodetic->lon > (self::pi)) {
+			$sat_geodetic->lon -= self::twopi;
+		}
+
+		$sat->az = Predict_Math::Degrees($obs_set->az);
+		$sat->el = Predict_Math::Degrees($obs_set->el);
+		$sat->range = $obs_set->range;
+		$sat->range_rate = $obs_set->range_rate;
+		$sat->ssplat = Predict_Math::Degrees($sat_geodetic->lat);
+		$sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
+		$sat->alt = $sat_geodetic->alt;
+		$sat->ma = Predict_Math::Degrees($sat->phase);
+		$sat->ma *= 256.0 / 360.0;
+		$sat->phase = Predict_Math::Degrees($sat->phase);
+
+		/* same formulas, but the one from predict is nicer */
+		//sat->footprint = 2.0 * xkmper * acos (xkmper/sat->pos.w);
+		$sat->footprint = 12756.33 * acos(self::xkmper / (self::xkmper + $sat->alt));
+		$age = $sat->jul_utc - $sat->jul_epoch;
+		$sat->orbit = floor(($sat->tle->xno * self::xmnpda / self::twopi +
+						$age * $sat->tle->bstar * self::ae) * $age +
+						$sat->tle->xmo / self::twopi) + $sat->tle->revnum - 1;
+	}
+
+	/** Find the LOS time of the next pass.
+	 *  @author Alexandru Csete, OZ9AEC
+	 *  @author John A. Magliacane, KD2BD
+	 *  @param Predict_Sat $sat The satellite data.
+	 *  @param Predict_QTH $qth The QTH observer location data.
+	 *  @param float       $start The time where calculation should start. (Julian Date)
+	 *  @param int         $maxdt The upper time limit in days (0.0 = no limit)
+	 *  @return The time (julian date) of the next LOS or 0.0 if the satellite has no LOS.
+	 *
+	 * This function finds the time of LOS for the first coming pass taking place
+	 * no earlier that start.
+	 * If the satellite is currently out of range, the function first calls
+	 * find_aos to get the next AOS time. Then the calculations are done using
+	 * the new start time.
+	 * The function has a built-in watchdog to ensure that we don't end up in
+	 * lengthy loops.
+	 *
+	 */
+	public function find_los(Predict_Sat $sat, Predict_QTH $qth, $start, $maxdt)
+	{
+		$t = $start;
+		$lostime = 0.0;
+
+
+		$this->predict_calc($sat, $qth, $start);
+
+		/* check whether satellite has aos */
+		if (($sat->otype == Predict_SGPSDP::ORBIT_TYPE_GEO) ||
+			($sat->otype == Predict_SGPSDP::ORBIT_TYPE_DECAYED) ||
+			!$this->has_aos ($sat, $qth)) {
+
+			return 0.0;
+		}
+
+		if ($sat->el < 0.0) {
+			$t = $this->find_aos($sat, $qth, $start, $maxdt) + 0.001; // +1.5 min
+		}
+
+		/* invalid time (potentially returned by find_aos) */
+		if ($t < 0.01) {
+			return 0.0;
+		}
+
+		/* update satellite data */
+		$this->predict_calc($sat, $qth, $t);
+
+		/* use upper time limit */
+		if ($maxdt > 0.0) {
+
+			/* coarse steps */
+			while (($sat->el >= 1.0) && ($t <= ($start + $maxdt))) {
+				$t += cos(($sat->el - 1.0) * self::de2ra) * sqrt($sat->alt) / 25000.0;
+				$this->predict_calc($sat, $qth, $t);
+			}
+
+			/* fine steps */
+			while (($lostime == 0.0) && ($t <= ($start + $maxdt)))  {
+
+				$t += $sat->el * sqrt($sat->alt) / 502500.0;
+				$this->predict_calc($sat, $qth, $t);
+
+				if (abs($sat->el) < 0.005) {
+					$lostime = $t;
+				}
+			}
+		} else {
+		/* don't use upper limit */
+
+			/* coarse steps */
+			while ($sat->el >= 1.0) {
+				$t += cos(($sat->el - 1.0) * self::de2ra) * sqrt($sat->alt) / 25000.0;
+				$this->predict_calc($sat, $qth, $t);
+			}
+
+			/* fine steps */
+			while ($lostime == 0.0) {
+
+				$t += $sat->el * sqrt($sat->alt) / 502500.0;
+				$this->predict_calc($sat, $qth, $t);
+
+				if (abs($sat->el) < 0.005)
+					$lostime = $t;
+			}
+		}
+
+		return $lostime;
+	}
+
+	/** Find AOS time of current pass.
+	 *  @param Predict_Sat $sat   The satellite to find AOS for.
+	 *  @param Predict_QTH $qth   The ground station.
+	 *  @param float       $start Start time, prefereably now.
+	 *  @return The time of the previous AOS or 0.0 if the satellite has no AOS.
+	 *
+	 * This function can be used to find the AOS time in the past of the
+	 * current pass.
+	 */
+	public function find_prev_aos(Predict_Sat $sat, Predict_QTH $qth, $start)
+	{
+		$aostime = $start;
+
+		/* make sure current sat values are
             in sync with the time
         */
-        $this->predict_calc($sat, $qth, $start);
-
-        /* check whether satellite has aos */
-        if (($sat->otype == Predict_SGPSDP::ORBIT_TYPE_GEO) ||
-            ($sat->otype == Predict_SGPSDP::ORBIT_TYPE_DECAYED) ||
-            !$this->has_aos($sat, $qth)) {
-
-            return 0.0;
-        }
-
-        while ($sat->el >= 0.0) {
-            $aostime -= 0.0005; // 0.75 min
-            $this->predict_calc($sat, $qth, $aostime);
-        }
-
-        return $aostime;
-    }
-
-    /** Determine whether satellite ever reaches AOS.
-     *  @author John A. Magliacane, KD2BD
-     *  @author Alexandru Csete, OZ9AEC
-     *  @param Predict_Sat $sat The satellite data.
-     *  @param Predict_QTH $qth The observer's location data
-     *  @return bool true if the satellite will reach AOS, false otherwise.
-     *
-     */
-    public function has_aos(Predict_Sat $sat, Predict_QTH $qth)
-    {
-         $retcode = false;
-
-         /* FIXME */
-         if ($sat->meanmo == 0.0) {
-              $retcode = false;
-         } else {
-
-            /* xincl is already in RAD by select_ephemeris */
-            $lin = $sat->tle->xincl;
-            if ($lin >= self::pio2) {
-                $lin = self::pi - $lin;
-            }
-
-            $sma = 331.25 * exp(log(1440.0 / $sat->meanmo) * (2.0 / 3.0));
-            $apogee = $sma * (1.0 + $sat->tle->eo) - self::xkmper;
-
-            if ((acos(self::xkmper / ($apogee + self::xkmper)) + ($lin)) > abs($qth->lat * self::de2ra)) {
-                $retcode = true;
-            } else {
-                $retcode = false;
-            }
-        }
-
-        return $retcode;
-    }
-
-    /** Predict passes after a certain time.
-     *
-     *
-     * This function calculates num upcoming passes with AOS no earlier
-     * than t = start and not later that t = (start+maxdt). The function will
-     *  repeatedly call get_pass until
-     * the number of predicted passes is equal to num, the time has reached
-     * limit or the get_pass function returns NULL.
-     *
-     * note For no time limit use maxdt = 0.0
-     *
-     * note the data in sat will be corrupt (future) and must be refreshed
-     *      by the caller, if the caller will need it later on (eg. if the caller
-     *      is GtkSatList).
-     *
-     * note Prepending to a singly linked list is much faster than appending.
-     *      Therefore, the elements are prepended whereafter the GSList is
-     *      reversed
-     *
-     *
-     * @param Predict_Sat  $sat The satellite data
-     * @param Predict_QTH  $qth The observer's location data
-     * @param float $start The start julian date
-     * @param int   $maxdt The max # of days to look
-     * @param int   $num   The max # of passes to get
-     * @return array of Predict_Pass instances if found, empty array otherwise
-     */
-    public function get_passes(Predict_Sat $sat, Predict_QTH $qth, $start, $maxdt, $num = 0)
-    {
-        $passes = array();
-
-        /* if no number has been specified
+		$this->predict_calc($sat, $qth, $start);
+
+		/* check whether satellite has aos */
+		if (($sat->otype == Predict_SGPSDP::ORBIT_TYPE_GEO) ||
+			($sat->otype == Predict_SGPSDP::ORBIT_TYPE_DECAYED) ||
+			!$this->has_aos($sat, $qth)) {
+
+			return 0.0;
+		}
+
+		while ($sat->el >= 0.0) {
+			$aostime -= 0.0005; // 0.75 min
+			$this->predict_calc($sat, $qth, $aostime);
+		}
+
+		return $aostime;
+	}
+
+	/** Determine whether satellite ever reaches AOS.
+	 *  @author John A. Magliacane, KD2BD
+	 *  @author Alexandru Csete, OZ9AEC
+	 *  @param Predict_Sat $sat The satellite data.
+	 *  @param Predict_QTH $qth The observer's location data
+	 *  @return bool true if the satellite will reach AOS, false otherwise.
+	 *
+	 */
+	public function has_aos(Predict_Sat $sat, Predict_QTH $qth)
+	{
+		 $retcode = false;
+
+		 /* FIXME */
+		 if ($sat->meanmo == 0.0) {
+			  $retcode = false;
+		 } else {
+
+			/* xincl is already in RAD by select_ephemeris */
+			$lin = $sat->tle->xincl;
+			if ($lin >= self::pio2) {
+				$lin = self::pi - $lin;
+			}
+
+			$sma = 331.25 * exp(log(1440.0 / $sat->meanmo) * (2.0 / 3.0));
+			$apogee = $sma * (1.0 + $sat->tle->eo) - self::xkmper;
+
+			if ((acos(self::xkmper / ($apogee + self::xkmper)) + ($lin)) > abs($qth->lat * self::de2ra)) {
+				$retcode = true;
+			} else {
+				$retcode = false;
+			}
+		}
+
+		return $retcode;
+	}
+
+	/** Predict passes after a certain time.
+	 *
+	 *
+	 * This function calculates num upcoming passes with AOS no earlier
+	 * than t = start and not later that t = (start+maxdt). The function will
+	 *  repeatedly call get_pass until
+	 * the number of predicted passes is equal to num, the time has reached
+	 * limit or the get_pass function returns NULL.
+	 *
+	 * note For no time limit use maxdt = 0.0
+	 *
+	 * note the data in sat will be corrupt (future) and must be refreshed
+	 *      by the caller, if the caller will need it later on (eg. if the caller
+	 *      is GtkSatList).
+	 *
+	 * note Prepending to a singly linked list is much faster than appending.
+	 *      Therefore, the elements are prepended whereafter the GSList is
+	 *      reversed
+	 *
+	 *
+	 * @param Predict_Sat  $sat The satellite data
+	 * @param Predict_QTH  $qth The observer's location data
+	 * @param float $start The start julian date
+	 * @param int   $maxdt The max # of days to look
+	 * @param int   $num   The max # of passes to get
+	 * @return array of Predict_Pass instances if found, empty array otherwise
+	 */
+	public function get_passes(Predict_Sat $sat, Predict_QTH $qth, $start, $maxdt, $num = 0)
+	{
+		$passes = array();
+
+		/* if no number has been specified
             set it to something big */
-        if ($num == 0) {
-            $num = 100;
-        }
+		if ($num == 0) {
+			$num = 100;
+		}
 
-        $t = $start;
+		$t = $start;
 
-        for ($i = 0; $i < $num; $i++) {
-            $pass = $this->get_pass($sat, $qth, $t, $maxdt);
+		for ($i = 0; $i < $num; $i++) {
+			$pass = $this->get_pass($sat, $qth, $t, $maxdt);
 
-            if ($pass != null) {
-                $passes[] = $pass;
-                $t = $pass->los + 0.014; // +20 min
+			if ($pass != null) {
+				$passes[] = $pass;
+				$t = $pass->los + 0.014; // +20 min
 
-                /* if maxdt > 0.0 check whether we have reached t = start+maxdt
+				/* if maxdt > 0.0 check whether we have reached t = start+maxdt
                     if yes finish predictions
                 */
-                if (($maxdt > 0.0) && ($t >= ($start + $maxdt))) {
-                    $i = $num;
-                }
-            } else {
-                /* we can't get any more passes */
-                $i = $num;
-            }
-        }
-
-        return $passes;
-    }
-
-    /**
-     * Filters out visible passes and adds the visible aos, tca, los, and
-     * corresponding az and ele for each.
-     *
-     * @param array $passes The passes returned from get_passes()
-     *
-     * @author Bill Shupp
-     * @return array
-     */
-    public function filterVisiblePasses(array $passes)
-    {
-        $filtered = array();
-
-        foreach ($passes as $result) {
-            // Dummy check
-            if ($result->vis[0] != 'V') {
-                continue;
-            }
-
-            $aos    = false;
-            $aos_az = false;
-            $aos    = false;
-            $tca    = false;
-            $los_az = false;
-            $max_el = 0;
-
-            foreach ($result->details as $detail) {
-                if ($detail->vis != Predict::SAT_VIS_VISIBLE) {
-                    continue;
-                }
-                if ($detail->el < $this->minEle) {
-                    continue;
-                }
-
-                if ($aos == false) {
-                    $aos       = $detail->time;
-                    $aos_az    = $detail->az;
-                    $aos_el    = $detail->el;
-                    $tca       = $detail->time;
-                    $los       = $detail->time;
-                    $los_az    = $detail->az;
-                    $los_el    = $detail->el;
-                    $max_el    = $detail->el;
-                    $max_el_az = $detail->el;
-                    continue;
-                }
-                $los    = $detail->time;
-                $los_az = $detail->az;
-                $los_el = $detail->el;
-
-                if ($detail->el > $max_el) {
-                    $tca       = $detail->time;
-                    $max_el    = $detail->el;
-                    $max_el_az = $detail->az;
-                }
-            }
-
-            if ($aos === false) {
-                // Does not reach minimum elevation, skip
-                continue;
-            }
-
-            $result->visible_aos       = $aos;
-            $result->visible_aos_az    = $aos_az;
-            $result->visible_aos_el    = $aos_el;
-            $result->visible_tca       = $tca;
-            $result->visible_max_el    = $max_el;
-            $result->visible_max_el_az = $max_el_az;
-            $result->visible_los       = $los;
-            $result->visible_los_az    = $los_az;
-            $result->visible_los_el    = $los_el;
-
-            $filtered[] = $result;
-        }
-
-        return $filtered;
-    }
-
-    /**
-     * Translates aziumuth degrees to compass direction:
-     *
-     * N (0°), NNE (22.5°), NE (45°), ENE (67.5°), E (90°), ESE (112.5°),
-     * SE (135°), SSE (157.5°), S (180°), SSW (202.5°), SW (225°),
-     * WSW (247.5°), W (270°), WNW (292.5°), NW (315°), NNW (337.5°)
-     *
-     * @param int $az The azimuth in degrees, defaults to 0
-     *
-     * @return string
-     */
-    public function azDegreesToDirection($az = 0)
-    {
-        $i = floor($az / 22.5);
-        $m = (22.5 * (2 * $i + 1)) / 2;
-        $i = ($az >= $m) ? $i + 1 : $i;
-
-        return trim(substr('N  NNENE ENEE  ESESE SSES  SSWSW WSWW  WNWNW NNWN  ', $i * 3, 3));
-    }
+				if (($maxdt > 0.0) && ($t >= ($start + $maxdt))) {
+					$i = $num;
+				}
+			} else {
+				/* we can't get any more passes */
+				$i = $num;
+			}
+		}
+
+		return $passes;
+	}
+
+	/**
+	 * Filters out visible passes and adds the visible aos, tca, los, and
+	 * corresponding az and ele for each.
+	 *
+	 * @param array $passes The passes returned from get_passes()
+	 *
+	 * @author Bill Shupp
+	 * @return array
+	 */
+	public function filterVisiblePasses(array $passes)
+	{
+		$filtered = array();
+
+		foreach ($passes as $result) {
+			// Dummy check
+			if ($result->vis[0] != 'V') {
+				continue;
+			}
+
+			$aos    = false;
+			$aos_az = false;
+			$aos    = false;
+			$tca    = false;
+			$los_az = false;
+			$max_el = 0;
+
+			foreach ($result->details as $detail) {
+				if ($detail->vis != Predict::SAT_VIS_VISIBLE) {
+					continue;
+				}
+				if ($detail->el < $this->minEle) {
+					continue;
+				}
+
+				if ($aos == false) {
+					$aos       = $detail->time;
+					$aos_az    = $detail->az;
+					$aos_el    = $detail->el;
+					$tca       = $detail->time;
+					$los       = $detail->time;
+					$los_az    = $detail->az;
+					$los_el    = $detail->el;
+					$max_el    = $detail->el;
+					$max_el_az = $detail->el;
+					continue;
+				}
+				$los    = $detail->time;
+				$los_az = $detail->az;
+				$los_el = $detail->el;
+
+				if ($detail->el > $max_el) {
+					$tca       = $detail->time;
+					$max_el    = $detail->el;
+					$max_el_az = $detail->az;
+				}
+			}
+
+			if ($aos === false) {
+				// Does not reach minimum elevation, skip
+				continue;
+			}
+
+			$result->visible_aos       = $aos;
+			$result->visible_aos_az    = $aos_az;
+			$result->visible_aos_el    = $aos_el;
+			$result->visible_tca       = $tca;
+			$result->visible_max_el    = $max_el;
+			$result->visible_max_el_az = $max_el_az;
+			$result->visible_los       = $los;
+			$result->visible_los_az    = $los_az;
+			$result->visible_los_el    = $los_el;
+
+			$filtered[] = $result;
+		}
+
+		return $filtered;
+	}
+
+	/**
+	 * Translates aziumuth degrees to compass direction:
+	 *
+	 * N (0°), NNE (22.5°), NE (45°), ENE (67.5°), E (90°), ESE (112.5°),
+	 * SE (135°), SSE (157.5°), S (180°), SSW (202.5°), SW (225°),
+	 * WSW (247.5°), W (270°), WNW (292.5°), NW (315°), NNW (337.5°)
+	 *
+	 * @param int $az The azimuth in degrees, defaults to 0
+	 *
+	 * @return string
+	 */
+	public function azDegreesToDirection($az = 0)
+	{
+		$i = floor($az / 22.5);
+		$m = (22.5 * (2 * $i + 1)) / 2;
+		$i = ($az >= $m) ? $i + 1 : $i;
+
+		return trim(substr('N  NNENE ENEE  ESESE SSES  SSWSW WSWW  WNWNW NNWN  ', $i * 3, 3));
+	}
 }

Spacing +93 added lines, -93 removed lines

@@ -52,60 +52,60 @@ 
  */
 class Predict
 {
-    const de2ra    =  1.74532925E-2;   /* Degrees to Radians */
-    const pi       =  3.1415926535898; /* Pi */
-    const pio2     =  1.5707963267949; /* Pi/2 */
-    const x3pio2   =  4.71238898;      /* 3*Pi/2 */
-    const twopi    =  6.2831853071796; /* 2*Pi  */
-    const e6a      =  1.0E-6;
-    const tothrd   =  6.6666667E-1;    /* 2/3 */
-    const xj2      =  1.0826158E-3;    /* J2 Harmonic */
-    const xj3      = -2.53881E-6;      /* J3 Harmonic */
-    const xj4      = -1.65597E-6;      /* J4 Harmonic */
-    const xke      =  7.43669161E-2;
-    const xkmper   =  6.378135E3;      /* Earth radius km */
-    const xmnpda   =  1.44E3;          /* Minutes per day */
-    const km2mi    =  0.621371;        /* Kilometers per Mile */
-    const ae       =  1.0;
-    const ck2      =  5.413079E-4;
-    const ck4      =  6.209887E-7;
-    const __f      =  3.352779E-3;
-    const ge       =  3.986008E5;
-    const __s__    =  1.012229;
-    const qoms2t   =  1.880279E-09;
-    const secday   =  8.6400E4;        /* Seconds per day */
-    const omega_E  =  1.0027379;
-    const omega_ER =  6.3003879;
-    const zns      =  1.19459E-5;
-    const c1ss     =  2.9864797E-6;
-    const zes      =  1.675E-2;
-    const znl      =  1.5835218E-4;
-    const c1l      =  4.7968065E-7;
-    const zel      =  5.490E-2;
-    const zcosis   =  9.1744867E-1;
-    const zsinis   =  3.9785416E-1;
+    const de2ra    = 1.74532925E-2; /* Degrees to Radians */
+    const pi       = 3.1415926535898; /* Pi */
+    const pio2     = 1.5707963267949; /* Pi/2 */
+    const x3pio2   = 4.71238898; /* 3*Pi/2 */
+    const twopi    = 6.2831853071796; /* 2*Pi  */
+    const e6a      = 1.0E-6;
+    const tothrd   = 6.6666667E-1; /* 2/3 */
+    const xj2      = 1.0826158E-3; /* J2 Harmonic */
+    const xj3      = -2.53881E-6; /* J3 Harmonic */
+    const xj4      = -1.65597E-6; /* J4 Harmonic */
+    const xke      = 7.43669161E-2;
+    const xkmper   = 6.378135E3; /* Earth radius km */
+    const xmnpda   = 1.44E3; /* Minutes per day */
+    const km2mi    = 0.621371; /* Kilometers per Mile */
+    const ae       = 1.0;
+    const ck2      = 5.413079E-4;
+    const ck4      = 6.209887E-7;
+    const __f      = 3.352779E-3;
+    const ge       = 3.986008E5;
+    const __s__    = 1.012229;
+    const qoms2t   = 1.880279E-09;
+    const secday   = 8.6400E4; /* Seconds per day */
+    const omega_E  = 1.0027379;
+    const omega_ER = 6.3003879;
+    const zns      = 1.19459E-5;
+    const c1ss     = 2.9864797E-6;
+    const zes      = 1.675E-2;
+    const znl      = 1.5835218E-4;
+    const c1l      = 4.7968065E-7;
+    const zel      = 5.490E-2;
+    const zcosis   = 9.1744867E-1;
+    const zsinis   = 3.9785416E-1;
     const zsings   = -9.8088458E-1;
-    const zcosgs   =  1.945905E-1;
-    const zcoshs   =  1;
-    const zsinhs   =  0;
-    const q22      =  1.7891679E-6;
-    const q31      =  2.1460748E-6;
-    const q33      =  2.2123015E-7;
-    const g22      =  5.7686396;
-    const g32      =  9.5240898E-1;
-    const g44      =  1.8014998;
-    const g52      =  1.0508330;
-    const g54      =  4.4108898;
-    const root22   =  1.7891679E-6;
-    const root32   =  3.7393792E-7;
-    const root44   =  7.3636953E-9;
-    const root52   =  1.1428639E-7;
-    const root54   =  2.1765803E-9;
-    const thdt     =  4.3752691E-3;
-    const rho      =  1.5696615E-1;
-    const mfactor  =  7.292115E-5;
-    const __sr__   =  6.96000E5;      /*Solar radius - kilometers (IAU 76)*/
-    const AU       =  1.49597870E8;   /*Astronomical unit - kilometers (IAU 76)*/
+    const zcosgs   = 1.945905E-1;
+    const zcoshs   = 1;
+    const zsinhs   = 0;
+    const q22      = 1.7891679E-6;
+    const q31      = 2.1460748E-6;
+    const q33      = 2.2123015E-7;
+    const g22      = 5.7686396;
+    const g32      = 9.5240898E-1;
+    const g44      = 1.8014998;
+    const g52      = 1.0508330;
+    const g54      = 4.4108898;
+    const root22   = 1.7891679E-6;
+    const root32   = 3.7393792E-7;
+    const root44   = 7.3636953E-9;
+    const root52   = 1.1428639E-7;
+    const root54   = 2.1765803E-9;
+    const thdt     = 4.3752691E-3;
+    const rho      = 1.5696615E-1;
+    const mfactor  = 7.292115E-5;
+    const __sr__   = 6.96000E5; /*Solar radius - kilometers (IAU 76)*/
+    const AU       = 1.49597870E8; /*Astronomical unit - kilometers (IAU 76)*/
 
     /* visibility constants */
     const SAT_VIS_NONE     = 0;
@@ -163,18 +163,18 @@ 
      */
     public function get_pass(Predict_Sat $sat_in, Predict_QTH $qth, $start, $maxdt)
     {
-        $aos = 0.0;    /* time of AOS */
-        $tca = 0.0;    /* time of TCA */
-        $los = 0.0;    /* time of LOS */
-        $dt = 0.0;     /* time diff */
-        $step = 0.0;   /* time step */
+        $aos = 0.0; /* time of AOS */
+        $tca = 0.0; /* time of TCA */
+        $los = 0.0; /* time of LOS */
+        $dt = 0.0; /* time diff */
+        $step = 0.0; /* time step */
         $t0 = $start;
-        $tres = 0.0;   /* required time resolution */
+        $tres = 0.0; /* required time resolution */
         $max_el = 0.0; /* maximum elevation */
         $pass = null;
         $detail = null;
         $done = false;
-        $iter = 0;      /* number of iterations */
+        $iter = 0; /* number of iterations */
         /* FIXME: watchdog */
 
         /*copy sat_in to a working structure*/
@@ -182,7 +182,7 @@ 
         $sat_working = clone $sat_in;
 
         /* get time resolution; sat-cfg stores it in seconds */
-        $tres = $this->timeRes / 86400.0;
+        $tres = $this->timeRes/86400.0;
 
         /* loop until we find a pass with elevation > SAT_CFG_INT_PRED_MIN_EL
             or we run out of time
@@ -207,7 +207,7 @@ 
             /* aos = 0.0 means no aos */
             if ($aos == 0.0) {
                 $done = true;
-            } else if (($maxdt > 0.0) && ($aos > ($start + $maxdt)) ) {
+            } else if (($maxdt > 0.0) && ($aos > ($start + $maxdt))) {
                 /* check whether we are within time limits;
                     maxdt = 0 mean no time limit.
                 */
@@ -217,7 +217,7 @@ 
                 $dt = $los - $aos;
 
                 /* get time step, which will give us the max number of entries */
-                $step = $dt / $this->numEntries;
+                $step = $dt/$this->numEntries;
 
                 /* but if this is smaller than the required resolution
                     we go with the resolution
@@ -369,9 +369,9 @@ 
         $solar_set = new Predict_ObsSet();
 
         /* FIXME: could be passed as parameter */
-        $obs_geodetic->lon   = $qth->lon * self::de2ra;
-        $obs_geodetic->lat   = $qth->lat * self::de2ra;
-        $obs_geodetic->alt   = $qth->alt / 1000.0;
+        $obs_geodetic->lon   = $qth->lon*self::de2ra;
+        $obs_geodetic->lat   = $qth->lat*self::de2ra;
+        $obs_geodetic->alt   = $qth->alt/1000.0;
         $obs_geodetic->theta = 0;
 
         Predict_Solar::Calculate_Solar_Position($jul_utc, $solar_vector);
@@ -452,7 +452,7 @@ 
 
             /* coarse time steps */
             while (($sat->el < -1.0) && ($t <= ($start + $maxdt))) {
-                $t -= 0.00035 * ($sat->el * (($sat->alt / 8400.0) + 0.46) - 2.0);
+                $t -= 0.00035*($sat->el*(($sat->alt/8400.0) + 0.46) - 2.0);
                 $this->predict_calc($sat, $qth, $t);
             }
 
@@ -462,7 +462,7 @@ 
                 if (abs($sat->el) < 0.005) {
                     $aostime = $t;
                 } else {
-                    $t -= $sat->el * sqrt($sat->alt) / 530000.0;
+                    $t -= $sat->el*sqrt($sat->alt)/530000.0;
                     $this->predict_calc($sat, $qth, $t);
                 }
             }
@@ -472,7 +472,7 @@ 
             /* coarse time steps */
             while ($sat->el < -1.0) {
 
-                $t -= 0.00035 * ($sat->el * (($sat->alt / 8400.0) + 0.46) - 2.0);
+                $t -= 0.00035*($sat->el*(($sat->alt/8400.0) + 0.46) - 2.0);
                 $this->predict_calc($sat, $qth, $t);
             }
 
@@ -482,7 +482,7 @@ 
                 if (abs($sat->el) < 0.005) {
                     $aostime = $t;
                 } else {
-                    $t -= $sat->el * sqrt($sat->alt) / 530000.0;
+                    $t -= $sat->el*sqrt($sat->alt)/530000.0;
                     $this->predict_calc($sat, $qth, $t);
                 }
 
@@ -504,17 +504,17 @@ 
         $sat_geodetic = new Predict_Geodetic();
         $obs_geodetic = new Predict_Geodetic();
 
-        $obs_geodetic->lon   = $qth->lon * self::de2ra;
-        $obs_geodetic->lat   = $qth->lat * self::de2ra;
-        $obs_geodetic->alt   = $qth->alt / 1000.0;
+        $obs_geodetic->lon   = $qth->lon*self::de2ra;
+        $obs_geodetic->lat   = $qth->lat*self::de2ra;
+        $obs_geodetic->alt   = $qth->alt/1000.0;
         $obs_geodetic->theta = 0;
 
         $sat->jul_utc = $t;
-        $sat->tsince = ($sat->jul_utc - $sat->jul_epoch) * self::xmnpda;
+        $sat->tsince = ($sat->jul_utc - $sat->jul_epoch)*self::xmnpda;
 
         /* call the norad routines according to the deep-space flag */
         $sgpsdp = Predict_SGPSDP::getInstance($sat);
-        if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
+        if ($sat->flags&Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
             $sgpsdp->SDP4($sat, $sat->tsince);
         } else {
             $sgpsdp->SGP4($sat, $sat->tsince);
@@ -523,7 +523,7 @@ 
         Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
 
         /* get the velocity of the satellite */
-        $sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
+        $sat->vel->w = sqrt($sat->vel->x*$sat->vel->x + $sat->vel->y*$sat->vel->y + $sat->vel->z*$sat->vel->z);
         $sat->velo = $sat->vel->w;
         Predict_SGPObs::Calculate_Obs($sat->jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
         Predict_SGPObs::Calculate_LatLonAlt($sat->jul_utc, $sat->pos, $sat_geodetic);
@@ -544,16 +544,16 @@ 
         $sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
         $sat->alt = $sat_geodetic->alt;
         $sat->ma = Predict_Math::Degrees($sat->phase);
-        $sat->ma *= 256.0 / 360.0;
+        $sat->ma *= 256.0/360.0;
         $sat->phase = Predict_Math::Degrees($sat->phase);
 
         /* same formulas, but the one from predict is nicer */
         //sat->footprint = 2.0 * xkmper * acos (xkmper/sat->pos.w);
-        $sat->footprint = 12756.33 * acos(self::xkmper / (self::xkmper + $sat->alt));
+        $sat->footprint = 12756.33*acos(self::xkmper/(self::xkmper + $sat->alt));
         $age = $sat->jul_utc - $sat->jul_epoch;
-        $sat->orbit = floor(($sat->tle->xno * self::xmnpda / self::twopi +
-                        $age * $sat->tle->bstar * self::ae) * $age +
-                        $sat->tle->xmo / self::twopi) + $sat->tle->revnum - 1;
+        $sat->orbit = floor(($sat->tle->xno*self::xmnpda/self::twopi +
+                        $age*$sat->tle->bstar*self::ae)*$age +
+                        $sat->tle->xmo/self::twopi) + $sat->tle->revnum - 1;
     }
 
     /** Find the LOS time of the next pass.
@@ -585,7 +585,7 @@ 
         /* check whether satellite has aos */
         if (($sat->otype == Predict_SGPSDP::ORBIT_TYPE_GEO) ||
             ($sat->otype == Predict_SGPSDP::ORBIT_TYPE_DECAYED) ||
-            !$this->has_aos ($sat, $qth)) {
+            !$this->has_aos($sat, $qth)) {
 
             return 0.0;
         }
@@ -607,14 +607,14 @@ 
 
             /* coarse steps */
             while (($sat->el >= 1.0) && ($t <= ($start + $maxdt))) {
-                $t += cos(($sat->el - 1.0) * self::de2ra) * sqrt($sat->alt) / 25000.0;
+                $t += cos(($sat->el - 1.0)*self::de2ra)*sqrt($sat->alt)/25000.0;
                 $this->predict_calc($sat, $qth, $t);
             }
 
             /* fine steps */
-            while (($lostime == 0.0) && ($t <= ($start + $maxdt)))  {
+            while (($lostime == 0.0) && ($t <= ($start + $maxdt))) {
 
-                $t += $sat->el * sqrt($sat->alt) / 502500.0;
+                $t += $sat->el*sqrt($sat->alt)/502500.0;
                 $this->predict_calc($sat, $qth, $t);
 
                 if (abs($sat->el) < 0.005) {
@@ -626,14 +626,14 @@ 
 
             /* coarse steps */
             while ($sat->el >= 1.0) {
-                $t += cos(($sat->el - 1.0) * self::de2ra) * sqrt($sat->alt) / 25000.0;
+                $t += cos(($sat->el - 1.0)*self::de2ra)*sqrt($sat->alt)/25000.0;
                 $this->predict_calc($sat, $qth, $t);
             }
 
             /* fine steps */
             while ($lostime == 0.0) {
 
-                $t += $sat->el * sqrt($sat->alt) / 502500.0;
+                $t += $sat->el*sqrt($sat->alt)/502500.0;
                 $this->predict_calc($sat, $qth, $t);
 
                 if (abs($sat->el) < 0.005)
@@ -701,10 +701,10 @@ 
                 $lin = self::pi - $lin;
             }
 
-            $sma = 331.25 * exp(log(1440.0 / $sat->meanmo) * (2.0 / 3.0));
-            $apogee = $sma * (1.0 + $sat->tle->eo) - self::xkmper;
+            $sma = 331.25*exp(log(1440.0/$sat->meanmo)*(2.0/3.0));
+            $apogee = $sma*(1.0 + $sat->tle->eo) - self::xkmper;
 
-            if ((acos(self::xkmper / ($apogee + self::xkmper)) + ($lin)) > abs($qth->lat * self::de2ra)) {
+            if ((acos(self::xkmper/($apogee + self::xkmper)) + ($lin)) > abs($qth->lat*self::de2ra)) {
                 $retcode = true;
             } else {
                 $retcode = false;
@@ -866,10 +866,10 @@ 
      */
     public function azDegreesToDirection($az = 0)
     {
-        $i = floor($az / 22.5);
-        $m = (22.5 * (2 * $i + 1)) / 2;
+        $i = floor($az/22.5);
+        $m = (22.5*(2*$i + 1))/2;
         $i = ($az >= $m) ? $i + 1 : $i;
 
-        return trim(substr('N  NNENE ENEE  ESESE SSES  SSWSW WSWW  WNWNW NNWN  ', $i * 3, 3));
+        return trim(substr('N  NNENE ENEE  ESESE SSES  SSWSW WSWW  WNWNW NNWN  ', $i*3, 3));
     }
 }

Braces +3 added lines, -2 removed lines

@@ -636,8 +636,9 @@
                 $t += $sat->el * sqrt($sat->alt) / 502500.0;
                 $this->predict_calc($sat, $qth, $t);
 
-                if (abs($sat->el) < 0.005)
-                    $lostime = $t;
+                if (abs($sat->el) < 0.005) {
+                                    $lostime = $t;
+                }
             }
         }
 

require/libs/Predict/Predict/Math.php

Doc Comments +16 added lines

@@ -42,6 +42,10 @@ 
     }
 
     /* Returns arccosine of rgument */
+
+    /**
+     * @param double $arg
+     */
     public static function ArcCos($arg)
     {
         return Predict::pio2 - self::ArcSin($arg);
@@ -68,6 +72,10 @@ 
     }
 
     /* Multiplies the vector v1 by the scalar k to produce the vector v2 */
+
+    /**
+     * @param integer $k
+     */
     public static function Scalar_Multiply($k, Predict_Vector $v1, Predict_Vector $v2)
     {
         $v2->x = $k * $v1->x;
@@ -155,6 +163,10 @@ 
     }
 
     /* Returns arg1 mod arg2 */
+
+    /**
+     * @param double $arg1
+     */
     public static function Modulus($arg1, $arg2)
     {
         $ret_val  = $arg1;
@@ -169,6 +181,10 @@ 
     }
 
     /* Returns fractional part of double argument */
+
+    /**
+     * @param double $arg
+     */
     public static function Frac($arg)
     {
         return $arg - floor($arg);

Indentation +174 added lines, -174 removed lines

@@ -19,178 +19,178 @@
  */
 class Predict_Math
 {
-    /* Returns sign of a float */
-    public static function Sign($arg)
-    {
-        if ($arg > 0 ) {
-            return 1;
-        } else if ($arg < 0 ) {
-            return -1;
-        } else {
-            return 0;
-        }
-    }
-
-    /* Returns the arcsine of the argument */
-    public static function ArcSin($arg)
-    {
-        if (abs($arg) >= 1 ) {
-            return (self::Sign($arg) * Predict::pio2);
-        } else {
-            return(atan($arg / sqrt(1 - $arg * $arg)));
-        }
-    }
-
-    /* Returns arccosine of rgument */
-    public static function ArcCos($arg)
-    {
-        return Predict::pio2 - self::ArcSin($arg);
-    }
-
-    /* Adds vectors v1 and v2 together to produce v3 */
-    public static function Vec_Add(Predict_Vector $v1, Predict_Vector $v2, Predict_Vector $v3)
-    {
-        $v3->x = $v1->x + $v2->x;
-        $v3->y = $v1->y + $v2->y;
-        $v3->z = $v1->z + $v2->z;
-
-        $v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
-    }
-
-    /* Subtracts vector v2 from v1 to produce v3 */
-    public static function Vec_Sub(Predict_Vector $v1, Predict_Vector $v2, Predict_Vector $v3)
-    {
-        $v3->x = $v1->x - $v2->x;
-        $v3->y = $v1->y - $v2->y;
-        $v3->z = $v1->z - $v2->z;
-
-        $v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
-    }
-
-    /* Multiplies the vector v1 by the scalar k to produce the vector v2 */
-    public static function Scalar_Multiply($k, Predict_Vector $v1, Predict_Vector $v2)
-    {
-        $v2->x = $k * $v1->x;
-        $v2->y = $k * $v1->y;
-        $v2->z = $k * $v1->z;
-        $v2->w = abs($k) * $v1->w;
-    }
-
-    /* Multiplies the vector v1 by the scalar k */
-    public static function Scale_Vector($k, Predict_Vector $v)
-    {
-        $v->x *= $k;
-        $v->y *= $k;
-        $v->z *= $k;
-
-        $v->w = sqrt($v->x * $v->x + $v->y * $v->y + $v->z * $v->z);
-    }
-
-    /* Returns the dot product of two vectors */
-    public static function Dot(Predict_Vector $v1, Predict_Vector $v2)
-    {
-        return ($v1->x * $v2->x + $v1->y * $v2->y + $v1->z * $v2->z);
-    }
-
-    /* Calculates the angle between vectors v1 and v2 */
-    public static function Angle(Predict_Vector $v1, Predict_Vector $v2)
-    {
-        $v1->w = sqrt($v1->x * $v1->x + $v1->y * $v1->y + $v1->z * $v1->z);
-        $v2->w = sqrt($v2->x * $v2->x + $v2->y * $v2->y + $v2->z * $v2->z);
-        return (self::ArcCos(self::Dot($v1, $v2) / ($v1->w * $v2->w)));
-    }
-
-    /* Produces cross product of v1 and v2, and returns in v3 */
-    public static function Cross(Predict_Vector $v1, Predict_Vector $v2 ,Predict_Vector $v3)
-    {
-        $v3->x = $v1->y * $v2->z - $v1->z * $v2->y;
-        $v3->y = $v1->z * $v2->x - $v1->x * $v2->z;
-        $v3->z = $v1->x * $v2->y - $v1->y * $v2->x;
-
-        $v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
-    }
-
-    /* Normalizes a vector */
-    public static function Normalize(Predict_Vector $v )
-    {
-        $v->x /= $v->w;
-        $v->y /= $v->w;
-        $v->z /= $v->w;
-    }
-
-    /* Four-quadrant arctan function */
-    public static function AcTan($sinx, $cosx)
-    {
-        if ($cosx == 0) {
-            if ($sinx > 0) {
-                return Predict::pio2;
-            } else {
-                return Predict::x3pio2;
-            }
-        } else {
-            if ($cosx > 0) {
-                if ($sinx > 0) {
-                    return atan($sinx / $cosx);
-                } else {
-                    return Predict::twopi + atan($sinx / $cosx);
-                }
-            } else {
-                return Predict::pi + atan($sinx / $cosx);
-            }
-        }
-    }
-
-    /* Returns mod 2pi of argument */
-    public static function FMod2p($x)
-    {
-        $ret_val  = $x;
-        $i        = (int) ($ret_val / Predict::twopi);
-        $ret_val -= $i * Predict::twopi;
-
-        if ($ret_val < 0) {
-            $ret_val += Predict::twopi;
-        }
-
-        return $ret_val;
-    }
-
-    /* Returns arg1 mod arg2 */
-    public static function Modulus($arg1, $arg2)
-    {
-        $ret_val  = $arg1;
-        $i        = (int) ($ret_val / $arg2);
-        $ret_val -= $i * $arg2;
-
-        if ($ret_val < 0) {
-            $ret_val += $arg2;
-        }
-
-        return $ret_val;
-    }
-
-    /* Returns fractional part of double argument */
-    public static function Frac($arg)
-    {
-        return $arg - floor($arg);
-    }
-
-    /* Converts the satellite's position and velocity  */
-    /* vectors from normalised values to km and km/sec */
-    public static function Convert_Sat_State(Predict_Vector $pos, Predict_Vector $vel)
-    {
-        self::Scale_Vector(Predict::xkmper, $pos);
-        self::Scale_Vector(Predict::xkmper * Predict::xmnpda / Predict::secday, $vel);
-    }
-
-    /* Returns angle in radians from arg in degrees */
-    public static function Radians($arg)
-    {
-        return $arg * Predict::de2ra;
-    }
-
-    /* Returns angle in degrees from arg in rads */
-    public static function Degrees($arg)
-    {
-      return $arg / Predict::de2ra;
-    }
+	/* Returns sign of a float */
+	public static function Sign($arg)
+	{
+		if ($arg > 0 ) {
+			return 1;
+		} else if ($arg < 0 ) {
+			return -1;
+		} else {
+			return 0;
+		}
+	}
+
+	/* Returns the arcsine of the argument */
+	public static function ArcSin($arg)
+	{
+		if (abs($arg) >= 1 ) {
+			return (self::Sign($arg) * Predict::pio2);
+		} else {
+			return(atan($arg / sqrt(1 - $arg * $arg)));
+		}
+	}
+
+	/* Returns arccosine of rgument */
+	public static function ArcCos($arg)
+	{
+		return Predict::pio2 - self::ArcSin($arg);
+	}
+
+	/* Adds vectors v1 and v2 together to produce v3 */
+	public static function Vec_Add(Predict_Vector $v1, Predict_Vector $v2, Predict_Vector $v3)
+	{
+		$v3->x = $v1->x + $v2->x;
+		$v3->y = $v1->y + $v2->y;
+		$v3->z = $v1->z + $v2->z;
+
+		$v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
+	}
+
+	/* Subtracts vector v2 from v1 to produce v3 */
+	public static function Vec_Sub(Predict_Vector $v1, Predict_Vector $v2, Predict_Vector $v3)
+	{
+		$v3->x = $v1->x - $v2->x;
+		$v3->y = $v1->y - $v2->y;
+		$v3->z = $v1->z - $v2->z;
+
+		$v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
+	}
+
+	/* Multiplies the vector v1 by the scalar k to produce the vector v2 */
+	public static function Scalar_Multiply($k, Predict_Vector $v1, Predict_Vector $v2)
+	{
+		$v2->x = $k * $v1->x;
+		$v2->y = $k * $v1->y;
+		$v2->z = $k * $v1->z;
+		$v2->w = abs($k) * $v1->w;
+	}
+
+	/* Multiplies the vector v1 by the scalar k */
+	public static function Scale_Vector($k, Predict_Vector $v)
+	{
+		$v->x *= $k;
+		$v->y *= $k;
+		$v->z *= $k;
+
+		$v->w = sqrt($v->x * $v->x + $v->y * $v->y + $v->z * $v->z);
+	}
+
+	/* Returns the dot product of two vectors */
+	public static function Dot(Predict_Vector $v1, Predict_Vector $v2)
+	{
+		return ($v1->x * $v2->x + $v1->y * $v2->y + $v1->z * $v2->z);
+	}
+
+	/* Calculates the angle between vectors v1 and v2 */
+	public static function Angle(Predict_Vector $v1, Predict_Vector $v2)
+	{
+		$v1->w = sqrt($v1->x * $v1->x + $v1->y * $v1->y + $v1->z * $v1->z);
+		$v2->w = sqrt($v2->x * $v2->x + $v2->y * $v2->y + $v2->z * $v2->z);
+		return (self::ArcCos(self::Dot($v1, $v2) / ($v1->w * $v2->w)));
+	}
+
+	/* Produces cross product of v1 and v2, and returns in v3 */
+	public static function Cross(Predict_Vector $v1, Predict_Vector $v2 ,Predict_Vector $v3)
+	{
+		$v3->x = $v1->y * $v2->z - $v1->z * $v2->y;
+		$v3->y = $v1->z * $v2->x - $v1->x * $v2->z;
+		$v3->z = $v1->x * $v2->y - $v1->y * $v2->x;
+
+		$v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
+	}
+
+	/* Normalizes a vector */
+	public static function Normalize(Predict_Vector $v )
+	{
+		$v->x /= $v->w;
+		$v->y /= $v->w;
+		$v->z /= $v->w;
+	}
+
+	/* Four-quadrant arctan function */
+	public static function AcTan($sinx, $cosx)
+	{
+		if ($cosx == 0) {
+			if ($sinx > 0) {
+				return Predict::pio2;
+			} else {
+				return Predict::x3pio2;
+			}
+		} else {
+			if ($cosx > 0) {
+				if ($sinx > 0) {
+					return atan($sinx / $cosx);
+				} else {
+					return Predict::twopi + atan($sinx / $cosx);
+				}
+			} else {
+				return Predict::pi + atan($sinx / $cosx);
+			}
+		}
+	}
+
+	/* Returns mod 2pi of argument */
+	public static function FMod2p($x)
+	{
+		$ret_val  = $x;
+		$i        = (int) ($ret_val / Predict::twopi);
+		$ret_val -= $i * Predict::twopi;
+
+		if ($ret_val < 0) {
+			$ret_val += Predict::twopi;
+		}
+
+		return $ret_val;
+	}
+
+	/* Returns arg1 mod arg2 */
+	public static function Modulus($arg1, $arg2)
+	{
+		$ret_val  = $arg1;
+		$i        = (int) ($ret_val / $arg2);
+		$ret_val -= $i * $arg2;
+
+		if ($ret_val < 0) {
+			$ret_val += $arg2;
+		}
+
+		return $ret_val;
+	}
+
+	/* Returns fractional part of double argument */
+	public static function Frac($arg)
+	{
+		return $arg - floor($arg);
+	}
+
+	/* Converts the satellite's position and velocity  */
+	/* vectors from normalised values to km and km/sec */
+	public static function Convert_Sat_State(Predict_Vector $pos, Predict_Vector $vel)
+	{
+		self::Scale_Vector(Predict::xkmper, $pos);
+		self::Scale_Vector(Predict::xkmper * Predict::xmnpda / Predict::secday, $vel);
+	}
+
+	/* Returns angle in radians from arg in degrees */
+	public static function Radians($arg)
+	{
+		return $arg * Predict::de2ra;
+	}
+
+	/* Returns angle in degrees from arg in rads */
+	public static function Degrees($arg)
+	{
+	  return $arg / Predict::de2ra;
+	}
 }

Spacing +32 added lines, -32 removed lines

@@ -22,9 +22,9 @@ 
     /* Returns sign of a float */
     public static function Sign($arg)
     {
-        if ($arg > 0 ) {
+        if ($arg > 0) {
             return 1;
-        } else if ($arg < 0 ) {
+        } else if ($arg < 0) {
             return -1;
         } else {
             return 0;
@@ -34,10 +34,10 @@ 
     /* Returns the arcsine of the argument */
     public static function ArcSin($arg)
     {
-        if (abs($arg) >= 1 ) {
-            return (self::Sign($arg) * Predict::pio2);
+        if (abs($arg) >= 1) {
+            return (self::Sign($arg)*Predict::pio2);
         } else {
-            return(atan($arg / sqrt(1 - $arg * $arg)));
+            return(atan($arg/sqrt(1 - $arg*$arg)));
         }
     }
 
@@ -54,7 +54,7 @@ 
         $v3->y = $v1->y + $v2->y;
         $v3->z = $v1->z + $v2->z;
 
-        $v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
+        $v3->w = sqrt($v3->x*$v3->x + $v3->y*$v3->y + $v3->z*$v3->z);
     }
 
     /* Subtracts vector v2 from v1 to produce v3 */
@@ -64,16 +64,16 @@ 
         $v3->y = $v1->y - $v2->y;
         $v3->z = $v1->z - $v2->z;
 
-        $v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
+        $v3->w = sqrt($v3->x*$v3->x + $v3->y*$v3->y + $v3->z*$v3->z);
     }
 
     /* Multiplies the vector v1 by the scalar k to produce the vector v2 */
     public static function Scalar_Multiply($k, Predict_Vector $v1, Predict_Vector $v2)
     {
-        $v2->x = $k * $v1->x;
-        $v2->y = $k * $v1->y;
-        $v2->z = $k * $v1->z;
-        $v2->w = abs($k) * $v1->w;
+        $v2->x = $k*$v1->x;
+        $v2->y = $k*$v1->y;
+        $v2->z = $k*$v1->z;
+        $v2->w = abs($k)*$v1->w;
     }
 
     /* Multiplies the vector v1 by the scalar k */
@@ -83,35 +83,35 @@ 
         $v->y *= $k;
         $v->z *= $k;
 
-        $v->w = sqrt($v->x * $v->x + $v->y * $v->y + $v->z * $v->z);
+        $v->w = sqrt($v->x*$v->x + $v->y*$v->y + $v->z*$v->z);
     }
 
     /* Returns the dot product of two vectors */
     public static function Dot(Predict_Vector $v1, Predict_Vector $v2)
     {
-        return ($v1->x * $v2->x + $v1->y * $v2->y + $v1->z * $v2->z);
+        return ($v1->x*$v2->x + $v1->y*$v2->y + $v1->z*$v2->z);
     }
 
     /* Calculates the angle between vectors v1 and v2 */
     public static function Angle(Predict_Vector $v1, Predict_Vector $v2)
     {
-        $v1->w = sqrt($v1->x * $v1->x + $v1->y * $v1->y + $v1->z * $v1->z);
-        $v2->w = sqrt($v2->x * $v2->x + $v2->y * $v2->y + $v2->z * $v2->z);
-        return (self::ArcCos(self::Dot($v1, $v2) / ($v1->w * $v2->w)));
+        $v1->w = sqrt($v1->x*$v1->x + $v1->y*$v1->y + $v1->z*$v1->z);
+        $v2->w = sqrt($v2->x*$v2->x + $v2->y*$v2->y + $v2->z*$v2->z);
+        return (self::ArcCos(self::Dot($v1, $v2)/($v1->w*$v2->w)));
     }
 
     /* Produces cross product of v1 and v2, and returns in v3 */
-    public static function Cross(Predict_Vector $v1, Predict_Vector $v2 ,Predict_Vector $v3)
+    public static function Cross(Predict_Vector $v1, Predict_Vector $v2, Predict_Vector $v3)
     {
-        $v3->x = $v1->y * $v2->z - $v1->z * $v2->y;
-        $v3->y = $v1->z * $v2->x - $v1->x * $v2->z;
-        $v3->z = $v1->x * $v2->y - $v1->y * $v2->x;
+        $v3->x = $v1->y*$v2->z - $v1->z*$v2->y;
+        $v3->y = $v1->z*$v2->x - $v1->x*$v2->z;
+        $v3->z = $v1->x*$v2->y - $v1->y*$v2->x;
 
-        $v3->w = sqrt($v3->x * $v3->x + $v3->y * $v3->y + $v3->z * $v3->z);
+        $v3->w = sqrt($v3->x*$v3->x + $v3->y*$v3->y + $v3->z*$v3->z);
     }
 
     /* Normalizes a vector */
-    public static function Normalize(Predict_Vector $v )
+    public static function Normalize(Predict_Vector $v)
     {
         $v->x /= $v->w;
         $v->y /= $v->w;
@@ -130,12 +130,12 @@ 
         } else {
             if ($cosx > 0) {
                 if ($sinx > 0) {
-                    return atan($sinx / $cosx);
+                    return atan($sinx/$cosx);
                 } else {
-                    return Predict::twopi + atan($sinx / $cosx);
+                    return Predict::twopi + atan($sinx/$cosx);
                 }
             } else {
-                return Predict::pi + atan($sinx / $cosx);
+                return Predict::pi + atan($sinx/$cosx);
             }
         }
     }
@@ -144,8 +144,8 @@ 
     public static function FMod2p($x)
     {
         $ret_val  = $x;
-        $i        = (int) ($ret_val / Predict::twopi);
-        $ret_val -= $i * Predict::twopi;
+        $i        = (int) ($ret_val/Predict::twopi);
+        $ret_val -= $i*Predict::twopi;
 
         if ($ret_val < 0) {
             $ret_val += Predict::twopi;
@@ -158,8 +158,8 @@ 
     public static function Modulus($arg1, $arg2)
     {
         $ret_val  = $arg1;
-        $i        = (int) ($ret_val / $arg2);
-        $ret_val -= $i * $arg2;
+        $i        = (int) ($ret_val/$arg2);
+        $ret_val -= $i*$arg2;
 
         if ($ret_val < 0) {
             $ret_val += $arg2;
@@ -179,18 +179,18 @@ 
     public static function Convert_Sat_State(Predict_Vector $pos, Predict_Vector $vel)
     {
         self::Scale_Vector(Predict::xkmper, $pos);
-        self::Scale_Vector(Predict::xkmper * Predict::xmnpda / Predict::secday, $vel);
+        self::Scale_Vector(Predict::xkmper*Predict::xmnpda/Predict::secday, $vel);
     }
 
     /* Returns angle in radians from arg in degrees */
     public static function Radians($arg)
     {
-        return $arg * Predict::de2ra;
+        return $arg*Predict::de2ra;
     }
 
     /* Returns angle in degrees from arg in rads */
     public static function Degrees($arg)
     {
-      return $arg / Predict::de2ra;
+      return $arg/Predict::de2ra;
     }
 }

require/libs/Predict/Predict/Sat.php

Doc Comments +7 added lines, -7 removed lines

@@ -124,8 +124,8 @@ 
     }
 
     /** Initialise satellite data.
-     *  @param sat The satellite to initialise.
-     *  @param qth Optional QTH info, use (0,0) if NULL.
+     *  @param sat Predict_Sat satellite to initialise.
+     *  @param qth Predict_QTH QTH info, use (0,0) if NULL.
      *
      * This function calculates the satellite data at t = 0, ie. epoch time
      * The function is called automatically by gtk_sat_data_read_sat.
@@ -216,8 +216,8 @@ 
 
     /** Determinte whether satellite is in geostationary orbit.
      *  @author John A. Magliacane, KD2BD
-     *  @param sat Pointer to satellite data.
-     *  @return TRUE if the satellite appears to be in geostationary orbit,
+     *  @param sat Predict_Sat to satellite data.
+     *  @return boolean if the satellite appears to be in geostationary orbit,
      *          FALSE otherwise.
      *
      * A satellite is in geostationary orbit if
@@ -239,8 +239,8 @@ 
     /** Determine whether satellite has decayed.
      *  @author John A. Magliacane, KD2BD
      *  @author Alexandru Csete, OZ9AEC
-     *  @param sat Pointer to satellite data.
-     *  @return TRUE if the satellite appears to have decayed, FALSE otherwise.
+     *  @param sat Predict_Sat to satellite data.
+     *  @return boolean if the satellite appears to have decayed, FALSE otherwise.
      *  @bug Modified version of the predict code but it is not tested.
      *
      * A satellite is decayed if
@@ -273,7 +273,7 @@ 
      * @param float       $time The daynum the satellite is calculated for
      * @param Predict_QTH $qth  The observer location
      *
-     * @return null on failure, float otherwise
+     * @return null|double on failure, float otherwise
      */
     public function calculateApparentMagnitude($time, Predict_QTH $qth)
     {

Indentation +296 added lines, -296 removed lines

@@ -22,304 +22,304 @@
  */
 class Predict_Sat
 {
-    // Fifth root of a hundred, used for magnitude calculation
-    const POGSONS_RATIO = 2.5118864315096;
-
-    public $name     = null;
-    public $nickname = null;
-    public $website  = null;
-
-    public $tle      = null;   /*!< Keplerian elements */
-    public $flags    = 0;      /*!< Flags for algo ctrl */
-    public $sgps     = null;
-    public $dps      = null;
-    public $deep_arg = null;
-    public $pos      = null;   /*!< Raw position and range */
-    public $vel      = null;   /*!< Raw velocity */
-
-    /*** FIXME: REMOVE */
-    public $bearing = null;   /*!< Az, El, range and vel */
-    public $astro   = null;   /*!< Ra and Decl */
-    /*** END */
-
-    /* time keeping fields */
-    public $jul_epoch = null;
-    public $jul_utc   = null;
-    public $tsince    = null;
-    public $aos       = null;    /*!< Next AOS. */
-    public $los       = null;    /*!< Next LOS */
-
-    public $az         = null;   /*!< Azimuth [deg] */
-    public $el         = null;   /*!< Elevation [deg] */
-    public $range      = null;   /*!< Range [km] */
-    public $range_rate = null;   /*!< Range Rate [km/sec] */
-    public $ra         = null;   /*!< Right Ascension [deg] */
-    public $dec        = null;   /*!< Declination [deg] */
-    public $ssplat     = null;   /*!< SSP latitude [deg] */
-    public $ssplon     = null;   /*!< SSP longitude [deg] */
-    public $alt        = null;   /*!< altitude [km] */
-    public $velo       = null;   /*!< velocity [km/s] */
-    public $ma         = null;   /*!< mean anomaly */
-    public $footprint  = null;   /*!< footprint */
-    public $phase      = null;   /*!< orbit phase */
-    public $meanmo     = null;   /*!< mean motion kept in rev/day */
-    public $orbit      = null;   /*!< orbit number */
-    public $otype      = null;   /*!< orbit type. */
-
-    public function __construct(Predict_TLE $tle)
-    {
-        $headerParts    = explode(' ', $tle->header);
-        $this->name     = $headerParts[0];
-        $this->nickname = $this->name;
-        $this->tle      = $tle;
-        $this->pos      = new Predict_Vector();
-        $this->vel      = new Predict_Vector();
-        $this->sgps     = new Predict_SGSDPStatic();
-        $this->deep_arg = new Predict_DeepArg();
-        $this->dps      = new Predict_DeepStatic();
-
-        $this->select_ephemeris();
-        $this->sat_data_init_sat($this);
-    }
-
-    /* Selects the apropriate ephemeris type to be used */
-    /* for predictions according to the data in the TLE */
-    /* It also processes values in the tle set so that  */
-    /* they are apropriate for the sgp4/sdp4 routines   */
-    public function select_ephemeris()
-    {
-        /* Preprocess tle set */
-        $this->tle->xnodeo *= Predict::de2ra;
-        $this->tle->omegao *= Predict::de2ra;
-        $this->tle->xmo    *= Predict::de2ra;
-        $this->tle->xincl  *= Predict::de2ra;
-        $temp = Predict::twopi / Predict::xmnpda / Predict::xmnpda;
-
-        /* store mean motion before conversion */
-        $this->meanmo       = $this->tle->xno;
-        $this->tle->xno     = $this->tle->xno * $temp * Predict::xmnpda;
-        $this->tle->xndt2o *= $temp;
-        $this->tle->xndd6o  = $this->tle->xndd6o * $temp / Predict::xmnpda;
-        $this->tle->bstar  /= Predict::ae;
-
-        /* Period > 225 minutes is deep space */
-        $dd1 = Predict::xke / $this->tle->xno;
-        $dd2 = Predict::tothrd;
-        $a1 = pow($dd1, $dd2);
-        $r1 = cos($this->tle->xincl);
-        $dd1 = 1.0 - $this->tle->eo * $this->tle->eo;
-        $temp = Predict::ck2 * 1.5 * ($r1 * $r1 * 3.0 - 1.0) / pow($dd1, 1.5);
-        $del1 = $temp / ($a1 * $a1);
-        $ao = $a1 * (1.0 - $del1 * (Predict::tothrd * 0.5 + $del1 *
-                                 ($del1 * 1.654320987654321 + 1.0)));
-        $delo = $temp / ($ao * $ao);
-        $xnodp = $this->tle->xno / ($delo + 1.0);
-
-        /* Select a deep-space/near-earth ephemeris */
-        if (Predict::twopi / $xnodp / Predict::xmnpda >= .15625) {
-            $this->flags |= Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
-        } else {
-            $this->flags &= ~Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
-        }
-    }
-
-    /** Initialise satellite data.
-     *  @param sat The satellite to initialise.
-     *  @param qth Optional QTH info, use (0,0) if NULL.
-     *
-     * This function calculates the satellite data at t = 0, ie. epoch time
-     * The function is called automatically by gtk_sat_data_read_sat.
-     */
-    public function sat_data_init_sat(Predict_Sat $sat, Predict_QTH $qth = null)
-    {
-        $obs_geodetic = new Predict_Geodetic();
-        $obs_set = new Predict_ObsSet();
-        $sat_geodetic = new Predict_Geodetic();
-        /* double jul_utc, age; */
-
-        $jul_utc = Predict_Time::Julian_Date_of_Epoch($sat->tle->epoch); // => tsince = 0.0
-        $sat->jul_epoch = $jul_utc;
-
-        /* initialise observer location */
-        if ($qth != null) {
-            $obs_geodetic->lon = $qth->lon * Predict::de2ra;
-            $obs_geodetic->lat = $qth->lat * Predict::de2ra;
-            $obs_geodetic->alt = $qth->alt / 1000.0;
-            $obs_geodetic->theta = 0;
-        }
-        else {
-            $obs_geodetic->lon = 0.0;
-            $obs_geodetic->lat = 0.0;
-            $obs_geodetic->alt = 0.0;
-            $obs_geodetic->theta = 0;
-        }
-
-        /* execute computations */
-        $sdpsgp = Predict_SGPSDP::getInstance($sat);
-        if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
-            $sdpsgp->SDP4($sat, 0.0);
-        } else {
-            $sdpsgp->SGP4($sat, 0.0);
-        }
-
-        /* scale position and velocity to km and km/sec */
-        Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
-
-        /* get the velocity of the satellite */
-        $sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
-        $sat->velo = $sat->vel->w;
-        Predict_SGPObs::Calculate_Obs($jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
-        Predict_SGPObs::Calculate_LatLonAlt($jul_utc, $sat->pos, $sat_geodetic);
-
-        while ($sat_geodetic->lon < -Predict::pi) {
-            $sat_geodetic->lon += Predict::twopi;
-        }
-
-        while ($sat_geodetic->lon > Predict::pi) {
-            $sat_geodetic->lon -= Predict::twopi;
-        }
-
-        $sat->az = Predict_Math::Degrees($obs_set->az);
-        $sat->el = Predict_Math::Degrees($obs_set->el);
-        $sat->range = $obs_set->range;
-        $sat->range_rate = $obs_set->range_rate;
-        $sat->ssplat = Predict_Math::Degrees($sat_geodetic->lat);
-        $sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
-        $sat->alt = $sat_geodetic->alt;
-        $sat->ma = Predict_Math::Degrees($sat->phase);
-        $sat->ma *= 256.0 / 360.0;
-        $sat->footprint = 2.0 * Predict::xkmper * acos (Predict::xkmper/$sat->pos->w);
-        $age = 0.0;
-        $sat->orbit = floor(($sat->tle->xno * Predict::xmnpda / Predict::twopi +
-                                   $age * $sat->tle->bstar * Predict::ae) * $age +
-                                  $sat->tle->xmo / Predict::twopi) + $sat->tle->revnum - 1;
-
-        /* orbit type */
-        $sat->otype = $sat->get_orbit_type($sat);
-    }
-
-    public function get_orbit_type(Predict_Sat $sat)
-    {
-         $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
-
-         if ($this->geostationary($sat)) {
-              $orbit = Predict_SGPSDP::ORBIT_TYPE_GEO;
-         } else if ($this->decayed($sat)) {
-              $orbit = Predict_SGPSDP::ORBIT_TYPE_DECAYED;
-         } else {
-              $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
-         }
-
-         return $orbit;
-    }
-
-
-    /** Determinte whether satellite is in geostationary orbit.
-     *  @author John A. Magliacane, KD2BD
-     *  @param sat Pointer to satellite data.
-     *  @return TRUE if the satellite appears to be in geostationary orbit,
-     *          FALSE otherwise.
-     *
-     * A satellite is in geostationary orbit if
-     *
-     *     fabs (sat.meanmotion - 1.0027) < 0.0002
-     *
-     * Note: Appearantly, the mean motion can deviate much more from 1.0027 than 0.0002
-     */
-    public function geostationary(Predict_Sat $sat)
-    {
-         if (abs($sat->meanmo - 1.0027) < 0.0002) {
-              return true;
-         } else {
-              return false;
-        }
-    }
-
-
-    /** Determine whether satellite has decayed.
-     *  @author John A. Magliacane, KD2BD
-     *  @author Alexandru Csete, OZ9AEC
-     *  @param sat Pointer to satellite data.
-     *  @return TRUE if the satellite appears to have decayed, FALSE otherwise.
-     *  @bug Modified version of the predict code but it is not tested.
-     *
-     * A satellite is decayed if
-     *
-     *    satepoch + ((16.666666 - sat.meanmo) / (10.0*fabs(sat.drag))) < "now"
-     *
-     */
-    public function decayed(Predict_Sat $sat)
-    {
-        /* tle.xndt2o/(twopi/xmnpda/xmnpda) is the value before converted the
+	// Fifth root of a hundred, used for magnitude calculation
+	const POGSONS_RATIO = 2.5118864315096;
+
+	public $name     = null;
+	public $nickname = null;
+	public $website  = null;
+
+	public $tle      = null;   /*!< Keplerian elements */
+	public $flags    = 0;      /*!< Flags for algo ctrl */
+	public $sgps     = null;
+	public $dps      = null;
+	public $deep_arg = null;
+	public $pos      = null;   /*!< Raw position and range */
+	public $vel      = null;   /*!< Raw velocity */
+
+	/*** FIXME: REMOVE */
+	public $bearing = null;   /*!< Az, El, range and vel */
+	public $astro   = null;   /*!< Ra and Decl */
+	/*** END */
+
+	/* time keeping fields */
+	public $jul_epoch = null;
+	public $jul_utc   = null;
+	public $tsince    = null;
+	public $aos       = null;    /*!< Next AOS. */
+	public $los       = null;    /*!< Next LOS */
+
+	public $az         = null;   /*!< Azimuth [deg] */
+	public $el         = null;   /*!< Elevation [deg] */
+	public $range      = null;   /*!< Range [km] */
+	public $range_rate = null;   /*!< Range Rate [km/sec] */
+	public $ra         = null;   /*!< Right Ascension [deg] */
+	public $dec        = null;   /*!< Declination [deg] */
+	public $ssplat     = null;   /*!< SSP latitude [deg] */
+	public $ssplon     = null;   /*!< SSP longitude [deg] */
+	public $alt        = null;   /*!< altitude [km] */
+	public $velo       = null;   /*!< velocity [km/s] */
+	public $ma         = null;   /*!< mean anomaly */
+	public $footprint  = null;   /*!< footprint */
+	public $phase      = null;   /*!< orbit phase */
+	public $meanmo     = null;   /*!< mean motion kept in rev/day */
+	public $orbit      = null;   /*!< orbit number */
+	public $otype      = null;   /*!< orbit type. */
+
+	public function __construct(Predict_TLE $tle)
+	{
+		$headerParts    = explode(' ', $tle->header);
+		$this->name     = $headerParts[0];
+		$this->nickname = $this->name;
+		$this->tle      = $tle;
+		$this->pos      = new Predict_Vector();
+		$this->vel      = new Predict_Vector();
+		$this->sgps     = new Predict_SGSDPStatic();
+		$this->deep_arg = new Predict_DeepArg();
+		$this->dps      = new Predict_DeepStatic();
+
+		$this->select_ephemeris();
+		$this->sat_data_init_sat($this);
+	}
+
+	/* Selects the apropriate ephemeris type to be used */
+	/* for predictions according to the data in the TLE */
+	/* It also processes values in the tle set so that  */
+	/* they are apropriate for the sgp4/sdp4 routines   */
+	public function select_ephemeris()
+	{
+		/* Preprocess tle set */
+		$this->tle->xnodeo *= Predict::de2ra;
+		$this->tle->omegao *= Predict::de2ra;
+		$this->tle->xmo    *= Predict::de2ra;
+		$this->tle->xincl  *= Predict::de2ra;
+		$temp = Predict::twopi / Predict::xmnpda / Predict::xmnpda;
+
+		/* store mean motion before conversion */
+		$this->meanmo       = $this->tle->xno;
+		$this->tle->xno     = $this->tle->xno * $temp * Predict::xmnpda;
+		$this->tle->xndt2o *= $temp;
+		$this->tle->xndd6o  = $this->tle->xndd6o * $temp / Predict::xmnpda;
+		$this->tle->bstar  /= Predict::ae;
+
+		/* Period > 225 minutes is deep space */
+		$dd1 = Predict::xke / $this->tle->xno;
+		$dd2 = Predict::tothrd;
+		$a1 = pow($dd1, $dd2);
+		$r1 = cos($this->tle->xincl);
+		$dd1 = 1.0 - $this->tle->eo * $this->tle->eo;
+		$temp = Predict::ck2 * 1.5 * ($r1 * $r1 * 3.0 - 1.0) / pow($dd1, 1.5);
+		$del1 = $temp / ($a1 * $a1);
+		$ao = $a1 * (1.0 - $del1 * (Predict::tothrd * 0.5 + $del1 *
+								 ($del1 * 1.654320987654321 + 1.0)));
+		$delo = $temp / ($ao * $ao);
+		$xnodp = $this->tle->xno / ($delo + 1.0);
+
+		/* Select a deep-space/near-earth ephemeris */
+		if (Predict::twopi / $xnodp / Predict::xmnpda >= .15625) {
+			$this->flags |= Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
+		} else {
+			$this->flags &= ~Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
+		}
+	}
+
+	/** Initialise satellite data.
+	 *  @param sat The satellite to initialise.
+	 *  @param qth Optional QTH info, use (0,0) if NULL.
+	 *
+	 * This function calculates the satellite data at t = 0, ie. epoch time
+	 * The function is called automatically by gtk_sat_data_read_sat.
+	 */
+	public function sat_data_init_sat(Predict_Sat $sat, Predict_QTH $qth = null)
+	{
+		$obs_geodetic = new Predict_Geodetic();
+		$obs_set = new Predict_ObsSet();
+		$sat_geodetic = new Predict_Geodetic();
+		/* double jul_utc, age; */
+
+		$jul_utc = Predict_Time::Julian_Date_of_Epoch($sat->tle->epoch); // => tsince = 0.0
+		$sat->jul_epoch = $jul_utc;
+
+		/* initialise observer location */
+		if ($qth != null) {
+			$obs_geodetic->lon = $qth->lon * Predict::de2ra;
+			$obs_geodetic->lat = $qth->lat * Predict::de2ra;
+			$obs_geodetic->alt = $qth->alt / 1000.0;
+			$obs_geodetic->theta = 0;
+		}
+		else {
+			$obs_geodetic->lon = 0.0;
+			$obs_geodetic->lat = 0.0;
+			$obs_geodetic->alt = 0.0;
+			$obs_geodetic->theta = 0;
+		}
+
+		/* execute computations */
+		$sdpsgp = Predict_SGPSDP::getInstance($sat);
+		if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
+			$sdpsgp->SDP4($sat, 0.0);
+		} else {
+			$sdpsgp->SGP4($sat, 0.0);
+		}
+
+		/* scale position and velocity to km and km/sec */
+		Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
+
+		/* get the velocity of the satellite */
+		$sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
+		$sat->velo = $sat->vel->w;
+		Predict_SGPObs::Calculate_Obs($jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
+		Predict_SGPObs::Calculate_LatLonAlt($jul_utc, $sat->pos, $sat_geodetic);
+
+		while ($sat_geodetic->lon < -Predict::pi) {
+			$sat_geodetic->lon += Predict::twopi;
+		}
+
+		while ($sat_geodetic->lon > Predict::pi) {
+			$sat_geodetic->lon -= Predict::twopi;
+		}
+
+		$sat->az = Predict_Math::Degrees($obs_set->az);
+		$sat->el = Predict_Math::Degrees($obs_set->el);
+		$sat->range = $obs_set->range;
+		$sat->range_rate = $obs_set->range_rate;
+		$sat->ssplat = Predict_Math::Degrees($sat_geodetic->lat);
+		$sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
+		$sat->alt = $sat_geodetic->alt;
+		$sat->ma = Predict_Math::Degrees($sat->phase);
+		$sat->ma *= 256.0 / 360.0;
+		$sat->footprint = 2.0 * Predict::xkmper * acos (Predict::xkmper/$sat->pos->w);
+		$age = 0.0;
+		$sat->orbit = floor(($sat->tle->xno * Predict::xmnpda / Predict::twopi +
+								   $age * $sat->tle->bstar * Predict::ae) * $age +
+								  $sat->tle->xmo / Predict::twopi) + $sat->tle->revnum - 1;
+
+		/* orbit type */
+		$sat->otype = $sat->get_orbit_type($sat);
+	}
+
+	public function get_orbit_type(Predict_Sat $sat)
+	{
+		 $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
+
+		 if ($this->geostationary($sat)) {
+			  $orbit = Predict_SGPSDP::ORBIT_TYPE_GEO;
+		 } else if ($this->decayed($sat)) {
+			  $orbit = Predict_SGPSDP::ORBIT_TYPE_DECAYED;
+		 } else {
+			  $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
+		 }
+
+		 return $orbit;
+	}
+
+
+	/** Determinte whether satellite is in geostationary orbit.
+	 *  @author John A. Magliacane, KD2BD
+	 *  @param sat Pointer to satellite data.
+	 *  @return TRUE if the satellite appears to be in geostationary orbit,
+	 *          FALSE otherwise.
+	 *
+	 * A satellite is in geostationary orbit if
+	 *
+	 *     fabs (sat.meanmotion - 1.0027) < 0.0002
+	 *
+	 * Note: Appearantly, the mean motion can deviate much more from 1.0027 than 0.0002
+	 */
+	public function geostationary(Predict_Sat $sat)
+	{
+		 if (abs($sat->meanmo - 1.0027) < 0.0002) {
+			  return true;
+		 } else {
+			  return false;
+		}
+	}
+
+
+	/** Determine whether satellite has decayed.
+	 *  @author John A. Magliacane, KD2BD
+	 *  @author Alexandru Csete, OZ9AEC
+	 *  @param sat Pointer to satellite data.
+	 *  @return TRUE if the satellite appears to have decayed, FALSE otherwise.
+	 *  @bug Modified version of the predict code but it is not tested.
+	 *
+	 * A satellite is decayed if
+	 *
+	 *    satepoch + ((16.666666 - sat.meanmo) / (10.0*fabs(sat.drag))) < "now"
+	 *
+	 */
+	public function decayed(Predict_Sat $sat)
+	{
+		/* tle.xndt2o/(twopi/xmnpda/xmnpda) is the value before converted the
            value matches up with the value in predict 2.2.3 */
-        /*** FIXME decayed is treated as a static quantity.
+		/*** FIXME decayed is treated as a static quantity.
              It is time dependent. Also sat->jul_utc is often zero
              when this function is called
         ***/
-        if ((10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda))) == 0) {
-    		return true;
-    	} elseif ($sat->jul_epoch + ((16.666666 - $sat->meanmo) /
-                               (10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda)))) < $sat->jul_utc) {
-              return true;
-        } else {
-              return false;
-        }
-    }
-
-    /**
-     * Experimental attempt at calculating apparent magnitude.  Known intrinsic
-     * magnitudes are listed inside the function for now.
-     *
-     * @param float       $time The daynum the satellite is calculated for
-     * @param Predict_QTH $qth  The observer location
-     *
-     * @return null on failure, float otherwise
-     */
-    public function calculateApparentMagnitude($time, Predict_QTH $qth)
-    {
-        // Recorded intrinsic magnitudes and their respective
-        // illumination and distance from heavens-above.com
-        static $intrinsicMagnitudes = array(
-            '25544' => array(
-                'mag'      => -1.3,
-                'illum'    => .5,
-                'distance' => 1000,
-            )
-        );
-
-        // Return null if we don't have a record of the intrinsic mag
-        if (!isset($intrinsicMagnitudes[$this->tle->catnr])) {
-            return null;
-        }
-        $imag = $intrinsicMagnitudes[$this->tle->catnr];
-
-        // Convert the observer's geodetic info to radians and km so
-        // we can compare vectors
-        $observerGeo      = new Predict_Geodetic();
-        $observerGeo->lat = Predict_Math::Radians($qth->lat);
-        $observerGeo->lon = Predict_Math::Radians($qth->lon);
-        $observerGeo->alt = $qth->alt * 1000;
-
-        // Now determine the sun and observer positions
-        $observerPos      = new Predict_Vector();
-        $observerVel      = new Predict_Vector();
-        $solarVector      = new Predict_Vector();
-        Predict_Solar::Calculate_Solar_Position($time, $solarVector);
-        Predict_SGPObs::Calculate_User_PosVel($time, $observerGeo, $observerPos, $observerVel);
-
-        // Determine the solar phase and and thus the percent illumination
-        $observerSatPos = new Predict_Vector();
-        Predict_Math::Vec_Sub($this->pos, $observerPos, $observerSatPos);
-        $phaseAngle = Predict_Math::Degrees(Predict_Math::Angle($solarVector, $observerSatPos));
-        $illum      = $phaseAngle / 180;
-
-        $illuminationChange            = $illum / $imag['illum'];
-        $inverseSquareOfDistanceChange = pow(($imag['distance'] / $this->range), 2);
-        $changeInMagnitude             = log(
-            $illuminationChange * $inverseSquareOfDistanceChange,
-            self::POGSONS_RATIO
-        );
-
-        return $imag['mag'] - $changeInMagnitude;
-    }
+		if ((10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda))) == 0) {
+			return true;
+		} elseif ($sat->jul_epoch + ((16.666666 - $sat->meanmo) /
+							   (10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda)))) < $sat->jul_utc) {
+			  return true;
+		} else {
+			  return false;
+		}
+	}
+
+	/**
+	 * Experimental attempt at calculating apparent magnitude.  Known intrinsic
+	 * magnitudes are listed inside the function for now.
+	 *
+	 * @param float       $time The daynum the satellite is calculated for
+	 * @param Predict_QTH $qth  The observer location
+	 *
+	 * @return null on failure, float otherwise
+	 */
+	public function calculateApparentMagnitude($time, Predict_QTH $qth)
+	{
+		// Recorded intrinsic magnitudes and their respective
+		// illumination and distance from heavens-above.com
+		static $intrinsicMagnitudes = array(
+			'25544' => array(
+				'mag'      => -1.3,
+				'illum'    => .5,
+				'distance' => 1000,
+			)
+		);
+
+		// Return null if we don't have a record of the intrinsic mag
+		if (!isset($intrinsicMagnitudes[$this->tle->catnr])) {
+			return null;
+		}
+		$imag = $intrinsicMagnitudes[$this->tle->catnr];
+
+		// Convert the observer's geodetic info to radians and km so
+		// we can compare vectors
+		$observerGeo      = new Predict_Geodetic();
+		$observerGeo->lat = Predict_Math::Radians($qth->lat);
+		$observerGeo->lon = Predict_Math::Radians($qth->lon);
+		$observerGeo->alt = $qth->alt * 1000;
+
+		// Now determine the sun and observer positions
+		$observerPos      = new Predict_Vector();
+		$observerVel      = new Predict_Vector();
+		$solarVector      = new Predict_Vector();
+		Predict_Solar::Calculate_Solar_Position($time, $solarVector);
+		Predict_SGPObs::Calculate_User_PosVel($time, $observerGeo, $observerPos, $observerVel);
+
+		// Determine the solar phase and and thus the percent illumination
+		$observerSatPos = new Predict_Vector();
+		Predict_Math::Vec_Sub($this->pos, $observerPos, $observerSatPos);
+		$phaseAngle = Predict_Math::Degrees(Predict_Math::Angle($solarVector, $observerSatPos));
+		$illum      = $phaseAngle / 180;
+
+		$illuminationChange            = $illum / $imag['illum'];
+		$inverseSquareOfDistanceChange = pow(($imag['distance'] / $this->range), 2);
+		$changeInMagnitude             = log(
+			$illuminationChange * $inverseSquareOfDistanceChange,
+			self::POGSONS_RATIO
+		);
+
+		return $imag['mag'] - $changeInMagnitude;
+	}
 }

Spacing +55 added lines, -55 removed lines

@@ -29,42 +29,42 @@ 
     public $nickname = null;
     public $website  = null;
 
-    public $tle      = null;   /*!< Keplerian elements */
-    public $flags    = 0;      /*!< Flags for algo ctrl */
+    public $tle      = null; /*!< Keplerian elements */
+    public $flags    = 0; /*!< Flags for algo ctrl */
     public $sgps     = null;
     public $dps      = null;
     public $deep_arg = null;
-    public $pos      = null;   /*!< Raw position and range */
-    public $vel      = null;   /*!< Raw velocity */
+    public $pos      = null; /*!< Raw position and range */
+    public $vel      = null; /*!< Raw velocity */
 
     /*** FIXME: REMOVE */
-    public $bearing = null;   /*!< Az, El, range and vel */
-    public $astro   = null;   /*!< Ra and Decl */
+    public $bearing = null; /*!< Az, El, range and vel */
+    public $astro   = null; /*!< Ra and Decl */
     /*** END */
 
     /* time keeping fields */
     public $jul_epoch = null;
     public $jul_utc   = null;
     public $tsince    = null;
-    public $aos       = null;    /*!< Next AOS. */
-    public $los       = null;    /*!< Next LOS */
-
-    public $az         = null;   /*!< Azimuth [deg] */
-    public $el         = null;   /*!< Elevation [deg] */
-    public $range      = null;   /*!< Range [km] */
-    public $range_rate = null;   /*!< Range Rate [km/sec] */
-    public $ra         = null;   /*!< Right Ascension [deg] */
-    public $dec        = null;   /*!< Declination [deg] */
-    public $ssplat     = null;   /*!< SSP latitude [deg] */
-    public $ssplon     = null;   /*!< SSP longitude [deg] */
-    public $alt        = null;   /*!< altitude [km] */
-    public $velo       = null;   /*!< velocity [km/s] */
-    public $ma         = null;   /*!< mean anomaly */
-    public $footprint  = null;   /*!< footprint */
-    public $phase      = null;   /*!< orbit phase */
-    public $meanmo     = null;   /*!< mean motion kept in rev/day */
-    public $orbit      = null;   /*!< orbit number */
-    public $otype      = null;   /*!< orbit type. */
+    public $aos       = null; /*!< Next AOS. */
+    public $los       = null; /*!< Next LOS */
+
+    public $az         = null; /*!< Azimuth [deg] */
+    public $el         = null; /*!< Elevation [deg] */
+    public $range      = null; /*!< Range [km] */
+    public $range_rate = null; /*!< Range Rate [km/sec] */
+    public $ra         = null; /*!< Right Ascension [deg] */
+    public $dec        = null; /*!< Declination [deg] */
+    public $ssplat     = null; /*!< SSP latitude [deg] */
+    public $ssplon     = null; /*!< SSP longitude [deg] */
+    public $alt        = null; /*!< altitude [km] */
+    public $velo       = null; /*!< velocity [km/s] */
+    public $ma         = null; /*!< mean anomaly */
+    public $footprint  = null; /*!< footprint */
+    public $phase      = null; /*!< orbit phase */
+    public $meanmo     = null; /*!< mean motion kept in rev/day */
+    public $orbit      = null; /*!< orbit number */
+    public $otype      = null; /*!< orbit type. */
 
     public function __construct(Predict_TLE $tle)
     {
@@ -93,30 +93,30 @@ 
         $this->tle->omegao *= Predict::de2ra;
         $this->tle->xmo    *= Predict::de2ra;
         $this->tle->xincl  *= Predict::de2ra;
-        $temp = Predict::twopi / Predict::xmnpda / Predict::xmnpda;
+        $temp = Predict::twopi/Predict::xmnpda/Predict::xmnpda;
 
         /* store mean motion before conversion */
         $this->meanmo       = $this->tle->xno;
-        $this->tle->xno     = $this->tle->xno * $temp * Predict::xmnpda;
+        $this->tle->xno     = $this->tle->xno*$temp*Predict::xmnpda;
         $this->tle->xndt2o *= $temp;
-        $this->tle->xndd6o  = $this->tle->xndd6o * $temp / Predict::xmnpda;
+        $this->tle->xndd6o  = $this->tle->xndd6o*$temp/Predict::xmnpda;
         $this->tle->bstar  /= Predict::ae;
 
         /* Period > 225 minutes is deep space */
-        $dd1 = Predict::xke / $this->tle->xno;
+        $dd1 = Predict::xke/$this->tle->xno;
         $dd2 = Predict::tothrd;
         $a1 = pow($dd1, $dd2);
         $r1 = cos($this->tle->xincl);
-        $dd1 = 1.0 - $this->tle->eo * $this->tle->eo;
-        $temp = Predict::ck2 * 1.5 * ($r1 * $r1 * 3.0 - 1.0) / pow($dd1, 1.5);
-        $del1 = $temp / ($a1 * $a1);
-        $ao = $a1 * (1.0 - $del1 * (Predict::tothrd * 0.5 + $del1 *
-                                 ($del1 * 1.654320987654321 + 1.0)));
-        $delo = $temp / ($ao * $ao);
-        $xnodp = $this->tle->xno / ($delo + 1.0);
+        $dd1 = 1.0 - $this->tle->eo*$this->tle->eo;
+        $temp = Predict::ck2*1.5*($r1*$r1*3.0 - 1.0)/pow($dd1, 1.5);
+        $del1 = $temp/($a1*$a1);
+        $ao = $a1*(1.0 - $del1*(Predict::tothrd*0.5 + $del1*
+                                 ($del1*1.654320987654321 + 1.0)));
+        $delo = $temp/($ao*$ao);
+        $xnodp = $this->tle->xno/($delo + 1.0);
 
         /* Select a deep-space/near-earth ephemeris */
-        if (Predict::twopi / $xnodp / Predict::xmnpda >= .15625) {
+        if (Predict::twopi/$xnodp/Predict::xmnpda >= .15625) {
             $this->flags |= Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
         } else {
             $this->flags &= ~Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
@@ -142,9 +142,9 @@ 
 
         /* initialise observer location */
         if ($qth != null) {
-            $obs_geodetic->lon = $qth->lon * Predict::de2ra;
-            $obs_geodetic->lat = $qth->lat * Predict::de2ra;
-            $obs_geodetic->alt = $qth->alt / 1000.0;
+            $obs_geodetic->lon = $qth->lon*Predict::de2ra;
+            $obs_geodetic->lat = $qth->lat*Predict::de2ra;
+            $obs_geodetic->alt = $qth->alt/1000.0;
             $obs_geodetic->theta = 0;
         }
         else {
@@ -156,7 +156,7 @@ 
 
         /* execute computations */
         $sdpsgp = Predict_SGPSDP::getInstance($sat);
-        if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
+        if ($sat->flags&Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
             $sdpsgp->SDP4($sat, 0.0);
         } else {
             $sdpsgp->SGP4($sat, 0.0);
@@ -166,7 +166,7 @@ 
         Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
 
         /* get the velocity of the satellite */
-        $sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
+        $sat->vel->w = sqrt($sat->vel->x*$sat->vel->x + $sat->vel->y*$sat->vel->y + $sat->vel->z*$sat->vel->z);
         $sat->velo = $sat->vel->w;
         Predict_SGPObs::Calculate_Obs($jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
         Predict_SGPObs::Calculate_LatLonAlt($jul_utc, $sat->pos, $sat_geodetic);
@@ -187,12 +187,12 @@ 
         $sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
         $sat->alt = $sat_geodetic->alt;
         $sat->ma = Predict_Math::Degrees($sat->phase);
-        $sat->ma *= 256.0 / 360.0;
-        $sat->footprint = 2.0 * Predict::xkmper * acos (Predict::xkmper/$sat->pos->w);
+        $sat->ma *= 256.0/360.0;
+        $sat->footprint = 2.0*Predict::xkmper*acos(Predict::xkmper/$sat->pos->w);
         $age = 0.0;
-        $sat->orbit = floor(($sat->tle->xno * Predict::xmnpda / Predict::twopi +
-                                   $age * $sat->tle->bstar * Predict::ae) * $age +
-                                  $sat->tle->xmo / Predict::twopi) + $sat->tle->revnum - 1;
+        $sat->orbit = floor(($sat->tle->xno*Predict::xmnpda/Predict::twopi +
+                                   $age*$sat->tle->bstar*Predict::ae)*$age +
+                                  $sat->tle->xmo/Predict::twopi) + $sat->tle->revnum - 1;
 
         /* orbit type */
         $sat->otype = $sat->get_orbit_type($sat);
@@ -256,10 +256,10 @@ 
              It is time dependent. Also sat->jul_utc is often zero
              when this function is called
         ***/
-        if ((10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda))) == 0) {
+        if ((10.0*abs($sat->tle->xndt2o/(Predict::twopi/Predict::xmnpda/Predict::xmnpda))) == 0) {
     		return true;
-    	} elseif ($sat->jul_epoch + ((16.666666 - $sat->meanmo) /
-                               (10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda)))) < $sat->jul_utc) {
+    	} elseif ($sat->jul_epoch + ((16.666666 - $sat->meanmo)/
+                               (10.0*abs($sat->tle->xndt2o/(Predict::twopi/Predict::xmnpda/Predict::xmnpda)))) < $sat->jul_utc) {
               return true;
         } else {
               return false;
@@ -298,7 +298,7 @@ 
         $observerGeo      = new Predict_Geodetic();
         $observerGeo->lat = Predict_Math::Radians($qth->lat);
         $observerGeo->lon = Predict_Math::Radians($qth->lon);
-        $observerGeo->alt = $qth->alt * 1000;
+        $observerGeo->alt = $qth->alt*1000;
 
         // Now determine the sun and observer positions
         $observerPos      = new Predict_Vector();
@@ -311,12 +311,12 @@ 
         $observerSatPos = new Predict_Vector();
         Predict_Math::Vec_Sub($this->pos, $observerPos, $observerSatPos);
         $phaseAngle = Predict_Math::Degrees(Predict_Math::Angle($solarVector, $observerSatPos));
-        $illum      = $phaseAngle / 180;
+        $illum      = $phaseAngle/180;
 
-        $illuminationChange            = $illum / $imag['illum'];
-        $inverseSquareOfDistanceChange = pow(($imag['distance'] / $this->range), 2);
+        $illuminationChange            = $illum/$imag['illum'];
+        $inverseSquareOfDistanceChange = pow(($imag['distance']/$this->range), 2);
         $changeInMagnitude             = log(
-            $illuminationChange * $inverseSquareOfDistanceChange,
+            $illuminationChange*$inverseSquareOfDistanceChange,
             self::POGSONS_RATIO
         );
 

Braces +1 added lines, -2 removed lines

@@ -146,8 +146,7 @@
             $obs_geodetic->lat = $qth->lat * Predict::de2ra;
             $obs_geodetic->alt = $qth->alt / 1000.0;
             $obs_geodetic->theta = 0;
-        }
-        else {
+        } else {
             $obs_geodetic->lon = 0.0;
             $obs_geodetic->lat = 0.0;
             $obs_geodetic->alt = 0.0;

require/libs/Predict/Predict/SGPObs.php

Doc Comments +4 added lines

@@ -56,6 +56,10 @@
     /* It is intended to be used to determine the ground track of */
     /* a satellite.  The calculations  assume the earth to be an  */
     /* oblate spheroid as defined in WGS '72.                     */
+
+    /**
+     * @param double $_time
+     */
     public static function Calculate_LatLonAlt($_time, Predict_Vector $pos,  Predict_Geodetic $geodetic)
     {
         /* Reference:  The 1992 Astronomical Almanac, page K12. */

Indentation +127 added lines, -127 removed lines

@@ -25,131 +25,131 @@
  */
 class Predict_SGPObs
 {
-    /* Procedure Calculate_User_PosVel passes the user's geodetic position */
-    /* and the time of interest and returns the ECI position and velocity  */
-    /* of the observer. The velocity calculation assumes the geodetic      */
-    /* position is stationary relative to the earth's surface.             */
-    public static function Calculate_User_PosVel(
-        $_time, Predict_Geodetic $geodetic, Predict_Vector $obs_pos, Predict_Vector $obs_vel
-    )
-    {
-        /* Reference:  The 1992 Astronomical Almanac, page K11. */
-
-        $sinGeodeticLat = sin($geodetic->lat); /* Only run sin($geodetic->lat) once */
-
-        $geodetic->theta = Predict_Math::FMod2p(Predict_Time::ThetaG_JD($_time) + $geodetic->lon);/*LMST*/
-        $c = 1 / sqrt(1 + Predict::__f * (Predict::__f - 2) * $sinGeodeticLat * $sinGeodeticLat);
-        $sq = (1 - Predict::__f) * (1 - Predict::__f) * $c;
-        $achcp = (Predict::xkmper * $c + $geodetic->alt) * cos($geodetic->lat);
-        $obs_pos->x = $achcp * cos($geodetic->theta); /*kilometers*/
-        $obs_pos->y = $achcp * sin($geodetic->theta);
-        $obs_pos->z = (Predict::xkmper * $sq + $geodetic->alt) * $sinGeodeticLat;
-        $obs_vel->x = -Predict::mfactor * $obs_pos->y; /*kilometers/second*/
-        $obs_vel->y =  Predict::mfactor * $obs_pos->x;
-        $obs_vel->z =  0;
-        $obs_pos->w = sqrt($obs_pos->x * $obs_pos->x + $obs_pos->y * $obs_pos->y + $obs_pos->z * $obs_pos->z);
-        $obs_vel->w = sqrt($obs_vel->x * $obs_vel->x + $obs_vel->y * $obs_vel->y + $obs_vel->z * $obs_vel->z);
-    }
-
-    /* Procedure Calculate_LatLonAlt will calculate the geodetic  */
-    /* position of an object given its ECI position pos and time. */
-    /* It is intended to be used to determine the ground track of */
-    /* a satellite.  The calculations  assume the earth to be an  */
-    /* oblate spheroid as defined in WGS '72.                     */
-    public static function Calculate_LatLonAlt($_time, Predict_Vector $pos,  Predict_Geodetic $geodetic)
-    {
-        /* Reference:  The 1992 Astronomical Almanac, page K12. */
-
-        /* double r,e2,phi,c; */
-
-        $geodetic->theta = Predict_Math::AcTan($pos->y, $pos->x); /*radians*/
-        $geodetic->lon = Predict_Math::FMod2p($geodetic->theta - Predict_Time::ThetaG_JD($_time)); /*radians*/
-        $r = sqrt(($pos->x * $pos->x) + ($pos->y * $pos->y));
-        $e2 = Predict::__f * (2 - Predict::__f);
-        $geodetic->lat = Predict_Math::AcTan($pos->z, $r); /*radians*/
-
-        do {
-            $phi    = $geodetic->lat;
-            $sinPhi = sin($phi);
-            $c      = 1 / sqrt(1 - $e2 * ($sinPhi * $sinPhi));
-            $geodetic->lat = Predict_Math::AcTan($pos->z + Predict::xkmper * $c * $e2 * $sinPhi, $r);
-        } while (abs($geodetic->lat - $phi) >= 1E-10);
-
-        $geodetic->alt = $r / cos($geodetic->lat) - Predict::xkmper * $c;/*kilometers*/
-
-        if ($geodetic->lat > Predict::pio2) {
-            $geodetic->lat -= Predict::twopi;
-        }
-    }
-
-    /* The procedures Calculate_Obs and Calculate_RADec calculate         */
-    /* the *topocentric* coordinates of the object with ECI position,     */
-    /* {pos}, and velocity, {vel}, from location {geodetic} at {time}.    */
-    /* The {obs_set} returned for Calculate_Obs consists of azimuth,      */
-    /* elevation, range, and range rate (in that order) with units of     */
-    /* radians, radians, kilometers, and kilometers/second, respectively. */
-    /* The WGS '72 geoid is used and the effect of atmospheric refraction */
-    /* (under standard temperature and pressure) is incorporated into the */
-    /* elevation calculation; the effect of atmospheric refraction on     */
-    /* range and range rate has not yet been quantified.                  */
-
-    /* The {obs_set} for Calculate_RADec consists of right ascension and  */
-    /* declination (in that order) in radians.  Again, calculations are   */
-    /* based on *topocentric* position using the WGS '72 geoid and        */
-    /* incorporating atmospheric refraction.                              */
-    public static function Calculate_Obs($_time, Predict_Vector $pos, Predict_Vector $vel, Predict_Geodetic $geodetic, Predict_ObsSet $obs_set)
-    {
-        $obs_pos = new Predict_Vector();
-        $obs_vel = new Predict_Vector();
-        $range   = new Predict_Vector();
-        $rgvel   = new Predict_Vector();
-
-        self::Calculate_User_PosVel($_time, $geodetic, $obs_pos, $obs_vel);
-
-        $range->x = $pos->x - $obs_pos->x;
-        $range->y = $pos->y - $obs_pos->y;
-        $range->z = $pos->z - $obs_pos->z;
-
-        $rgvel->x = $vel->x - $obs_vel->x;
-        $rgvel->y = $vel->y - $obs_vel->y;
-        $rgvel->z = $vel->z - $obs_vel->z;
-
-        $range->w = sqrt($range->x * $range->x + $range->y * $range->y + $range->z * $range->z);
-
-        $sin_lat   = sin($geodetic->lat);
-        $cos_lat   = cos($geodetic->lat);
-        $sin_theta = sin($geodetic->theta);
-        $cos_theta = cos($geodetic->theta);
-        $top_s = $sin_lat * $cos_theta * $range->x
-            + $sin_lat * $sin_theta * $range->y
-            - $cos_lat * $range->z;
-        $top_e = -$sin_theta * $range->x
-            + $cos_theta * $range->y;
-        $top_z = $cos_lat * $cos_theta * $range->x
-            + $cos_lat * $sin_theta * $range->y
-            + $sin_lat * $range->z;
-        $azim = atan(-$top_e / $top_s); /*Azimuth*/
-        if ($top_s > 0) {
-            $azim = $azim + Predict::pi;
-        }
-        if ($azim < 0 ) {
-            $azim = $azim + Predict::twopi;
-        }
-        $el = Predict_Math::ArcSin($top_z / $range->w);
-        $obs_set->az = $azim;        /* Azimuth (radians)  */
-        $obs_set->el = $el;          /* Elevation (radians)*/
-        $obs_set->range = $range->w; /* Range (kilometers) */
-
-        /* Range Rate (kilometers/second)*/
-        $obs_set->range_rate = Predict_Math::Dot($range, $rgvel) / $range->w;
-
-        /* Corrections for atmospheric refraction */
-        /* Reference:  Astronomical Algorithms by Jean Meeus, pp. 101-104    */
-        /* Correction is meaningless when apparent elevation is below horizon */
-        //	obs_set->el = obs_set->el + Radians((1.02/tan(Radians(Degrees(el)+
-        //							      10.3/(Degrees(el)+5.11))))/60);
-        if ($obs_set->el < 0) {
-            $obs_set->el = $el;  /*Reset to true elevation*/
-        }
-    }
+	/* Procedure Calculate_User_PosVel passes the user's geodetic position */
+	/* and the time of interest and returns the ECI position and velocity  */
+	/* of the observer. The velocity calculation assumes the geodetic      */
+	/* position is stationary relative to the earth's surface.             */
+	public static function Calculate_User_PosVel(
+		$_time, Predict_Geodetic $geodetic, Predict_Vector $obs_pos, Predict_Vector $obs_vel
+	)
+	{
+		/* Reference:  The 1992 Astronomical Almanac, page K11. */
+
+		$sinGeodeticLat = sin($geodetic->lat); /* Only run sin($geodetic->lat) once */
+
+		$geodetic->theta = Predict_Math::FMod2p(Predict_Time::ThetaG_JD($_time) + $geodetic->lon);/*LMST*/
+		$c = 1 / sqrt(1 + Predict::__f * (Predict::__f - 2) * $sinGeodeticLat * $sinGeodeticLat);
+		$sq = (1 - Predict::__f) * (1 - Predict::__f) * $c;
+		$achcp = (Predict::xkmper * $c + $geodetic->alt) * cos($geodetic->lat);
+		$obs_pos->x = $achcp * cos($geodetic->theta); /*kilometers*/
+		$obs_pos->y = $achcp * sin($geodetic->theta);
+		$obs_pos->z = (Predict::xkmper * $sq + $geodetic->alt) * $sinGeodeticLat;
+		$obs_vel->x = -Predict::mfactor * $obs_pos->y; /*kilometers/second*/
+		$obs_vel->y =  Predict::mfactor * $obs_pos->x;
+		$obs_vel->z =  0;
+		$obs_pos->w = sqrt($obs_pos->x * $obs_pos->x + $obs_pos->y * $obs_pos->y + $obs_pos->z * $obs_pos->z);
+		$obs_vel->w = sqrt($obs_vel->x * $obs_vel->x + $obs_vel->y * $obs_vel->y + $obs_vel->z * $obs_vel->z);
+	}
+
+	/* Procedure Calculate_LatLonAlt will calculate the geodetic  */
+	/* position of an object given its ECI position pos and time. */
+	/* It is intended to be used to determine the ground track of */
+	/* a satellite.  The calculations  assume the earth to be an  */
+	/* oblate spheroid as defined in WGS '72.                     */
+	public static function Calculate_LatLonAlt($_time, Predict_Vector $pos,  Predict_Geodetic $geodetic)
+	{
+		/* Reference:  The 1992 Astronomical Almanac, page K12. */
+
+		/* double r,e2,phi,c; */
+
+		$geodetic->theta = Predict_Math::AcTan($pos->y, $pos->x); /*radians*/
+		$geodetic->lon = Predict_Math::FMod2p($geodetic->theta - Predict_Time::ThetaG_JD($_time)); /*radians*/
+		$r = sqrt(($pos->x * $pos->x) + ($pos->y * $pos->y));
+		$e2 = Predict::__f * (2 - Predict::__f);
+		$geodetic->lat = Predict_Math::AcTan($pos->z, $r); /*radians*/
+
+		do {
+			$phi    = $geodetic->lat;
+			$sinPhi = sin($phi);
+			$c      = 1 / sqrt(1 - $e2 * ($sinPhi * $sinPhi));
+			$geodetic->lat = Predict_Math::AcTan($pos->z + Predict::xkmper * $c * $e2 * $sinPhi, $r);
+		} while (abs($geodetic->lat - $phi) >= 1E-10);
+
+		$geodetic->alt = $r / cos($geodetic->lat) - Predict::xkmper * $c;/*kilometers*/
+
+		if ($geodetic->lat > Predict::pio2) {
+			$geodetic->lat -= Predict::twopi;
+		}
+	}
+
+	/* The procedures Calculate_Obs and Calculate_RADec calculate         */
+	/* the *topocentric* coordinates of the object with ECI position,     */
+	/* {pos}, and velocity, {vel}, from location {geodetic} at {time}.    */
+	/* The {obs_set} returned for Calculate_Obs consists of azimuth,      */
+	/* elevation, range, and range rate (in that order) with units of     */
+	/* radians, radians, kilometers, and kilometers/second, respectively. */
+	/* The WGS '72 geoid is used and the effect of atmospheric refraction */
+	/* (under standard temperature and pressure) is incorporated into the */
+	/* elevation calculation; the effect of atmospheric refraction on     */
+	/* range and range rate has not yet been quantified.                  */
+
+	/* The {obs_set} for Calculate_RADec consists of right ascension and  */
+	/* declination (in that order) in radians.  Again, calculations are   */
+	/* based on *topocentric* position using the WGS '72 geoid and        */
+	/* incorporating atmospheric refraction.                              */
+	public static function Calculate_Obs($_time, Predict_Vector $pos, Predict_Vector $vel, Predict_Geodetic $geodetic, Predict_ObsSet $obs_set)
+	{
+		$obs_pos = new Predict_Vector();
+		$obs_vel = new Predict_Vector();
+		$range   = new Predict_Vector();
+		$rgvel   = new Predict_Vector();
+
+		self::Calculate_User_PosVel($_time, $geodetic, $obs_pos, $obs_vel);
+
+		$range->x = $pos->x - $obs_pos->x;
+		$range->y = $pos->y - $obs_pos->y;
+		$range->z = $pos->z - $obs_pos->z;
+
+		$rgvel->x = $vel->x - $obs_vel->x;
+		$rgvel->y = $vel->y - $obs_vel->y;
+		$rgvel->z = $vel->z - $obs_vel->z;
+
+		$range->w = sqrt($range->x * $range->x + $range->y * $range->y + $range->z * $range->z);
+
+		$sin_lat   = sin($geodetic->lat);
+		$cos_lat   = cos($geodetic->lat);
+		$sin_theta = sin($geodetic->theta);
+		$cos_theta = cos($geodetic->theta);
+		$top_s = $sin_lat * $cos_theta * $range->x
+			+ $sin_lat * $sin_theta * $range->y
+			- $cos_lat * $range->z;
+		$top_e = -$sin_theta * $range->x
+			+ $cos_theta * $range->y;
+		$top_z = $cos_lat * $cos_theta * $range->x
+			+ $cos_lat * $sin_theta * $range->y
+			+ $sin_lat * $range->z;
+		$azim = atan(-$top_e / $top_s); /*Azimuth*/
+		if ($top_s > 0) {
+			$azim = $azim + Predict::pi;
+		}
+		if ($azim < 0 ) {
+			$azim = $azim + Predict::twopi;
+		}
+		$el = Predict_Math::ArcSin($top_z / $range->w);
+		$obs_set->az = $azim;        /* Azimuth (radians)  */
+		$obs_set->el = $el;          /* Elevation (radians)*/
+		$obs_set->range = $range->w; /* Range (kilometers) */
+
+		/* Range Rate (kilometers/second)*/
+		$obs_set->range_rate = Predict_Math::Dot($range, $rgvel) / $range->w;
+
+		/* Corrections for atmospheric refraction */
+		/* Reference:  Astronomical Algorithms by Jean Meeus, pp. 101-104    */
+		/* Correction is meaningless when apparent elevation is below horizon */
+		//	obs_set->el = obs_set->el + Radians((1.02/tan(Radians(Degrees(el)+
+		//							      10.3/(Degrees(el)+5.11))))/60);
+		if ($obs_set->el < 0) {
+			$obs_set->el = $el;  /*Reset to true elevation*/
+		}
+	}
 }

Spacing +34 added lines, -34 removed lines

@@ -37,18 +37,18 @@ 
 
         $sinGeodeticLat = sin($geodetic->lat); /* Only run sin($geodetic->lat) once */
 
-        $geodetic->theta = Predict_Math::FMod2p(Predict_Time::ThetaG_JD($_time) + $geodetic->lon);/*LMST*/
-        $c = 1 / sqrt(1 + Predict::__f * (Predict::__f - 2) * $sinGeodeticLat * $sinGeodeticLat);
-        $sq = (1 - Predict::__f) * (1 - Predict::__f) * $c;
-        $achcp = (Predict::xkmper * $c + $geodetic->alt) * cos($geodetic->lat);
-        $obs_pos->x = $achcp * cos($geodetic->theta); /*kilometers*/
-        $obs_pos->y = $achcp * sin($geodetic->theta);
-        $obs_pos->z = (Predict::xkmper * $sq + $geodetic->alt) * $sinGeodeticLat;
-        $obs_vel->x = -Predict::mfactor * $obs_pos->y; /*kilometers/second*/
-        $obs_vel->y =  Predict::mfactor * $obs_pos->x;
-        $obs_vel->z =  0;
-        $obs_pos->w = sqrt($obs_pos->x * $obs_pos->x + $obs_pos->y * $obs_pos->y + $obs_pos->z * $obs_pos->z);
-        $obs_vel->w = sqrt($obs_vel->x * $obs_vel->x + $obs_vel->y * $obs_vel->y + $obs_vel->z * $obs_vel->z);
+        $geodetic->theta = Predict_Math::FMod2p(Predict_Time::ThetaG_JD($_time) + $geodetic->lon); /*LMST*/
+        $c = 1/sqrt(1 + Predict::__f*(Predict::__f - 2)*$sinGeodeticLat*$sinGeodeticLat);
+        $sq = (1 - Predict::__f)*(1 - Predict::__f)*$c;
+        $achcp = (Predict::xkmper*$c + $geodetic->alt)*cos($geodetic->lat);
+        $obs_pos->x = $achcp*cos($geodetic->theta); /*kilometers*/
+        $obs_pos->y = $achcp*sin($geodetic->theta);
+        $obs_pos->z = (Predict::xkmper*$sq + $geodetic->alt)*$sinGeodeticLat;
+        $obs_vel->x = -Predict::mfactor*$obs_pos->y; /*kilometers/second*/
+        $obs_vel->y = Predict::mfactor*$obs_pos->x;
+        $obs_vel->z = 0;
+        $obs_pos->w = sqrt($obs_pos->x*$obs_pos->x + $obs_pos->y*$obs_pos->y + $obs_pos->z*$obs_pos->z);
+        $obs_vel->w = sqrt($obs_vel->x*$obs_vel->x + $obs_vel->y*$obs_vel->y + $obs_vel->z*$obs_vel->z);
     }
 
     /* Procedure Calculate_LatLonAlt will calculate the geodetic  */
@@ -56,7 +56,7 @@ 
     /* It is intended to be used to determine the ground track of */
     /* a satellite.  The calculations  assume the earth to be an  */
     /* oblate spheroid as defined in WGS '72.                     */
-    public static function Calculate_LatLonAlt($_time, Predict_Vector $pos,  Predict_Geodetic $geodetic)
+    public static function Calculate_LatLonAlt($_time, Predict_Vector $pos, Predict_Geodetic $geodetic)
     {
         /* Reference:  The 1992 Astronomical Almanac, page K12. */
 
@@ -64,18 +64,18 @@ 
 
         $geodetic->theta = Predict_Math::AcTan($pos->y, $pos->x); /*radians*/
         $geodetic->lon = Predict_Math::FMod2p($geodetic->theta - Predict_Time::ThetaG_JD($_time)); /*radians*/
-        $r = sqrt(($pos->x * $pos->x) + ($pos->y * $pos->y));
-        $e2 = Predict::__f * (2 - Predict::__f);
+        $r = sqrt(($pos->x*$pos->x) + ($pos->y*$pos->y));
+        $e2 = Predict::__f*(2 - Predict::__f);
         $geodetic->lat = Predict_Math::AcTan($pos->z, $r); /*radians*/
 
         do {
             $phi    = $geodetic->lat;
             $sinPhi = sin($phi);
-            $c      = 1 / sqrt(1 - $e2 * ($sinPhi * $sinPhi));
-            $geodetic->lat = Predict_Math::AcTan($pos->z + Predict::xkmper * $c * $e2 * $sinPhi, $r);
+            $c      = 1/sqrt(1 - $e2*($sinPhi*$sinPhi));
+            $geodetic->lat = Predict_Math::AcTan($pos->z + Predict::xkmper*$c*$e2*$sinPhi, $r);
         } while (abs($geodetic->lat - $phi) >= 1E-10);
 
-        $geodetic->alt = $r / cos($geodetic->lat) - Predict::xkmper * $c;/*kilometers*/
+        $geodetic->alt = $r/cos($geodetic->lat) - Predict::xkmper*$c; /*kilometers*/
 
         if ($geodetic->lat > Predict::pio2) {
             $geodetic->lat -= Predict::twopi;
@@ -114,34 +114,34 @@ 
         $rgvel->y = $vel->y - $obs_vel->y;
         $rgvel->z = $vel->z - $obs_vel->z;
 
-        $range->w = sqrt($range->x * $range->x + $range->y * $range->y + $range->z * $range->z);
+        $range->w = sqrt($range->x*$range->x + $range->y*$range->y + $range->z*$range->z);
 
         $sin_lat   = sin($geodetic->lat);
         $cos_lat   = cos($geodetic->lat);
         $sin_theta = sin($geodetic->theta);
         $cos_theta = cos($geodetic->theta);
-        $top_s = $sin_lat * $cos_theta * $range->x
-            + $sin_lat * $sin_theta * $range->y
-            - $cos_lat * $range->z;
-        $top_e = -$sin_theta * $range->x
-            + $cos_theta * $range->y;
-        $top_z = $cos_lat * $cos_theta * $range->x
-            + $cos_lat * $sin_theta * $range->y
-            + $sin_lat * $range->z;
-        $azim = atan(-$top_e / $top_s); /*Azimuth*/
+        $top_s = $sin_lat*$cos_theta*$range->x
+            + $sin_lat*$sin_theta*$range->y
+            - $cos_lat*$range->z;
+        $top_e = -$sin_theta*$range->x
+            + $cos_theta*$range->y;
+        $top_z = $cos_lat*$cos_theta*$range->x
+            + $cos_lat*$sin_theta*$range->y
+            + $sin_lat*$range->z;
+        $azim = atan(-$top_e/$top_s); /*Azimuth*/
         if ($top_s > 0) {
             $azim = $azim + Predict::pi;
         }
-        if ($azim < 0 ) {
+        if ($azim < 0) {
             $azim = $azim + Predict::twopi;
         }
-        $el = Predict_Math::ArcSin($top_z / $range->w);
-        $obs_set->az = $azim;        /* Azimuth (radians)  */
-        $obs_set->el = $el;          /* Elevation (radians)*/
+        $el = Predict_Math::ArcSin($top_z/$range->w);
+        $obs_set->az = $azim; /* Azimuth (radians)  */
+        $obs_set->el = $el; /* Elevation (radians)*/
         $obs_set->range = $range->w; /* Range (kilometers) */
 
         /* Range Rate (kilometers/second)*/
-        $obs_set->range_rate = Predict_Math::Dot($range, $rgvel) / $range->w;
+        $obs_set->range_rate = Predict_Math::Dot($range, $rgvel)/$range->w;
 
         /* Corrections for atmospheric refraction */
         /* Reference:  Astronomical Algorithms by Jean Meeus, pp. 101-104    */
@@ -149,7 +149,7 @@ 
         //	obs_set->el = obs_set->el + Radians((1.02/tan(Radians(Degrees(el)+
         //							      10.3/(Degrees(el)+5.11))))/60);
         if ($obs_set->el < 0) {
-            $obs_set->el = $el;  /*Reset to true elevation*/
+            $obs_set->el = $el; /*Reset to true elevation*/
         }
     }
 }