| 查看: 813 | 回复: 4 | ||
| 【奖励】 本帖被评价4次，作者rorob增加金币 3.2 个 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[资源]
NAMD Pressure profile calculation 输出来计算溶液表面张力的shell脚本
|
||
|
因为一直在用namd跑溶液体系，最近翻namd ug 13.2看到关于 Pressure profile calculation （侧向压计算），于是决定试试用这个来计算表面张力。 数据量比较多，于是写了这个shell脚本 首先在linux系统中打开控制台 chmod +x interface_tension.awk #首先要给权限 ./interface_tension.awk *.log #输入文件为 namd的输出log文件 具体见 NAMD ug 13.2 这个是计算侧向压的输出log跟通常的格式不一样 Input the thickness of one slab(angstorm) #提示输入slab的厚度 单位为1埃米 10 ”Enter“ #输入10车 最后输出： slabs tension # 默认输出两列， 沿Z轴的不同slab的张力值 类似密度分布 两边为零中间也趋近零 具体结果与模拟参数以及体系有很大关系 N (mN/m) 1 0 2 0 3 0 4 0 5 0 6 0 7 0 8 -0.97385 9 35.9439 10 -12.935 11 0.993465 12 1.97719 13 0.80234 14 1.93347 15 0.790208 16 1.10239 17 -13.2716 18 37.2449 19 13.0549 20 0 21 0 22 0 23 0 24 0 25 0 26 0 27 0 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : interface_tension.awk
- 附件 2 : Tension8PF2EONVTX3.log
2013-12-28 12:42:56, 708 bytes
2013-12-28 12:53:25, 7.22 M
» 猜你喜欢
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有9人回复
-
版面费该交吗
已经有9人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
面上可以超过30页吧？
已经有4人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有5人回复
-
什么是人一生最重要的？
已经有4人回复
3楼2016-12-19 13:05:51
简单回复
dummy12楼
2015-04-02 04:55
回复
五星好评 顶一下，感谢分享！
carldxl4楼
2016-12-22 18:48
回复
五星好评 顶一下，感谢分享！
2017-07-16 08:55
回复
五星好评 顶一下，感谢分享！