get_mesh_corner.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_corner (UG, &
                                  KLUOUT,PCORNER_LAT,PCORNER_LON)
!     ##############################################################
!
!!**** *GET_MESH_CORNER* get the grid cell corner for each (lat,lon)
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    B. Decharme         Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    10/2013
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
!
USE modi_get_mesh_corner_cartesian
USE modi_get_mesh_corner_conf_proj
USE modi_get_mesh_corner_gauss
USE modi_get_mesh_corner_ign
USE modi_get_mesh_corner_lonlat_reg
USE modi_get_mesh_corner_lonlatval
!
USE modi_abor1_sfx
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
!
INTEGER,                         INTENT(IN)    :: KLUOUT      ! output listing
REAL,    DIMENSION(:,:),         INTENT(OUT)   :: PCORNER_LAT ! Grid corner Latitude
REAL,    DIMENSION(:,:),         INTENT(OUT)   :: PCORNER_LON ! Grid corner Longitude
!
!*    0.2    Declaration of other local variables
!            ------------------------------------
!
INTEGER :: IL, IC
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!----------------------------------------------------------------------------
!
!*    1.     Get position
!            ------------
!
IF (lhook) CALL dr_hook('GET_MESH_CORNER',0,zhook_handle)
!
il = SIZE(pcorner_lat,1)
ic = SIZE(pcorner_lat,2)
!
SELECT CASE (ug%G%CGRID)
!     
  CASE("LONLAT REG")
      CALL get_mesh_corner_lonlat_reg(ug%G%NGRID_PAR,il,ic,ug%G%XGRID_PAR,pcorner_lat,pcorner_lon) 
  CASE("CARTESIAN")
      CALL get_mesh_corner_cartesian(ug%G%NGRID_PAR,il,ic,ug%G%XGRID_PAR,pcorner_lat,pcorner_lon)        
  CASE("CONF PROJ")
      CALL get_mesh_corner_conf_proj(ug%G%NGRID_PAR,il,ic,ug%G%XGRID_PAR,pcorner_lat,pcorner_lon)      
  CASE("GAUSS     ")
      CALL get_mesh_corner_gauss(ug%G%NGRID_PAR,il,ic,ug%G%XGRID_PAR,pcorner_lat,pcorner_lon)
  CASE("IGN       ")
      CALL get_mesh_corner_ign(ug%G%NGRID_PAR,il,ic,ug%G%XGRID_PAR,pcorner_lat,pcorner_lon)  
  CASE("LONLATVAL ")
      CALL get_mesh_corner_lonlatval(ug%G%NGRID_PAR,il,ic,ug%G%XGRID_PAR,pcorner_lat,pcorner_lon)  
  CASE DEFAULT
    WRITE(kluout,*) 'error in grid cell corner computations (routine GET_MESH_CORNER)'
    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_CORNER: IMPOSSIBLE TO CALCULATE GRID CELL CORNER')
!    
END SELECT
!
IF (lhook) CALL dr_hook('GET_MESH_CORNER',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE get_mesh_corner