SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
modd_isba_par.F90
Go to the documentation of this file.
1 !SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
3 !SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
4 !SFX_LIC for details. version 1.
5 ! ######################
6  MODULE modd_isba_par
7 ! ######################
8 !
9 !!**** *MODD_ISBA_PAR* - declaration of ISBA parameters
10 !!
11 !! PURPOSE
12 !! -------
13 ! The purpose of this declarative module is to specify the
14 ! parameters related to the surface parameterization ISBA.
15 !
16 !!
17 !!
18 !!
19 !! AUTHOR
20 !! ------
21 !! S. Belair *Meteo France*
22 !!
23 !! MODIFICATIONS
24 !! -------------
25 !! Original 29/04/95
26 !! (V.Masson) 05/10/98+ add XCDZ0EFF, XRHOSMIN, XRHOSMAX
27 !! (V.Masson) 15/03/99 add number of layers
28 !! (A.Boone) 02/05/02 add ISBA-ES parameters
29 !! (A.Boone) 21/11/11 add Rsmax
30 !! (S.Gollvik) 20/02/12 add XFLXMAX
31 !! (A.Boone) 20/02/12 add ISBA-MEB parameters
32 !! (B. Decharme) 07/15 Add numerical adjustement for F2 soilstress function
33 !-------------------------------------------------------------------------------
34 !
35 !* 0. DECLARATIONS
36 ! ------------
37 !
38 IMPLICIT NONE
39 !
40 !------------------------------------------------------------------------------
41 ! Vegetation: (for additional parameters related to vegetation, see
42 ! modd_co2v_par.f90)
43 !------------------------------------------------------------------------------
44 !
45 ! vegetation emissivity
46 !
47 REAL, PARAMETER :: XEMISVEG = 0.97
48 !
49 ! drag coefficient in z0eff computation
50 !
51 REAL, PARAMETER :: XCDZ0EFF = 0.8
52 !
53 ! minimum vegetation fraction (for C3 grassland: for the case with large
54 ! VEG and low LAI, such as wintertime)
55 !
56 REAL, PARAMETER :: XVEGMIN = 0.95
57 !
58 ! Maximum stomatal resistance (s m-1)
59 !
60 REAL, PARAMETER :: XRS_MAX = 5000.
61 !
62 ! Factor to restore explicit Cv value (DIF option)
63 !
64 REAL, PARAMETER :: XCVHEATF = 0.20
65 !
66 ! Numerical factor to prevent division by 0 for F2 soilstress function
67 !
68 REAL, PARAMETER :: XDENOM_MIN = 1.E-12
69 !
70 !--------------------------------------------------------------------------------
71 ! Soil:
72 !--------------------------------------------------------------------------------
73 !
74 ! Caracteristic time for ice in force-restore (s)
75 !
76 REAL, PARAMETER :: XTAU_ICE = 3300.
77 !
78 ! Bare soil emissivity
79 !
80 REAL, PARAMETER :: XEMISSOIL = 0.94
81 !
82 ! Minimum allowable volumetric liquid water content of soil
83 !
84 REAL, PARAMETER :: XWGMIN = 0.001 ! (m3 m-3)
85 !
86 ! Peters-Lidard et al. (JAS, 1998) from method of Johanssen (1975)
87 ! thermal conductivity (option) parameters:
88 !
89 REAL, PARAMETER :: XSPHSOIL = 733. ! J/(kg K) Soil specific heat
90 REAL, PARAMETER :: XDRYWGHT = 2700.0 ! kg/m3 Soil solids dry weight
91 REAL, PARAMETER :: XCONDQRTZ = 7.7 ! W/(m K) Quartz thermal conductivity
92 REAL, PARAMETER :: XCONDOTH1 = 2.0 ! W/(m K) Other thermal conductivity
93 REAL, PARAMETER :: XCONDOTH2 = 3.0 ! W/(m K) Other thermal conductivity
94 REAL, PARAMETER :: XCONDWTR = 0.57 ! W/(m K) Water thermal conductivity
95 !
96 REAL, PARAMETER :: XOMRHO = 1300. !Organic mater density (kg.m-3)
97 REAL, PARAMETER :: XOMSPH = 1926. !Organic mater specific heat (J/(kg K))
98 REAL, PARAMETER :: XOMCONDDRY = 0.05 !Organic mater dry thermal conductivity (W.m–1.K–1)
99 REAL, PARAMETER :: XOMCONDSLD = 0.25 !Organic mater solid thermal conductivity (W.m–1.K–1)
100 !
101 ! Maximum depth of the water table for soil thermal computation
102 !
103 REAL, PARAMETER :: XWTD_MAXDEPTH = 100. !m
104 !
105 ! Minimun depth of permafrost and limit area
106 !
107 REAL, PARAMETER :: XPERMFRAC = 0.25 ! permafrost limit area (fraction)
108 !
109 REAL, PARAMETER :: XPERMDEPTH = 12.0 ! permafrost depth (m)
110 !
111 !--------------------------------------------------------------------------------
112 ! Vegetation radiative properties
113 !--------------------------------------------------------------------------------
114 !
115 ! Wavelength between near-infra-red and visible parts of the solar spectra
116 !
117 REAL, PARAMETER :: XRED_EDGE = 0.0000007 ! (m) 0.7 micro-m
118 !
119 !
120 ! Wavelength between visible and UV parts of the solar spectra
121 !
122 REAL, PARAMETER :: XUV_EDGE = 0.0000002 ! (m) 0.1 micro-m
123 !
124 !--------------------------------------------------------------------------------
125 ! MEB: Multiple energy balance parameters
126 !--------------------------------------------------------------------------------
127 !
128 REAL, PARAMETER :: XFLXMAX = 5000. ! [kg/(m**2*s)]
129 ! Maximum value of exchange coeffient
130 ! (should go to infinity, for some cases, i.e. when lai=>0)
131 !
132 REAL, PARAMETER :: XLIMH = 2.0 ! m
133 ! Minimum forcing height above vegetation top (turbulence computations)
134 !
135 ! Soil geometry if DF option
136 !--------------------------------------------------------------------------------
137 !
138 INTEGER, PARAMETER :: NOPTIMLAYER=14
139 REAL, DIMENSION(NOPTIMLAYER), PARAMETER :: XOPTIMGRID = &
140  (/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/)
141 !
142 !--------------------------------------------------------------------------------
143 !
144 END MODULE modd_isba_par
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156