c/*@@
c  @file      AHFinder_find3.F
c  @date      March 1999
c  @author    Lars Nerger
c  @desc 
c             Muliplicate gridfunctions for final grid function
c  @enddesc 
c@@*/

c Note that including cactus.h will also include AHFinder.h
!#include "cactus.h"
#include "cctk.h"
#include "cctk_parameters.h"
#include "cctk_arguments.h"

!     To output single gridfunctions when searching for 3 horizons
!     the gridfunctions are here multiplicated.

      subroutine AHFinder_find3(CCTK_FARGUMENTS,mtype)

      use AHFinder_dat

      implicit none

      DECLARE_CCTK_FARGUMENTS
      DECLARE_CCTK_PARAMETERS

      integer mtype

      REAL one

!     Description of variables
!
!     mtype     Type of surface: see AHFinder.F


!     **********************
!     ***   DEFINE one   ***
!     **********************

      one  = 1.0D0


!     *************************
!     ***   START ROUTINE   ***
!     *************************

      if ((mfind.eq.0).and.(mtype.eq.1)) then
         ahfgrid3 = ahfgrid
         ahf_exp3 = ahf_exp
      else if ((mfind.eq.0).and.(mtype.ne.1)) then
         ahfgrid3 = one
         ahf_exp3 = one
      else if ((mfind.eq.1).and.(mtype.eq.1)) then
         ahfgrid3 = ahfgrid*ahfgrid3
         ahf_exp3 = ahf_exp*ahf_exp3
      else if ((mfind.eq.2).and.(mtype.eq.1)) then
         ahfgrid3 = ahfgrid*ahfgrid3
         ahf_exp3 = ahf_exp*ahf_exp3
      end if


!     ***************
!     ***   END   ***
!     ***************

      end subroutine AHFinder_find3