Possible fix to broken large tables reported by Alec -- switching to 64-bit indexing of variable length arrays.

Author: Cameron-Van-Eck

polspectra.py

@@ -456,7 +456,7 @@ def write_FITS(self,filename,overwrite=False):
                 col_format=pf.column._convert_record2fits(tabcol.dtype)
             else: #Channelized columns
                 subtype=np.result_type(np.array(tabcol[0])) #get the type of each element in 2D array
-                col_format='P'+pf.column._convert_record2fits(subtype)+'()'
+                col_format='Q'+pf.column._convert_record2fits(subtype)+'()'
             if tabcol.unit != None:
                 unit=tabcol.unit.to_string()
             else: