read_nam_pgd_topd.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_pgd_topd(HPROGRAM,OCOUPL_TOPD,HCAT,PF_PARAM_BV,PC_DEPTH_RATIO_BV)
!     ##############################################################
!
!!**** *READ_NAM_TOPD_n* reads namelist NAM_TOPD
!!
!!    PURPOSE
!!    -------
!!    NAM_TOPD is a namelist used only for Topmodel coupling
!!    It permits to define the different catchments studied.
!!    This routine aims at reading and initialising those names.
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    B. Vincendon        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original   11/2006
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
!
USE modd_topd_par, ONLY : jpcat
USE modd_topodyn, ONLY : nncat
!
USE mode_pos_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
 CHARACTER(LEN=6),                   INTENT(IN)   :: hprogram     ! Type of program
LOGICAL, INTENT(OUT) :: ocoupl_topd
 CHARACTER(LEN=15), DIMENSION(JPCAT),INTENT(OUT)  :: hcat         ! Names of catchments         
REAL, DIMENSION(JPCAT),INTENT(OUT)               :: pf_param_bv
REAL, DIMENSION(JPCAT),INTENT(OUT)               :: pc_depth_ratio_bv 
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
 CHARACTER(LEN=15), DIMENSION(JPCAT) :: ccat
LOGICAL                           :: lcoupl_topd
REAL, DIMENSION(JPCAT)            :: xf_param_bv
REAL, DIMENSION(JPCAT)            :: xc_depth_ratio_bv
!
INTEGER                           :: iluout    ! output listing logical unit
INTEGER                           :: ilunam    ! namelist file logical unit
LOGICAL                           :: gfound    ! flag when namelist is present
REAL(KIND=JPRB) :: zhook_handle
!
!*    0.3    Declaration of namelists
!  
!
namelist/nam_pgd_topd/ccat, lcoupl_topd, xf_param_bv, xc_depth_ratio_bv   
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('READ_NAM_PGD_TOPD',0,zhook_handle)
!
!*    1.      Initializations of defaults
!             ---------------------------
!
lcoupl_topd = .false.
 ccat(:) = '   '
xf_param_bv(:) = 2.5
xc_depth_ratio_bv(:) = 1.
!
 CALL get_luout(hprogram,iluout)
!
!-------------------------------------------------------------------------------
!
!*    2.      Reading of namelist
!             -------------------
!
 CALL open_namelist(hprogram,ilunam)
!CALL OPEN_NAMELIST(HPROGRAM,'SURF  ',ILUNAM)
!
 CALL posnam(ilunam,'NAM_PGD_TOPD',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_pgd_topd)
!
 CALL close_namelist(hprogram,ilunam)
!
!         2.   Initialises number of catchments and time step variables
!              -------------------------------------------------------
!
nncat=count(ccat(:)/='   ')
!
!-------------------------------------------------------------------------------
!
!*    3.      Fills output arguments
!             ----------------------
!
ocoupl_topd = lcoupl_topd
hcat(1:nncat) = ccat(1:nncat)
pf_param_bv(1:nncat) = xf_param_bv(1:nncat)
pc_depth_ratio_bv(1:nncat) = xc_depth_ratio_bv(1:nncat)
!
IF (lhook) CALL dr_hook('READ_NAM_PGD_TOPD',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_nam_pgd_topd