!!----------------------------------------------------------------------
   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
   !! $Id: obs_level_search.h90 2287 2010-10-18 07:53:52Z smasson $
   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
   !!----------------------------------------------------------------------

   SUBROUTINE obs_level_search( kgrd, pgrddep, kobs, pobsdep, kobsk )
      !!----------------------------------------------------------------------
      !!                    ***  ROUTINE obs_level_search ***
      !!
      !! ** Purpose : Search levels to find matching level to observed depth
      !!
      !! ** Method  : Straightforward search
      !!
      !! ** Action  : 
      !!
      !! History :
      !!        !  2001-11  (N. Daget, A. Weaver)
      !!        !  2006-03  (A. Weaver) NEMOVAR migration.
      !!        !  2006-05  (K. Mogensen) Moved to to separate routine.
      !!        !  2006-10  (A. Weaver) Cleanup
      !!        !  2008-10  (K. Mogensen) Remove assumptions on grid.
      !!----------------------------------------------------------------------

      !! * Arguments
      INTEGER, INTENT(IN) :: kgrd     ! Number of gridpoints
      REAL(KIND=wp), DIMENSION(kgrd), INTENT(INOUT) :: &
         &   pgrddep  ! Depths of gridpoints
      INTEGER, INTENT(IN) :: &
         &   kobs     ! Number of observations
      REAL(KIND=wp), DIMENSION(kobs), INTENT(INOUT) :: &
         &   pobsdep  ! Depths of observations
      INTEGER ,DIMENSION(kobs), INTENT(OUT) :: &
         &   kobsk    ! Level indices of observations
  
      !! * Local declarations
      INTEGER :: ji
      INTEGER :: jk

      !------------------------------------------------------------------------
      ! Search levels for each observations to find matching level
      !------------------------------------------------------------------------
      DO ji = 1, kobs 
         kobsk(ji) = 1
         depk: DO jk = 2, kgrd
            IF ( pgrddep(jk) >= pobsdep(ji) ) EXIT depk
         END DO depk
         kobsk(ji) = jk
      END DO

   END SUBROUTINE obs_level_search