| 6 | 1/1 | 返回列表 |
| 查看: 1611 | 回復(fù): 5 | ||
趙紅霞鐵桿木蟲 (著名寫手)
|
[求助]
求VMD腳本 計(jì)算軌跡文件的溶劑可及表面sasa 已有1人參與
|
|
我的體系包含polymer和solvent,我需要計(jì)算polymer的溶劑可及表面,搜到的一些命令,如: set all [atomselect top "all"] set some [atomselect top "resid 1 to 5"] measure sasa 1.4 $all -restrict $some 是可以計(jì)算一幀的sasa,但是我要計(jì)算軌跡文件的sasa隨模擬時(shí)間的變化,也就是把每一幀的sasa都輸出,求腳本。 |
木蟲 (正式寫手)

|
去 |
木蟲 (正式寫手)

木蟲 (正式寫手)
|
set numframe [molinfo top get numframes] set out [open zhx.txt w] for {set i 0 } {$i<numframe} {incr i} { set some [atomselect top "resid 1 to 5"] set zhx [ measure sasa 1.4 $all -restrict $some] $some frame $i $zhx frame $i $zhx update $some update puts $out $zhx } close $out |

金蟲 (正式寫手)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|