Dcdreporter subset

Addresses issue #18

Here's the use case:

```
[...]
atomIndices = [atom.index for atom in simulation.topology.atoms() if atom.residue.name != 'HOH']
simulation.reporters.append(DCDReporter('dry.dcd', 1000, atomIndices=atomIndices))
[...]
```