Implement command axial_avg for averaging over the cylindrical phi direction