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(UG,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_surf_atm_grid_n, ONLY : surf_atm_grid_t
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
!            ------------------------
!
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
!
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 (ug%G%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 (ug%G%CGRID=="CONF PROJ ") THEN
      CALL get_mesh_index_conf_proj(isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy)  
    ENDIF
    IF (ug%G%CGRID=="LONLAT REG") THEN
      IF (PRESENT(pvalue) .AND. PRESENT(pnodata)) THEN
        CALL get_mesh_index_lonlat_reg(isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy, &
                                     pvalue,pnodata)
      ELSE            
        CALL get_mesh_index_lonlat_reg(isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy)  
      ENDIF
    ENDIF
    IF (ug%G%CGRID=="GAUSS     ") THEN
      IF (PRESENT(pvalue) .AND. PRESENT(pnodata)) THEN
        CALL get_mesh_index_gauss(knblines,isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy, &
                                     pvalue,pnodata)
      ELSE
        CALL get_mesh_index_gauss(knblines,isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy)
      ENDIF              
    ENDIF
    IF (ug%G%CGRID=="IGN       ") THEN
      IF (PRESENT(pvalue) .AND. PRESENT(pnodata)) THEN
        CALL get_mesh_index_ign(isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy, &
                                     pvalue,pnodata)
      ELSE       
        CALL get_mesh_index_ign(isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy)
      ENDIF  
    ENDIF
    IF (ug%G%CGRID=="LONLATVAL ") &
      CALL get_mesh_index_lonlatval(isso,ug%XGRID_FULL_PAR,plat,plon,kindex,iissox,iissoy)  
    IF (ug%G%CGRID=="LONLAT ROT") THEN
      CALL get_mesh_index_lonlat_rot(SIZE(plat),ug%XGRID_FULL_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 = ', ug%G%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