SURFEX v7.3: /home/dasprezs/EXPORT_v7_3/src/SURFEX/get_grid_coord

Go to the documentation of this file.
00001 !     ###############################################
00002       SUBROUTINE GET_GRID_COORD_GAUSS(KGRID_PAR,KL,PGRID_PAR,PLON_XY,PLAT_XY)
00003 !     ###############################################
00004 !
00005 !!****  *GET_GRID_COORD_GAUSS* - gets pseudo longitude LON_XY and 
00006 !                                pseudo-latitude LAT_XY of the grid points
00007 !!
00008 !!    PURPOSE
00009 !!    -------
00010 !!
00011 !!**  METHOD
00012 !!    ------
00013 !!
00014 !!    EXTERNAL
00015 !!    --------
00016 !!
00017 !!
00018 !!    IMPLICIT ARGUMENTS
00019 !!    ------------------
00020 !!
00021 !!    REFERENCE
00022 !!    ---------
00023 !!
00024 !!
00025 !!    AUTHOR
00026 !!    ------
00027 !!
00028 !!    MODIFICATIONS
00029 !!    -------------
00030 !-------------------------------------------------------------------------------
00031 !
00032 !*       0.    DECLARATIONS
00033 !              ------------
00034 !
00035 USE MODE_GRIDTYPE_GAUSS
00036 !
00037 !
00038 USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
00039 USE PARKIND1  ,ONLY : JPRB
00040 !
00041 IMPLICIT NONE
00042 !
00043 !*       0.1   Declarations of arguments
00044 !              -------------------------
00045 !
00046 INTEGER,                    INTENT(IN)  :: KGRID_PAR  ! size of PGRID_PAR
00047 INTEGER,                    INTENT(IN)  :: KL         ! number of points
00048 REAL, DIMENSION(KGRID_PAR), INTENT(IN)  :: PGRID_PAR  ! parameters defining this grid
00049 REAL, DIMENSION(KL),        INTENT(OUT) :: PLON_XY    ! pseudo-latitude  (deg)
00050 REAL, DIMENSION(KL),        INTENT(OUT) :: PLAT_XY    ! pseudo-longitude (deg)
00051 !
00052 !*       0.2   Declarations of local variables
00053 !              -------------------------------
00054 !
00055 REAL    :: ZLAPO   ! latitude  of the rotated pole (deg)
00056 REAL    :: ZLOPO   ! longitude of the rotated pole (deg)
00057 REAL    :: ZCODIL  ! stretching factor (must be greater than or equal to 1)
00058 REAL(KIND=JPRB) :: ZHOOK_HANDLE
00059 !
00060 !---------------------------------------------------------------------------
00061 !
00062 !*       1.    Gets latitudes and longitudes
00063 !              -----------------------------
00064 !
00065 IF (LHOOK) CALL DR_HOOK('GET_GRID_COORD_GAUSS',0,ZHOOK_HANDLE)
00066  CALL GET_GRIDTYPE_GAUSS(PGRID_PAR,PLAPO=ZLAPO,PLOPO=ZLOPO,PCODIL=ZCODIL, &
00067                           PLAT_XY=PLAT_XY,PLON_XY=PLON_XY)  
00068 IF (LHOOK) CALL DR_HOOK('GET_GRID_COORD_GAUSS',1,ZHOOK_HANDLE)
00069 !
00070 !---------------------------------------------------------------------------
00071 !
00072 END SUBROUTINE GET_GRID_COORD_GAUSS