| 查看: 1472 | 回复: 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都输出，求脚本。 |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
VMD 中调用脚本出现 domain error 的问题怎么解决？
已经有2人回复- gromac跑完动力学以后，用VMD计算每一帧构象的溶剂可及表面积，求TCL脚本
已经有0人回复
- 使用AdNDP方法以及ELF/LOL、多中心键级研究多中心键
已经有14人回复
声梦奇缘001
木虫 (正式写手)
- 应助: 56 (初中生)
- 金币: 5597.4
- 红花: 5
- 帖子: 614
- 在线: 218.7小时
- 虫号: 2406187
- 注册: 2013-04-07
- 专业: 化学工程及工业化学
2楼2015-08-26 18:09:16
|
去 |
3楼2015-08-26 19:07:50
声梦奇缘001
木虫 (正式写手)
- 应助: 56 (初中生)
- 金币: 5597.4
- 红花: 5
- 帖子: 614
- 在线: 218.7小时
- 虫号: 2406187
- 注册: 2013-04-07
- 专业: 化学工程及工业化学
4楼2015-08-26 19:52:37
声梦奇缘001
木虫 (正式写手)
- 应助: 56 (初中生)
- 金币: 5597.4
- 红花: 5
- 帖子: 614
- 在线: 218.7小时
- 虫号: 2406187
- 注册: 2013-04-07
- 专业: 化学工程及工业化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
赵红霞: 金币+30, ★★★★★最佳答案, good 2015-08-27 09:13:03
赵红霞: 金币+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
lao7
金虫 (正式写手)
- 应助: 22 (小学生)
- 金币: 1305.3
- 散金: 30
- 红花: 5
- 帖子: 787
- 在线: 176.6小时
- 虫号: 1041612
- 注册: 2010-06-14
- 性别: MM
- 专业: 电磁场与电路
6楼2020-01-14 16:39:12