get_mesh_dim_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_mesh_dim_gauss(KGRID_PAR,KL,PGRID_PAR,PMESHSIZE,PDX,PDY)
!     ##############################################################
!
!!**** *GET_MESH_DIM_GAUSS* get the grid mesh dimensions
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson         Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    03/2004
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_csts, ONLY : xpi, xradius
!
USE mode_gridtype_gauss
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration 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 ! grid parameters
REAL,    DIMENSION(KL),          INTENT(IN)    :: pmeshsize ! mesh size (m2)
REAL,    DIMENSION(KL),          INTENT(OUT)   :: pdx       ! mean dimension in x dir. (meters)
REAL,    DIMENSION(KL),          INTENT(OUT)   :: pdy       ! mean dimension in y dir. (meters)
!
!*    0.2    Declaration of other local variables
!            ------------------------------------
!
INTEGER :: inlati    ! number of pseudo-latitudes
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)
INTEGER, DIMENSION(:), ALLOCATABLE :: inlopa ! number of pseudo-longitudes on each
                                             ! pseudo-latitude circle
REAL,  DIMENSION(KL) :: zlat  ! latitudes
REAL,  DIMENSION(KL) :: zlon  ! longitudes
REAL,  DIMENSION(KL) :: zxinf ! minimum pseudo longitude of the grid point (deg)
REAL,  DIMENSION(KL) :: zxsup ! maximum pseudo longitude of the grid point (deg)
REAL,  DIMENSION(KL) :: zyinf ! minimum pseudo latitude  of the grid point (deg)
REAL,  DIMENSION(KL) :: zysup ! maximum pseudo latitude  of the grid point (deg)
REAL,  DIMENSION(KL) :: zmap  ! map factor
REAL(KIND=JPRB) :: zhook_handle
!
!----------------------------------------------------------------------------
!
!*       1.    Gets grid definition
!              --------------------
!
IF (lhook) CALL dr_hook('GET_MESH_DIM_GAUSS',0,zhook_handle)
 CALL get_gridtype_gauss(pgrid_par,inlati)
!
ALLOCATE(inlopa(inlati))
 CALL get_gridtype_gauss(pgrid_par,inlati,zlapo,zlopo,zcodil,inlopa(:),plat=zlat,plon=zlon)
!
!---------------------------------------------------------------------------
!
!*       2.    Computation of pseudo latitude and pseudo longitude limits for each mesh
!              ------------------------------------------------------------------------
!
 CALL gauss_grid_limits(inlati,inlopa,zxinf,zxsup,zyinf,zysup)
!
!-----------------------------------------------------------------------------
!
!*       3.    Map factor
!              ----------
!
 CALL map_factor_gauss(zlapo,zlopo,zcodil,zlat,zlon,zmap)
!
!-----------------------------------------------------------------------------
!
!*       4.    Compute mean grid dimension in each direction
!              ---------------------------------------------
!
pdy(:) = xradius * xpi/180. * (zysup(:)-zyinf(:)) * zmap(:)
pdx(:) = pmeshsize(:) / pdy(:)
!
!-------------------------------------------------------------------------------
!
DEALLOCATE(inlopa)
IF (lhook) CALL dr_hook('GET_MESH_DIM_GAUSS',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE get_mesh_dim_gauss