| 查看: 1128 | 回复: 1 | |||
congmingbing木虫 (正式写手)
|
[交流]
【求助】lammps计算heat/flux出错 已有1人参与
|
|
在lammps手册上看到了compute heat/flux command，用提供的输入文件执行，出错了，错误信息是 ERROR: Illegal compute heat/flux command，不知道这个错误怎么解决，感觉像是版本的问题。我用的lammps是在主页上下载的lammps-15Jan10.tar.gz，这个应该是最新版了吧。 刚开始学lammps,还请高手指点！万分感谢。 下面是in文件。 ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ atom_style atomic units real dimension 3 boundary p p p lattice fcc 5.376 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region box block 0 4 0 4 0 4 create_box 1 box create_atoms 1 box mass 1 39.948 pair_style lj/cut 13.0 pair_coeff * * 0.2381 3.405 group every region box velocity all create 70 102486 mom yes rot yes dist gaussian timestep 4.0 thermo 10 # ------------- Equilibration and thermalisation ---------------- fix NPT all npt 70 70 10 xyz 0.0 0.0 100.0 drag 0.2 run 8000 unfix NPT # --------------- Equilibration in nve ----------------- fix NVE all nve run 8000 # -------------- Flux calculation in nve --------------- reset_timestep 0 compute myKE all ke/atom compute myPE all pe/atom compute myStress all stress/atom virial compute flux all heat/flux myKE myPE myStress log flux.log variable J equal c_flux[1]/vol thermo_style custom step temp v_J run 100000 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
面上意见出来了
已经有17人回复
-
找导师
已经有5人回复
-
国社科又开始会评了，不知道这次命运如何
已经有3人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
congmingbing
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2089.9
- 散金: 106
- 帖子: 333
- 在线: 412.3小时
- 虫号: 446335
- 注册: 2007-10-29
- 专业: 凝聚态物性I:结构、力学和
2楼2010-06-17 19:32:50
