<div><div><div>I'm modeling a water system which contains 27 molecules on QE V 6.1 and i need to get a center of centers of orbitals for each molecule. </div><div>So cp.x outputs 4 centers of orbitals for each H2O molecule in .wfc file. If i take some kind of 'center of mass' of those i get a Wannier centroid. But the problem is that i cant find logic in what order wfcs are printing.</div><div>I counted manually which molecule possesses orbitals and here is what i got</div><div>I:    41, 76, 86, 100</div><div>II:    21, 42, 77, 99</div><div>III:    14, 39, 61, 85</div><div>IV:    30, 38, 68, 101</div><div>V:    44, 59, 69, 79</div><div>VI:    52, 62, 74, 92</div><div>VII:    15, 34, 80, 88</div><div>VIII:    11, 13, 36, 96</div><div>IX:    55, 67, 71, 78</div><div>X:    7, 32, 47, 81</div><div>XI:    12, 46, 98, 108</div><div>XII:    43, 54, 89, 91</div><div>XIII:    4, 23, 83, 106</div><div>XIV:    19, 50, 84, 87</div><div>XV:    40, 56, 63, 73</div><div>XVI:    10, 18, 65, 93</div><div>XVII:    58, 66, 70, 72</div><div>XVIII:    3, 25, 51, 103</div><div>XIX:    16, 33, 37, 104</div><div>XX:    8, 31, 64, 95</div><div>XXI:    2, 6, 22, 107</div><div>XXII:    27, 35, 60, 90</div><div>XXIII:    1, 5, 9, 28</div><div>XXIV:    26, 57, 97, 102</div><div>XV:    29, 45, 75, 105</div><div>XXVI:    17, 20, 49, 53</div><div>XXVII:    24, 48, 82, 94</div><div> </div><div>One should read that as "Nth molecule: lines in output file with corresponding wfcs".</div><div> </div><div>Here is the question. How can i get output in a such convenient way, so i don't need to search manually for which line with wfc corresponds to which molecule in a system with arbitrary number of molecules or some completely different system? Or if there is a logic in the .wfc output please point it out.</div><div> </div><div>(Don't mind my previous message. I just misclicked.)</div></div></div>