| 5 | 1/1 | 返回列表 |
| 查看: 1770 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
小小鑫000鐵蟲 (小有名氣)
|
[求助]
我自己寫了個(gè)vasp的shell腳本,高手幫我看看哪錯(cuò)了?
|
|||
|
各位高手,大家好: 我想寫一個(gè)腳本,目的是把TOTAL FORCE下的力的數(shù)據(jù)呈現(xiàn)到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 這三行是我這個(gè)機(jī)器的運(yùn)行指令,這個(gè)絕對(duì)沒錯(cuò),但這之外的就是我瞎寫的了。 最后的結(jié)果就是comment文件空白,機(jī)器出現(xiàn)如下提示: fatal: not enough arguments to satisfy format string `%12.6f ' ^ ran out for this one 求各位前輩幫我改正一下,謝謝! |

鐵蟲 (小有名氣)
|
前輩,謝謝您的指導(dǎo),最后comment的結(jié)果如下呈現(xiàn): 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 但我又想讓它完全呈現(xiàn)下面的全部?jī)?nèi)容及寫法和格式: 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) 這一行也要加上啊!前輩該怎嗎改? ![]() |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
#!/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 |
鐵蟲 (小有名氣)
|
前輩,這樣寫的話還對(duì)嗎? 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 哪? |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +13 | 阿斯蒂芬2004 2026-03-28 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 材料專碩調(diào)劑 +3 | 椰椰。 2026-03-29 | 3/150 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +7 | cccchenso 2026-03-29 | 7/350 |
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 289求調(diào)劑 +13 | 新時(shí)代材料 2026-03-27 | 13/650 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|