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