| 查看: 3514 | 回复: 24 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
homo和lumo批量提取的脚本
|
||
最近在弄一批分子，一万个左右 ，想提取其中的homo和lumo能量值，但是发现很难实现，因为才接触脚本半年，还是自学的，所以搞不定了，有哪位大侠给个小脚本？我只会一点shell，提取个Energy还行，homo和lumo不好提取，若可以分享一下不胜感激！ |
回复此楼» 猜你喜欢
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有2人回复
-
分子动力学模拟合作需求
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有158人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
招收博士一名-宁波大学
已经有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
