Dear Sabine, 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:
#!/bin/bash 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 end eof 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. Cheers, Tim On Fri, Feb 18, 2011 at 03:54:58PM +0100, Sabine Schneider wrote: > Hello, > > 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! > > Sabine -- -- Tim Gruene Institut fuer anorganische Chemie Tammannstr. 4 D-37077 Goettingen phone: +49 (0)551 39 22149 GPG Key ID = A46BEE1A
signature.asc
Description: Digital signature