24小时热门版块排行榜

返回列表

【奖励】本帖被评价20次，作者tephoon78增加金币 17.3 个

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tephoon78

木虫 (正式写手)

应助: 4 (幼儿园)
金币: 1682
帖子: 677
在线: 130.3小时
虫号: 330598

[资源] 【原创】成功编译NBO 5.G on GAMESS (January 12, 2009 R3.)

NBO 5的源码下载
http://www.mdbbs.org/viewthread.php?tid=12591&highlight=nbo
2009年1月的源码。
另外在厦门大学的量子化学论坛上找到的NBO 5是02年的，在链接时出错。

1.  首先确保自己会、且已编译Gamess，可以参考
http://hi.baidu.com/sobereva/blo ... 8dd04fd6887db6.html
很详细，但是有点乱。

2. 先在windows下解压，copy到linux下
3、vi nbo_5g.src
:26638  ! 到达文件的26639行
dG
:w NBO_5GA.SRC !将nbo_5g.src的剩余内容写到一个新文件内
:q!  !不保存退出。
vi nbo_5g.src
:1 !将光标定位在第一行
:d 26638 ! 从文件的第1行到26638行全删掉
:w NBO_5GB.SRC !将nbo_5g.src的剩余内容写到一个新文件内
:q!  !不保存退出。
注：这一步的工作是将nbo_5g.src这个文件分成2个文件NBO_5GA.SRC和NBO_5GB.SRC
4、$ifort –o enable enable.f    !编译enable.f，我用的是intel fortran 11.1
$ ./enable                   !运行enable下面是过程。
Currently supported environments:
(1) 32-bit Unix UNIX/LINUX
(2) 64-bit Unix UNIX/LINUX
(3) DOS          PC DOS
(4) VMS          VAX VMS

Environment to enable (1-4)? 1
Currently supported NBO versions:
(1) GEN    GENNBO  (standalone)
(2) G03    Gaussian 03
(3) G98    Gaussian 98
(4) G94    Gaussian 94
(5) GMS    GAMESS
(6) NWC    NWChem
(7) AM1    MOPAC or AMPAC

NBO Program version to enable (1-7)? 5

Module  gmsnbo.src  will be created; OK to proceed [Y]/N? Y

NBO source (NBO_5G.SRC ) contains 52194 lines of code.
A total of 44308 lines were written to  gmsnbo.src
including 272 utility and  2 operating system dependent lines.
－－－－－－
结束后产生一个gmsnbo.src文件。

5、vi gmsnbo.src
8427行 IF(MORB.EQ.0.AND.NAMEAT(IATNO(NCTR.NE.'gh') THEN
修改成 IF(MORB.EQ.0.AND.NAMEAT(IATNO(NCTR.NE.'gh'))) THEN

15421行  IF(MORB.EQ.0.AND.NAMEAT(IATNO(NCTR).NE.'gh') THEN
修改成 IF(MORB.EQ.0.AND.NAMEAT(IATNO(NCTR).NE.'gh')) THEN

6、cp gmsnbo.src ~gamess/source/
cd ~gamess/source/
7、vi prppop.src
将1111行修改成
IF(NPROC.EQ.1) CALL RUNNBO
8、vi rhfuhf.src
1676行插入一行：
         COMMON /ENEDA / E0
4723行在ETOT = EHF+EN行插入1行
         IF(ITER.EQ.1) E0 = ETOT

9、vi  iolib.src:
269行    IF(LEN.LE.0          ) GO TO 810
修改成    IF(LEN.LT.0          ) GO TO 810

281行       CALL DARD(V(IS),LENW,IDAF,NSP,IDTYP)
修改成       IF(LENW .GT. 0) CALL DARD(V(IS),LENW,IDAF,NSP,IDTYP)

549行       IF (N .GT. 0 .AND. LEN .NE. IFILEN(NREC)) GO TO 800
修改成       IF (N .GT. 0 .AND. LEN .GT. IFILEN(NREC)) GO TO 800

10、cd ..
vi compall
在 265行，插入一行
         ./comp gmsnbo

11. vi lked
在 801行最后加入 gmsnbo.o

12、./compall >& compall.log&       !重新编译gamess
13、./lked gamess 01 >& lnked.log& ！重新链接程序  01可以是其它数值
14、vi rungms
修改set VERNO=01 或你自己设置的数值。
完成。注意nbo 不能并行。

http://hi.baidu.com/quchem/blog/ ... 8c40d6a3cc2b00.html

[ Last edited by tephoon78 on 2010-1-21 at 14:40 ]

回复此楼