get_near_meshes_lonlatval.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
      SUBROUTINE get_near_meshes_lonlatval(KGRID_PAR,KL,PGRID_PAR,KNEAR_NBR,KNEAR)
!     ##############################################################
!
!!**** *GET_NEAR_MESHES_LONLATVAL* get the near grid mesh indices
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson         Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    03/2004
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE mode_gridtype_lonlatval
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
INTEGER,                         INTENT(IN)    :: kgrid_par ! size of PGRID_PAR
INTEGER,                         INTENT(IN)    :: kl        ! number of points
INTEGER,                         INTENT(IN)    :: knear_nbr ! number of nearest points wanted
REAL,    DIMENSION(KGRID_PAR),   INTENT(IN)    :: pgrid_par ! grid parameters
INTEGER, DIMENSION(:,:),POINTER  :: knear     ! near mesh indices
!
!*    0.2    Declaration of other local variables
!            ------------------------------------
!
REAL, DIMENSION(KL) :: zdis
REAL,DIMENSION(KL)    :: zx
REAL,DIMENSION(KL)    :: zy
REAL,DIMENSION(KL)    :: zdx
REAL,DIMENSION(KL)    :: zdy
REAL :: zdmax
INTEGER :: id0
INTEGER :: jp1, jp2, jn, il
REAL(KIND=JPRB) :: zhook_handle
!
!----------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('GET_NEAR_MESHES_LONLATVAL',0,zhook_handle)
!
 CALL get_gridtype_lonlatval(pgrid_par,il,zx,zy,zdx,zdy)
!
knear(:,:) = 0
!
! calcul de la distance de tous les points 2 à 2
!
zdis = 1.e20
!
DO jp1=1,kl
  !
  DO jp2=1,kl
    zdis(jp2) = sqrt((zx(jp1)-zx(jp2))**2+(zy(jp1)-zy(jp2))**2)
  ENDDO
  zdmax = maxval(zdis(:)) + 1.
  zdis(jp1) = zdmax
  !
  ! on prend les knear_nbr premiers, pour chaque
  !
  DO jn=1,min(kl-1,knear_nbr)
    !
    id0 = maxval(minloc(zdis(:)))       
    !
    knear(jp1,jn) = id0
    zdis(id0) = zdmax
    !
  ENDDO
  !
ENDDO
!
IF (lhook) CALL dr_hook('GET_NEAR_MESHES_LONLATVAL',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE get_near_meshes_lonlatval