Read before you write. Corrected dimensions of work vector, now tb01pd is
functioning