24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4107  |  回复: 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的回帖

wangzg9292

木虫 (正式写手)


★★★ 三星级,支持鼓励

good!顶了
2楼2010-01-24 12:13:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)


★★★ 三星级,支持鼓励

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

bingmou

金虫 (著名写手)


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

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

dsagfds

银虫 (小有名气)


★★★ 三星级,支持鼓励

顶 值得推荐
6楼2010-02-06 21:35:50
已阅   回复此楼   关注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的回帖

valenhou001

至尊木虫 (职业作家)


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

如果做自旋极化的计算,还需考虑spin density= charge density (spin-up) - charge density (spin-down)。

如果楼主有兴趣,可以将你的处理vasp输出结果的c代码整合成一个类似lev00的工具。
10楼2010-05-22 09:01:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bigcontinent

金虫 (正式写手)


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

楼主勇猛
11楼2010-07-05 22:13:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youzhizhe

荣誉版主 (职业作家)


帖子真精彩!
已经收录到淘贴专辑《材料计算模拟实用技巧
13楼2011-08-30 14:38:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stefsui

木虫 (小有名气)


引用回帖:
5楼: Originally posted by bingmou at 2010-01-26 15:41:50:
谢了,用linux的shell在这一点还是不行,awk只能进行单个文件的处理(或者我不懂),所以这个软件太棒了!一定要学习一下,以后自己可以写。希望看不懂问lz,lz可以回复。谢谢

awk可以处理,单纯相减的话
awk很强大.....
14楼2011-11-26 00:52:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)


awk可以处理多个文件的
18楼2012-02-06 03:27:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuzhaonan

铁杆木虫 (著名写手)


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

楼主您好,编译过程中出现,“undefined reference to `sqrt'”应该如何处理呢,谢谢您!
20楼2012-05-20 20:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuzhaonan

铁杆木虫 (著名写手)


已解决谢谢楼主
21楼2012-05-20 20:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinye87

禁虫 (小有名气)

本帖内容被屏蔽

26楼2012-12-29 21:56:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppshuang

新虫 (初入文坛)


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

太赞了!!!
32楼2014-04-25 18:09:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
fox-fox4楼
2010-01-26 12:30   回复  
 
qzhaosdu8楼
2010-05-06 10:27   回复  
 
xh51212楼
2011-08-29 09:06   回复  
五星好评  顶一下,感谢分享!
thinkllh15楼
2011-11-26 10:50   回复  
五星好评  顶一下,感谢分享!
huyuxiav16楼
2011-12-30 10:17   回复  
五星好评  顶一下,感谢分享!
senlinwolf17楼
2012-02-05 08:40   回复  
三星好评  顶一下,感谢分享!
seaso19楼
2012-02-08 11:35   回复  
五星好评  顶一下,感谢分享!
sjzxbe22楼
2012-06-05 14:51   回复  
五星好评  顶一下,感谢分享!
dufuming23楼
2012-08-28 15:35   回复  
五星好评  顶一下,感谢分享!
sunnsi24楼
2012-09-13 09:10   回复  
五星好评  顶一下,感谢分享!
qinye8725楼
2012-12-29 21:55   回复  
五星好评  顶一下,感谢分享!
2013-03-04 16:38   回复  
五星好评  顶一下,感谢分享!
cenchao61328楼
2013-09-13 21:33   回复  
一般  顶一下,感谢分享!
wz540902529楼
2013-11-01 19:12   回复  
五星好评  顶一下,感谢分享!
baiyujie30楼
2013-11-08 16:40   回复  
五星好评  顶一下,感谢分享!
xingeruc31楼
2014-01-10 19:00   回复  
五星好评  顶一下,感谢分享!
wucx123433楼
2014-04-30 15:19   回复  
五星好评  顶一下,感谢分享!
2014-06-13 17:16   回复  
五星好评  顶一下,感谢分享!
2008089035楼
2014-07-02 11:00   回复  
五星好评  顶一下,感谢分享!
50418372336楼
2015-03-03 15:09   回复  
三星好评  
wanby7637楼
2015-03-04 15:56   回复  
五星好评  顶一下,感谢分享!
hczheng89938楼
2015-08-11 17:17   回复  
五星好评  顶一下,感谢分享!
hczheng89939楼
2015-12-26 21:36   回复  
顶一下,感谢分享!
sportbotao40楼
2016-01-11 12:54   回复  
五星好评  顶一下,感谢分享!
叹奇41楼
2016-10-27 18:53   回复  
五星好评  顶一下,感谢分享!
aa560749842楼
2018-03-13 23:07   回复  
五星好评  顶一下,感谢分享!
2018-05-30 09:41   回复  
五星好评  顶一下,感谢分享!
yfwanhy44楼
2021-03-14 19:00   回复  
五星好评  牛
相关版块跳转 我要订阅楼主 neweroica 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见