A numerical exploration of hyporheic zone solute transport behavior estimated from electrical resistivity inversions

1). We run the Comsol model (ER3D_cs300ws6ds20_Zxl1.m) with an electrical push-pull loop, and then the voltage at every electrode of each time is printed to model_name.out and resistivity to model_name.ro. 

2). superposition.m is used to create final voltage results for later use in R2 inversions.

3). postproc_c_dual.m is used to create the pixel breakthrough curves of the simulated results.

4). plot_c_comsol.m read the resistivity result from step (1), convert the resistivity to bulk EC, and plot the bulk EC contour map. 

5). protocal_pc.m is used to take the voltages, add some random noise, and make the input files and batch scripts for R2.

6). runr2.bat and runr2_TL.bat are used to generate inverted bulk EC.

7). plot_C_r2_dual.m is used to plot the inversed bulk EC.