204
VI = heave;
205
VIorg = R.V(VI/mmpdTomps);
206
VI = VIorg;
207
VB = pene;
208
VBorg = R.V(VB/mmpdTomps);
209
VB = VBorg;
210
GRADTB = gradtb;
211
GRADTBorg = R.GRAD(R.RTHETA(gradtb));
212
GRADTB = GRADTBorg;
213
}
214
215
void Tbnds::InitQuasi(double press, double pene, double gradtb){
216
217
UIorg = R.P(press*kP2P);
218
VBorg = R.V(pene/mmpdTomps);
219
GRADTBorg = R.GRAD(R.RTHETA(gradtb));
220
}
221
222
Ttol::Ttol(){
223
FirstTol(0.1,0.1,100);
224
}
225
226
void Ttol::FirstTol(double prec, double resol, double lcmax){
227
orgprecision = prec;
228
orgresolution = resol;
229
LCMAX= lcmax;
230
InitTol();
231
}
232
void Ttol::InitTol(){
233
precision = orgprecision;
234
resolution = orgresolution;
235
}
236
237
238
Trigidice::Trigidice(){
239
// specific gravity of ice?
240
//
Y=rhoI/rhoW;
241
//
lamda=1.0e-6;
242
//
eta=0.01;
243
//
FW = 1.0;
244
//
FW = F(FW*(-rhoW*grav));
245
// reduced latent heat of fusion, per unit volume of melt
246
//
H =rhoW*(lamda/gammaIA)*hIW;
247
InitTCALC();
248
}
249
250
void Trigidice::InitTCALC(){
251
QW1=0;
252
QW2=0;// reduced mass flux
253
DW=0; // change in water content across layer
254
GRADUW1=0;
255
GRADUW2=0;// reduced water pressure gradient
256
QH1=0;
257
22
Previous Page
