#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"

#include "cctk_Functions.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