Cut l<m for power spectrum estimation.

Now the first part, i.e., runprod,  works fine.