| 查看: 1697 | 回复: 5 | ||
[求助]
我自己写了个vasp的shell脚本，高手帮我看看哪错了？
|
|
各位高手，大家好： 我想写一个脚本，目的是把TOTAL FORCE下的力的数据呈现到comment文件中。也就是我想把OUTCAR里的一下部分通过脚本写到comment文件中。如下： POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 0.00000 0.00000 4.05312 0.000010 -0.000006 -0.005155 0.00000 0.00000 10.13242 0.000022 -0.000013 -0.030471 1.24096 0.71646 2.02643 -0.000046 0.000027 -0.050123 1.24096 0.71646 8.10599 -0.000039 0.000023 0.059639 -0.00001 1.43293 0.00000 -0.000014 0.000008 0.029145 -0.00001 1.43293 6.07955 0.000067 -0.000039 -0.003035 我自己写的脚本如下： #!/bin/sh rm WAVECAR -rf source /opt/openmpi/openmpi-intel.sh source /etc/profile mpirun -np 16 -machinefile ma vasp E=`grep -A 7 TOTAL OUTCAR | tail -6 | awk '{printf "%12.6f \n"}'` echo $E >>comment 其中source /opt/openmpi/openmpi-intel.sh source /etc/profile mpirun -np 16 -machinefile ma vasp 这三行是我这个机器的运行指令，这个绝对没错，但这之外的就是我瞎写的了。 最后的结果就是comment文件空白，机器出现如下提示： fatal: not enough arguments to satisfy format string `%12.6f ' ^ ran out for this one 求各位前辈帮我改正一下，谢谢！ |
回复此楼» 本帖@通知
» 猜你喜欢
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有6人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于linux下的任务管理软件比如PBS
已经有10人回复
- 提取VASP计算前后的坐标转换脚本（shell）
已经有9人回复
- python之蟒开启理论计算化学的新时代
已经有22人回复
- 计算费米面时遇到的问题
已经有26人回复
- BASH脚本的一个小问题
已经有3人回复
- Linux下怎样使用批处理让任务一个接一个算？
已经有13人回复
2楼2014-02-21 22:59:53
|
前辈，谢谢您的指导，最后comment的结果如下呈现： 0.000010 -0.000006 -0.005155 0.000022 -0.000013 -0.030471 -0.000046 0.000027 -0.050123 -0.000039 0.000023 0.059639 -0.000014 0.000008 0.029145 0.000067 -0.000039 -0.003035 但我又想让它完全呈现下面的全部内容及写法和格式： POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 0.00000 0.00000 4.05312 0.000010 -0.000006 -0.005155 0.00000 0.00000 10.13242 0.000022 -0.000013 -0.030471 1.24096 0.71646 2.02643 -0.000046 0.000027 -0.050123 1.24096 0.71646 8.10599 -0.000039 0.000023 0.059639 -0.00001 1.43293 0.00000 -0.000014 0.000008 0.029145 -0.00001 1.43293 6.07955 0.000067 -0.000039 -0.003035 不要像0.000010 -0.000006 -0.005155 0.000022 -0.000013 -0.030471 -0.000046 0.000027 -0.050123 -0.000039 0.000023 0.059639 -0.000014 0.000008 0.029145 0.000067 -0.000039 -0.003035 这样杂乱，并且有POSITION TOTAL-FORCE (eV/Angst) 这一行也要加上啊!前辈该怎吗改？  |
3楼2014-02-22 11:49:34
★ ★
identation: 金币+2, 谢谢交流 2014-02-25 15:53:51
identation: 金币+2, 谢谢交流 2014-02-25 15:53:51
|
#!/bin/sh rm WAVECAR -rf source /opt/openmpi/openmpi-intel.sh source /etc/profile mpirun -np 16 -machinefile ma vasp grep -A 1 TOTAL OUTCAR > comment grep -A 7 TOTAL OUTCAR | tail -6 | \ awk '{printf "%12.6f %12.6f %12.6f %12.6f %12.6f %12.6f \n",$1,$2,$3,$4,$5,$6}' \ >> comment |
4楼2014-02-23 18:42:32
|
前辈，这样写的话还对吗？ grep TOTAL OUTCAR |tail -1 >>comment grep -A 7 TOTAL OUTCAR | tail -6 | \ awk '{printf "%12.6f %12.6f %12.6f %12.6f %12.6f %12.6f \n",$1,$2,$3,$4,$5,$6}' \ >> comment 或grep -A 1 TOTAL OUTCAR |tail -1 >>comment grep -A 7 TOTAL OUTCAR | tail -6 | \ awk '{printf "%12.6f %12.6f %12.6f %12.6f %12.6f %12.6f \n",$1,$2,$3,$4,$5,$6}' \ >> comment 哪？ |
5楼2014-02-25 12:43:20
6楼2014-02-25 14:14:26