What I want to do is to draw cluster structures and write them to cor files, but if I write like this,
imag byre sele all end xcen ?XAVE ycen ?YAVE zcen ?ZAVE
upda cutim 15
dele atom sele .byres. (.not.(ires @aimres .around. 10)) end
write coor card name ../../res-cor/@email@example.com
I will lose residues across the range 10, which is not expected.
How should I do to select?
Or can I just write them use `ires @aimres .around. 10` ?
also, the step `upda cutim 15` seems very time consuming, so are there any other ways more efficient or elegant to achieve my goal ?