24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1659  |  回复: 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的回帖

声梦奇缘001

木虫 (正式写手)

【答案】应助回帖

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

声梦奇缘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的回帖

lao7

金虫 (正式写手)

引用回帖:
5楼: Originally posted by 声梦奇缘001 at 2015-08-26 19:56:57
set numframe
set out
for {set i 0 } {$i&lt;numframe} {incr i} {
  set some  
  set zhx
  $some frame $i
  $zhx frame $i
  $zhx update
  $some update
  puts $out $zhx
}
close $out...

运行起来总是报错,肯定是某些语句错误 但找不到原因
6楼2020-01-14 16:39:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Andsharia 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[考研] 化学070300 求调剂 +21 哈哈哈^_^ 2026-04-12 21/1050 2026-04-13 21:36 by nxybio2007
[考研] 食品与营养(0955)271求调剂 +9 升格阿达 2026-04-12 9/450 2026-04-13 14:12 by maddjdld
[考研] 一志愿厦大生物学332求调剂 +11 池池池池池池 2026-04-08 11/550 2026-04-13 14:10 by 科研论
[考研] 2本,初试303,0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 电子信息270求调剂 +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv风雷
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 287求调剂 +15 Fnhc 2026-04-07 21/1050 2026-04-10 19:09 by chemisry
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 复试调剂,一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
信息提示
请填处理意见