24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 1471  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

赵红霞

铁杆木虫 (著名写手)

[求助] 求VMD脚本 计算轨迹文件的溶剂可及表面sasa已有1人参与

我的体系包含polymer和solvent,我需要计算polymer的溶剂可及表面,搜到的一些命令,如:
set all [atomselect top "all"]
set some [atomselect top "resid 1 to 5"]
measure sasa 1.4 $all -restrict $some
是可以计算一帧的sasa,但是我要计算轨迹文件的sasa随模拟时间的变化,也就是把每一帧的sasa都输出,求脚本。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵红霞

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by 声梦奇缘001 at 2015-08-26 18:09:16
VMD 的measure 命令可以实现这个功能,具体的可以参考 VMD 手册

3楼2015-08-26 19:07:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

声梦奇缘001

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
VMD 的measure 命令可以实现这个功能,具体的可以参考 VMD 手册
只做不说,明知没地位,坚信有机会
2楼2015-08-26 18:09:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

声梦奇缘001

木虫 (正式写手)

红姐说..................
只做不说,明知没地位,坚信有机会
4楼2015-08-26 19:52:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

声梦奇缘001

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
赵红霞: 金币+30, ★★★★★最佳答案, good 2015-08-27 09:13:03
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
只做不说,明知没地位,坚信有机会
5楼2015-08-26 19:56:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见