if you are using bash (or probably any other POSIX compliant shell more
sophisticated than sh), rename your PDB-file to myfile001.pdb and you should be
able to use something like this:
for i in $(seq 100); do
pdbset xyzin myfile$(printf "%03d" $i).pdb \
xyzout myfile$(printf "%03d" $((i + 1)) ) << eof
shift 10 0 0
rotate matrix 0.9848 -0.1736 0 0.1736 0.9848 0 0 0 1
Where the entries of the rotation matrix are explained e.g. on wikipedia (this
example rotates about the z-axis and shifts along x).
You could also use moleman2 instead.
On Fri, Feb 18, 2011 at 03:54:58PM +0100, Sabine Schneider wrote:
> I want to apply to a small structure a rotation and translation (i.e
> 10A translation and 10dg rotation), write out the coordinates and
> apply the same rotation/translation to the new coordinates and
> repeat this lets say 100 times.
> I guess the program of choice would be pdbset but I am not sure how
> to script it to do the repetition?
> Thanks a lot for your help!
Institut fuer anorganische Chemie
phone: +49 (0)551 39 22149
GPG Key ID = A46BEE1A