double
DIV(double divSI);
39
double
T(double tSI);
40
double
C(double cSI);
41
double
KH(double khSI);
42
double
QH(double qhSI);
43
};
44
45
46
class Tsoil{
47
public:
48
49
50
51
52
53
54
55
56
57
58
REDUCE R;
59
// saturated volumetric water content
60
double WSAT;
61
// lower limit of drying
62
double WD;
63
64
double KWSAT;
65
66
public:
67
// Constructor
68
Tsoil();
69
void InitWATER(double wsat, double wd, double ksat);
70
// calculates water content from degree of saturation
71
double W (double DSAT);
72
73
private:
74
// volume fraction of ice
75
double I;
76
// volume fraction of soil
77
double G;
78
// thermal conductivity of water
79
double KHW;
80
// thermal conductivity of ice
81
double KHI;
82
// thermal conductivity of soil grains
83
double KHG;
84
public:
85
void InitTherm(double khw, double khi, double khg);
86
87
double RKH (double W);
88
89
double RKHI();
90
91
protected:
92
14