root/cached_images/trunk/CGI_wrappers/wms_seacoos_in_situ

#!/usr/bin/perl -w

use CGI qw(:standard);
my $query = new CGI;

# look at an 1.5 hours ago as default time
my $this_time_stamp = gmtime(time-60*60*1.5);

# see if an offset from now was passed
$this_time_offset_hours = $query->param('time_offset_hours');
if (length($this_time_offset_hours) > 0) {
 $this_time_offset_hours = $query->param('time_offset_hours');
 $this_time_stamp = gmtime(time-60*60*$this_time_offset_hours);
}

# get the passed time_stamp
if ($query->param('time')) {
  $this_time_stamp = $query->param('time');
}
elsif ($query->param('TIME')) {
  $this_time_stamp = $query->param('TIME');
}
elsif ($query->param('Time')) {
  $this_time_stamp = $query->param('Time');
}
if ($this_time_stamp =~ /current/) {
  $this_time_stamp = gmtime(time-60*60*1.5);
}

# any time zone (only for time formatting) -- assume GMT
my $time_zone_for_formatting_only = 'GMT';
if ($query->param('time_zone_for_formatting_only')) {
  $time_zone_for_formatting_only = $query->param('time_zone_for_formatting_only');
}

# get current time
$cmd = 'date -d \''.$this_time_stamp.' UTC\' +%s';
my $this_time_stamp_sec = `$cmd`;
chop($this_time_stamp_sec);

# get bbox
my $this_bbox_param = $query->param('bbox');
if (length($this_bbox_param) < 1) {
  $this_bbox_param = $query->param('BBOX');
}
if (length($this_bbox_param) < 1) {
  $this_bbox_param = '0,0,0,0';
}
my @this_bbox = split(/,/,$this_bbox_param);

# logging stuff
my $log_map_name = 'wms_seacoos_in_situ';
my $log_map_file = '/opt/ogc/logs/wms_in_situ_log.txt';
my $log_map_time_stamp = "timestamp with time zone '$this_time_stamp UTC'";
my $log_remote_ip = $query->remote_addr();
my $log_extents = 'GeometryFromText(\'POLYGON(('
  .$this_bbox[0]
  .' '.$this_bbox[1]
  .', '.$this_bbox[0]
  .' '.$this_bbox[3]
  .', '.$this_bbox[2]
  .' '.$this_bbox[3]
  .', '.$this_bbox[2]
  .' '.$this_bbox[1]
  .', '.$this_bbox[0]
  .' '.$this_bbox[1]
  .'))\',-1)';
my $log_center_point = 'GeometryFromText(\'POINT('
  .(($this_bbox[2] + $this_bbox[0])/2)
  .' '.(($this_bbox[3] + $this_bbox[1])/2)
  .', '.$this_bbox[0]
  .')\',-1)';
  
sub log_layer {
  $this_layer_name = shift @_;
  $log_sql = 'insert into map_log ('
    .'access_time_stamp'
    .',map_time_stamp'
    .',source_ip'
    .',referer'
    .',layer_name'
    .',map_name'
    .',extents'
    .',center_point'
    .')'
    .' select'
    .' timestamp with time zone \''
    .gmtime.' UTC'
    .'\''
    .','.$log_map_time_stamp
    .',\''.$log_remote_ip.'\''
    .',\''.$log_referer.'\''
    .',\''.$this_layer_name.'\''
    .",'$log_map_name'"
    .','.$log_extents
    .','.$log_center_point
   .';';
  $log_cmd = "echo \"$log_sql\" >> $log_map_file";
  `$log_cmd`;
}

# check to see if $this_time_stamp should require activation of the archive layers
# < 2 weeks
my $archive_flag = 0;
if ($this_time_stamp_sec < time-60*60*24*14) {
  $archive_flag = 1;
}

# station_id filter is optional as a parameter
my $station_id = '';
if ($query->param('station_id')) {
  $station_id = $query->param('station_id');
}

# water_level_reference is optional as a parameter but required come showtime
# valid values: MSL (default), MLLW, NAVD88
my $water_level_reference = 'MSL';
if ($query->param('water_level_reference')) {
  $water_level_reference = $query->param('water_level_reference');
}

# set default units, and get any if passed
# valid values: mps (default), mph, knots
my $wind_label_char_column = 'label_char';
if ($query->param('wind_speed_units')) {
  $wind_label_char_column .= '_'.$query->param('wind_speed_units');
  if ($wind_label_char_column =~ /mps/) {
    $wind_label_char_column = 'label_char';
  }
}

# get all the layers
my $this_layers = $query->param('layers');
if (length($this_layers) < 1) {
  $this_layers = $query->param('LAYERS');
}

# produce layer-specific params
my $layers_params = '';
if ($this_layers =~ /sst_obs_hourly/) {
  if ($archive_flag == 1) {
    $layers_params .= 'sst_obs_hourly_archive,';
  }
  else {
    $layers_params .= 'sst_obs_hourly_recent,';
  }
  log_layer('sst_obs_hourly');
}
if ($this_layers =~ /wind_obs_hourly/) {
  if ($archive_flag == 1) {
    $layers_params .= 'wind_obs_hourly_archive,';
  }
  else {
    $layers_params .= 'wind_obs_hourly_recent,';
  }
  log_layer('wind_obs_hourly');
}
if ($this_layers =~ /wind_obs_normalized_hourly/) {
  if ($archive_flag == 1) {
    $layers_params .= 'wind_obs_normalized_hourly_archive,';
  }
  else {
    $layers_params .= 'wind_obs_normalized_hourly_recent,';
  }
  log_layer('wind_obs_normalized_hourly');
}
if ($this_layers =~ /water_level_obs_hourly/) {
  if ($archive_flag == 1) {
    $layers_params .= 'water_level_obs_hourly_archive,';
  }
  else {
    $layers_params .= 'water_level_obs_hourly_recent,';
  }
  log_layer('water_level_obs_hourly');
}
if ($this_layers =~ /timestamp_label/) {
  $layers_params .= 'timestamp_label,';
}

# get rid of last comma
chop($layers_params);

# get all the query_layers
my $this_query_layers = $query->param('query_layers');
if (length($this_query_layers) < 1) {
  $this_query_layers = $query->param('QUERY_LAYERS');
}

# produce query_layer-specific params
my $query_layers_params = '';
if ($this_query_layers =~ /sst_obs_hourly/) {
  if ($archive_flag == 1) {
    $query_layers_params .= 'sst_obs_hourly_archive,';
  }
  else {
    $query_layers_params .= 'sst_obs_hourly_recent,';
  }
}
if ($this_query_layers =~ /wind_obs_hourly/) {
  if ($archive_flag == 1) {
    $query_layers_params .= 'wind_obs_hourly_archive,';
  }
  else {
    $query_layers_params .= 'wind_obs_hourly_recent,';
  }
}
if ($this_query_layers =~ /wind_obs_normalized_hourly/) {
  if ($archive_flag == 1) {
    $query_layers_params .= 'wind_obs_normalized_hourly_archive,';
  }
  else {
    $query_layers_params .= 'wind_obs_normalized_hourly_recent,';
  }
}
if ($this_query_layers =~ /water_level_obs_hourly/) {
  if ($archive_flag == 1) {
    $query_layers_params .= 'water_level_obs_hourly_archive,';
  }
  else {
    $query_layers_params .= 'water_level_obs_hourly_recent,';
  }
}

# get rid of last comma
chop($query_layers_params);

# get only the params (subtract the URL)
my $orig_url_params = substr($query->url(-path_info=>0,-query=>1),length($query->url()));

# hide original layer request
$orig_url_params =~ s/layers=/layers_hide=/g;
$orig_url_params =~ s/LAYERS=/layers_hide=/g;

# build complete target URL
my $target_url = 'http://maury.marine.unc.edu/cgi-bin/in_situ'
  .$orig_url_params
  .'&formatted_time_stamp='.$this_time_stamp
  .'&layers='.$layers_params
  .'&query_layers='.$query_layers_params
  .'&wind_label_char_column='.$wind_label_char_column
  .'&water_level_reference='.$water_level_reference
  .'&station_id='.$station_id
  .'&time_zone_for_formatting_only='.$time_zone_for_formatting_only;
  
# replace ; w/ & (don't know why ;'s are done by default)
$target_url =~ s/;/&/g;
# url-encode (sort of) a space
$target_url =~ s/ /%20/g;

# redirect to the target URL w/ the completed raster filenames
print redirect($target_url);