24小时热门版块排行榜    

查看: 1998  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 求调剂! +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
信息提示
请填处理意见