read_gridtype.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 read_gridtype (&
                                hprogram,hgrid,kgrid_par,klu,oread,pgrid_par,kresp,hdir)
!     #########################################
!
!!****  *READ_GRID* - routine to initialise the horizontal grid of a scheme
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2003 
!!                  10/2007  E. Martin  IGN Grids
!!                  12/2012  P. Samuelsson  Rotated lonlat
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_read_gridtype_cartesian
!
USE modi_read_gridtype_conf_proj
!
USE modi_read_gridtype_gauss
!
USE modi_read_gridtype_ign
!
USE modi_read_gridtype_lonlat_reg
!
USE modi_read_gridtype_lonlatval
!
USE modi_read_gridtype_lonlat_rot
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
!
 CHARACTER(LEN=6),             INTENT(IN)    :: hprogram   ! calling program
 CHARACTER(LEN=10),            INTENT(IN)    :: hgrid      ! type of horizontal grid
INTEGER,                      INTENT(INOUT) :: kgrid_par  ! size of PGRID_PAR
INTEGER,                      INTENT(IN)    :: klu        ! number of points
LOGICAL,                      INTENT(IN)    :: oread      ! flag to read the grid
REAL, DIMENSION(:), OPTIONAL, INTENT(OUT)   :: pgrid_par  ! parameters defining this grid
INTEGER,            OPTIONAL, INTENT(OUT)   :: kresp      ! error return code
 CHARACTER(LEN=1),   OPTIONAL, INTENT(IN)    :: hdir       ! reading directive ('A': all field;'H': this proc)
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL, DIMENSION(:), ALLOCATABLE :: zgrid_par
INTEGER                         :: igrid_par
INTEGER                         :: iresp
 CHARACTER(LEN=1)                :: ydir
REAL(KIND=JPRB) :: zhook_handle
!
!---------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('READ_GRIDTYPE',0,zhook_handle)
IF (oread) THEN
  igrid_par = kgrid_par
ELSE
  igrid_par = 0
END IF
!
IF (present(hdir)) THEN
  ydir = hdir
ELSE
  ydir = 'H'
END IF
!
ALLOCATE(zgrid_par(igrid_par))
!
SELECT CASE (hgrid)
  CASE("NONE      ")
    DEALLOCATE(zgrid_par)
    IF (lhook) CALL dr_hook('READ_GRIDTYPE',1,zhook_handle)
    RETURN

  CASE ("CONF PROJ ")
    CALL read_gridtype_conf_proj(&
                                 hprogram,kgrid_par,klu,oread,igrid_par,zgrid_par,iresp,ydir)

  CASE ("CARTESIAN ")
    CALL read_gridtype_cartesian(&
                                 hprogram,kgrid_par,klu,oread,igrid_par,zgrid_par,iresp,ydir)

  CASE ("LONLAT REG")
    CALL read_gridtype_lonlat_reg(&
                                  hprogram,kgrid_par,klu,oread,igrid_par,zgrid_par,iresp,ydir)

  CASE ("GAUSS     ")
    CALL read_gridtype_gauss(&
                             hprogram,kgrid_par,klu,oread,igrid_par,zgrid_par,iresp,ydir)

  CASE ("IGN       ")
    CALL read_gridtype_ign(&
                           hprogram,kgrid_par,klu,oread,igrid_par,zgrid_par,iresp,ydir)

  CASE ("LONLATVAL ")
    CALL read_gridtype_lonlatval(&
                                 hprogram,kgrid_par,klu,oread,igrid_par,zgrid_par,iresp,ydir)

  CASE ("LONLAT ROT")
    CALL read_gridtype_lonlat_rot(&
                                  hprogram,kgrid_par,klu,oread,igrid_par,zgrid_par,iresp,ydir)

END SELECT
!
IF (oread) pgrid_par = zgrid_par
!
IF (present(kresp)) kresp = iresp
!
DEALLOCATE(zgrid_par)
IF (lhook) CALL dr_hook('READ_GRIDTYPE',1,zhook_handle)
!---------------------------------------------------------------------------
!
END SUBROUTINE read_gridtype