#!/usr/bin/perl
#

$d = 5.82;

$other_file = "./drag_2015.x.dat";
open(OTHER_FILE, ">$other_file") || die("can't open file $other_file");


while (<STDIN>) {
  $line = $_;
  if ($line =~ /^#/) { 
    next;
  }
  @words = split(/\s+/, " " . $line);
  $numclips = $words[1];
  $mass = 0.0016 + 0.00038*$numclips;
  $fair = $mass*9.8;

  $num = 0;
  $sum = 0;
  $sumsq = 0;
  for ($i = 2; $i <= $#words; $i++) {
    $time = $words[$i];
    $vel = $d/$time;
    printf OTHER_FILE " %10.5f %6.3f \n", $fair, $vel;

    if ($debug > 0) {
      printf " i = %02d  time_array[%2d] = %lf \n", 
             $i, $words[$i];
    }
    $sum += $time;
    $sumsq += ($time*$time);
    $count ++;
  }
  $mean = $sum / $count;
  if ($count == 1) {
    $stdev = 0.0;
  } 
  else {
    $stdev = sqrt( ($sumsq - $count*$mean*$mean) / ($count - 1.0) );
  }

  $avg_time = $mean;
  $stdev_time = $stdev;

  $vel = $d/$avg_time;
  $maxvel = $d/($avg_time - $stdev_time);
  $minvel = $d/($avg_time + $stdev_time);
  $velsig = ($maxvel - $minvel)/2.0;


  printf " %3d  %8.5f  %4.2f +/- %4.2f  %5.3f +/- %5.3f   %7.4f \n",
         $numclips, $mass, $avg_time, $stdev_time, $vel, $velsig,
         $fair;



}

exit 0;