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.