Date: Tue, 24 Sep 2002 12:28:09 -0700 (PDT)
From: Thomas Cheatham
Subject: Re: watershell in ptraj


> I have a question about the watershell command in ptraj. The manual
> says: "...this option will count the number of waters within a certain
> distance of the atoms in the mask in order to represent the first and
> second solvation shells..." Now my question is when is a "water"
> molecule counted as beeing in the cutoff region? Is the center of mass,
> the oxygen or something else taken.

This is a good question since I failed to properly document this. It
looks like what I do is mark the solvent as being in the shell if any
atom of the solvent is within the shell. There is no center of mass
option. However, there is a solvent atom selection, so with water if you
only wanted to know if the oxygen was within the shell, then you would
specify :WAT@O as the solventmask. This routine ASSUMES that you have
solvent in the prmtop or have specified which residues are solvent with
the solvent command in ptraj. Let me know if this does not properly
clarify the functionality.