Hi Robert,
> select x > x1 and x < x2 and y > y1 and y < y2 and z > z1 and z < z2
This is how I did it.
> boundbox @{point(x1, y1, z1)} @{point(x2, y2, z2)}
> select within(boundbox)
> boundbox {0 0 0} {10 10 10}
> select within(boundbox)
> select within(group, within(boundbox))
> select not within(group, not within(boundbox))
I never used the within() function before. I am so glad that you introduced it
to me!
> Usually when people label a residue, they just label the alpha carbon:
> select *.ca and not within(group, not within(boundbox))
I see. But let's say there's a residue which is partly in the box, but its
alpha carbon is not.
How can I make sure that it's labeled? I tried this:
[javascript]
...... #only those atoms inside the box are visible.
no1 = {visible}.resno.min;
no2 = {visible}.resno.max;
for (i=no1; i<no2+1; i++) {
select {visible and resno=i}[1]; label %r;
}
[/javascript]
It turns out that most residues are labeled correctly, but some partial
residues are not
Do you know why? Or is there a better way of doing it?
> delete file=2
I see.
> You are asking great questions -- and look -- you found a bug! Thank you!
Thank you!
Yingjie
------------------------------------------------------------------------------
Enable your software for Intel(R) Active Management Technology to meet the
growing manageability and security demands of your customers. Businesses
are taking advantage of Intel(R) vPro (TM) technology - will your software
be a part of the solution? Download the Intel(R) Manageability Checker
today! http://p.sf.net/sfu/intel-dev2devmar
_______________________________________________
Jmol-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/jmol-users
