# Set up lattice
vector_a = [2.7411436080585054, 0.0, 0.0]*Angstrom
vector_b = [-0.913714536019502, 2.584374977952439, 0.0]*Angstrom
vector_c = [0.0, 0.0, 67.1433438381166]*Angstrom
lattice = UnitCell(vector_a, vector_b, vector_c)

# Define elements
elements = [Tungsten, Tungsten, Tungsten, Tungsten, Tungsten, Tungsten,
            Tungsten, Tungsten, Iron, Iron]

# Define coordinates
fractional_coordinates = [[ 0.5           ,  0.5           ,  0.366663811725],
                          [ 0.            ,  0.            ,  0.399997486456],
                          [ 0.5           ,  0.5           ,  0.433331161186],
                          [ 0.            ,  0.            ,  0.466664835917],
                          [ 0.5           ,  0.500000000001,  0.500502858333],
                          [ 0.000000000001,  0.000000000002,  0.534153685787],
                          [ 0.5           ,  0.500000000001,  0.567805820039],
                          [-0.000000000001, -0.000000000001,  0.601344961788],
                          [ 0.5           ,  0.499999999999,  0.631395384651],
                          [ 0.000000000001,  0.            ,  0.657049211718]]

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