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

39. 39. Initialize parameters.

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


 * 1) define OH 0
 * 2) define Cl 1
 * 3) define S04 2
 * 4) define CO3 3
 * 5) define H 0
 * 6) define Ca 1
 * 7) define Na 2
 * 8) define K 3
 * 9) define TIME_LIMIT 0
 * 10) define STRUCT_LIMIT 1
 * 11) 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 */