<?php
include ("/usr/share/jpgraph/jpgraph.php");
include ("/usr/share/jpgraph/jpgraph_scatter.php");


import_request_variables("gP", "rvar_");  // define all $rvar_varnames

// figure out some timeing issues
// default : JD_U=today - (span=7days)
// Use the UNIX Julian Day as the variable 

$curUJD = time()/24/3600;  // Unix Julian Date
$span = 30;
$var = "Wind_spd";
$utcOffset = 0;

if ( $rvar_JD_U )  $curUJD = (float)$rvar_JD_U;
if ( $rvar_span )  $span = (float)$rvar_span;
if ( $rvar_utcOffset )  $utcOffset = (float)$rvar_utcOffset;
//if ( $rvar_var )  $var = $rvar_var;

$lookupVar = array("Wind_dir" , "Wind_spd" , "Wind_spd_std" , "Wind_dir_std" , "Watertemp" , "Salinity" , "D_O" , "depth" , "pH" , "Turbidity" , "Chlorophyll" , "Fluoresence" , "Batteryvolt", "water_depth" );
if ( $rvar_var ) {
   $var = array_values(preg_grep("/$rvar_var/",$lookupVar)); 
   $var=$var[0];
}

//print_r(preg_grep('/.*/',$lookupVar));

// $curUJD = (float)$rvar_JD_U;

 $endISO=    date('Y-m-d\TH:i:sO', $curUJD   *3600*24);
 $startISO = date('Y-m-d\TH:i:sO',($curUJD-$span)*3600*24);

// matlab's datenum(1970,1,1) = 719529 # = days from 0000-1-1 to 1970-1-1

$conn=odbc_connect("surapgtest",'webserver','webserver');
// The 'surapgtest' ODBC DSN points to:
//     server = sura-vims-pe6600-1.vims.edu
//     port   = 5432
// ODBC ipostgres drivers for many platforms ar available at 
//     http://pgfoundry.org/projects/psqlodbc/

if (! $conn )   echo "didn't connect to the surapgtest database\n";

$datatable="VB_H2O_active_raw";
if ($var == 'water_depth') $datatable = 'vb_acu_active_raw';

$sql="select matlab_datenum, year,yearday, $var , time_obs
  from $datatable
  where CAST (yearday* 24*60 as int8) % 1 = 0 
  and time_obs between '$startISO' and '$endISO' 
  order by time_obs;";
$rs = odbc_exec($conn,$sql);


// find the greatest month before $curUJD
$thismonthUJD = mktime(0, 0, 0, date('m',$curUJD*24*3600),1 ) /(3600*24) +1;

while (odbc_fetch_row($rs))
{
//$xdata[] =date('z', odbc_result($rs,"matlab_datenum") * 24 * 3600);
//$xdata[] =odbc_result($rs,"matlab_datenum") - $thismonthUJD; 

$xdata[] =mktime(0,0,odbc_result($rs,'yearday')*24*3600,1,
	1,
	odbc_result($rs,'year')) /(24*3600)-$thismonthUJD-$utcOffset/24;

$tmp = odbc_result($rs,"$var");
if ( $tmp > 1e6 ) // mask the missing values from the graph 
  $tmp = '.'; 

$ydata[] = $tmp;
}

// Create the graph. 
$graph = new Graph(450,350,"auto");	
$graph->SetScale("linlin");
//$graph->img->SetMargin(30,90,40,50);
$graph->xaxis->SetFont(FF_ARIAL);
//$graph->xscale->ticks->Set(7,1);
$graph->xaxis->SetLabelAngle(25); // 45 deg upwards

//echo $graph->SetTickDensity(TICKD_NORMAL,TICKD_SPARSE);
//$graph->title->Set($sql);
$graph->title->Set("$rvar_title UJD=$curUJD UJD_m=$thismonthUJD
 $startISO $endISO");
$graph->xaxis->title->Set(sprintf("fractional day of %s",date("Y-M",$thismonthUJD*24*3600)));
$graph->yaxis->title->Set("$var (units?)");
// Create the linear plot
$plot=new ScatterPlot($ydata,$xdata);
//$plot->SetLegend("mph");
$plot->mark->SetType(0);
$plot->mark->SetWidth(5);

// $sp1->SetWeight(3);
//$plot->mark->SetWidth(4);
$plot->SetLinkPoints();


// Add the plot to the graph
$graph->Add($plot);

// Display the graph
$graph->Stroke();
?>