Subprogram modesroutinesmodule::modalmasses (sups, modes, sam, eigVectors, normFactors, modeOrder, ierr): It works now for LANCZ2 and for the -undamped solvers (although I encountered something that looks like a bug in the -undamped calculation of the eigenvectors). Furthermore, local orientation of triads are not (yet) taken into account, but should normally not be critical as long as the local triad masses are the same in all coordinate directions (for inertias it will induce errors, however). Component modes seems to be treated correctly, but has not been tested too much. Some more elegance could also be put into the creation of unit displacement vectors, just to make sure that all kinds of joints and boundary conditions are treated correctly. In fact, there might be a better and easier solution to this than the approach I have used here.