admin 管理员组

文章数量: 1086019

步骤一：准备gnuplot 画图脚本potential.gnu

#potential.gnu内容

move= 
set border lw 2
set arrow   from  0,0  to  40,0     nohead dashty 2 lw 1
set xlabel 'Distance (Å)'

set ylabel 'E (eV)'
set grid 
set xtics nomirror  
set ytics nomirror
#set yrange [6:16]

set xlabel font "Times New Roman,25" offset 0,-0.5
set ylabel font "Times New Roman,25" offset 0.1,0
set xtics font "Times New Roman,18"
set ytics font "Times New Roman,18"
set key font "Times New Roman, 18"
set key outside ###将key 放在外面
set key right top 


plot 'PLANAR_AVERAGE.dat' u 1:($2-move) w l linewidth 2 lc rgb 'red' title ''

步骤二：画图。我将potential.gnu 位置设置为 ~/script/potential.gnu

这里我利用了for循环，画出了0.0~(-0.1)文件下scf的势能图。并将每个图的title设置为了a.$i

#用gnulot 画静电势能(利用了vaspkit)

E=$(seq 0.0 -0.1 -1.3)
for i in $E
do
cd a.$i;  #打开a.$i文件
rm scf/potential.gnu
cp ~/script/potential.gnu scf;        
path="$(pwd)"; 
sed -i "1i set title \"`basename "$path"`\" font \"Times New Roman,25\"" scf/potential.gnu;
cd scf;    #进入a.$i/scf文件中
echo -e "426\n3\n"|vaspkit
ss=$(echo "911" |vaspkit |grep -m 1 "Fermi Energy (eV):"| awk '{print $5}');
sed -i "s/^move=.*/move=$ss/g" potential.gnu;
unset ss;
gnuplot -p  potential.gnu;
cd ../..;   #退出a.$i.scf
done





#拓展
#path="$(dirname "$(pwd)")";ss=$(basename "$path")
#sed -i "s/^\#BSUB -J.*/\#BSUB -J $ss/g" sub;

本文标签： 势能 并将 画出 静电 文件

版权声明：本文标题：在linux 系统下将~.scricppotential 复制到a.$iscf下，并设置title 为a.$i；并将每个文件的静电势能图用gnuplot 画出 内容由网友自发贡献，该文观点仅代表作者本人， 转载请联系作者并注明出处：http://roclinux.cn/b/1738258930a1952299.html， 本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容，一经查实，本站将立刻删除。