| 查看: 4604 | 回复: 98 | |||
| 本帖产生 1 个 模拟EPI ,点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
【求助】DL_POLY的config文件如何生成?? 已有9人参与
|
|||
|
我想模拟聚合物在金属团簇表面的吸附,想请教大家这样的体系的config文件如何生成? 另外我把手册已经看了一遍,但是希望运行一下.有么有人根据手册的算例给写一个模拟过程,给传一个谢谢!! |
» 猜你喜欢
退学或坚持读
已经有25人回复
遇见不省心的家人很难过
已经有13人回复
博士延得我,科研能力直往上蹿
已经有4人回复
免疫学博士有名额,速联系
已经有14人回复
面上基金申报没有其他的参与者成吗
已经有4人回复
多组分精馏求助
已经有6人回复

★ ★ ★
lei0736(金币+3):谢谢 2010-03-23 08:54
lei0736(金币+3):谢谢 2010-03-23 08:54
|
program dlpolyinput implicit none integer i, numcoordinate, pbc, numatom real, PARAMETER :: K=0 character *20 atomname, systemname real x(200),y(200),z(200), lbox open(20,file='CONFIG',status='unknown') !打开一个输出文件 open(1,file='input.txt',status='unknown')!打开文件,并读取 read(1,*) systemname write(20,"(A80)" systemnameread(1,*) numcoordinate, pbc, numatom write(20,"(I10,I10,I10)" numcoordinate, pbc, numatomread(1,*) lbox write(20,"(F20.3,F20.3,F20.3)" lbox, K, Kwrite(20,"(F20.3,F20.3,F20.3)" K, lbox, Kwrite(20,"(F20.3,F20.3,F20.3)" K, K, lbox do i=1, 200 read(1,*)atomname,x(i),y(i),z(i) write(20,"(A8,I10)" atomname, iwrite(20,"(F20.3,F20.3,F20.3)" x(i), y(i) ,z(i)end do end |

49楼2010-03-22 22:18:53
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2010-01-24 16:53:03

3楼2010-01-24 20:21:22
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2010-01-25 09:45:24













回复此楼
systemname