93
// unfrozen water content exponent
94
double ALPHA;
95
private:
96
// frozen hydraulic conductivity exponent
97
double BETA;
98
protected:
99
// reduced ice intrusion pressure
100
double PHIB;
101
public:
102
void InitBC(double alpha, double beta, double phib);
103
double DSAT(double PHI);
104
double RKW(double PHI);
105
106
// running sum
107
private:
108
double PROD;
109
int FlagSny;
110
public:
111
// initialize starting parameters
112
void InitStress();
113
// calculates Snyder's factor
114
double SNYDER (double PHI, double DSAT, double deltaDSAT);
115
};
116
117
118
// boundary conditions
119
class Tbnds
120
{
121
public: REDUCE R;
122
protected:
123
// heave rate
124
double VI, VIorg,
125
// penetration rate
126
VB, VBorg,
127
// unfrozen temperature gradient
128
GRADTB, GRADTBorg,
129
// sought after heave pressure
130
UIorg;
131
132
public:
133
void InitQuasi(double press, double pene, double gradtb);
134
void InitContinuum(double heave, double pene, double gradtb);
135
Tbnds();
136
};
137
138
139
// variables used for numeric precision
140
class Ttol
141
{
142
protected:
143
double
orgresolution,
144
orgprecision;
145
146
public:
15
Previous Page
