Only implemented for spherical1d.
spherical1d
Works somehow, but results are not comparable to FV+TVB method.