应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » homo和lumo批量提取的脚本
253/3返回列表上一页123
查看: 2890  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

秋雨子

金虫 (小有名气)
引用回帖:
18楼: Originally posted by asaka at 2011-12-10 16:33:17:
以上作者只考虑闭壳层计算的情况，你有开壳层的计算吗？那时候你需要总的HOMO,LUMO还是分别alpha以及beta的HOMO和LUMO?
他们的程序都只有一点小问题，相信LZ可以搞定。
mchen10 的程序少了一个 |，多了一个奇怪 ...

非常感谢高手的赐教，我的分子群全是闭壳层的，当然也有许多别的分子信息需要提取，那些都相对简单一些，自己可以学着搞定！非常感谢！
赞一下 回复此楼
To the time to life, rather than to life in time
21楼2011-12-11 20:10:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋雨子

金虫 (小有名气)
送鲜花一朵
引用回帖:
19楼: Originally posted by mchen10 at 2011-12-11 04:06:32:
[code]#!/bin/bash
# print homo lumo energy of gaussian output
# usage: print_homo_lumo filename(s)
for file in "$@"
do
echo $file
tac $file | tr "\n" "@" | egre ...

谢谢不吝赐教！我会努力的学习的，才接触脚本，总是觉得深奥难懂，心理的排斥应该克服！总是想绕过难题去走捷径是不行的！非常感谢！
赞一下 回复此楼
To the time to life, rather than to life in time
22楼2011-12-11 20:13:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)
`ls -l *.log | awk '{print $9}' | sed 's/.log//g' > file_list`
这一行中的9可能是8，需要自己数一下ls -l *.log出现多少列。
赞一下 回复此楼
23楼2011-12-12 13:28:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)
再来个Fortran
PROGRAM HOLU2
!
! Program to extract HOMO and LUMO data from multiple Gaussian output files
! SRK 2011/06/30
! SRK 2011/10/18 Modified to produce 3 output files, bug fix
!
! Input : A file containing a list of filenames, one per line
! Output:  One line per input filename, 3 columns, holding the filename, the HOMO and the LUMO
!          Input filename, HOMO and LUMO are written to file 'HOMOandLUMO.txt'
!          HOMO is written to file 'HOMO.txt'
!          LUMO is written to file 'LUMO.txt'
!
IMPLICIT NONE
! Variable declarations
INTEGER, PARAMETER :: inunit=22, listunit=23 ! Unit numbers for input files and list file
INTEGER, PARAMETER :: hlunit=24, hunit=25, lunit=26 ! Unit numbers for output files
INTEGER :: i          ! Loop variable
INTEGER :: n          ! Number of characters in current line
INTEGER :: stat       ! Holds the current state of the input file, 0=OK
INTEGER :: numfiles   ! Number of filenames detected in the input
INTEGER :: blockflag  ! Flag variable to find the correct HOMO-LUMO data
CHARACTER (LEN=100) :: listname ! File containing the list of filenames to be processed
CHARACTER (LEN=100) :: inpfile  ! Current input filename
CHARACTER (LEN=100) :: line     ! Current line in the input file
CHARACTER (LEN=100) :: prevline ! Previous line in the input file
CHARACTER (LEN=9) :: homostring, lumostring ! Strings containing the HOMO/LUMO eigenvalues

!
! Executable code

WRITE (*,*) 'HOLU - a program to extract HOMO and LUMO data'
WRITE (*,*) 'Enter the name of an input file containing a list of filenames (no blank lines!)'
READ (*,'(A)') listname
OPEN(listunit,FILE=listname)
!
! Find out how many files are in the list
numfiles = 0
stat = 0
DO
   READ(listunit,'(A)',iostat=stat) line
   IF (stat /= 0) EXIT
   numfiles = numfiles + 1
END DO
CLOSE(listunit)
WRITE (*,*) 'Detected ',numfiles,' lines in list file'

!
! Open output files
OPEN(hlunit,FILE='HOMOandLUMO.txt')
OPEN(hunit,FILE='HOMO.txt')
OPEN(lunit,FILE='LUMO.txt')

! Now loop over these files
OPEN(listunit,FILE=listname)
DO i = 1, numfiles
   READ(listunit,'(A)') inpfile
   OPEN(inunit,file=TRIM(inpfile))
   blockflag = 0
!
! Search for the needed information in the input file
   prevline = ''
   DO
      READ(inunit,FMT='(A)',IOSTAT=stat) line ! Read in the next line of the file using '(A)' format
      IF (stat /= 0 ) EXIT
      IF (index(line,'Alpha virt.') > 0 ) THEN   ! found a line listing virtual states
        IF (blockflag == 0) blockflag = 1
        IF (blockflag == 2) THEN                    ! we have found the second block
           lumostring = line(30:40)                 ! first number on this line
           n = len_trim(prevline)
           homostring = prevline(n-8:n)             ! last number on previous line
           WRITE(*,*) TRIM(inpfile),' ', 'HOMO', ' ', homostring,' ', 'LUMO', ' ', lumostring
           WRITE(hlunit,'(A),(A),(A),(A)') TRIM(inpfile),' ', 'HOMO', ' ', homostring,' ', 'LUMO', ' ', lumostring
           WRITE(hunit,'(A)') TRIM(inpfile),' ', 'HOMO', ' ', homostring
           WRITE(lunit,'(A)') TRIM(inpfile),' ', 'LUMO', ' ', lumostring
           EXIT
        END IF
      ELSE
        IF (blockflag == 1) blockflag = 2      
      END IF
      prevline = line
   END DO
   CLOSE(inunit)
END DO


! Close all files that are still open
CLOSE(listunit)
CLOSE(hlunit)
CLOSE(hunit)
CLOSE(lunit)

!
END PROGRAM HOLU2
这是source code,需要编译。
编译后，`ls -l *.log | awk '{print $9}' > file_list
再运行编译后的程序，输入文件就是file_list。所有的输出结果都显示在屏幕上，可以拷入txt文档，进行后续的计算。
赞一下 回复此楼
24楼2011-12-12 13:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋雨子

金虫 (小有名气)
高手如云啊！学习了学习了！可以成为经典贴了！
赞一下 回复此楼
To the time to life, rather than to life in time
25楼2011-12-13 23:43:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 秋雨子 的主题更新
253/3返回列表上一页123
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 5/250 2026-03-17 17:27 by Little-xue
[考研] 344求调剂 +4 knight344 2026-03-16 4/200 2026-03-17 17:27 by ruiyingmiao
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 268求调剂 +6 好运连绵不绝 2026-03-12 7/350 2026-03-17 14:56 by 呦呦忧郁
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 一志愿中科院，化学方向，295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 328化工专硕求调剂 +4 。，。，。，。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭