Sorry, I misread that somehow.  Then you should do pq.Quantity(subdata[:, n], unit)