Page:4SIGHT manual- a computer program for modelling degradation of underground low level waste concrete vaults (IA 4sightmanualcomp5612snyd).pdf/111

39. Initialize parameters.

The scale and physical parameters must be set based upon the information from the parsing and ion initialization routines.

#define OH 0
#define Cl 1
#define S04 2
#define CO3 3
#define H 0 
#define Ca 1
#define Na 2
#define K 3
#define TIME_LIMIT 0
#define STRUCT_LIMIT 1
#define PRINTOUTS 20.0

⟨Function declarations 23⟩ $\equiv$
 void initialize-parameters ();

40.

⟨Global variables 4⟩ $+\equiv$
 real sample_length = 1.0;    /* length of sample in meters */
 real sulfate_failure_year
 real chloride_failure_year
 real k[NUM_CELLS];     /* perm, coefficient due to cracking */
 real sulfate_rate;    /* dimensionless rate */
 real sulfate_rate_ms     /* m/sec */
 real sulfate_depth = 0.0;   /* depth of sulfate penetration */
 real chloride_depth = 0.0;   /* depth of chloride penetration */
 real rebar_depth = 1.0;
 boolean termination = FALSE;     /* simulation termination */
 int termination_type= TIME.LIMIT;     /* limiting condition for termination */
 MATERIAL MaxDay = {10000.0,TRUE};     /* time limit (day) */
 real MaxDepth = 100.0;      /* penetration depth limit (m) */
 MATERIAL crack.width =  {0.000100, TRUE};
 real crackspacing=100.0;
 real crack_depth= 0.0;
 real crack_permeability;
 real joint_width = 0.0;
 real joint_spacing = 100.0;
 real joint_lifetime = 20000.0;
 boolean joint_change_flag = FALSE;
 boolean joint_is_specified = FALSE;
 real joint_permeability = 0.0;
 real Pout    /* dimensionless ext. pressure */
 real dP     /* dimensionless rho g h pressure */
 real v[NUM_CELLS];     /* Darcy velocity */
 real $\mu$ = 0.001;    /* pore fluid viscosity */
 real g = 9.8;    /* gravitational constant(m/sec/sec) */
 real $\rho$ = 1000;    /* density of pore fluid (kg/m^3) */ 
 real $\kappa _{B}$    /* bulk permeability */
 real $D_{B}$    /* bulk diffusivity */