Adding file that was left out of the previous commit.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@37 f75ba9e5-694f-0410-ac2c-87ea7ce7132b

1 files changed, 29 insertions, 0 deletions

diff --git a/src/fixedsphere.F90 b/src/fixedsphere.F90
new file mode 100644
index 0000000..9b4fba6
--- /dev/null
+++ b/src/fixedsphere.F90
@@ -0,0 +1,29 @@
+
+#include "cctk.h"
+#include "cctk_Parameters.h"
+#include "cctk_Arguments.h"
+
+subroutine Lego_FixedSphere(CCTK_ARGUMENTS)
+
+  implicit none
+
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_FUNCTIONS
+
+  if (CCTK_EQUALS(fixed_excision,"sphere")) then
+    where (r < fixed_size)
+      emask = 0.d0
+    elsewhere
+      emask = 1.d0
+    end where
+  else if (CCTK_EQUALS(fixed_excision,"cube")) then
+    where ((abs(x) < fixed_size).and.(abs(z) < fixed_size).and.&
+         (abs(y) < fixed_size))
+      emask = 0.d0
+    elsewhere
+      emask = 1.d0
+    end where
+  end if
+  
+end subroutine Lego_FixedSphere