24小时热门版块排行榜

返回列表

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

秋雨子

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1578.8
散金: 13
帖子: 163
在线: 78.2小时
虫号: 1027042
注册: 2010-05-24
专业: 理论环境化学

[求助] homo和lumo批量提取的脚本

最近在弄一批分子，一万个左右，想提取其中的homo和lumo能量值，但是发现很难实现，因为才接触脚本半年，还是自学的，所以搞不定了，有哪位大侠给个小脚本？我只会一点shell，提取个Energy还行，homo和lumo不好提取，若可以分享一下不胜感激！

回复此楼

» 猜你喜欢

半透性容器进行高湿或者低湿研究已经有0人回复
电化学基础知识与资源网站已经有3人回复
物理化学论文润色/翻译怎么收费? 已经有249人回复
中国科大电池方向任晓迪课题组招收2026级博士生-电解液/电池安全性/人工智能方向已经有20人回复
26年秋季博士申请已经有0人回复
推荐给英语教学者的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有59人回复
推荐给教师的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有30人回复
核磁分析软件MestReNova打开文件时报错已经有0人回复

» 本主题相关商家推荐: (我也要在这里推广)

To the time to life, rather than to life in time

1楼 2011-12-05 22:47:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 1882
红花: 7
帖子: 736
在线: 228.7小时
虫号: 655151
注册: 2008-11-15
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★
zhou2009(金币+5): 2011-12-11 20:12:18
秋雨子(金币+40): 没有错误，不用修改可以直接搞定！也不用再学习新的语言，非常感谢！ 2011-12-11 20:26:28

直接shell:
#!/bin/sh
# ZDB, HNU, 2011/11/15.
# ZDB, HNU, 2011/11/23 modified to be more pratical.
`ls -l *.log | awk '{print $9}' | sed 's/.log//g' > file_list`
cat file_list | while read each_file
do
file1=./$each_file".log"
HOMO=`grep 'Alpha  occ.' $file1 | tail -n 1 | awk '{print $NF}'`
`grep '^.Alpha virt.' $file1 > tmp`
total_line=`awk '{print $1}' tmp | wc -l`
blockflag=`grep 'Population analysis' $file1 | wc -l`
if [ $blockflag == "1" ]
then
      start_n=`expr ${total_line} \* 1 / 2 + 1`
elif [ $blockflag == "2" ]
then
      start_n=`expr ${total_line} \* 1 / 2 + 1`
else
      start_n=`expr ${total_line} \* 2 / 3 + 1`
fi
LUMO=`sed -n ''$start_n'p' tmp | awk '{print $5}'`
echo $each_file $HOMO $LUMO >> extract.txt
done