24小时热门版块排行榜    

查看: 714  |  回复: 4
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长

[交流] 【求助】gromacs中程序grompp警告太多的错误?

运行grompp之后报错如下:

NOTE 2 [file HMG.top, line 38]:
System has non-zero total charge: -1.999994e+00



processing coordinates...
double-checking input for internal consistency...
Velocities were taken from a Maxwell distribution at 310 K
Reading position restraint coords from mini_water.gro
renumbering atomtypes...
converting bonded parameters...
initialising group options...
processing index file...
Analysing residue names:
Opening library file /usr/local/gromacs-4.0.2/share/gromacs/top/aminoacids.dat
There are: 162095      OTHER residues
There are:   842    PROTEIN residues
There are:     0        DNA residues
Analysing Protein...
Analysing Other...
Making dummy/rest group for Acceleration containing 494183 elements
Making dummy/rest group for Freeze containing 494183 elements
Making dummy/rest group for Energy Mon. containing 494183 elements
Making dummy/rest group for VCM containing 494183 elements
Number of degrees of freedom in T-Coupling group Protein is 15703.95
Number of degrees of freedom in T-Coupling group Non-Protein is 972567.06
Making dummy/rest group for User1 containing 494183 elements
Making dummy/rest group for User2 containing 494183 elements
Making dummy/rest group for XTC containing 494183 elements
Making dummy/rest group for Or. Res. Fit containing 494183 elements
Making dummy/rest group for QMMM containing 494183 elements
T-Coupling       has 2 element(s): Protein Non-Protein
Energy Mon.      has 1 element(s): rest
Acceleration     has 1 element(s): rest
Freeze           has 1 element(s): rest
User1            has 1 element(s): rest
User2            has 1 element(s): rest
VCM              has 1 element(s): rest
XTC              has 1 element(s): rest
Or. Res. Fit     has 1 element(s): rest
QMMM             has 1 element(s): rest
Checking consistency between energy and charge groups...
Calculating fourier grid dimensions for X Y Z
Using a fourier grid of 60x60x60, spacing 0.284 0.284 0.284
Estimate for the relative computational load of the PME mesh part: 0.20
This run will generate roughly 476 Mb of data
writing run input file...

There were 2 notes

There was 1 warning

-------------------------------------------------------
Program grompp, VERSION 4.0.2
Source code file: gmx_fatal.c, line: 481

Fatal error:
Too many warnings (1), grompp terminated.
If you are sure all warnings are harmless, use the -maxwarn option.
-------------------------------------------------------


请问这是什么问题

[ Last edited by lei0736 on 2009-11-13 at 21:36 ]
回复此楼

» 猜你喜欢

个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuefei06

木虫 (正式写手)

★ ★ ★
gwdavid(金币+3,VIP+0):非常感谢!
你写出的警告是体系的电荷总数不为零,即整个体系不是电中性的,你需要重新检查一下你的top文件哩的原子的电荷有没有写错,别的警告没有给全,你可以看下aftermd0文件里面有详细的警告内容
4楼2008-12-05 23:49:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

suny2005

木虫 (小有名气)

★ ★ ★
gwdavid(金币+3,VIP+0):非常感谢!
就是违反常规的东西太多!需要你注意!如果你确定这些警告对你的结果没有影响就在grommp命令中加上 -maxwarn 选项就可以了
3楼2008-12-05 21:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiaobf


gwdavid(金币+1,VIP+0):谢谢
grompp里面的每一个warning都是要注意的,除非你真的认为没有问题
5楼2008-12-15 11:32:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见