for i in {1..8};do cp 400 $((400+$i*50)) -r ;sed -i "10s/400/$((400+$i*50))/g" $((400+$i*50))/INCAR;done
上述为创造400到800的截断能的命令,前提需要创造一个400的截断能的文件
下面为配套的提交命令
for i in {1..9};do cd $((350+$i*50));sh run_vasp.sh;cd $OLDPWD;done
for i in *; do cd $i ; qsub sub4; cd $OLDPWD; done
输出截断能,能量,时间的提交命令
for i in *0; do echo -e $i "\t" $(grep User $i/OUTCAR |awk '{print $4}')"\t" $(grep ' without' $i/OUTCAR | tail -n 1 | awk '{print $4}'); done
画图的简要脚本
#!/usr/bin/python3
import matplotlib.pyplot as plt
import numpy as np
x,y,z=np.loadtxt('date1.dat',delimiter=',',usecols=(0,1,2),unpack=True)
plt.xlabel('ENCUT/eV')
plt.ylabel('TIME/s')
plt.plot(x,y,'rs-',linewidth=2.0)
plt.plot(x,z,'bs-',linewidth=2.5)
plt.show()
处理POSCAR文件晶胞大小快捷命令
for i in $(seq 10 2 20);do cp 8 $i -r;sed -i "3,5s/8.0/$i.0/g" $i/POSCAR;done