24小时热门版块排行榜    

查看: 2515  |  回复: 12
【奖励】 本帖被评价8次,作者cxyuan09增加金币 6.7 个
本帖产生 1 个 QC强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cxyuan09

金虫 (正式写手)


[资源] 【cxyuan09个人文集】一个在linux下批量产生cube文件的小程序

由于linux生成cube(可以用GV读入.cube文件直接画分子轨道)的速度比在windows下快得多,如果要画大量的分子轨道,可以用本程序在linux批量产生。

程序内容如下:

#!/bin/csh
foreach j ( $argv )
@ i=27
while ( $i <= 39 )
cubegen 0 mo=$i $j.fchk ${j}_$i.cube 0 h
@ i++
end
end

使用方法:
1、在linux下创建程序文件: 比如命名为mos.csh
2、把上述程序内容粘贴进去,保存退出
3、为程序增加可执行属性:  chmod +x mos.csh
4、程序执行方式:
./mos.csh filename

注意:
1、filename不用加扩展名.fchk,
2、需要事先将chk文件转化为fchk文件,批量转化chk文件可参考我之前一个帖子中我回复的内容
http://muchong.com/bbs/viewthread.php?tid=2082151&fpage=1&view=&highlight=&page=2


本例中,将产生某分子从第21低能量的轨道到第39低能量的轨道,大家可以根据需要修改。

如果有很多同分异构体,需要看它们的价轨道,可以多输入几个文件名,使用方式:
./mos.csh filename1 filename2 filename3 ......

注意:filename不用加扩展名“.fchk”

如果觉得有用,别忘了五星支持啊!

[ Last edited by yjcmwgk on 2010-9-25 at 15:09 ]
回复此楼

» 收录本帖的淘帖专辑推荐

量化软件学习

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

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

oyljw

至尊木虫 (职业作家)


你们是学编程的? 这么厉害呀!! 我什么时候能达到你们徒弟的水平啊!
5楼2010-09-23 18:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

bitgreen

金虫 (小有名气)


★★★★★ 五星级,优秀推荐

★ ★
erylingjet(金币+2):感谢交流! 2010-09-23 18:28:47
建议把
cubegen 0 mo=$i $j.fchk ${j}_$i.cube 0 h
中的$j.fchk 改成$j
这样既可执行
./mos.csh *.fchk来批量处理目录下的所有fchk
4楼2010-09-23 17:48:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)


引用回帖:
Originally posted by bitgreen at 2010-09-23 17:48:32:
建议把
cubegen 0 mo=$i $j.fchk ${j}_$i.cube 0 h
中的$j.fchk 改成$j
这样既可执行
./mos.csh *.fchk来批量处理目录下的所有fchk

我看了看程序
如果把$j.fchk 改成$j的话,后面生成的cube的文件名就需要额外的语句来处理了,否则不利于资料整理
6楼2010-09-24 14:15:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bitgreen

金虫 (小有名气)


引用回帖:
Originally posted by hairan at 2010-09-24 14:15:42:


我看了看程序
如果把$j.fchk 改成$j的话,后面生成的cube的文件名就需要额外的语句来处理了,否则不利于资料整理

简单,把$j用cut切前半部分输出
7楼2010-09-24 14:40:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:29 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:56 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:54 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:13 by Z7qfCsdIiD1T
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 20:07 by Z7qfCsdIiD1T
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
信息提示
请填处理意见