r = pybamm.SpatialVariable("r", ["current collector"], coord_sys="cylindrical polar") theta = pybamm.SpatialVariable( "theta", ["current collector"], coord_sys ...