def METHOD_NAME(self):
    command = self.command + ' -s ca --pdb {1:s} {0:s} {0:s}'.format(
                            self.dcdpath, self.pdbpath)
    namespace = prody_parser.parse_args(shlex.split(command))
    namespace.func(namespace)
    select = self.ag.ca
    assert_equal(select.numAtoms(), 10)
    coords = self.dcd[:]
    coords.setAtoms(select)
    coords = coords._getCoordsets()
    concat = parseDCD(self.output)
    assert_equal(concat.numAtoms(), select.numAtoms())
    concat = concat._getCoordsets()
    assert_equal(select.numAtoms(), coords.shape[1])
    assert_equal(select.numAtoms(), concat.shape[1])
    assert_equal(coords, concat[:3])
    assert_equal(coords, concat[3:])