版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

shixha

银虫 (小有名气)

应助: 4 (幼儿园)
金币: 194.4
散金: 539
红花: 2
帖子: 119
在线: 138.2小时
虫号: 1114260
注册: 2010-10-05
性别: GG
专业: 半导体晶体与薄膜材料

[交流] chgdiff.pl 处理电荷差分密度已有1人参与

本来想用chgdiff.pl处理电荷差分密度, 但处理三个文件时有个小bug的, 浪费半天时间  OMG  写下来供借鉴
1. 因为原子信息保留了file1的信息, 应该倒过来,  file1-file2
2. printf  OUT输出前加个空格, 保持与VASP格式一致

修改:
81c81
<       $line1[$j] = $line1[$j]-$line2[$j];
---
>       $line1[$j] = $line2[$j]-$line1[$j];
83c83
<    printf OUT " %18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5];
---
>    printf OUT "%18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5];

chgdiff.pl CHGCAR1 CHGCAR2
mv CHGCAR_diff  tmp
chgdiff.pl tmp CHGCAR3
得到差分密度:CHGCAR1-CHGCAR2-CHGCAR3
用网上找的fortran处理得到plane-averaged difference electronic density,  注意删除原子符号那行
齐活!

回复此楼

» 本帖已获得的红花（最新10朵）

朱全喜

» 猜你喜欢

292求调剂已经有6人回复
化学工程085602 305分求调剂已经有7人回复
332求调剂已经有5人回复
材料求调剂已经有4人回复
303求调剂已经有5人回复
334分一志愿武理-080500 材料求调剂已经有4人回复
求调剂已经有3人回复
299求调剂已经有3人回复
299求调剂已经有8人回复
316求调剂已经有8人回复

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

1楼 2015-01-09 16:54:28

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

朱全喜

新虫 (小有名气)

应助: 19 (小学生)
金币: 2977.9
红花: 12
帖子: 205
在线: 132.6小时
虫号: 1972785
注册: 2012-09-04
性别: GG
专业: 理论和计算化学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

送红花一朵

非常感谢你的这篇帖子，希望以后还可以交流。

赞一下