subroutine hunt(xx ,n ,x ,jlo )
!----- GPL ---------------------------------------------------------------------
!
! Copyright (C) Stichting Deltares, 2011-2020.
!
! This program is free software: you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by
! the Free Software Foundation version 3.
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with this program. If not, see .
!
! contact: delft3d.support@deltares.nl
! Stichting Deltares
! P.O. Box 177
! 2600 MH Delft, The Netherlands
!
! All indications and logos of, and references to, "Delft3D" and "Deltares"
! are registered trademarks of Stichting Deltares, and remain the property of
! Stichting Deltares. All rights reserved.
!
!-------------------------------------------------------------------------------
! $Id$
! $HeadURL$
!!--description-----------------------------------------------------------------
! NONE
!!--pseudo code and references--------------------------------------------------
! NONE
!!--declarations----------------------------------------------------------------
use precision
implicit none
!
! Global variables
!
integer :: jlo
integer , intent(in) :: n
real(hp), intent(in) :: x
real(hp), dimension(n), intent(in) :: xx
!
!
! Local variables
!
integer :: inc
integer :: jhi
integer :: jm
logical :: ascnd
!
!
!! executable statements -------------------------------------------------------
!
!
!
!
ascnd = xx(n)>=xx(1)
if (jlo<=0 .or. jlo>n) then
jlo = 0
jhi = n + 1
goto 3
endif
inc = 1
if (x>=xx(jlo) .eqv. ascnd) then
1 continue
jhi = jlo + inc
if (jhi>n) then
jhi = n + 1
elseif (x>=xx(jhi) .eqv. ascnd) then
jlo = jhi
inc = inc + inc
goto 1
else
endif
else
jhi = jlo
2 continue
jlo = jhi - inc
if (jlo<1) then
jlo = 0
elseif (xxx(jm) .eqv. ascnd) then
jlo = jm
else
jhi = jm
endif
goto 3
end subroutine hunt