The vibronic structure of the B system in the two-photon ionization spectrum of Na3 is investigated theoretically, based on an earlier coupling scheme with three interacting potential energy surfaces (pseudo Jahn-Teller coupling). This is extended in the present work to allow for additional ab initio data in the modeling, to include the totally symmetric vibrational mode and to treat a higher excited state in the electronic manifold. Important features of the experimental recording can thus be reproduced, although some unsatisfactory aspects remain. The implications for the nonadiabatic nature of the underlying nuclear motion, as well as directions for future work, are discussed.