[8] | 1 | !$Header: /home/ssipsl/CVSREP/ORCHIDEE/src_parameters/constantes_soil.f90,v 1.11 2010/04/06 14:38:48 ssipsl Exp $ |
---|
| 2 | !IPSL (2006) |
---|
| 3 | ! This software is governed by the CeCILL licence see ORCHIDEE/ORCHIDEE_CeCILL.LIC |
---|
| 4 | !- |
---|
| 5 | MODULE constantes_soil |
---|
| 6 | !!-------------------------------------------------------------------- |
---|
| 7 | !! "constantes_soil" module contains public data for the soils |
---|
| 8 | !!-------------------------------------------------------------------- |
---|
| 9 | USE constantes |
---|
| 10 | !- |
---|
| 11 | IMPLICIT NONE |
---|
| 12 | !- |
---|
| 13 | ! Dimensioning parameters |
---|
| 14 | !- |
---|
| 15 | ! Number of soil level |
---|
| 16 | INTEGER(i_std),PARAMETER :: ngrnd=7 |
---|
| 17 | ! Number of diagnostic levels in the soil |
---|
| 18 | INTEGER(i_std),PARAMETER :: nbdl=11 |
---|
| 19 | !MM : if you want to compare hydrology variables with old TAG 1.6 and lower, |
---|
| 20 | ! you must set the Number of diagnostic levels in the soil to 6 : |
---|
| 21 | ! INTEGER(i_std),PARAMETER :: nbdl=6 |
---|
| 22 | ! Number of levels in CWRR |
---|
| 23 | INTEGER(i_std),PARAMETER :: nslm=11 |
---|
| 24 | ! Number of soil types |
---|
| 25 | INTEGER(i_std),PARAMETER :: nstm=3 |
---|
| 26 | !- |
---|
| 27 | !- Parameters for soil thermodynamics |
---|
| 28 | !- |
---|
| 29 | ! Average Thermal Conductivity of soils |
---|
| 30 | REAL(r_std),PARAMETER :: so_cond = 1.5396 |
---|
| 31 | ! Average Heat capacity of soils |
---|
| 32 | REAL(r_std),PARAMETER :: so_capa = 2.0514e+6 |
---|
| 33 | !- |
---|
| 34 | ! Values taken from : PIELKE,'MESOSCALE METEOROLOGICAL MODELING',P.384 |
---|
| 35 | ! Dry soil heat capacity was decreased and conductivity increased. |
---|
| 36 | !- |
---|
| 37 | ! Dry soil Heat capacity of soils |
---|
| 38 | !*REAL(r_std),PARAMETER :: so_capa_dry = 1.35e+6 |
---|
| 39 | REAL(r_std),PARAMETER :: so_capa_dry = 1.80e+6 |
---|
| 40 | ! Dry soil Thermal Conductivity of soils |
---|
| 41 | !*REAL(r_std),PARAMETER :: so_cond_dry = 0.28 |
---|
| 42 | REAL(r_std),PARAMETER :: so_cond_dry = 0.40 |
---|
| 43 | !- |
---|
| 44 | ! Wet soil Heat capacity of soils |
---|
| 45 | REAL(r_std),PARAMETER :: so_capa_wet = 3.03e+6 |
---|
| 46 | ! Wet soil Thermal Conductivity of soils |
---|
| 47 | REAL(r_std),PARAMETER :: so_cond_wet = 1.89 |
---|
| 48 | !- |
---|
| 49 | ! Thermal Conductivity of snow |
---|
| 50 | REAL(r_std),PARAMETER :: sn_cond = 0.3 |
---|
| 51 | ! Snow density for the soil thermodynamics |
---|
| 52 | REAL(r_std),PARAMETER :: sn_dens = 330.0 |
---|
| 53 | ! Heat capacity for snow |
---|
| 54 | REAL(r_std),PARAMETER :: sn_capa = 2100.0_r_std*sn_dens |
---|
| 55 | !- |
---|
| 56 | ! Constantes from the Choisnel hydrology |
---|
| 57 | !- |
---|
| 58 | ! Wilting point (Has a numerical role for the moment) |
---|
| 59 | REAL(r_std),PARAMETER :: qwilt = 5.0 |
---|
| 60 | ! Total depth of soil reservoir (for hydrolc) |
---|
| 61 | REAL(r_std),SAVE :: dpu_cste=2.0_r_std |
---|
| 62 | ! The minimal size we allow for the upper reservoir (m) |
---|
| 63 | REAL(r_std),PARAMETER :: min_resdis = 2.e-5 |
---|
| 64 | ! Diffusion constant for the slow regime |
---|
| 65 | ! (This is for the diffusion between reservoirs) |
---|
| 66 | REAL(r_std),PARAMETER :: min_drain = 0.001 |
---|
| 67 | ! Diffusion constant for the fast regime |
---|
| 68 | REAL(r_std),PARAMETER :: max_drain = 0.1 |
---|
| 69 | ! The exponential in the diffusion law |
---|
| 70 | REAL(r_std),PARAMETER :: exp_drain = 1.5 |
---|
| 71 | ! Transforms leaf area index into size of interception reservoir |
---|
| 72 | REAL(r_std),SAVE :: qsintcst = 0.1 |
---|
| 73 | ! Maximum quantity of water (Kg/M3) |
---|
| 74 | REAL(r_std),PARAMETER :: mx_eau_eau = 150. |
---|
| 75 | !- |
---|
| 76 | ! Constant in the computation of resistance for bare soil evaporation |
---|
| 77 | REAL(r_std),PARAMETER :: rsol_cste = 33.E3 |
---|
| 78 | ! Scaling depth for litter humidity (m) |
---|
| 79 | !SZ changed this according to SP from 0.03 to 0.08, 080806 |
---|
| 80 | REAL(r_std),PARAMETER :: hcrit_litter=0.08_r_std |
---|
| 81 | !- |
---|
| 82 | ! Parameters for soil type distribution |
---|
| 83 | !- |
---|
| 84 | ! Default soil texture distribution in the following order : |
---|
| 85 | ! sand, loam and clay |
---|
| 86 | REAL(r_std),DIMENSION(nstm),SAVE :: soiltype_default = & |
---|
| 87 | & (/ 0.0, 1.0, 0.0 /) |
---|
| 88 | !- |
---|
| 89 | ! Parameters specific for the CWRR hydrology. |
---|
| 90 | !- |
---|
| 91 | ! Van genuchten coefficient n |
---|
| 92 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: nvan = & |
---|
| 93 | & (/ 1.89_r_std, 1.56_r_std, 1.31_r_std /) |
---|
| 94 | !!$! Van genuchten coefficient a (cm^{-1}) |
---|
| 95 | !!$ REAL(r_std),PARAMETER,DIMENSION(nstm) :: avan = & |
---|
| 96 | !!$ & (/ 0.036_r_std, 0.036_r_std, 0.036_r_std /) |
---|
| 97 | !TdO |
---|
| 98 | ! Van genuchten coefficient a (mm^{-1}) |
---|
| 99 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: avan = & |
---|
| 100 | & (/ 0.0075_r_std, 0.0036_r_std, 0.0019_r_std /) |
---|
| 101 | ! CWRR linearisation |
---|
| 102 | INTEGER(i_std),PARAMETER :: imin = 1 |
---|
| 103 | ! number of interval for CWRR |
---|
| 104 | INTEGER(i_std),PARAMETER :: nbint = 100 |
---|
| 105 | ! number of points for CWRR |
---|
| 106 | INTEGER(i_std),PARAMETER :: imax = nbint+1 |
---|
| 107 | ! Residual soil water content |
---|
| 108 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: mcr = & |
---|
| 109 | & (/ 0.065_r_std, 0.078_r_std, 0.095_r_std /) |
---|
| 110 | ! Saturated soil water content |
---|
| 111 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: mcs = & |
---|
| 112 | & (/ 0.41_r_std, 0.43_r_std, 0.41_r_std /) |
---|
| 113 | ! Total depth of soil reservoir (m) |
---|
| 114 | REAL(r_std),SAVE,DIMENSION(nstm) :: dpu = & |
---|
| 115 | & (/ 2.0_r_std, 2.0_r_std, 2.0_r_std /) |
---|
| 116 | ! dpu must be constant over the different soil types |
---|
| 117 | ! Hydraulic conductivity Saturation (mm/d) |
---|
| 118 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: ks = & |
---|
| 119 | & (/ 1060.8_r_std, 249.6_r_std, 62.4_r_std /) |
---|
| 120 | ! Soil moisture above which transpir is max |
---|
| 121 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: pcent = & |
---|
| 122 | & (/ 0.5_r_std, 0.5_r_std, 0.5_r_std /) |
---|
| 123 | ! Max value of the permeability coeff at the bottom of the soil |
---|
| 124 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: free_drain_max = & |
---|
| 125 | & (/ 1.0_r_std, 1.0_r_std, 1.0_r_std /) |
---|
| 126 | ! Volumetric water content field capacity |
---|
| 127 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: mcf = & |
---|
| 128 | & (/ 0.32_r_std, 0.32_r_std, 0.32_r_std /) |
---|
| 129 | ! Volumetric water content Wilting pt |
---|
| 130 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: mcw = & |
---|
| 131 | & (/ 0.10_r_std, 0.10_r_std, 0.10_r_std /) |
---|
| 132 | ! Vol. wat. cont. above which albedo is cst |
---|
| 133 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: mc_awet = & |
---|
| 134 | & (/ 0.25_r_std, 0.25_r_std, 0.25_r_std /) |
---|
| 135 | ! Vol. wat. cont. below which albedo is cst |
---|
| 136 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: mc_adry = & |
---|
| 137 | & (/ 0.1_r_std, 0.1_r_std, 0.1_r_std /) |
---|
| 138 | ! Matrix potential at saturation (mm) |
---|
| 139 | REAL(r_std),PARAMETER,DIMENSION(nstm) :: psis = & |
---|
| 140 | & (/ -300.0_r_std, -300.0_r_std, -300.0_r_std /) |
---|
| 141 | ! Time weighting for discretisation |
---|
| 142 | REAL(r_std),PARAMETER :: w_time = 1.0_r_std |
---|
| 143 | !- |
---|
| 144 | ! Diagnostic variables |
---|
| 145 | !- |
---|
| 146 | ! The lower limit of the layer on which soil moisture (relative) |
---|
| 147 | ! and temperature are going to be diagnosed. |
---|
| 148 | ! These variables are made for transfering the information |
---|
| 149 | ! to the biogeophyical processes modelled in STOMATE. |
---|
| 150 | !- |
---|
| 151 | REAL(r_std),DIMENSION(nbdl),SAVE :: diaglev |
---|
| 152 | !------------------------- |
---|
| 153 | END MODULE constantes_soil |
---|