RE: AMBER: angle between two vectors. .

From: Lwin, ThuZar <ThuZar.Lwin.stjude.org>
Date: Mon, 16 Oct 2006 15:51:28 -0500

Hi,
   Thank you very much to all who have responded.
Thanks,
ThuZar


You can use tcl script in VMD.
You may want to modify the following tcl script.
Enjoy it.


+++++++++++++++++++++++++++++++++++++++
set outfile [open cn-m2-vmd.dat w]
set nf [molinfo top get numframes]
set g34 [atomselect top "resname GLY and name CA"]
set v27 [atomselect top "resid 7 34 61 88 and name CA"]
set c5 [atomselect top "resname HMT and name C5"]
set n1 [atomselect top "resname HMT and name N1"]
set all [atomselect top all]
set conv [expr 180.0/acos(-1.0)]

# angle calculation loop
for { set i 0 } { $i < $nf } { incr i } {
  $g34 frame $i
  $v27 frame $i
  set com1 [measure center $g34]
  set com2 [measure center $v27]
  set m2 [vecsub $com2 $com1]
  set m2mag [veclength $m2]

  $c5 frame $i
  $n1 frame $i
  set com3 [measure center $c5]
  set com4 [measure center $n1]
  set cn [vecsub $com4 $com3]
  set cnmag [veclength $cn]

  set dotp [vecdot $m2 $cn]
  set theta [expr $conv * acos($dotp / ($m2mag * $cnmag))]
  if { $theta > 90.0 } {
    set theta [expr 180.0 - $theta]
  }
  puts $outfile "[expr ($i + 1) * 10] $theta"
}
close $outfile
++++++++++++++++++++++++++++++++++++++++++




On 10/12/06, David A. Case <case.scripps.edu> wrote:
> On Thu, Oct 12, 2006, Lwin, ThuZar wrote:
>
> > Do you have any suggestion on how to do the following using
amber
> > ptraj?
> > Vector 1: Assign a vector between atom A and B.
> > Vector 2: Assign another vector between atom C and D.
> > Then calculate angle between vector 1 and vector 2.
>
> My best guess is to use the vector command to write out files with the
two
> vectors as a function of time, then use a short awk or perl script to
compute
> the angles.
>
> (This would be a good functionality to have inside ptraj, but as far
as I
> know, it is not there now.)
>

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Oct 18 2006 - 06:07:19 PDT
Custom Search