# Load the monomers into memory
methyl_acrylate = nlread('Methyl_Acrylate.hdf5', class_type=MoleculeConfiguration)[-1]
propylene = nlread('Propylene.hdf5', class_type=MoleculeConfiguration)[-1]

# Create the polymer sequence with two polymer blocks
# This is done by setting the monomer_additions to enforce two different blocks.
# To make an A-B block copolymer the sum of the block sizes need to be as long as
# the polymer chain
sequence = PolymerSequence(
    [methyl_acrylate, propylene],
    number_of_monomers=40,
    number_of_chains=20,
    monomer_additions=[[0, 1], [1, 0]],
    monomer_block_sizes=[10, 30],
)