default_isba.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 default_isba(PTSTEP, POUT_TSTEP,                        &
                              hrough, hrunoff, halbedo, hscond,          &
                              hc1dry, hsoilfrz, hdifsfcond, hsnowres,    &
                              hcpsurf, pcgmax, pcdrag, hksat, osoc,      &
                              hrain, hhort, oglacier, ocanopy_drag,      &
                              ovegupd, ospinupcarbs, ospinupcarbw,       &
                              pspinmaxs, pspinmaxw, pco2_start, pco2_end,&
                              knbyearspins, knbyearspinw,                &
                              onitro_dilu                                )
!     ########################################################################
!
!!****  *DEFAULT_ISBA* - routine to set default values for the configuration for ISBA
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004 
!!      B.Decharme  04/2013 delete HTOPREG (never used)
!!                          water table / surface coupling 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_surf_par,   ONLY : xundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
REAL,              INTENT(OUT) :: ptstep     ! time-step for run
REAL,              INTENT(OUT) :: pout_tstep ! time-step for writing
 CHARACTER(LEN=4),  INTENT(OUT) :: hrough   ! type of roughness length
 CHARACTER(LEN=4),  INTENT(OUT) :: halbedo  ! albedo type
!                                          ! 'DRY ' 
!                                          ! 'EVOL' 
!                                          ! 'WET ' 
!                                          ! 'USER' 
 CHARACTER(LEN=4),  INTENT(OUT) :: hscond   ! Thermal conductivity
!                                          ! 'DEF ' = DEFault: NP89 implicit method
!                                          ! 'PL98' = Peters-Lidard et al. 1998 used
!                                          ! for explicit computation of CG
 CHARACTER(LEN=4),  INTENT(OUT) :: hc1dry   ! C1 formulation for dry soils
!                                          ! 'DEF ' = DEFault: Giard-Bazile formulation
!                                          ! 'GB93' = Giordani 1993, Braud 1993 
!                                          !discontinuous at WILT
 CHARACTER(LEN=3),  INTENT(OUT) :: hsoilfrz ! soil freezing-physics option
!                                          ! 'DEF' = Default (Boone et al. 2000; 
!                                          !        Giard and Bazile 2000)
!                                          ! 'LWT' = Phase changes as above,
!                                          !         but relation between unfrozen 
!                                          !         water and temperature considered
!                            NOTE that when using the YISBA='DIF' multi-layer soil option,
!                            the 'LWT' method is used. It is only an option
!                            when using the force-restore soil method ('2-L' or '3-L')
!
 CHARACTER(LEN=4),  INTENT(OUT) :: hdifsfcond ! Mulch effects
!                                          ! 'MLCH' = include the insulating effect of
!                                          ! leaf litter/mulch on the surf. thermal cond.
!                                          ! 'DEF ' = no mulch effect
!                           NOTE: Only used when YISBA = DIF
!
 CHARACTER(LEN=3), INTENT(OUT) :: hsnowres  ! Turbulent exchanges over snow
!                                          ! 'DEF' = Default: Louis (ISBA)
!                                          ! 'RIL' = Maximum Richardson number limit
!                                          !         for stable conditions ISBA-SNOW3L
!                                          !         turbulent exchange option
 CHARACTER(LEN=3), INTENT(OUT) :: hcpsurf   ! SPECIFIC HEAT
!                                          ! 'DRY' = dry Cp
!                                          ! 'HUM' = Cp fct of qs
REAL,              INTENT(OUT) :: pcgmax   ! maximum soil heat capacity
!
REAL,              INTENT(OUT) :: pcdrag   ! drag coefficient in canopy
!
 CHARACTER(LEN=4),  INTENT(OUT) :: hrunoff  ! surface runoff formulation
!                                          ! 'WSAT'
!                                          ! 'DT92'
!                                          ! 'SGH ' Topmodel
!
 CHARACTER(LEN=3), INTENT(OUT) :: hksat     ! SOIL HYDRAULIC CONDUCTIVITY PROFILE OPTION
!                                          ! 'DEF'  = ISBA homogenous soil
!                                          ! 'SGH'  = ksat exponential decay
!
LOGICAL, INTENT(OUT) ::          osoc      ! SOIL ORGANIC CARBON PROFILE OPTION
!                                          ! False  = ISBA homogenous soil
!                                          ! True   = SOC profile effect
!
 CHARACTER(LEN=3), INTENT(OUT) :: hrain     ! Rainfall spatial distribution
                                           ! 'DEF' = No rainfall spatial distribution
                                           ! 'SGH' = Rainfall exponential spatial distribution
                                           ! 
! 
 CHARACTER(LEN=3), INTENT(OUT) :: hhort     ! Horton runoff
                                           ! 'DEF' = no Horton runoff
                                           ! 'SGH' = Horton runoff
!                                         
LOGICAL, INTENT(OUT)          :: oglacier  ! True = Over permanent snow and ice, 
!                                                   initialise WGI=WSAT, 
!                                                   Hsnow>=3.3m and allow 0.8<SNOALB<0.85
                                           ! False = No specific treatment
LOGICAL, INTENT(OUT)          :: ocanopy_drag ! T: drag activated in SBL scheme within the canopy
!
LOGICAL, INTENT(OUT)          :: ovegupd   ! T: update vegetation parameters 
                                           !    every decade
                                           ! F: keep vegetation parameters
                                           !    constant in time
!
LOGICAL, INTENT(OUT)          :: ospinupcarbs ! T: carbon spinup soil
LOGICAL, INTENT(OUT)          :: ospinupcarbw ! T: carbon spinup wood
REAL,    INTENT(OUT)          :: pspinmaxs    ! max number of times CARBON_SOIL subroutine is called
REAL,    INTENT(OUT)          :: pspinmaxw    ! max number of times the wood is accelerated 
REAL,    INTENT(OUT)          :: pco2_start   ! Pre-industrial CO2 concentration
REAL,    INTENT(OUT)          :: pco2_end     ! Begin-transient CO2 concentration
INTEGER, INTENT(OUT)          :: knbyearspins ! nbr years needed to reaches soil equilibrium
INTEGER, INTENT(OUT)          :: knbyearspinw ! nbr years needed to reaches wood equilibrium
!
LOGICAL, INTENT(OUT)          :: onitro_dilu ! nitrogen dilution fct of CO2 (Calvet et al. 2008)
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('DEFAULT_ISBA',0,zhook_handle)
!
ptstep     = xundef
pout_tstep = xundef
!!!!!do not phased!!!!!!!
!HROUGH  = "NONE"
!!!!!do not phased!!!!!!!
hrough  = "UNDE"  ! undefined. Needs further information on canopy scheme use to set default
hscond  = "PL98"
halbedo = "DRY "
!
hc1dry     = 'DEF '
hsoilfrz   = 'DEF'
hdifsfcond = 'DEF '
hsnowres   = 'DEF'
hcpsurf    = 'DRY'
!
hrunoff    = "WSAT"
hksat      = 'DEF'
osoc       = .false.
hrain      = 'DEF'
hhort      = 'DEF'
!
pcgmax   = 2.0e-5
!
pcdrag   = 0.15
!
oglacier = .false.
!
ocanopy_drag = .false.
!
ovegupd = .true.
!
ospinupcarbs = .false.
ospinupcarbw = .false.
!
pspinmaxs = 0.
pspinmaxw = 0.
pco2_start= xundef
pco2_end  = xundef
knbyearspins = 0
knbyearspinw = 0
!
onitro_dilu = .false.
!
IF (lhook) CALL dr_hook('DEFAULT_ISBA',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE default_isba