blob: e4680a893061d61f14ebfbff8df046d7c343a128 (
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
60
61
|
#! /usr/bin/perl -w
#/*@@
# @file mergeCarpetIOScalar.pl
# @date Tue 15 August 2006
# @author Thomas Radke
# @desc
# Perl script to merge CarpetIOScalar output files,
# eliminating duplicate timesteps.
#
# Source code comments also by Thomas Radke :-)
# @enddesc
#@@*/
my $help = $#ARGV < 0;
for (my $arg = 0; $arg <= $#ARGV; $arg++) {
$help |= ($ARGV[$arg] eq '-h' or
$ARGV[$arg] eq '-help' or
$ARGV[$arg] eq '--help');
}
if ($help) {
print << "EOF";
Usage: $0 [-h | -help | --help] <list of files>
This script can be used to merge CarpetIOScalar output written before
and after recovery. It reads one or more files in CarpetIOScalar format
and writes their contents to STDOUT, eliminating duplicate timesteps
(all but the last occurance are discarded).
Example: $0 alp.norm1.asc > alp.norm1.asc.merged
EOF
exit;
}
# Rauten-Feld zum Merken der Anzahl und Haeufigkeit vorhandener Datensaetze
my %timesteps = ();
# Liste aller Eingabe-Dateien
my @filelist = @ARGV;
# lies zeilenweise alle Eingabe-Dateien
while (<>) {
# falls diese Zeile eine Datenzeile ist:
# vermerke den Datensatz mit seiner Iterationsnummer
++$timesteps{$1} if (/^(\d+)\s/);
}
# stelle die Liste aller Eingabe-Dateien wieder her
@ARGV = @filelist;
# lies zeilenweise alle Eingabe-Dateien
while (<>) {
# falls diese Zeile eine Datenzeile ist:
# ueberspringe alle Datensaetze mit dieser Iterationsnummer
# bis auf den letzten
print unless (/^(\d+)\s/ and --$timesteps{$1} > 0);
}
|
