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(G,KL,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 modd_sfx_grid_n, ONLY : grid_t
!
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
!              -------------------------
!
TYPE(grid_t), INTENT(INOUT) :: G
INTEGER,            INTENT(IN)  :: KL         ! number of points
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 (g%CGRID)
!
!*    1.      Conformal projection grid
!             -------------------------
!
  CASE ('CONF PROJ ')
    CALL latlon_gridtype_conf_proj(g,kl,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(g,kl,zdir)
    IF (PRESENT(pdir)) pdir = zdir

!
!*    3.      Cartesian grid
!             --------------
!
  CASE ('CARTESIAN ')
    CALL latlon_gridtype_cartesian(g,kl,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(g,kl,zdir)
    IF (PRESENT(pdir)) pdir = zdir
!
!*    5.      IGN grid
!             --------
!
  CASE ('IGN       ')
    CALL latlon_gridtype_ign(g,kl,zdir)
    IF (PRESENT(pdir)) pdir = zdir
!
!*    6.      lonlatval grid
!             --------
!
  CASE ('LONLATVAL ')
    CALL latlon_gridtype_lonlatval(g,kl,zdir)
    IF (PRESENT(pdir)) pdir = zdir
!
!*    7.      Rotated lonlat grid
!             -------------------
!
  CASE ('LONLAT ROT')
    CALL latlon_gridtype_lonlat_rot(g,kl,zdir)
    IF (PRESENT(pdir)) pdir = zdir
!
!
  CASE DEFAULT
    CALL abor1_sfx('LATLON_GRID: GRID TYPE NOT SUPPORTED '//g%CGRID)

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