diff options
| author | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2002-05-20 12:43:50 +0000 |
|---|---|---|
| committer | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2002-05-20 12:43:50 +0000 |
| commit | f4c43aa510c36db424ab465e3b3894aa9463f25b (patch) | |
| tree | d67f42a2b40aeadec5da6261257171a2c4e28fc9 /src/EHFinder_mod.F90 | |
| parent | d15331dcc9cb3e73e829fe9f4145baf09ed37ed3 (diff) | |
First partial working version. Only re-parametrization with pde using
second order upwinded differences are working. Haven't checked octant or
any other symmetries yet. Approximate re-parametrization was broken when
I tried to add support for octant symmetry. Not yet converted to Einstein 2.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@5 2a26948c-0e4f-0410-aee8-f1d3e353619c
Diffstat (limited to 'src/EHFinder_mod.F90')
| -rw-r--r-- | src/EHFinder_mod.F90 | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/src/EHFinder_mod.F90 b/src/EHFinder_mod.F90 new file mode 100644 index 0000000..2c88ee5 --- /dev/null +++ b/src/EHFinder_mod.F90 @@ -0,0 +1,24 @@ +#include "cctk.h" + +module EHFinder_mod + + use EHFinder_Constants + + CCTK_INT, dimension(0:5), parameter :: ix = (/ -1, 1, 0, 0, 0, 0 /) + CCTK_INT, dimension(0:5), parameter :: jy = (/ 0, 0, -1, 1, 0, 0 /) + CCTK_INT, dimension(0:5), parameter :: kz = (/ 0, 0, 0, 0, -1, 1 /) + CCTK_INT, dimension(0:5), parameter :: ll = (/ 1, 2, 4, 8, 16, 32 /) +! CCTK_INT, parameter :: ixl = 1 +! CCTK_INT, parameter :: ixr = 2 +! CCTK_INT, parameter :: iyl = 4 +! CCTK_INT, parameter :: iyr = 8 +! CCTK_INT, parameter :: izl = 16 +! CCTK_INT, parameter :: izr = 32 + CCTK_REAL :: hfac + CCTK_INT :: nx, ny, nz, ngx, ngy, ngz, ierr + CCTK_INT :: max_handle, min_handle, sum_handle + CCTK_INT :: ncalls, rep_count, rep_total + CCTK_REAL :: delta + CCTK_REAL :: h, hlocal + logical :: mask_first = .true. +end module EHFinder_mod |