147
double precision,
// check on Wd
148
resolution,
// layer thickness
149
LCMAX;
// cutoff number of layers
150
151
double DPHI;
152
void InitTol();
153
void FirstTol(double prec, double resol, double lcmax);
154
Ttol();
155
};
156
157
158 // variables used in calculations
159 class Trigidice:
160
public Tbnds,
161
public Ttol,
162
public REDUCE,
163
public Tsoil
164
{
165
public:
166
double
QW1, QW2,// reduced mass flux
167
DW,// reduced change in water content across layer
168
GRADUW1, GRADUW2,// reduced water pressure gradient
169
QH1, QH2,// reduced thermal flux
170
GRADTHETA1, GRADTHETA2,// reduced thermal gradient
171
PHI,// reduced phi
172
GRADPHI1, GRADPHI2,// reduced gradient of PHI from Clapeyron
173
Z,// reduced thickness
174
ZN,// reduce location of new lens
175
ZI,// reduce location of old lens
176
DZ,// reduced layer thickness
177
UW,// reduced water pressure
178
DUW,// reduced change in water pressure across layer
179
UI,// reduced ice pressure
180
CHI,//reduced stress partition factor
181
UN,// reduced neutral stress
182
UNold,// reduced previous neutral stress
183
UNmax,// reduced maximum neutral stress
184
THETA,// reduced temperature
185
DTHETA,// reduced temperature change across layer
186
WF,// unfrozen water in mature frozen zone
187
FSI;
188
189
Trigidice();
190
void LAYER(double rphi, double rdphi);
191
void PROFILE();
192
void GLOBAL();
193
194
double oldPHI;
195
double sum;
196
int iteration,
197
// error signal for GRADPHI < 0
198
GRADPROB,
199
// not enough layers
16