changed to match 2d_axisym better and to plot rho=0.5 isosurface.  Dimensional splitting is much faster