read_nam_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_topd(HPROGRAM,&
                                obudget_topd,knb_topd,&
                                ostock_topd,&
                                knb_stock,knb_restart,&
                                kfreq_maps_wg,kfreq_maps_asat,kfreq_maps_runoff,&
                                pspeedr,pspeedg,pspeedh,pqinit,prtop_d2)
!     ##############################################################
!
!!**** *READ_NAM TOPD* reads namelist NAM_TOPD
!!
!!    PURPOSE
!!    -------
!!
!!    NAM_TOPD is a namelist used to define whether Topmodel coupling
!!    is performed or not and the time step ratio between hydrological 
!!    model and ISBA.
!!    This routine aims at reading and initialising those variables.
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    B. Vincendon        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    11/2006
!!    B. Vincendon 02/2014 : adding possibility to choose the speed of water on hillslopes and
!!                           to write runoff maps on watersheds
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_topd_par, ONLY : jpcat
USE modd_topodyn, ONLY : nncat
!
USE mode_pos_surf
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
!
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)  :: obudget_topd ! budget computation
INTEGER,                            INTENT(OUT)  :: knb_topd     ! Ratio between Topmodel time step and ISBA time step
LOGICAL,                            INTENT(OUT)  :: ostock_topd  ! T if use of stock from previous simulation
INTEGER,                            INTENT(OUT)  :: knb_stock    ! number of time step to read in previous simulation
INTEGER,                            INTENT(OUT)  :: knb_restart  ! number of time step to write for next simulation
INTEGER,                            INTENT(OUT)  :: kfreq_maps_wg! 
INTEGER,                            INTENT(OUT)  :: kfreq_maps_asat! 
INTEGER,                            INTENT(OUT)  :: kfreq_maps_runoff
REAL, DIMENSION(JPCAT),INTENT(OUT)               :: pspeedr ! River speed
REAL, DIMENSION(JPCAT),INTENT(OUT)               :: pspeedg ! Ground speed
REAL, DIMENSION(JPCAT),INTENT(OUT)               :: pspeedh ! Hillslope speed
REAL, DIMENSION(JPCAT),INTENT(OUT)               :: pqinit  ! Initial discharge at catchments outlet
REAL, DIMENSION(JPCAT),INTENT(OUT)               :: prtop_d2
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
LOGICAL                           :: lbudget_topd
LOGICAL                           :: lstock_topd
INTEGER                           :: nnb_topd
INTEGER                           :: nfreq_maps_wg
INTEGER                           :: nfreq_maps_asat
INTEGER                           :: nfreq_maps_runoff
INTEGER                           :: nnb_stp_stock
INTEGER                           :: nnb_stp_restart
REAL, DIMENSION(JPCAT)            :: xspeedr
REAL, DIMENSION(JPCAT)            :: xspeedg
REAL, DIMENSION(JPCAT)            :: xspeedh
REAL, DIMENSION(JPCAT)            :: xqinit
REAL, DIMENSION(JPCAT)            :: xrtop_d2
!
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_topd/lbudget_topd, lstock_topd, nnb_topd, &
                  nfreq_maps_wg, nfreq_maps_asat, nfreq_maps_runoff,&
                  nnb_stp_stock, nnb_stp_restart, &
                  xspeedr, xspeedg, xspeedh, xqinit, xrtop_d2
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('READ_NAM_TOPD',0,zhook_handle)
!
!*    1.      Initializations of defaults
!             ---------------------------
!
lbudget_topd = .false.
lstock_topd = .false.
nnb_topd = 1
nfreq_maps_wg = 0
nfreq_maps_asat = 0
nfreq_maps_runoff = 0
nnb_stp_stock = 1
nnb_stp_restart = 1
xspeedr(:) = 3.0 ! default value of river speed (adapted for Cevennes zone)
xspeedg(:) = 0.3 ! default value of speed in the ground
xspeedh(:) = 0.3 ! default value of hillspeed
xqinit(:) = 0.
xrtop_d2(:) = 1.
!
 CALL get_luout(hprogram,iluout)
!
!-------------------------------------------------------------------------------
!
!*    2.      Reading of namelist
!             -------------------
!
 CALL open_namelist(hprogram,ilunam)
!
 CALL posnam(ilunam,'NAM_TOPD',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_topd)
!
 CALL close_namelist(hprogram,ilunam)
!
!-------------------------------------------------------------------------------
!
!*    3.      Fills output arguments
!             ----------------------
!
obudget_topd = lbudget_topd
ostock_topd = lstock_topd
knb_topd = nnb_topd
knb_stock = nnb_stp_stock
knb_restart = nnb_stp_restart
kfreq_maps_wg = nfreq_maps_wg
kfreq_maps_asat = nfreq_maps_asat
kfreq_maps_runoff = nfreq_maps_runoff
pspeedr(1:nncat) = xspeedr(1:nncat)
pspeedg(1:nncat) = xspeedg(1:nncat)
WHERE(xspeedh(1:nncat)/=0.3)
 pspeedh(1:nncat) = xspeedh(1:nncat)
ELSEWHERE
 pspeedh(1:nncat) = xspeedr(1:nncat)/10.
ENDWHERE
pqinit(1:nncat) = xqinit(1:nncat)
prtop_d2(1:nncat) = xrtop_d2(1:nncat)
!
WRITE(iluout,*) 'NAM_TOPD:'
WRITE(iluout,*) 'LBUDGET ',lbudget_topd
WRITE(iluout,*) 'NNB_TOP',nnb_topd
WRITE(iluout,*) 'LSTOCK',lstock_topd
WRITE(iluout,*) 'NNB_RESTART,NNB_STOCK',nnb_stp_restart,nnb_stp_stock
WRITE(iluout,*) 'NFREQ_MAPS_WG,NFREQ_MAPS_ASAT',nfreq_maps_wg,nfreq_maps_asat
WRITE(iluout,*) 'NFREQ_MAPS_RUNOFF',nfreq_maps_runoff
!
IF (lhook) CALL dr_hook('READ_NAM_TOPD',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_nam_topd