7 Post processing

Post-processing steps can be called in order to automate the treatment and display of the results. The pdepe function returns the solution as the array sol, where ui = sol(:,:,i) contains the concentrations of chemical species i. The elements ui(j,k) = sol(j,k,i) represents the numerical solution ui at (t,x)= (tspan(j),xmesh(k)). MATSEDLAB can plot the computed concentration depth profiles of the pore water and solid-bound chemical species against measured values. The user must save an Excel® file named ‘FIELD_DATA.xls’ containing the measured concentrations and depths in the working directory. This file should consist of L individual worksheets, each corresponding to one of the chemical species in the same order as in MATSEDLAB. If field data are present, the function result_plot.m , called in the MATSEDLAB script, will automatically produce depth-concentration plots on which measured concentrations are displayed together with the computed profiles at the end of the simulation. In addition, result_plot.m will create a series of associated figures in the working directory. The post-processing capabilities are further details in the manual and illustrated in the applications presented below.