| 查看: 2027 | 回复: 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 求各位前辈帮我改正一下，谢谢！ |
回复此楼» 本帖@通知
» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于linux下的任务管理软件比如PBS
已经有10人回复
- 提取VASP计算前后的坐标转换脚本（shell）
已经有9人回复
- python之蟒开启理论计算化学的新时代
已经有22人回复
- 计算费米面时遇到的问题
已经有26人回复
- BASH脚本的一个小问题
已经有3人回复
- Linux下怎样使用批处理让任务一个接一个算？
已经有13人回复
★ ★
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
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
|
前辈，这样写的话还对吗？ 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
200