96
// lower limit of drying
97
//
WD = 0.02;
98
WD = wd;
99
// saturated hydraulic conductivity
100
//
KWSAT=1.0e-8;
// m/s
101
KWSAT = ksat;
102
KWSAT=KWSAT/(rhoW*grav);
103
KWSAT = R.KW(KWSAT);
104
}
105
106
// given degree of saturation
107
double Tsoil::W(double DSAT){
108
return( DSAT*(WSAT-WD)+WD);
109
}
110
//
111
//
112
//
Thermal conductivity
113
114
void Tsoil::InitTherm(double khw, double khi, double khg){
115
// volume fraction of soil
116
G=1-WSAT;
117
// volume fraction of ice
118
I=1-G-WD;
119
120
//
KHW =0.52; // W/m
121
KHW = khw;
122
KHW = R.KH(KHW);
123
124
//
KHI =2.32; // W/m
125
KHI = khi;
126
KHI = R.KH(KHI);
127
128
//
KHG=3.42; // W/m
129
KHG = khg;
130
KHG = R.KH(KHG);
131
}
132
133
134
double Tsoil::RKH (double W){
135
I = 1 - G - W;
136
return ( pow(KHG,G)*pow(KHW,W)*pow(KHI,I));
137
}
138
139
// to get KHI
140
double Tsoil::RKHI(){
141
return (KHI);
142
}
143
//
144
145
//
146
//
Brooks & Corey equations
147
//
148
void Tsoil::InitBC(double alpha, double beta, double phib) {
149
//
ALPHA = 0.36;
20