APPENDIX A: RIGIDICE MODEL
RIGIDICE.H
#define grav ((const double) 9.800)
//
m/sec**2
// density of liquid water
#define rhoW ((const double) 1.000E03)
//
kg/m**3
// density of ice
#define rhoI ((const double) 0.917E03)
//
kg/m**3
// viscosity of water
#define nu ((const double) 1.787E-3)
//
kg/m*s
// specific surface free energy of a ice-air interface
0
0 #define gammaIA ((const double) 0.100)
//
n/m
1
// Latent heat of fusion per unit volume of melt
2
#define hIW ((const double) 3.335E05)
//
j/m**3
3
// standard melting point of bulk ice, or its analogue, when exposed
4
// to air at standard atmospheric pressure
5
#define theta0 ((const double) 273.15)
//
kelvins
6
// conversions
7
// mm/day / (m/s)
8
#define mmpdTomps ((const double) 8.64e7)
9
// kPa to Pa
10
#define kP2P ((const double) 1.00e3)
11
12
13
class REDUCE{
14
protected:
15
double
orglamda,
16
orgeta,
17
orgFW;
18
public:
19
// Constructor
20
REDUCE();
21
void InitScales();
22
void FirstScales(double micro, double macro, double body);
23
double
lamda;
24
double
eta;
25
//
double
body;
26
double FW;
27
double Y;
28
double L(double rSI);
29
double P(double pSI) ;
30
// reduced latent heat of fusion, per unit volume of melt
31
double H;
32
double RTHETA(double thetaSI);
33
double SIGMA(double sigmaSI);
34
double GRAD(double gradSI);
35
double F(double fSI);
36
double KW(double kwSI);
37
double V(double vSI);
38
13