From 33479637cddfa24c0150cb691ee0bc1695f42866 Mon Sep 17 00:00:00 2001
From: Erik Schnetter <schnetter@cct.lsu.edu>
Date: Sun, 14 Mar 2010 14:11:03 -0500
Subject: Correct boundary conditions

---
 m/fixup-ADMConstraints.pl | 36 ++++++++++++++++++++++++++++++++++++
 1 file changed, 36 insertions(+)
 create mode 100755 m/fixup-ADMConstraints.pl

(limited to 'm/fixup-ADMConstraints.pl')

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;
-- 
cgit v1.2.3