modd_isba_par.F90 Source File


Contents

Source Code


Source Code

!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 version 2.1
!SFX_LIC version 1. See LICENSE, Licence_CeCILL_V2.1-en.txt and Licence_CeCILL_V2.1-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