Correct boundary conditions

1 files changed, 36 insertions, 0 deletions

diff --git a/m/fixup-ADMConstraints.pl b/m/fixup-ADMConstraints.pl
new file mode 100755
index 0000000..faf68b6
--- /dev/null
+++ b/m/fixup-ADMConstraints.pl
@@ -0,0 +1,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;