aboutsummaryrefslogtreecommitdiff
path: root/src/AHFinder.F
diff options
context:
space:
mode:
Diffstat (limited to 'src/AHFinder.F')
-rw-r--r--src/AHFinder.F96
1 files changed, 1 insertions, 95 deletions
diff --git a/src/AHFinder.F b/src/AHFinder.F
index 0d646a6..8b25e30 100644
--- a/src/AHFinder.F
+++ b/src/AHFinder.F
@@ -444,98 +444,8 @@
nonaxi = .false.
end if
-! Find out reflection symmetries depending on
-! the values of the parameters.
- if ((CCTK_Equals(ahf_octant,"yes").eq.1).or.
- . (CCTK_Equals(ahf_octant,"high").eq.1)) then
-
- refx = .true.
- refy = .true.
- refz = .true.
-
- else
-
- if (ahf_refx.ne.0) then
- refx = .true.
- else
- refx = .false.
- end if
-
- if (ahf_refy.ne.0) then
- refy = .true.
- else
- refy = .false.
- end if
-
- if (ahf_refz.ne.0) then
- refz = .true.
- else
- refz = .false.
- end if
-
- end if
-
- if (ahf_phi.eq.0) then
- refx = .true.
- refy = .true.
- end if
-
-! Force grid symmetries. If the grid has some
-! symmetry (octant, quadrant, bitant), and the
-! horizon is centered on one of the symmetry
-! planes, then the grid symmetries are forced
-! on it regardless of the values of the symmetry
-! parameters.
-
- if (CCTK_Equals(domain,"octant").eq.1) then
-
- if (xc.eq.zero) refx = .true.
- if (yc.eq.zero) refy = .true.
- if (zc.eq.zero) refz = .true.
-
- else if (CCTK_Equals(domain,"quadrant").eq.1) then
-
- if (xc.eq.zero) refx = .true.
- if (yc.eq.zero) refy = .true.
-
- else if (CCTK_Equals(domain,"bitant").eq.1) then
-
- if (CCTK_Equals(bitant_plane,"xy").eq.1) then
- if (zc.eq.zero) refz = .true.
- else if (CCTK_Equals(bitant_plane,"xz").eq.1) then
- if (yc.eq.zero) refy = .true.
- else
- if (xc.eq.zero) refx = .true.
- end if
-
- end if
-
-! Check if we are using cartoon. If we are, then the
-! corresponding symmetries are forced on the horizon.
-
- if (ahf_cartoon.ne.0) then
- cartoon = .true.
- else
- cartoon = .false.
- end if
-
- if (cartoon) then
- nonaxi = .false.
- refx = .true.
- refy = .true.
- end if
-
-! If desired, write values of symmetry flags to screen.
-
- if (veryver) then
- write(*,*)
- write(*,*) 'Symmetries used for horizon:'
- write(*,*) 'refx = ',refx
- write(*,*) 'refy = ',refy
- write(*,*) 'refz = ',refz
- write(*,*)
- end if
+ call AHFinder_SetReflections(CCTK_ARGUMENTS)
! Find the values of {stepx,stepy,stepz} depending
! on the symmetries.
@@ -1620,7 +1530,3 @@
! ***************
end subroutine AHFinder
-
-
-
-