get_near_meshes_conf_proj.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_conf_proj(KGRID_PAR,KL,PGRID_PAR,KNEAR_NBR,KNEAR)
!     ##############################################################
!
!!**** *GET_NEAR_MESHES_CONF_PROJ* get the near grid mesh indices
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson         Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    03/2004
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE mode_gridtype_conf_proj
!
!
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
!            ------------------------------------
!
INTEGER                            :: iimax, ijmax
INTEGER                            :: ji, jj
INTEGER                            :: jx, jy
INTEGER                            :: jl
INTEGER                            :: idist
INTEGER                            :: icount
REAL(KIND=JPRB) :: zhook_handle
!----------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('GET_NEAR_MESHES_CONF_PROJ',0,zhook_handle)
 CALL get_gridtype_conf_proj(pgrid_par,kimax=iimax,kjmax=ijmax)
!
knear(:,:) = 0
!
idist = int(sqrt(float(knear_nbr)))
!
IF (iimax*ijmax==kl) THEN
  DO jj=1,ijmax
    DO ji=1,iimax
      icount = 0
      jl = ji + iimax * (jj-1)
      knear(jl,:) = 0      
      DO jx=-(idist-1)/2,idist/2
        DO jy=-(idist-1)/2,idist/2
          IF (ji+jx>0 .AND. ji+jx<iimax+1 .AND. jj+jy>0 .AND. jj+jy<ijmax+1) THEN
            icount = icount + 1
            knear(jl,icount) = (ji+jx) + iimax * (jj+jy-1)
          END IF
        END DO
      END DO
    END DO
  END DO
END IF
IF (lhook) CALL dr_hook('GET_NEAR_MESHES_CONF_PROJ',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE get_near_meshes_conf_proj