get_grid_coord_gauss.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_grid_coord_gauss(KGRID_PAR,KL,PGRID_PAR,PLON_XY,PLAT_XY)
!     ###############################################
!
!!****  *GET_GRID_COORD_GAUSS* - gets pseudo longitude LON_XY and 
!                                pseudo-latitude LAT_XY of the grid points
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!
!!    MODIFICATIONS
!!    -------------
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE mode_gridtype_gauss
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
INTEGER,                    INTENT(IN)  :: KGRID_PAR  ! size of PGRID_PAR
INTEGER,                    INTENT(IN)  :: KL         ! number of points
REAL, DIMENSION(KGRID_PAR), INTENT(IN)  :: PGRID_PAR  ! parameters defining this grid
REAL, DIMENSION(KL),        INTENT(OUT) :: PLON_XY    ! pseudo-latitude  (deg)
REAL, DIMENSION(KL),        INTENT(OUT) :: PLAT_XY    ! pseudo-longitude (deg)
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL    :: ZLAPO   ! latitude  of the rotated pole (deg)
REAL    :: ZLOPO   ! longitude of the rotated pole (deg)
REAL    :: ZCODIL  ! stretching factor (must be greater than or equal to 1)
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!---------------------------------------------------------------------------
!
!*       1.    Gets latitudes and longitudes
!              -----------------------------
!
IF (lhook) CALL dr_hook('GET_GRID_COORD_GAUSS',0,zhook_handle)
 CALL get_gridtype_gauss(pgrid_par,plapo=zlapo,plopo=zlopo,pcodil=zcodil, &
                          plat_xy=plat_xy,plon_xy=plon_xy)  
IF (lhook) CALL dr_hook('GET_GRID_COORD_GAUSS',1,zhook_handle)
!
!---------------------------------------------------------------------------
!
END SUBROUTINE get_grid_coord_gauss