A new method for the treatment of data from multidimensional solid-state NMR investigations is described. It approximates the theoretical NMR chemical shifts from the chemical shielding values obtained by first-principles calculations and subsequently treats these results to quantify the similarity between predicted and experimental chemical shift correlations. The test case of this approach is performed for the measured and several sets of computed C-13-H-1 heteronuclear correlations in the polymorphic form I of metergoline, which is relatively large, pharmaceutically active system. The proposed protocol is general, however, and it can be immediately applied to study other compounds and nuclei. (C) 2013 Elsevier B.V. All rights reserved.