blob: 5c312d66fc2b3b471fdabe890f8d0b318a9b6027 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
|
#! /usr/bin/perl -s
use FileHandle;
if(@ARGV != 3)
{
print "Usage: $0 direction <Inputfile> <Outputfile> \n";
exit;
}
open(CARPETFILE, "<$ARGV[1]") || die "Unable to find file \"$ARGV[1]\".";
$file = $ARGV[2].".xg";
my $fh = new FileHandle(">$file") || die "Unable to open file \"$file\".";
my %data;
my $time = -1;
my $new = 0;
while (<CARPETFILE>)
{
$line = $_;
if ($line =~ /^#/)
{
#Do nothing for headers!
}
elsif ($line =~ /([0-9+-ed.]+)*/)
{
@dataline = split(/ +/,$line);
if ($dataline[8] - $time > 1.e-10)
{
if ($new)
{
print $fh "\n\n\"Time = ",$time,"\n";
my @sortedcoords = sort numerically (keys %data);
my $localcoord;
foreach $localcoord (@sortedcoords)
{
print $fh $localcoord." ".$data{$localcoord};
}
}
$new++;
$time = $dataline[8];
}
my $coord = $dataline[9+$ARGV[0]];
my $val = $dataline[12];
$data{$coord} = $val;
}
}
print $fh "\n\"Time = ",$time,"\n";
my @sortedcoords = sort numerically (keys %data);
my $localcoord;
foreach $localcoord (@sortedcoords)
{
print $fh $localcoord." ".$data{$localcoord};
}
sub numerically {$a <=> $b;}
|