diff options
author | pollney <pollney@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2003-09-10 10:45:46 +0000 |
---|---|---|
committer | pollney <pollney@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2003-09-10 10:45:46 +0000 |
commit | 43df4d2ecd98e3112845c28b22c0b85d4370bec7 (patch) | |
tree | d60f92242c9038ea66794914c365bfec25a5c98e /src/fixedsphere.F90 | |
parent | 33ec7d1cb74c47a0082682c881cf9aa4f8eb3985 (diff) |
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
Diffstat (limited to 'src/fixedsphere.F90')
-rw-r--r-- | src/fixedsphere.F90 | 29 |
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 |