#! /usr/bin/env python

import numpy
import sasmol.sasmol as sasmol

m = sasmol.SasMol(0)
m.read_pdb('3FAD.pdb')

mask = numpy.array([1]*m.natoms())

copy_m = sasmol.SasMol(0)
error = m.copy_molecule_using_mask(copy_m,mask,0)

double_m = sasmol.SasMol(0)
error = double_m.merge_two_molecules(m,copy_m)

double_m.segnames()