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
|
#! /usr/bin/perl -w
use strict;
my $thorn = $ARGV[0];
print "Post-processing thorn $thorn\n";
my $filename = "$thorn/src/Boundaries.c";
open (FILE, "< $filename") or die "Cannot read file \"$filename\"";
my @lines = <FILE>;
close FILE or die;
#my $n=0;
#foreach my $line (@lines) {
# print "$n $line";
# ++$n;
#}
$lines[33] = <<EOF;
ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, 1, -1, "${thorn}::ML_Ham", "scalar");
if (ierr<0) CCTK_WARN (CCTK_WARN_ABORT, "Failed to select boundary condition for ${thorn}::ML_Ham");
ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, 1, -1, "${thorn}::ML_mom", "scalar");
if (ierr<0) CCTK_WARN (CCTK_WARN_ABORT, "Failed to select boundary condition for ${thorn}::ML_mom");
EOF
#my $n=0;
#foreach my $line (@lines) {
# print "$n $line";
# ++$n;
#}
open (FILE, "> $filename") or die "Cannot write file \"$filename\"";
print FILE join '', @lines;
close FILE or die;
|
