Use mask without using excision calls, this may need some more tweaking

depending on how excision labels points in the mask, but this will be
sorted out in the near future.

Converted to F77, then realised that this might be a bad idea since people
might be including F90 in the CalcTmunu stuff, so left the files labeled
as F90 for now.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@69 b7a48df3-cbbf-4440-997f-b4b717c9f7fc

17 files changed, 126 insertions, 0 deletions

diff --git a/test/ADMConstraints_Mask/ADMConstraints_Mask.par b/test/ADMConstraints_Mask/ADMConstraints_Mask.par
new file mode 100644
index 0000000..9aa0995
--- /dev/null
+++ b/test/ADMConstraints_Mask/ADMConstraints_Mask.par
@@ -0,0 +1,46 @@
+!DESC "Constraints for Schwarzschild with mask activated"
+
+# $Header$
+#
+# /*@@
+#   @file      test_ADM_1.par
+#   @date      Sunday 24th July
+#   @author    Gabrielle Allen
+#   @desc
+#   ADMConstraints test file
+#   @enddesc
+# @@*/           
+
+# Required thorns
+ActiveThorns = "Boundary Time ADMconstraints pugh pughreduce cartgrid3d ADMBase StaticConformal admmacros spacemask idanalyticbh ioutil iobasic"
+
+# GENERAL
+
+driver::global_nx  =  20
+driver::global_ny  =  20
+driver::global_nz  =  20
+
+grid::type         =  "byspacing"
+grid::dxyz         =  0.2
+grid::domain     =  "octant"
+
+cactus::cctk_itlast        =  0
+
+
+# Einstein
+
+SpaceMask::use_mask       = "yes"
+ADMBase::initial_data     = "schwarzschild"
+ADMBase::metric_type      = "static conformal"
+
+# Output
+
+IO::out_dir         = "ADMConstraints_Mask"
+IO::out_fileinfo   = "none"
+
+IOBasic::outScalar_every = 1 
+IOBasic::outScalar_vars  = "admconstraints::hamiltonian admconstraints::momentum"
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_vars = "admconstraints::ham admbase::gxx"
+
diff --git a/test/ADMConstraints_Mask/ham_maximum.xg b/test/ADMConstraints_Mask/ham_maximum.xg
new file mode 100644
index 0000000..ad2b029
--- /dev/null
+++ b/test/ADMConstraints_Mask/ham_maximum.xg
@@ -0,0 +1,5 @@
+"ham v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0002044852710
+0.2000000000000	0.0007982173981
+0.3000000000000	0.0017258133454
diff --git a/test/ADMConstraints_Mask/ham_minimum.xg b/test/ADMConstraints_Mask/ham_minimum.xg
new file mode 100644
index 0000000..f342654
--- /dev/null
+++ b/test/ADMConstraints_Mask/ham_minimum.xg
@@ -0,0 +1,5 @@
+"ham v time
+0.0000000000000	-0.0000000000000
+0.1000000000000	-0.0001124590969
+0.2000000000000	-0.0004379400554
+0.3000000000000	-0.0009426106873
diff --git a/test/ADMConstraints_Mask/ham_norm1.xg b/test/ADMConstraints_Mask/ham_norm1.xg
new file mode 100644
index 0000000..d923144
--- /dev/null
+++ b/test/ADMConstraints_Mask/ham_norm1.xg
@@ -0,0 +1,5 @@
+"ham v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0000159379883
+0.2000000000000	0.0000628594942
+0.3000000000000	0.0001383184346
diff --git a/test/ADMConstraints_Mask/ham_norm2.xg b/test/ADMConstraints_Mask/ham_norm2.xg
new file mode 100644
index 0000000..74a3289
--- /dev/null
+++ b/test/ADMConstraints_Mask/ham_norm2.xg
@@ -0,0 +1,5 @@
+"ham v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0000382103057
+0.2000000000000	0.0001490771389
+0.3000000000000	0.0003222790671
diff --git a/test/ADMConstraints_Mask/momx_maximum.xg b/test/ADMConstraints_Mask/momx_maximum.xg
new file mode 100644
index 0000000..e5594a2
--- /dev/null
+++ b/test/ADMConstraints_Mask/momx_maximum.xg
@@ -0,0 +1,5 @@
+"momx v time
+0.0000000000000	-0.0000000000000
+0.1000000000000	0.0055634120450
+0.2000000000000	0.0111123908075
+0.3000000000000	0.0166325230605
diff --git a/test/ADMConstraints_Mask/momx_minimum.xg b/test/ADMConstraints_Mask/momx_minimum.xg
new file mode 100644
index 0000000..dc5df80
--- /dev/null
+++ b/test/ADMConstraints_Mask/momx_minimum.xg
@@ -0,0 +1,5 @@
+"momx v time
+0.0000000000000	-0.0000000000000
+0.1000000000000	-0.0055634120450
+0.2000000000000	-0.0111123908075
+0.3000000000000	-0.0166325230605
diff --git a/test/ADMConstraints_Mask/momx_norm1.xg b/test/ADMConstraints_Mask/momx_norm1.xg
new file mode 100644
index 0000000..d9c62fe
--- /dev/null
+++ b/test/ADMConstraints_Mask/momx_norm1.xg
@@ -0,0 +1,5 @@
+"momx v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0000998180019
+0.2000000000000	0.0002004327396
+0.3000000000000	0.0003032309573
diff --git a/test/ADMConstraints_Mask/momx_norm2.xg b/test/ADMConstraints_Mask/momx_norm2.xg
new file mode 100644
index 0000000..ffbb440
--- /dev/null
+++ b/test/ADMConstraints_Mask/momx_norm2.xg
@@ -0,0 +1,5 @@
+"momx v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0003573831141
+0.2000000000000	0.0007153714352
+0.3000000000000	0.0010746503210
diff --git a/test/ADMConstraints_Mask/momy_maximum.xg b/test/ADMConstraints_Mask/momy_maximum.xg
new file mode 100644
index 0000000..fc20965
--- /dev/null
+++ b/test/ADMConstraints_Mask/momy_maximum.xg
@@ -0,0 +1,5 @@
+"momy v time
+0.0000000000000	-0.0000000000000
+0.1000000000000	0.0055634120450
+0.2000000000000	0.0111123908075
+0.3000000000000	0.0166325230605
diff --git a/test/ADMConstraints_Mask/momy_minimum.xg b/test/ADMConstraints_Mask/momy_minimum.xg
new file mode 100644
index 0000000..7b818df
--- /dev/null
+++ b/test/ADMConstraints_Mask/momy_minimum.xg
@@ -0,0 +1,5 @@
+"momy v time
+0.0000000000000	-0.0000000000000
+0.1000000000000	-0.0055634120450
+0.2000000000000	-0.0111123908075
+0.3000000000000	-0.0166325230605
diff --git a/test/ADMConstraints_Mask/momy_norm1.xg b/test/ADMConstraints_Mask/momy_norm1.xg
new file mode 100644
index 0000000..734505a
--- /dev/null
+++ b/test/ADMConstraints_Mask/momy_norm1.xg
@@ -0,0 +1,5 @@
+"momy v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0000998180019
+0.2000000000000	0.0002004327396
+0.3000000000000	0.0003032309573
diff --git a/test/ADMConstraints_Mask/momy_norm2.xg b/test/ADMConstraints_Mask/momy_norm2.xg
new file mode 100644
index 0000000..479b003
--- /dev/null
+++ b/test/ADMConstraints_Mask/momy_norm2.xg
@@ -0,0 +1,5 @@
+"momy v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0003573831141
+0.2000000000000	0.0007153714352
+0.3000000000000	0.0010746503210
diff --git a/test/ADMConstraints_Mask/momz_maximum.xg b/test/ADMConstraints_Mask/momz_maximum.xg
new file mode 100644
index 0000000..52bbbf3
--- /dev/null
+++ b/test/ADMConstraints_Mask/momz_maximum.xg
@@ -0,0 +1,5 @@
+"momz v time
+0.0000000000000	-0.0000000000000
+0.1000000000000	0.0055634120450
+0.2000000000000	0.0111123908075
+0.3000000000000	0.0166325230605
diff --git a/test/ADMConstraints_Mask/momz_minimum.xg b/test/ADMConstraints_Mask/momz_minimum.xg
new file mode 100644
index 0000000..51ac7bb
--- /dev/null
+++ b/test/ADMConstraints_Mask/momz_minimum.xg
@@ -0,0 +1,5 @@
+"momz v time
+0.0000000000000	-0.0000000000000
+0.1000000000000	-0.0055634120450
+0.2000000000000	-0.0111123908075
+0.3000000000000	-0.0166325230605
diff --git a/test/ADMConstraints_Mask/momz_norm1.xg b/test/ADMConstraints_Mask/momz_norm1.xg
new file mode 100644
index 0000000..4600896
--- /dev/null
+++ b/test/ADMConstraints_Mask/momz_norm1.xg
@@ -0,0 +1,5 @@
+"momz v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0000998180019
+0.2000000000000	0.0002004327396
+0.3000000000000	0.0003032309573
diff --git a/test/ADMConstraints_Mask/momz_norm2.xg b/test/ADMConstraints_Mask/momz_norm2.xg
new file mode 100644
index 0000000..11f584a
--- /dev/null
+++ b/test/ADMConstraints_Mask/momz_norm2.xg
@@ -0,0 +1,5 @@
+"momz v time
+0.0000000000000	0.0000000000000
+0.1000000000000	0.0003573831141
+0.2000000000000	0.0007153714352
+0.3000000000000	0.0010746503210