vmdを開き、File --> New Molecule ... よりpsf, dcdを読ませる。
Extensions --> Tk Console
にて、
set output [open rmsd.dat w]
set numf [molinfo top get numframes]
set ref [atomselect top "protein and noh" frame 0]
set sel [atomselect top "protein and noh"]
for { set i 1 } { $i < $numf } { incr i } {
$sel frame $i
$sel move [measure fit $sel $ref]
puts $output "[measure rmsd $sel $ref]"
}
close $output
でrmsd.datができる。
Friday, July 03, 2009
Wednesday, July 01, 2009
Tuesday, June 02, 2009
Monday, June 01, 2009
ssh/config覚え書き
普段サーバ上でscreenしているので気づかなかったけど、sshで繋いだ時にタイムアウトしないように信号を送る設定。
ServerAliveInterval 120
みたいにすると二分に一回信号を送るので切れない。
ServerAliveInterval 120
みたいにすると二分に一回信号を送るので切れない。
Subscribe to:
Posts (Atom)