#!/usr/bin/perl
use strict;

=comment

This program assumes creates a netcdf file to be pulled by the seacoos data loader/aggregation scripts

=cut

use XML::LibXML;

my ($operator_name,$platform_name) = @ARGV;

`unzip -p -qq /var/www/html/obskml/feeds/scdnr/scdnr_metadata_latest.kmz | cat > latest_obskml.kml`;

my $content;
my $measurement_date;

my $file_ref = "obskml_$operator_name\_$platform_name";
my $netcdf_header = './netcdf/obskml_header.txt';
my $netcdf_final = './netcdf/obskml_final.txt'; 

open (NETCDF_FILE, ">$netcdf_final");
$content = `cat $netcdf_header`;

my $xp = XML::LibXML->new->parse_file('latest_obskml.kml');

#my $xp_xml = $xp->serialize; #debug
#print $xp_xml;

foreach my $platform ($xp->findnodes('//Placemark')) {
        #print "platform:$platform\n"; #debug

	my $lon_lat_string = sprintf("%s",$platform->find('Point/coordinates'));
	my ($lon,$lat) = split(/,/,$lon_lat_string);
	#print "$lon $lat \n";

	$content =~ s/<LON>/$lon/g;
	$content =~ s/<LAT>/$lat/g;

	my $timestamp = sprintf("%s",$platform->find('TimeStamp'));
	#print "$time \n";
	$timestamp =~ s/T/ /g;

        $measurement_date = substr($timestamp,0,19);
        my $time_zone = substr($timestamp,19,3).substr($timestamp,23,2);
	#print "$measurement_date $time_zone \n";
        my $time_sec = `date --date='$measurement_date $time_zone' +%s`;
        chomp($time_sec);

        $content =~ s/<DATETIME>/$measurement_date$time_zone/g;
        $content =~ s/<TIME>/$time_sec/g;

        $content =~ s/<DEPTH_SURFACE>/0/g;

	my $operator_url = sprintf("%s",$platform->find('Metadata/obsList/operatorURL'));
        $content =~ s/<OPERATOR_URL>/$operator_url/g;

        $content =~ s/<OPERATOR_NAME>/$operator_name/g;
        $content =~ s/<PLATFORM_NAME>/$platform_name/g;

	foreach my $obs ($platform->findnodes('Metadata/obsList/obs')) {
		my $value = sprintf("%s",$obs->find('value'));
		my $elev = sprintf("%s",$obs->find('elev'));
		if ($elev <= 0) { $content =~ s/<DEPTH>/$elev/g; } #assuming one positive depth represents all depths
		my $obs_type = sprintf("%s",$obs->find('obsType'));
		#print "$obs_type $value $elev \n";

		#swap variable macros as values present
		if ($obs_type eq 'air_temperature') { $content =~ s/<AIR_TEMPERATURE>/$value/g; }
		if ($obs_type eq 'air_pressure') { $content =~ s/<AIR_PRESSURE>/$value/g; }
		if ($obs_type eq 'relative_humidity') { $content =~ s/<RELATIVE_HUMIDITY>/$value/g; }
	
		if ($obs_type eq 'wind_speed') { $content =~ s/<WIND_SPEED>/$value/g; }
		if ($obs_type eq 'wind_from_direction') { $content =~ s/<WIND_FROM_DIRECTION>/$value/g; }
		if ($obs_type eq 'wind_gust') { $content =~ s/<WIND_GUST>/$value/g; }

		if ($obs_type eq 'water_temperature') { $content =~ s/<WATER_TEMPERATURE>/$value/g; }
		if ($obs_type eq 'salinity') { $content =~ s/<SALINITY>/$value/g; }
			
	}

	#replace empty macros with missing value
	my @obs_array = qw(<WIND_SPEED> <WIND_FROM_DIRECTION> <WIND_GUST> <AIR_TEMPERATURE> <AIR_PRESSURE> <RELATIVE_HUMIDITY> <SOLAR_RADIATION> <SALINITY> <SEA_BOTTOM_SALINITY> <WATER_TEMPERATURE> <SEA_BOTTOM_TEMPERATURE> <SEA_BOTTOM_CURRENT_SPEED> <SEA_BOTTOM_CURRENT_TO_DIRECTION> <SEA_BOTTOM_EASTWARD_CURRENT> <SEA_BOTTOM_NORTHWARD_CURRENT> <SEA_SURFACE_CURRENT_SPEED> <SEA_SURFACE_CURRENT_TO_DIRECTION> <SEA_SURFACE_EASTWARD_CURRENT> <SEA_SURFACE_NORTHWARD_CURRENT> <CURRENT_SPEED> <CURRENT_TO_DIRECTION> <EASTWARD_CURRENT> <NORTHWARD_CURRENT>);

	foreach my $element (@obs_array) {
		$content =~ s/$element/-99999/g;
	}
}

print NETCDF_FILE $content;
close (NETCDF_FILE);

my $timestamp_prefix = `date --date='$measurement_date'  +"%Y%m%d%H%M%S"`;
chomp($timestamp_prefix);

my $filename_old = $file_ref."_".$timestamp_prefix."_old.nc";
my $filename_latest = $file_ref."_".$timestamp_prefix."_latest.nc";

my $netcdf_path = '/usr2/prod/buoys/processed/netcdf';
`/usr2/home/jcothran/storm_surge/ncgen -o $netcdf_path/$filename_old $netcdf_final`;
`cp -f $netcdf_path/$filename_old $netcdf_path/$filename_latest`;

exit 0;