modd_isba_par.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.
!     ######################
      MODULE modd_isba_par
!     ######################
!
!!****  *MODD_ISBA_PAR* - declaration of ISBA parameters
!!
!!    PURPOSE
!!    -------
!       The purpose of this declarative module is to specify  the 
!     parameters related to the surface parameterization ISBA. 
!
!!
!!      
!!
!!    AUTHOR
!!    ------
!!      S. Belair   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       29/04/95                      
!!      (V.Masson)     05/10/98+ add XCDZ0EFF, XRHOSMIN, XRHOSMAX
!!      (V.Masson)     15/03/99 add number of layers
!!      (A.Boone)      02/05/02 add ISBA-ES parameters
!!      (A.Boone)      21/11/11 add Rsmax
!!      (S.Gollvik)    20/02/12 add XFLXMAX
!!      (A.Boone)      20/02/12 add ISBA-MEB parameters
!!     (B. Decharme)      07/15 Add numerical adjustement for F2 soilstress function
!-------------------------------------------------------------------------------
!
!*       0.   DECLARATIONS
!             ------------
!
IMPLICIT NONE
!
!------------------------------------------------------------------------------
! Vegetation: (for additional parameters related to vegetation, see
! modd_co2v_par.f90)
!------------------------------------------------------------------------------
!
! vegetation emissivity
!
REAL, PARAMETER       :: xemisveg = 0.97
!
! drag coefficient in z0eff computation
!
REAL, PARAMETER       :: xcdz0eff = 0.8
!
! minimum vegetation fraction (for C3 grassland: for the case with large 
! VEG and low LAI, such as wintertime)
!
REAL, PARAMETER       :: xvegmin   = 0.95
!
! Maximum stomatal resistance (s m-1)
!
REAL, PARAMETER       :: xrs_max   = 5000. 
!
! Factor to restore explicit Cv value (DIF option)
!
REAL, PARAMETER       :: xcvheatf  = 0.20 
!
! Numerical factor to prevent division by 0 for F2 soilstress function
!
REAL, PARAMETER       :: xdenom_min  = 1.e-12 
!
!--------------------------------------------------------------------------------
! Soil:
!--------------------------------------------------------------------------------
!                        
! Caracteristic time for ice in force-restore (s)
!
REAL, PARAMETER       :: xtau_ice = 3300.
!                        
! Bare soil emissivity
!
REAL, PARAMETER       :: xemissoil = 0.94
!                        
! Minimum allowable volumetric liquid water content of soil
!
REAL, PARAMETER       :: xwgmin   = 0.001   ! (m3 m-3)
!
! Peters-Lidard et al. (JAS, 1998) from method of Johanssen (1975)
! thermal conductivity (option) parameters:
!
REAL, PARAMETER       :: xsphsoil  = 733.   ! J/(kg K) Soil specific heat
REAL, PARAMETER       :: xdrywght  = 2700.0 ! kg/m3    Soil solids dry weight
REAL, PARAMETER       :: xcondqrtz = 7.7    ! W/(m K)  Quartz thermal conductivity
REAL, PARAMETER       :: xcondoth1 = 2.0    ! W/(m K)  Other thermal conductivity
REAL, PARAMETER       :: xcondoth2 = 3.0    ! W/(m K)  Other thermal conductivity
REAL, PARAMETER       :: xcondwtr  = 0.57   ! W/(m K)  Water thermal conductivity
!
REAL, PARAMETER       :: xomrho     = 1300.   !Organic mater density (kg.m-3)
REAL, PARAMETER       :: xomsph     = 1926.   !Organic mater specific heat              (J/(kg K))
REAL, PARAMETER       :: xomconddry = 0.05    !Organic mater dry thermal conductivity   (W.m–1.K–1)
REAL, PARAMETER       :: xomcondsld = 0.25    !Organic mater solid thermal conductivity (W.m–1.K–1)
!                        
! Maximum depth of the water table for soil thermal computation
!
REAL, PARAMETER       :: xwtd_maxdepth = 100. !m
!                        
! Minimun depth of permafrost and limit area
!
REAL, PARAMETER :: xpermfrac  = 0.25   ! permafrost limit area (fraction)
!
REAL, PARAMETER :: xpermdepth = 12.0   ! permafrost depth (m)
!
!--------------------------------------------------------------------------------
! Vegetation radiative properties
!--------------------------------------------------------------------------------
!                        
! Wavelength between near-infra-red and visible parts of the solar spectra
!
REAL, PARAMETER       :: xred_edge = 0.0000007  ! (m)   0.7 micro-m
!
!                        
! Wavelength between visible and UV parts of the solar spectra
!
REAL, PARAMETER       :: xuv_edge  = 0.0000002 ! (m)   0.1 micro-m
!
!--------------------------------------------------------------------------------
! MEB: Multiple energy balance  parameters
!--------------------------------------------------------------------------------
!                        
REAL, PARAMETER       :: xflxmax = 5000.   ! [kg/(m**2*s)]
!                        Maximum value of exchange coeffient
!                        (should go to infinity, for some cases, i.e. when lai=>0) 
!                        
REAL, PARAMETER       :: xlimh       = 2.0 ! m
!                        Minimum forcing height above vegetation top (turbulence computations)
!
! Soil geometry if DF option
!--------------------------------------------------------------------------------
!
INTEGER,                      PARAMETER :: noptimlayer=14
REAL, DIMENSION(NOPTIMLAYER), PARAMETER :: xoptimgrid = & 
      (/0.01,0.04,0.10,0.20,0.40,0.60,0.80,1.00,1.50,2.00,3.00,5.00,8.00,12.0/)
!
!--------------------------------------------------------------------------------
!
END MODULE modd_isba_par