E=`seq -0.1 0.01 0.1` # 生成从 -0.1 到 0.1 步长为 0.01 的序列
for i in $E
do
mkdir $i
cp INCAR KPOINTS POTCAR POSCAR $i/ # 明确指定需要复制的文件
awk '{if(NR==2) print $1 * '$i'; else print $0}' POSCAR > $i/POSCAR.tmp
awk '{if(NR==5) print $1,$2,$3/(1+'$i'); else print $0}' $i/POSCAR.tmp > $i/POSCAR
rm $i/POSCAR.tmp
done
awk 双轴应变
于 2024-03-28 10:58:52 首次发布