prep_grid_extern.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 prep_grid_extern (&
                                   hfiletype,kluout,hgridtype,hinterp_type,kni)
!     ##########################################################################
!
!!****  *PREP_GRID_EXTERN* - reads EXTERNALIZED Surface grid.
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!
!!      V. Masson
!!
!!    MODIFICATIONS
!!    -------------
!!      Original   06/2003
!-------------------------------------------------------------------------------
!
!*      0. DECLARATIONS
!          ------------
!
!
!
!
USE modi_read_surf
USE modi_prep_grid_conf_proj
USE modi_prep_grid_cartesian
USE modi_prep_grid_gauss
USE modi_prep_grid_lonlat_reg
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
IMPLICIT NONE
!
!* 0.1. Declaration of arguments
!       ------------------------
!
!
!
 CHARACTER(LEN=6),  INTENT(IN)    :: hfiletype    ! file type
INTEGER,           INTENT(IN)    :: kluout       ! logical unit of output listing
 CHARACTER(LEN=10),  INTENT(OUT)  :: hgridtype    ! Grid type
 CHARACTER(LEN=6),  INTENT(OUT)   :: hinterp_type ! Grid type
INTEGER,           INTENT(OUT)   :: kni          ! number of points
!
!* 0.2 Declaration of local variables
!      ------------------------------
!
INTEGER :: iresp
REAL(KIND=JPRB) :: zhook_handle
!
!-----------------------------------------------------------------------
!
!*   1 Type of grid
!      ------------
!
IF (lhook) CALL dr_hook('PREP_GRID_EXTERN',0,zhook_handle)
 CALL read_surf(&
                hfiletype,'GRID_TYPE',hgridtype,iresp)
!
!-----------------------------------------------------------------------
!
!*   2 Reading of grid
!      ---------------
!
IF (hgridtype=='CONF PROJ ') THEN
  CALL prep_grid_conf_proj(&
                           hfiletype,hinterp_type,kni)
ELSE IF (hgridtype=='CARTESIAN ') THEN
  CALL prep_grid_cartesian(&
                           hfiletype,hinterp_type,kni)
ELSE IF (hgridtype=='GAUSS     ') THEN
  CALL prep_grid_gauss(&
                       hfiletype,hinterp_type,kni)
ELSE IF (hgridtype=='LONLAT REG') THEN
  hgridtype = 'LATLON    '
  CALL prep_grid_lonlat_reg(&
                            hfiletype,hinterp_type,kni)
ELSE
  WRITE(kluout,*) 'GRIDTYPE "',hgridtype,'" NOT ACCEPTED AS INPUT FILE FOR FIELD PREPARATION'
  CALL abor1_sfx('GRIDTYPE NOT ACCEPTED AS INPUT FILE FOR FIELD PREPARATION, '//hgridtype)
END IF
IF (lhook) CALL dr_hook('PREP_GRID_EXTERN',1,zhook_handle)
!
!-----------------------------------------------------------------------
!
END SUBROUTINE prep_grid_extern