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

61.ADVANCE SULFATE FRONT '''61. Advance sulfate front.'''

⟨Function declarations 23)⟩$$+\equiv$$ void advance.sulfate-front

62. void advance_sulfate_front ; {   int i, j, k, new_first_cell; sulfate_depth+=sulfate_rate * $$\Delta T$$; if (sulfate_depth>(real) (FIRST_CELL+1) $$\lor\Delta X$$[FIRST_CELL] <0.10) { new_first-cell=(int)(sulfate-depth +0.4);    /* penetration beyond next cell */ /* copy info from old first cell into new first cell */ for(k = FIRST_CELL + 1; k $$\le$$ new_first_cell; k++) { for(i = 0; i < num_cations; i++) cation[i].c[k] = cation[i].c[FIRST_CELL]; for(j = 0; j < num_anions; j++) anion[j].c[k] = anion[j].c[FIRST_CELL]; for (i = 0; i < num_cations; i++) for (j = 0; j < num_anions; j++) sol_array[i][j].s[k] = sol_array[i][j].s[FIRST_CELL]; }    $$\Psi$$[new_first_cell]=$$\Psi$$[FIRST_CELL]; FIRST_CELL = new_first_cell; $$\Delta X$$[FIRST_CELL] = (real) (FIRST_CELL+1)-sulfate_depth; } }