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,  注意删除原子符号那行
齐活!

回复此楼