Added functionality to smooth outputted rho
Changed computation of r1 and r2 to standard function
Fixed a bug in charge scaling
Added pair creation condition of B field strength
Fixed flux function computation bug
Fixed Dipole initial condition