import pulsar as psr
def load_ref_system():
    """ Returns d-sorbose as found in the IQMol fragment library.
        All credit to https://github.com/nutjunkie/IQmol
    """
    return psr.make_system("""
        C     -0.6212     -0.4667     -0.1735
        C     -2.0028      0.0205      0.2912
        C     -3.1049     -1.0304      0.2866
        H     -0.7389     -0.7757     -1.2437
        O     -0.3134     -1.6108      0.5918
        C      0.4916      0.6223     -0.1063
        O     -2.2064      1.1707      0.6305
        O     -4.4026     -0.5162      0.2122
        H     -2.9817     -1.6755      1.1803
        H     -3.0372     -1.6745     -0.6082
        C      1.7947      0.0731     -0.7576
        O      0.8195      0.9935      1.2141
        H      0.1528      1.5221     -0.6756
        H      1.5565     -0.3364     -1.7684
        C      2.8610      1.1864     -0.8634
        O      2.2501     -1.0784     -0.0783
        O      4.0789      0.7420     -1.3884
        H      3.1509      1.5835      0.1266
        H      2.4597      2.0260     -1.4637
        H      3.9035      0.3345     -2.2278
        H      2.7239     -0.7862      0.6928
        H     -4.5040      0.1298      0.9011
        H      0.6288     -1.7580      0.4801
        H      0.0271      1.3312      1.6105
        """)
