read_pgd_schemes.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_pgd_schemes(HPROGRAM,HNATURE,HSEA,HTOWN,HWATER)
!     ######################################
!!
!!    PURPOSE
!!    -------
!!   initializes the surface SCHEMES.
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson                   Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original     13/10/03
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modi_default_schemes
USE modi_test_nam_var_surf
USE modi_open_namelist
USE modi_close_namelist
USE modi_get_luout
!
USE modn_pgd_schemes
!
USE mode_pos_surf
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of dummy arguments
!            ------------------------------
!
 CHARACTER(LEN=6),  INTENT(IN)  :: hprogram ! program calling READ_PGD
 CHARACTER(LEN=6),  INTENT(OUT) :: hnature  ! scheme for natural surfaces
 CHARACTER(LEN=6),  INTENT(OUT) :: hsea     ! scheme for sea
 CHARACTER(LEN=6),  INTENT(OUT) :: htown    ! scheme for towns
 CHARACTER(LEN=6),  INTENT(OUT) :: hwater   ! scheme for inland water
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER           :: iluout     ! logical unit
INTEGER           :: ilunam     ! logical unit
!
LOGICAL           :: gfound
REAL(KIND=JPRB) :: zhook_handle
!
!------------------------------------------------------------------------------
!
!*       1.    defaults
! 
IF (lhook) CALL dr_hook('READ_PGD_SCHEMES',0,zhook_handle)
 CALL default_schemes(hprogram,cnature,csea,ctown,cwater)
!
!------------------------------------------------------------------------------
!
!*       2.    opening of namelist
! 
 CALL get_luout(hprogram,iluout)
!
 CALL open_namelist(hprogram,ilunam)
!
!-------------------------------------------------------------------------------
!
!*       3.    reading of namelist
! 
 CALL posnam(ilunam,'NAM_PGD_SCHEMES',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_pgd_schemes)
!
!-------------------------------------------------------------------------------
!
!*       4.    check of file type
! 
 CALL test_nam_var_surf(iluout,'CNATURE',cnature,'NONE  ','ISBA  ','TSZ0  ','FLUX  ')
 CALL test_nam_var_surf(iluout,'CSEA   ',csea   ,'NONE  ','SEAFLX','FLUX  ')
 CALL test_nam_var_surf(iluout,'CTOWN  ',ctown  ,'NONE  ','TEB   ','FLUX  ')
 CALL test_nam_var_surf(iluout,'CWATER ',cwater ,'NONE  ','WATFLX','FLUX  ','FLAKE ')
!
hnature = cnature
hsea    = csea
htown   = ctown
hwater  = cwater
!------------------------------------------------------------------------------
!
!*       5.    close namelist file
! 
 CALL close_namelist(hprogram,ilunam)
IF (lhook) CALL dr_hook('READ_PGD_SCHEMES',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_pgd_schemes