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

73. OUTPUT.

The following are output routines.

74. Print the state of the system for the first 4 concentraction surfaces.

(Print initial system state 74 ) =
 printf("\n\nInitial state of system:\n\n");
 printf("\n" I0N "\t\"EXTERNAL\"\t\" INTERNAL\"\n");
 for (i = 0; i < num.cations i++) {
   if (cation [i].c[FIRST_CELL] != 0 || cation [i].c[FIRST_CELL+ 1] != 0) { 
     printf ("\"%4s:\"\t",cation[i] cation[i].name)
     for (k = FIRST.CELL; k < FIRST_CELL + 2; k++) printf ("%8 5l\t", cation[i].c[k]);
     printf (" n");
   }
  }
  for (j=0; j < num.amons j++) {
    if (anion [j].c[FIRST_CELL]!=0 || anion [j].c[FIRST_CELL +1]!=0) {
      printf("\"%4s:\"\t", anion[j].name)
    for (k = FIRST_CELL; k < FIRST_CELL+2; k++) printf("%8.5lf\t",anion[j].c[k]);
    printf ("\n");
    }
  }
  printf/("\"  pH:\"\t\");
  for (k = FIRST_CELL; k < FIRST_CELL + 2; k++) print/ ("%8.5lf\t",-log10(cation[H].c[k}));
  printf("\n");
  printf("\n");

This  code is used in section 7.

75. Print header for permeability, sulfate penetration, and chloride penetration depths.

⟨Print header 75⟩
  printf("\"Day\"\t\"\L\"\t\"K  \"\t\"\D    "\t\"S04\"\t\"Cl\"\t\"Flux\"\t  \"pH\"\n");
  printf ("\n"   "\t\"\m"\t\"m/s\"\t\"\m^2/s\"\t\"\m  \"\t\"\m \"\t\"ml/dy/m2\"\n");
  printf ("\n");

This code is used in section 2.

Print intermediate results for permeability, sulfate and chloride penetration depths. The bulk values of permeability and diffusivity must be calculated here.

⟨Print intermediate results 76⟩ $\equiv$
 if ( Day = 0.0 V Day > print_day ) {
 print-day += print_day_interval
 printf ("%ld t%6.3lf",(long) Day ($\delta X$[FIRST_CELL] + NUM_CELLS - FIRST.CELL -  1) * L)
 $k_{B}$= $\delta X$[FIRST_CELL] * CUB($\vartheta _{0}<math>/<math>\xi$[FIRST_CELL]/$\kappa$[FIRST_CELL]
 $D_{b}$ = $\delta X$ [FIRST.CELL] * ($\vartheta _{0}<$/$\xi$[FIRST_CELL]);
 for (k = FIRST.CELL + 1; k < NUM.CELLS; k++) {
   $k_{B}$ += CUB($\vartheta _{0}$/$\xi$[k])/$\kappa$[k];
   $d_{B}$+= $\vartheta _{0}$//$\xi$[k];
 }
 $k_{B}$ = ($\delta X$[FIRST_CELL] + NUM.CELLS - FIRST.CELL - 1) * $k_{0}$/<math>k_B;