read_nam_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_nam_gridtype(HPROGRAM,HGRID,KGRID_PAR,PGRID_PAR,KL)
!     ##########################################################
!!
!!    PURPOSE
!!    -------
!!   Reads in namelist the grid type and parameters.
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson                   Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original     01/2004
!!    E. Martin    10/2007  IGN Grids
!!    P. Samuelsson SMHI 12/2012  Rotated lonlat
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
USE modi_read_nam_grid_cartesian
!
USE modi_read_nam_grid_conf_proj
!
USE modi_read_nam_grid_gauss
!
USE modi_read_nam_grid_ign
!
USE modi_read_nam_grid_lonlat_reg
!
USE modi_read_nam_grid_lonlatval
!
USE modi_read_nam_grid_lonlat_rot
IMPLICIT NONE
!
!*    0.1    Declaration of dummy arguments
!            ------------------------------
!
 CHARACTER(LEN=6),  INTENT(IN)   :: hprogram   ! program calling the surface
 CHARACTER(LEN=10), INTENT(IN)   :: hgrid      ! grid type
INTEGER,           INTENT(OUT)  :: kgrid_par  ! size of PGRID_PAR
REAL, DIMENSION(:), POINTER     :: pgrid_par  ! parameters defining this grid
INTEGER,           INTENT(OUT)  :: kl         ! number of points
REAL(KIND=JPRB) :: zhook_handle
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
!------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('READ_NAM_GRIDTYPE',0,zhook_handle)
SELECT CASE (hgrid)

!*    1.      Conformal projection grid and regular lat/lon
!             ---------------------------------------------
!
  CASE ('CONF PROJ ','LONLAT REG','CARTESIAN ','GAUSS     ','IGN       ','LONLATVAL ','LONLAT ROT')
    kgrid_par = 0
    ALLOCATE(pgrid_par(0))
    IF (hgrid=='CONF PROJ ')&
      CALL read_nam_grid_conf_proj(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='CARTESIAN ')&
      CALL read_nam_grid_cartesian(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='LONLAT REG')&
      CALL read_nam_grid_lonlat_reg(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='GAUSS     ')&
      CALL read_nam_grid_gauss(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='IGN       ')&
      CALL read_nam_grid_ign(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='LONLATVAL ')&
      CALL read_nam_grid_lonlatval(hprogram,kgrid_par,kl,pgrid_par)      
    IF (hgrid=='LONLAT ROT')&
      CALL read_nam_grid_lonlat_rot(hprogram,kgrid_par,kl,pgrid_par)  
    DEALLOCATE(pgrid_par)
    ALLOCATE(pgrid_par(kgrid_par))
    IF (hgrid=='CONF PROJ ')&
      CALL read_nam_grid_conf_proj(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='CARTESIAN ')&
      CALL read_nam_grid_cartesian(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='LONLAT REG')&
      CALL read_nam_grid_lonlat_reg(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='GAUSS     ')&
      CALL read_nam_grid_gauss(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='IGN       ')&
      CALL read_nam_grid_ign(hprogram,kgrid_par,kl,pgrid_par)  
    IF (hgrid=='LONLATVAL ')&
      CALL read_nam_grid_lonlatval(hprogram,kgrid_par,kl,pgrid_par)   
    IF (hgrid=='LONLAT ROT')&
      CALL read_nam_grid_lonlat_rot(hprogram,kgrid_par,kl,pgrid_par)  
        ! 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.      Other cases
!             -----------
!
  CASE default
    CALL abor1_sfx('READ_NAM_GRIDTYPE: GRID TYPE NOT SUPPORTED, '//hgrid)

END SELECT
IF (lhook) CALL dr_hook('READ_NAM_GRIDTYPE',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_nam_gridtype