#!/opt/perl/bin/perl
use Smart::Comments;
use File::Basename;
use File::Path;
use Cwd;
my $here = cwd;
my $progname=basename $0;
my $usage = qq(
  $progname dir 

);
my $dir = shift || die $usage;
my $series_dir = "$here/series";
if ( ! -d $series_dir ) {
   mkpath($series_dir) or die "couldn't make $series_dir $!\n";
}

my %gds=();
my ($sta,$t,@env);
#foreach my $f ( sort glob("$dir/*.envcov") ) { ### working on envcov $cal [====|   ] % done
opendir(my $dh,$dir) or die "Couldn't open dir $dir $!\n";

foreach my $f ( sort grep { /envcov/ } readdir($dh)  ) { ### working on envcov $cal [====|   ] % done
   $cal = unpack('a10',$f);
   open(my $fh,'<',"$dir/$f" ) or die "Couldn't open $f $!\n";
   while( <$fh> ) {
     chomp;
     next if /^#/;
     ($sta,$t,@env) = split(' ');
     push @{$gds{$sta}{lon}},"$t $env[0] $env[3] $sta LON $cal\n";
     push @{$gds{$sta}{lat}},"$t $env[1] $env[4] $sta LAT $cal\n";
     push @{$gds{$sta}{rad}},"$t $env[2] $env[5] $sta RAD $cal\n";
    }
}
my $c;
foreach  $c ( qw( lat lon rad ) ) {
  foreach $sta ( sort keys %gds ) { ### working on $sta $c [===|  ] % done
     open( my $fh,'>',"$series_dir/$sta.$c") or die "Couldn't open $sta.$c $!\n";
     print {$fh} @{$gds{$sta}{$c}};
     close $fh;
  }
}


