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

Go to the documentation of this file.
00001 !     #########
00002       SUBROUTINE GET_MESH_DIM_IGN(KGRID_PAR,KL,PGRID_PAR,PDX,PDY)
00003 !     ##############################################################
00004 !
00005 !!**** *GET_MESH_DIM_IGN* get the grid mesh dimensions
00006 !!
00007 !!    PURPOSE
00008 !!    -------
00009 !!
00010 !!    METHOD
00011 !!    ------
00012 !!   
00013 !!    REFERENCE
00014 !!    ---------
00015 !!
00016 !!    AUTHOR
00017 !!    ------
00018 !!
00019 !!    E. Martin         Meteo-France
00020 !!
00021 !!    MODIFICATION
00022 !!    ------------
00023 !!
00024 !!    Original    10/2007
00025 !!
00026 !----------------------------------------------------------------------------
00027 !
00028 !*    0.     DECLARATION
00029 !            -----------
00030 !
00031 USE MODE_GRIDTYPE_IGN
00032 !
00033 !
00034 USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
00035 USE PARKIND1  ,ONLY : JPRB
00036 !
00037 IMPLICIT NONE
00038 !
00039 !*    0.1    Declaration of arguments
00040 !            ------------------------
00041 !
00042 INTEGER,                         INTENT(IN)    :: KGRID_PAR ! size of PGRID_PAR
00043 INTEGER,                         INTENT(IN)    :: KL        ! number of points
00044 REAL,    DIMENSION(KGRID_PAR),   INTENT(IN)    :: PGRID_PAR ! grid parameters
00045 REAL,    DIMENSION(KL),          INTENT(OUT)   :: PDX       ! dimension in x dir. (meters)
00046 REAL,    DIMENSION(KL),          INTENT(OUT)   :: PDY       ! dimension in y dir. (meters)
00047 !
00048 !*    0.2    Declaration of other local variables
00049 !            ------------------------------------
00050 !
00051 INTEGER             :: ILAMBERT ! Lambert type
00052 REAL, DIMENSION(KL) :: ZX       ! X Lambert   coordinate
00053 REAL, DIMENSION(KL) :: ZY       ! Y Lambert   coordinate
00054 REAL, DIMENSION(KL) :: ZLAT     ! latitude
00055 REAL, DIMENSION(KL) :: ZLON     ! longitude
00056 REAL, DIMENSION(KL) :: ZMAP     ! map factor
00057 REAL(KIND=JPRB) :: ZHOOK_HANDLE
00058 
00059 !----------------------------------------------------------------------------
00060 !
00061 IF (LHOOK) CALL DR_HOOK('GET_MESH_DIM_IGN',0,ZHOOK_HANDLE)
00062  CALL GET_GRIDTYPE_IGN(PGRID_PAR,KLAMBERT=ILAMBERT, &
00063                               PX=ZX,PY=ZY,PDX=PDX,PDY=PDY      )  
00064 !
00065 !---------------------------------------------------------------------------
00066 !
00067 !*       2.    Computation of latitude and longitude
00068 !              -------------------------------------
00069 !
00070  CALL LATLON_IGN(ILAMBERT,ZX,ZY,ZLAT,ZLON)
00071 !
00072 !-----------------------------------------------------------------------------
00073 !
00074 !*       3.    Compute grif size
00075 !              -----------------
00076 !
00077  CALL MAP_FACTOR_IGN(ILAMBERT,ZX,ZY,ZMAP)
00078 !
00079 PDX(:) = PDX(:) / ZMAP(:)
00080 PDY(:) = PDY(:) / ZMAP(:)
00081 IF (LHOOK) CALL DR_HOOK('GET_MESH_DIM_IGN',1,ZHOOK_HANDLE)
00082 !
00083 !-------------------------------------------------------------------------------
00084 !
00085 END SUBROUTINE GET_MESH_DIM_IGN