for rxn in list(submodel.reactions):
                    submodel.reactions.remove(rxn)
                    merged_submodel.reactions.append(rxn)