| 查看: 3360 | 回复: 24 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
homo和lumo批量提取的脚本
|
||
最近在弄一批分子，一万个左右 ，想提取其中的homo和lumo能量值，但是发现很难实现，因为才接触脚本半年，还是自学的，所以搞不定了，有哪位大侠给个小脚本？我只会一点shell，提取个Energy还行，homo和lumo不好提取，若可以分享一下不胜感激！ |
回复此楼» 猜你喜欢
申请2027年材料/化学类博士
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
物理化学论文润色/翻译怎么收费?
已经有71人回复
-
求相关合成条件
已经有14人回复
-
2027春季申博自荐
已经有6人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有68人回复
-
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）
已经有78人回复
-
申博
已经有2人回复
-
核磁分析求助
已经有4人回复
-
有偿核磁分析求助
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
【答案】应助回帖
★ ★ ★ ★ ★
zhou2009(金币+5): 2011-12-11 20:12:18
秋雨子(金币+40): 没有错误，不用修改可以直接搞定！也不用再学习新的语言，非常感谢！ 2011-12-11 20:26:28
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 |
秋雨子16楼2011-12-10 13:31:18
2楼2011-12-06 14:11:26
   |
3楼2011-12-07 11:23:16
abbott
金虫 (著名写手)
不要用QQ问我东西
- 应助: 16 (小学生)
- 金币: 1071.6
- 散金: 2787
- 红花: 10
- 帖子: 1015
- 在线: 105.2小时
- 虫号: 452267
- 注册: 2007-11-05
- 性别: GG
- 专业: 计算机硬件技术
4楼2011-12-08 10:57:14
