/home/dasprezs/EXPORT_v7_3/src/SURFEX/default_isba.F90

Go to the documentation of this file.

!     #########
      SUBROUTINE DEFAULT_ISBA(PTSTEP, POUT_TSTEP,                     &
                              HROUGH, HRUNOFF, HALBEDO, HSCOND,       &
                              HC1DRY, HSOILFRZ, HDIFSFCOND, HSNOWRES, &
                              HCPSURF, PCGMAX, PCDRAG, HKSAT, HSOC,   &
                              HTOPREG, HRAIN, HHORT, OFLOOD, OTRIP,   &
                              OGLACIER, OCANOPY_DRAG, OVEGUPD,        &
                              OSPINUPCARBS, OSPINUPCARBW              )  
!     ########################################################################
!
!!****  *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 
!-------------------------------------------------------------------------------
!
!*       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) :: HTOPREG   ! Wolock and McCabe (2000) linear regression for Topmodel
                                           ! 'DEF' = reg
                                           ! 'NON' = no reg
!
 CHARACTER(LEN=3), INTENT(OUT) :: HKSAT     ! SOIL HYDRAULIC CONDUCTIVITY PROFILE OPTION
!                                          ! 'DEF'  = ISBA homogenous soil
!                                          ! 'SGH'  = ksat exponential decay
!
 CHARACTER(LEN=3), INTENT(OUT) :: HSOC      ! SOIL ORGANIC CARBON PROFILE OPTION
!                                          ! 'DEF'  = ISBA homogenous soil
!                                          ! 'SGH'  = 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)          :: OTRIP     ! True  = ISBA-TRIP coupling
                                           ! False = No ISBA-TRIP coupling
!                                          
LOGICAL, INTENT(OUT)          :: OFLOOD    ! True  = Flooding scheme
                                           ! False = No flooding scheme
!
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(KIND=JPRB) :: ZHOOK_HANDLE
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
!-------------------------------------------------------------------------------
!
IF (LHOOK) CALL DR_HOOK('DEFAULT_ISBA',0,ZHOOK_HANDLE)
PTSTEP     = XUNDEF
POUT_TSTEP = XUNDEF
!
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"
HTOPREG    = 'DEF'
HKSAT      = 'DEF'
HSOC       = 'DEF'
HRAIN      = 'DEF'
HHORT      = 'DEF'
!
OTRIP = .FALSE.
OFLOOD = .FALSE.
!
PCGMAX   = 2.0E-5
!
PCDRAG   = 0.15
!
OGLACIER = .FALSE.
!
OCANOPY_DRAG = .FALSE.
!
OVEGUPD = .TRUE.
!
OSPINUPCARBS = .FALSE.
OSPINUPCARBW = .FALSE.
!
IF (LHOOK) CALL DR_HOOK('DEFAULT_ISBA',1,ZHOOK_HANDLE)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE DEFAULT_ISBA