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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » homo和lumo批量提取的脚本
253/3返回列表上一页123
查看: 3250  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 重要消息，中午系统在维护 +6 yuleib84 2026-08-18 7/350 2026-08-18 15:02 by sytp
[教师之家] 两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 5/250 2026-08-18 14:55 by LNP@mRNA
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +46 医学老男孩 2026-08-13 102/5100 2026-08-18 14:41 by qinwangshen
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
关闭
请填处理意见
关闭