";
echo "gtime ";
print_r ($gtime);
echo "sunpos ";
print_r ($sunpos);
echo "sunazi ";
print_r ($sunazi);
echo "moonpos ";
print_r ($moonpos);
echo "moonazi ";
print_r ($moonazi);

##########################################################################################################################
Since the denominator in ArcTan [ +0.147 / (-0.320)] is negative, we must add 180° to the final answer:

                     A = -24.7° + 180°
                     A = 155.3°
folowing this I have a correct azimuth HALF the time, don't know how to handle the other half
*/
for ($i=0;$i<25;$i++){
if ($moonazi[$i] <= 0){$moonazi[$i] = $moonazi[$i]+180;}
else if ($moonazi[$i] >> 0){$moonazi[$i] = $moonazi[$i];}
}
##########################################################################################################################
for ($i=0;$i<=24;$i++){$tm11[$i]=$sunazi[$i];}
$saz[0] = $sunazi[25];
for ($i=0;$i<=24;$i++){$tm21[$i]=$moonazi[$i];}
$maz[0] = round($moonazi[25]);
$he5[0] = round($moonpos[25]);
for ($i=0;$i<=24;$i++){$he61[$i]=$moonpos[$i];}
$he7[0] = round($sunpos[25]);
for ($i=0;$i<=24;$i++){$he81[$i]=$sunpos[$i];}
#
##END FUNCTIONS###############################################################################
#
/* CREATE SUN PATH */

if (version_compare(PHP_VERSION, '5.2.1', '>')) {                             
	$sun_info = date_sun_info(strtotime(date("Y/n/j",time())), $lat, $lon);
} else {                                                                      
	$sun_info = date_sun_info(strtotime(date("Y/n/j",time())), $lon, $lat);  // A bug in the earlier versions.
}

$day=$sun_info["sunset"]-$sun_info["sunrise"];

######################################################################
$hours = $day / 3600; // 3600 seconds in an hour
#
$minutes = ($hours - floor($hours)) * 60;
#
$final_minutes = round($minutes);
if ($final_minutes <= 9){$final_minutes = "0" . $final_minutes;}
$dl = $final_hours = floor($hours) . ":" . $final_minutes;
######################################################################

$sr = date("h:i a",$sun_info["sunrise"]);
$zen = date("h:i a",$sun_info["transit"]);
$ss = date("h:i a",$sun_info["sunset"]);

$b = array();
  $a = ($sun_info["transit"]-36000);
for ($i=0;$i<=25;$i++){
$a=$a+3600;
$b[] = $a;
}
foreach ($b as $v) {
  $pos = sun_pos($v,$lat,$lon);
  $az[] = $pos[0];
  $he[] = $pos[1];
}

//Hi position
$sun_info = date_sun_info(strtotime("21 June".date("Y",time())), $lat,$lon);
$b = array();
  $a = ($sun_info["transit"]-36000);
for ($i=0;$i<=49;$i++){
$a=$a+1800;
$b[] = $a;
}
foreach ($b as $v) {
  $pos = sun_pos($v,$lat,$lon);
  $az0[] = $pos[0];
  $he0[] = $pos[1];
}

// Autumnal Eq
$sun_info = date_sun_info(strtotime("21 September".date("Y",time())),$lat,$lon);
$b = array();
  $a = ($sun_info["transit"]-36000);
for ($i=0;$i<=25;$i++){
$a=$a+3600;
$b[] = $a;
}
foreach ($b as $v) {
  $pos = sun_pos($v,$lat,$lon);
  $az3[] = $pos[0];
  $he3[] = $pos[1];
}

//Low position
$sun_info = date_sun_info(strtotime("21 December".date("Y",time())), $lat,$lon);
$b = array();
  $a = ($sun_info["transit"]-32400);
for ($i=0;$i<=25;$i++){
$a=$a+3600;
$b[] = $a;
}
foreach ($b as $v) {
  $pos = sun_pos($v,$lat,$lon);
  $az1[] = $pos[0];
  $he1[] = $pos[1];
}

function getTimestamp ($val) {
    list ($hour,$min) = preg_split("/:/",$val);
    $ttl = ($min * 60) + ($hour * 60 * 60);
    return($ttl);
}


function out_time ( $seconds , $mode = false) {
    $uday = (3600 * 24);
    $uhr = 3600;
    $umin = 60;
    
    if ($seconds < 0 ) {
        $neg = "-" ;
        $seconds = $seconds * -1 ;
        //echo "Seconds = " . $seconds . "
";
        
    } else {
        $neg = "" ;
    }
    
    // Calculate days
    $dd = intval($seconds / $uday );
    $mmremain = ($seconds - ($dd * $uday));
    
    // Calculate hours
    $hh = intval($mmremain / $uhr);    
    $ssremain = ($mmremain - ($hh * $uhr));  
    
    // Calculate minutes
    $mm = intval($ssremain / $umin);
    $ss = ($ssremain - ($mm * $umin));
    
    // If day or days
    if ($dd == 1) { $days = 'day';  }
    if ($dd > 1 ) { $days = "days"; }
    
    if (!$mode) {
        // String for if there are days
        if ( $dd > 0 ) {
            $out_string = sprintf("%d %s %02d hrs %02d mins", 
                $dd, $days, $hh, $mm);
        }
        
        // String if there are hours
        if ( $dd == 0 && $hh > 0 ) {
            $out_string = $neg . sprintf("%d hrs %02d mins", 
                $hh, $mm);
        }
        
        // String if there are minutes
        if ( $dd == 0 && $hh == 0 && $mm > 0 ) {
            $out_string = $neg . sprintf("%d mins", 
                $mm);
        }
        
        // Only output Seconds
        if ( $hh == 0 && $mm == 0 ) {
            $out_string = $neg . sprintf("%d secs", $ss) ;
        }
    } else {
        // String for if there are days
        if ( $dd > 0 ) {
            $out_string = $neg . sprintf("%d:%s %02d:%02d", 
                $dd, $days, $hh, $mm);
        }
        
        // String if there are hours
        if ( $dd == 0 && $hh > 0 ) {
            $out_string = $neg . sprintf("%d:%02d", 
                $hh, $mm);
        }
        
        // String if there are minutes
        if ( $dd == 0 && $hh == 0 && $mm > 0 ) {
            $out_string = $neg . sprintf("00:%02d", 
                $mm);
        }
        
        // Only output Seconds
        if ( $hh == 0 && $mm == 0 ) {
            $out_string = $neg . sprintf("%d secs", $ss) ;
        }       
        
    }
    
    // Return value
    return ($out_string);
}

function time24to12 ($val) {
    if (strlen($val) == 5 ) {
        $val .= ":00";
    }
    
    return( date("g:i", strtotime($val)) );
} 

$mrise = $clientrawextra[558];
$millum = $clientrawextra[560];

/* CREATE ACTUAL SUN POSITION */
$time = time();

// For testing using http://aa.usno.navy.mil/data/docs/AltAz.php & PHP function date_sun_info() results - (for Jim 90° 4' & 43° 48' & 5 hrs W)
/*
$testHr = 19;
$testMin = 25;
$time = mktime($testHr,$testMin,0,date('m'),date('j'),date('Y'));
*/

$pos = sun_pos($time, $lat, $lon);
$az2[0] = $pos[0];
$he2[0] = $pos[1];

// Create the graph. These two calls are always required
    $graph  = new Graph(500,300,"auto",1);
    $graph->SetMargin(25,20,25,55);

	$graph->SetBackgroundGradient('blue',$daycolor,GRAD_HOR,BGRAD_PLOT);                        
	$suncolor = "yellow";  
	                
	if ($az2[0] == 180) { 
		$skyword = 'Solar Noon';
	} else if ($az2[0] == 360) { 
		$skyword = 'Solar Midnight';
	$graph->SetBackgroundGradient('black',$nightcolor,GRAD_HOR,BGRAD_PLOT);    
		
	} else if ($he2[0] <= 3 && $he2[0] > -1) {     // Seems to work out about right
		$suncolor = "darkorange";
		$graph->SetBackgroundGradient('brown',$dawncolor,GRAD_HOR,BGRAD_PLOT); 
		if ($he2[0] >= 0) {                        // Not sunrise or sunset
			if ($az2[0] > 180) {
				$skyword = 'Dusk';
			} else {
				$skyword = 'Dawn';
			}		
		} 
	} else if ($he2[0] == -1) {
		$suncolor = "darkorange";
		$graph->SetBackgroundGradient('brown',$dawncolor,GRAD_HOR,BGRAD_PLOT); 
		if ($az2[0] > 180) {
			$skyword = 'Sunset';
		} else {
			$skyword = 'Sunrise';
		}
	} else if ($he2[0] > 0) {
		if ($az2[0] > 180) {
			$skyword = 'Afternoon';
		} else {
		 	$skyword = 'Morning';
		}
	} else if ($he2[0] >= -6) {          // Civil Twilight
		$graph->SetBackgroundGradient('darkorange',$ctlcolor,GRAD_HOR,BGRAD_PLOT);  
		$skyword = 'Civil Twilight';
	} else if ($he2[0] >= -12) {         // Nautical Twilight
		$graph->SetBackgroundGradient('black',$ntlcolor,GRAD_HOR,BGRAD_PLOT);
		$skyword = 'Nautical Twilight';
	} else if ($he2[0] >= -18) {         // Astronomical Twilight
		$graph->SetBackgroundGradient('black',$atlcolor,GRAD_HOR,BGRAD_PLOT);
		$skyword = 'Astronomical Twilight';
	} else {                             // Must be night
		$graph->SetBackgroundGradient('black',$nightcolor,GRAD_HOR,BGRAD_PLOT);	
		$skyword = 'Nighttime';
	}

  $graph->SetMarginColor("black");
	$fecha = date($dtstring, time()) . " local time";
  $graph->SetClipping();
	$graph->SetScale( "linlin",0,90,40,320);
	$graph ->xgrid->Show( true,true);
	$graph ->xgrid->SetColor('#99999');
	$graph->ygrid->Show(false, false);
//	$graph ->ygrid->Show( true,true); // Prevents background color from displaying!!!!
	$graph ->ygrid->SetColor('#99999');
	$graph->setTickDensity(TICKD_NORMAL,TICKD_DENSE);
	$graph->xaxis->SetPos("min");
	$graph->legend->Hide();
	$graph->xaxis->title->SetFont(FF_ARIAL,FS_BOLD,8);
	$graph->xaxis->title->SetColor("orange");
	$graph->xaxis->title->Set("Sun AZ: $az2[0]°                                       Zenith: $zen                                         Sun EL: $he2[0]°");
	$graph->xaxis->SetTitleMargin(10);
	$graph->yaxis->SetFont(FF_ARIAL,FS_BOLD,7);
	$graph->yaxis->SetColor('white');
	$graph->xaxis->SetFont(FF_ARIAL,FS_BOLD,7);
	$graph->xaxis->SetColor('white');
	$graph->footer->left->SetFont(FF_ARIAL,FS_BOLD,9);
	$graph->footer->left->Set("Sunrise: $sr");
	$graph->footer->left->SetColor("darkorange");
	$graph->footer->center->SetFont(FF_ARIAL,FS_BOLD,9);
	$graph->footer->center->Set("Sunlight: $dl hrs");
	$graph->footer->center->SetColor("darkorange");
	$graph->footer->right->SetFont(FF_ARIAL,FS_BOLD,9);
	$graph->footer->right->Set("Sunset: $ss");
	$graph->footer->right->SetColor("darkorange");
	
// Create the line plot
    $plot = new LinePlot($he, $az);
    $plot->SetColor("orange");
    $plot->SetWeight(1);
    $plot->SetStyle('dashed');

// Create the line plot
    $plot0 = new LinePlot($he0, $az0);
    $plot0->SetColor("green");
    $plot0->SetWeight(1);

// Create the line plot
    $plot1 = new LinePlot($he1, $az1);
    $plot1->SetColor("magenta");
    $plot1->SetWeight(1);

// Create the line plot
    $plot2 = new LinePlot($he3, $az3);
    $plot2->SetColor("cyan");
    $plot2->SetWeight(1);

// Create the line plot for Moon path
    $plot3 = new LinePlot($he61,$tm21);
    $plot3->SetColor("white");
    $plot3->SetWeight(1);
    $plot3->SetStyle('dashed');

// Top Title	
	$graph->title->Set('Sun/Moon position: '.$fecha.' - '.$skyword);      
	$graph->title->SetFont(FF_ARIAL,FS_BOLD,9);
	$graph->title->SetColor('orange','black','white');

// Create the scatter plot
	if ($he2[0] == -1) $he2[0]++;             // Push the sun back up on the horizon when at -1 which is sunrise/set
	$sp1 = new ScatterPlot($he2,$az2);
        $sp1->mark->SetType(MARK_IMG,'Sun48.png',0.8);
//	$sp1->mark->SetType(MARK_FILLEDCIRCLE);
	$sp1->mark->SetWidth(8);
	$sp1->SetImpuls();
	$sp1->mark->SetFillColor($suncolor);
	$sp1->mark->SetColor($suncolor);
	$sp2 = new ScatterPlot($he5,$maz);
    	$sp2->mark->SetType(MARK_IMG,'jpmoon.png',0.6);	

// Add plots to graph
	$graph->Add($plot0);
	$graph->Add($plot1);
	$graph->Add($plot2);
	$graph->Add($plot3); //Moon
	$graph->Add($plot);
	if($he2[0]>-1)$graph->Add($sp1);
	$graph->Add($sp2);

// Add labels for equinox peaks	                
	$txt =new Text("Jun 21");
	$txt->SetFont(FF_ARIAL,FS_BOLD,7);
	$txt->SetPos( 237,40);
	$txt->SetColor( "green");
	$graph->AddText($txt);	
	$txt1 =new Text("Dec 21");
	$txt1->SetFont(FF_ARIAL,FS_BOLD,7);
	$txt1->SetPos( 237,173);
	$txt1->SetColor( "magenta");
	$graph->AddText($txt1);	
	$txt2 =new Text("Equinox");
	$txt2->SetFont(FF_ARIAL,FS_BOLD,7);
	$txt2->SetPos( 233,96);
	$txt2->SetColor( "cyan");
	$graph->AddText($txt2);
	$txt3 =new Text("Moon EL: $he5[0]°");
	$txt3->SetFont(FF_ARIAL,FS_BOLD,8);
	$txt3->SetPos(408,276);
	$txt3->SetColor( "white");
	$graph->AddText($txt3);	
	$txt4 =new Text("Moonrise: $mrise");
	$txt4->SetFont(FF_ARIAL,FS_BOLD,8);
	$txt4->SetPos(15,276);
	$txt4->SetColor( "white");
	$graph->AddText($txt4);	
	$txt5 =new Text("Illumination: $millum%");
	$txt5->SetFont(FF_ARIAL,FS_BOLD,8);
	$txt5->SetPos(208,276);
	$txt5->SetColor( "white");
	$graph->AddText($txt5);	

// Save the graph
  $graph->Stroke();

?>