24小时热门版块排行榜    

查看: 2408  |  回复: 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文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +8 zju2000 2026-09-26 9/450 2026-09-29 13:47 by ahsoarli
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见