diff options
Diffstat (limited to 'src/AHFinder.F')
-rw-r--r-- | src/AHFinder.F | 96 |
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 - - - - |