Changed r=0 boundary back to zero-gradient as reflecting needs a bit of work. Need to ensure density along the z-axis behaves physically.