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

⟨Propagate ion concentration in concrete 15⟩ $\equiv$

 for (i = 1; i < num_cations; i++} {
 cation [i].moles [0] = cation [i].c[0] * litre [0];
 cation [i].moles[1] = cation[i].c[1] * litre[1];
 for (k = 2; k < NUM_SURFACES; k++) {
   cation [i].c[k] =cation [i].c[1];
   cation [i].moles[k] = cation [i].c[k] + litre [k];
 }
 }
 for(j = 1 ; j < num_anions; j++) {
   anion [j].moles [0] = anion [j].c[0]* litre [0];
    anion [j].moles [1] = anion [i].c[1] * litre [1];
   for (k = 2; k< NUM_SURFACES; k++) {
     anion [j].c[k]=anion [j].c[1];
     anion [j].moles[k] = anion [j].c[k] * litre[k];
   }
 }

This code is used in section 41,

46. Report material parameters to stdout.

⟨Preprocessor definitions 11⟩ $+\equiv$
#define GET_STAT(a) ((a.is_default = TRUE) ? "DEFAULT": "USER")

47.

 ⟨Print parameters to stdout 47⟩ $\equiv$
 printf("\"THICKNESS\"\t%8.51f\t\" (m) uuuuuu\"\t\"%s\"\n", sample_length, GET_STAT(thickness));
 printf("\"DIFFuuuuu\"\t%8. 11e\t\" (m^2/sec)\"\t\"%s\"\n", Do, GET_STAT(D));
 printf("\"PERMUuuuuu\"\t%8.11e\t\"(m/sec) uu\"\t\"%s\"\n", 9.8 - 10+06 * ko, GET_STAT(ko));
 printf("\"WCuuuuuuuu\"\t%8.51f\t\"Uuuuuuuuu\"\t\"%s\"\n", wc.value, GET_STAT(wc));
 printf("\"HEADUuuuuu\"\t%8.51f\t\"(m)uuuuuu\"\t\"%s\"\n", GET_STAT(head));
 printf("Sulfate Attack Parameters: \n");  
 printf("\"YOUNGSµ\"\t%8.11e\t\"(N/m^2}uu\"\t\"%s\"\n",  Youngs.value, GET_STAT(Youngs));
 printf("\"ROUGENESS\"\t%8.51f\t\"uuuuuuuuu\"\t\"%s\"\n", roughness.value, GET_STAT(roughness));
 printf("\"GAMMAuuuuu\"\t%8.51f\t\" (J/m^2) uu\"\t\"%s\"\n", $\gamma$.value, GET_STAT($\gamma$);
 printf("\"BETAuuuuuu\"\t%8. 11e\t\"uuuuuuuuu\"\t\"%s\"\n", $\beta$.value, GET_STAT($\beta$));
 printf("\"CEUuuuu\"\t%8.51f\t\" (Mol/m^3}\"\t\"%s\"\n", CE. value, GET_STAT(CE));
 printf("\"POISSON\"\t%8.51f\t\"uuuuuuuuu\"\t\"%s\"\n", $\nu$ .value, GET STAT($\nu$ ));
 printf("\n");
 if (MaxDepth < sample_length}
    printf("\"DEPTHuuuuu\"\t%8.5lf\t\"(day)    \"\t\"%s\"\n", Max Depth, "USER");
 if (rebar.depth < sample_length)
    printf("\"REBARuuuuu\"\t%8.51f\t\"(m)uuuuuu\"\t\"%s\"\n", rebar_depth, "USER");
 printf("\"TIMEUuuuuu\"\tu%8.01f\t\"(day) uuuu\"\t\"%s\"\n", MaxDay.value, GET_STAT(Max Day));
 if (crack_width.is_default = TRUE A neutral_axis_depth.is_default = TRUE) {
    crack_depth = 0.25 * sample_length;