Evaluate the spherical solution on uniform grid

[liu-ziyuan-math]

Dear authors,

I have come up with a need for solutions on uniform longitude and colatitude. I wonder if a function is already available in dedalus or scipy, to evaluate the solution on uniform grids with/without poles, let's say, for the S2 shallow-water example. It seems pretty like a scipy.special.eval_jacobi function.

Thank you!