| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051 |
- !!----------------------------------------------------------------------
- !! 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
|
