24小时热门版块排行榜    

查看: 5557  |  回复: 43
【奖励】 本帖被评价37次,作者neweroica增加金币 28.3 个
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

neweroica

木虫 (著名写手)


[资源] 【neweroica个人文集】操作VASP的CHGCAR文件的程序:chgcarman(v0.0)

在前面的帖子中,我发布将CHGCAR转换为cube文件的程序:chgcar2cube。见:http://muchong.com/bbs/viewthread.php?tid=1809216

有时候,我们需要对多个CHGCAR文件进行加减乘除等操作(比如计算电荷密度的变化)。Gaussian提供了cubeman程序用于操作cube文件,但我发现有时会出现bug。所以,我干脆编写了这个chgcarman程序,用于直接对CHGCAR文件进行进行操作。

目前chgcarman的主要功能有:

1. 输出CHGCAR文件的属性,包括:最大值,最小值,平均值,RMS值,正负值的RMS值。这样在绘制等值面图时,有助于选取合适的等值面的数值。
2. 将两个CHGCAR文件相加。
3. 将两个CHGCAR文件相减。
4. 将CHGCAR文件乘以一个缩放因子。
5. 将CHGCAR文件平方。
6. 将CHGCAR文件取常用对数(以10为底)。


建议在Linux下进行编译和使用,则可用标准C编译器进行编译,如:
gcc -o chgcarman -lm chgcarman.c

如果你非要用Windows的话,请先设置源文件中的MS_Window值为1,然后再用C编译器进行编译。这很重要!!!否则程序虽然编译通过,但转换结果会发生严重错误。

注意:目前该程序只在Linux下测试通过,Windows下我还没有测试过,也欢迎大家测试后向我反馈结果。


程序使用帮助:
        chgcarman [option] [ file ...]
        Options:
            -a
                   Add two CHGCAR files to produce a new one.
            -h
                   Show information and usage of this program and exit.
            -l
                   Compute the logarithm of a CHGCAR file to produce a new one.
            -p
                   Get properties of a CHGCAR file  without producing a new one.
            -q
                   Compute the square of a CHGCAR file to produce a new one.
            -s
                   Substract two CHGCAR files to produce a new one.
            -x factor
                   Scale a CHGCAR file by a scaling factor to produce a new one.
        Input files:
            VASP's CHGCAR or LOCPOT files

[ Last edited by yjcmwgk on 2010-6-16 at 21:41 ]
回复此楼

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

材料计算模拟实用技巧 第一性原理相关文档 vasp计算资料 Photochemistry
量化经验 淘淘

» 猜你喜欢

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

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

bingmou

金虫 (著名写手)


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

谢了,用linux的shell在这一点还是不行,awk只能进行单个文件的处理(或者我不懂),所以这个软件太棒了!一定要学习一下,以后自己可以写。希望看不懂问lz,lz可以回复。谢谢
5楼2010-01-26 15:41:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 44 个回答

hookah

木虫 (正式写手)


★★★ 三星级,支持鼓励

鼓励这样的小程序,试算一下哦!
3楼2010-01-24 19:55:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happyjwx

木虫 (小有名气)


★★★ 三星级,支持鼓励

下载试试看,先支持一下
7楼2010-02-28 09:41:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

greatqi

木虫 (正式写手)


★★★ 三星级,支持鼓励

谢谢,下来看看。
9楼2010-05-21 15:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
xh51212楼
2011-08-29 09:06   回复  
五星好评  顶一下,感谢分享!
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-10 01:40 by FChKX8meu1ds
[考博] 售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
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售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一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 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.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.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 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
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 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.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 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
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见