The aim of this paper is to address computational aspects of the total energy singular vector analysis of atmospheric chemical transport models. We discuss the symmetry of the tangent-linear/adjoint operator for stiff systems. Results for a 3D simulation with real data reveal that the total energy singular vectors depend on the target domain, simulation window, chemical reactions, and meteorological data.