latlon_grid.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 latlon_grid(HGRID,KGRID_PAR,KL,KLUOUT,PGRID_PAR,PLAT,PLON,PMESH_SIZE,PDIR)
!     #########################################################################
!
!!****  *LATLON_GRID* - routine to compute the horizontal geographic fields
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004 
!!                  10/2007 (E. Martin) IGN grids
!!                  12/2012 (P. Samuelsson SMHI) Rotated lonlat
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
USE modi_latlon_gridtype_cartesian
!
USE modi_latlon_gridtype_conf_proj
!
USE modi_latlon_gridtype_gauss
!
USE modi_latlon_gridtype_ign
!
USE modi_latlon_gridtype_lonlat_reg
!
USE modi_latlon_gridtype_lonlatval
!
USE modi_latlon_gridtype_lonlat_rot
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
 CHARACTER(LEN=10),  INTENT(IN)  :: hgrid      ! grid type
INTEGER,            INTENT(IN)  :: kgrid_par  ! size of PGRID_PAR
INTEGER,            INTENT(IN)  :: kl         ! number of points
INTEGER,            INTENT(IN)  :: kluout     ! output listing logical unit
REAL, DIMENSION(:), INTENT(IN)  :: pgrid_par  ! parameters defining this grid
REAL, DIMENSION(:), INTENT(OUT) :: plat       ! latitude  (degrees)
REAL, DIMENSION(:), INTENT(OUT) :: plon       ! longitude (degrees)
REAL, DIMENSION(:), INTENT(OUT) :: pmesh_size ! mesh size (m2)
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL :: pdir ! direction of main grid Y axis (deg. from N, clockwise)
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL, DIMENSION(KL) :: zdir
REAL(KIND=JPRB) :: zhook_handle
!---------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('LATLON_GRID',0,zhook_handle)
SELECT CASE (hgrid)
!
!*    1.      Conformal projection grid
!             -------------------------
!
  CASE ('CONF PROJ ')
    CALL latlon_gridtype_conf_proj(kgrid_par,kl,pgrid_par,plat,plon,pmesh_size,zdir)
    IF (present(pdir)) pdir = zdir
    ! note that all points of the grid will be kept, whatever the surface
    ! type under consideration (e.g. sea points will be kept even for
    ! initialization of continents)
    !

!*    2.      latitude/longitude grid
!             -----------------------
!
  CASE ('LONLAT REG')
    CALL latlon_gridtype_lonlat_reg(kgrid_par,kl,pgrid_par,plat,plon,pmesh_size,zdir)
    IF (present(pdir)) pdir = zdir

!
!*    3.      Cartesian grid
!             --------------
!
  CASE ('CARTESIAN ')
    CALL latlon_gridtype_cartesian(kgrid_par,kl,pgrid_par,plat,plon,pmesh_size,zdir)
    IF (present(pdir)) pdir = zdir
    ! note that all points of the grid will be kept, whatever the surface
    ! type under consideration (e.g. sea points will be kept even for
    ! initialization of continents)
    !
!*    4.      gaussian grid
!             -------------
!
  CASE ('GAUSS     ')
    CALL latlon_gridtype_gauss(kgrid_par,kl,pgrid_par,plat,plon,pmesh_size,zdir)
    IF (present(pdir)) pdir = zdir
!
!*    5.      IGN grid
!             --------
!
  CASE ('IGN       ')
    CALL latlon_gridtype_ign(kgrid_par,kl,pgrid_par,plat,plon,pmesh_size,zdir)
    IF (present(pdir)) pdir = zdir
!
!*    6.      lonlatval grid
!             --------
!
  CASE ('LONLATVAL ')
    CALL latlon_gridtype_lonlatval(kgrid_par,kl,pgrid_par,plat,plon,pmesh_size,zdir)
    IF (present(pdir)) pdir = zdir
!
!*    7.      Rotated lonlat grid
!             -------------------
!
  CASE ('LONLAT ROT')
    CALL latlon_gridtype_lonlat_rot(kgrid_par,kl,pgrid_par,plat,plon,pmesh_size,zdir)
    IF (present(pdir)) pdir = zdir
!
!
  CASE default
    CALL abor1_sfx('LATLON_GRID: GRID TYPE NOT SUPPORTED '//hgrid)

END SELECT
IF (lhook) CALL dr_hook('LATLON_GRID',1,zhook_handle)
!
!---------------------------------------------------------------------------
!
END SUBROUTINE latlon_grid