We should probably check when ifheat=T to make sure that p8 .ne. 0

Otherwise, it messes up the velocity field... :(

I'm unclear what exactly the pmtm function is outputting? It returns three np.arrays and the documentation doesn't say what it returns.

I'm not sure what to do to those output arrays to get the actual multitapered output and produce a plot similar to what is returned from the show=True pa

I am looking at the docs and see that Function and Array also needs to be imported for the example to execute.

from shenfun import Basis, Function, Array
N = 8
T = Basis(N, 'Chebyshev', bc=None)

Also, later on in the code right before Operators you need to import Dx for this code to work.

From the Multipdimensional part, one needs to import TensorProductSpace, define comm and also i