get_mesh_index.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_mesh_index(KLUOUT,KNBLINES,PLAT,PLON,KINDEX,PVALUE,PNODATA,KSSO,KISSOX,KISSOY)
!     ##############################################################
!
!!**** *GET_MESH_INDEX* get the grid mesh where point (lat,lon) is located
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson         Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    12/09/95
!!    P. Samuelsson  SMHI  10/2014   Rotated lonlat
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_pgd_grid, ONLY : cgrid, xgrid_par, ngrid_par
USE modd_point_overlay
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
USE modi_get_mesh_index_conf_proj
!
USE modi_get_mesh_index_gauss
!
USE modi_get_mesh_index_ign
!
USE modi_get_mesh_index_lonlat_reg
!
USE modi_get_mesh_index_lonlatval
!
USE modi_get_mesh_index_lonlat_rot
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
INTEGER,                         INTENT(IN)    :: kluout  ! output listing
INTEGER,                         INTENT(IN)    :: knblines
REAL,    DIMENSION(:),           INTENT(IN)    :: plat    ! latitude of the point
REAL,    DIMENSION(:),           INTENT(IN)    :: plon    ! longitude of the point
INTEGER, DIMENSION(:,:),         INTENT(OUT)   :: kindex  ! index of the grid mesh where the point is
!
REAL, DIMENSION(:), OPTIONAL,     INTENT(IN)   :: pvalue  ! value of the point to add
REAL, OPTIONAL,                   INTENT(IN)   :: pnodata
!
INTEGER,               OPTIONAL, INTENT(IN)    :: ksso    ! number of subgrid mesh in each direction
INTEGER, DIMENSION(:,:), OPTIONAL, INTENT(OUT) :: kissox  ! X index of the subgrid mesh where the point is
INTEGER, DIMENSION(:,:), OPTIONAL, INTENT(OUT) :: kissoy  ! Y index of the subgrid mesh where the point is
!
!*    0.2    Declaration of other local variables
!            ------------------------------------
!
INTEGER                        :: isso
INTEGER, DIMENSION(NOVMX,SIZE(PLAT)) :: iissox
INTEGER, DIMENSION(NOVMX,SIZE(PLAT)) :: iissoy
REAL(KIND=JPRB) :: zhook_handle
!----------------------------------------------------------------------------
!
!*    1.     Get position
!            ------------
!
IF (lhook) CALL dr_hook('GET_MESH_INDEX',0,zhook_handle)
SELECT CASE (cgrid)
!     
  CASE("CONF PROJ ","LONLAT REG","GAUSS     ","IGN      ","LONLATVAL ","LONLAT ROT")
    IF (present(ksso) .AND. present(kissox) .AND. present(kissoy)) THEN
      isso = ksso
    ELSE
      isso = 0
    ENDIF
    !
    IF (cgrid=="CONF PROJ ") THEN
      CALL get_mesh_index_conf_proj(ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy)  
    ENDIF
    IF (cgrid=="LONLAT REG") THEN
      IF (present(pvalue) .AND. present(pnodata)) THEN
        CALL get_mesh_index_lonlat_reg(ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy, &
                                     pvalue,pnodata)
      ELSE            
        CALL get_mesh_index_lonlat_reg(ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy)  
      ENDIF
    ENDIF
    IF (cgrid=="GAUSS     ") THEN
      IF (present(pvalue) .AND. present(pnodata)) THEN
        CALL get_mesh_index_gauss(knblines,ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy, &
                                     pvalue,pnodata)
      ELSE
        CALL get_mesh_index_gauss(knblines,ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy)
      ENDIF              
    ENDIF
    IF (cgrid=="IGN       ") THEN
      IF (present(pvalue) .AND. present(pnodata)) THEN
        CALL get_mesh_index_ign(ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy, &
                                     pvalue,pnodata)
      ELSE       
        CALL get_mesh_index_ign(ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy)
      ENDIF  
    ENDIF
    IF (cgrid=="LONLATVAL ") &
      CALL get_mesh_index_lonlatval(ngrid_par,isso,xgrid_par,plat,plon,kindex,iissox,iissoy)  
    IF (cgrid=="LONLAT ROT") THEN
      CALL get_mesh_index_lonlat_rot(ngrid_par,SIZE(plat),xgrid_par,plat,plon,kindex,isso,iissox,iissoy)  
    ENDIF
    !
    IF (present(ksso) .AND. present(kissox) .AND. present(kissoy)) THEN
      kissox = iissox
      kissoy = iissoy
    ENDIF

  CASE default
    WRITE(kluout,*) 'error in physiographic fields computations (routine GET_MESH_INDEX)'
    WRITE(kluout,*) 'It is impossible to retrieve geographical coordinates (latitude, longitude)'
    WRITE(kluout,*) 'for the following grid type: CGRID = ', cgrid
    CALL abor1_sfx('GET_MESH_INDEX: IMPOSSIBLE TO RETRIEVE GEOGRAPHICAL COORDINATES')
END SELECT
IF (lhook) CALL dr_hook('GET_MESH_INDEX',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE get_mesh_index