# -------------------------------------------------------------
# Bulk Configuration
# -------------------------------------------------------------

# Set up lattice
vector_a = [2.866, 0.0, 0.0]*Angstrom
vector_b = [0.0, 2.866, 0.0]*Angstrom
vector_c = [0.0, 0.0, 31.13975]*Angstrom
bravais_lattice = UnitCell(vector_a, vector_b, vector_c)

# Define elements
elements = [Iron, Iron, Iron, Iron, Iron, Iron, Magnesium, Oxygen, Oxygen,
            Magnesium, Magnesium, Oxygen, Oxygen, Magnesium, Magnesium, Oxygen,
            Oxygen, Magnesium, Iron, Iron, Iron, Iron, Iron, Iron]

# Define coordinates
cartesian_coordinates = [[  0.716500000003,   0.716500000003,   0.716499999992],
                         [  2.149500000003,   2.149500000003,   2.149500000008],
                         [  0.716500000003,   0.716500000003,   3.582499999992],
                         [  2.149500000003,   2.149500000003,   5.015500000008],
                         [  0.7165        ,   0.7165        ,   6.461176916299],
                         [  2.1495        ,   2.1495        ,   7.883348084818],
                         [  0.7165        ,   0.7165        ,  10.046913166035],
                         [  2.1495        ,   2.1495        ,  10.125854994682],
                         [  0.7165        ,   0.7165        ,  12.2742460054  ],
                         [  2.1495        ,   2.1495        ,  12.287565384199],
                         [  0.7165        ,   0.7165        ,  14.47251361208 ],
                         [  2.1495        ,   2.1495        ,  14.47501183558 ],
                         [  0.7165        ,   0.7165        ,  16.667197114   ],
                         [  2.1495        ,   2.1495        ,  16.669729913304],
                         [  0.7165        ,   0.7165        ,  18.854675912015],
                         [  2.1495        ,   2.1495        ,  18.867967064779],
                         [  0.7165        ,   0.7165        ,  21.016369375396],
                         [  2.1495        ,   2.1495        ,  21.095331583047],
                         [  0.7165        ,   0.7165        ,  23.258901803911],
                         [  2.1495        ,   2.1495        ,  24.681132523465],
                         [  0.716500000003,   0.716500000003,  26.124249999992],
                         [  2.149500000003,   2.149500000003,  27.557250000008],
                         [  0.716500000003,   0.716500000003,  28.990249999992],
                         [  2.149500000003,   2.149500000003,  30.423250000008]]*Angstrom

# Set up configuration
bulk_configuration = BulkConfiguration(
    bravais_lattice=bravais_lattice,
    elements=elements,
    cartesian_coordinates=cartesian_coordinates,
    )