24小时热门版块排行榜    

查看: 2510  |  回复: 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的回帖

abbott

金虫 (著名写手)


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

技术支持
2楼2010-09-22 07:48:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)


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


不错
3楼2010-09-22 10:04:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

oyljw

至尊木虫 (职业作家)


你们是学编程的? 这么厉害呀!! 我什么时候能达到你们徒弟的水平啊!
5楼2010-09-23 18:34:20
已阅   回复此楼   关注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的回帖

suosuosky

铁杆木虫 (著名写手)


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

非常感谢分享!!
9楼2011-05-25 20:32:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)


楼主麻烦注释一下,好方便我等菜鸟学习。
11楼2011-06-30 02:26:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

向日葵的香气

新虫 (初入文坛)


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

真好啊,太实用了
13楼2014-10-17 11:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
czyzsu8楼
2010-09-25 13:14   回复  
 谢谢分享!
82454868110楼
2011-06-26 20:52   回复  
五星好评  
senen12楼
2012-04-11 13:09   回复  
五星好评  
相关版块跳转 我要订阅楼主 cxyuan09 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:40 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:23 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 by Z7qfCsdIiD1T
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 04:12 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文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:54 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 01:21 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.551.O54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 00:59 by Z7qfCsdIiD1T
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 22:39 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:12 by Z7qfCsdIiD1T
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考研] 售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一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 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文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
信息提示
请填处理意见