| 6 | 1/1 | 返回列表 |
| 查看: 1767 | 回復(fù): 5 | ||||
小小鑫000鐵蟲 (小有名氣)
|
[求助]
我自己寫了個vasp的shell腳本,高手幫我看看哪錯了?
|
|
各位高手,大家好: 我想寫一個腳本,目的是把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 這三行是我這個機(jī)器的運(yùn)行指令,這個絕對沒錯,但這之外的就是我瞎寫的了。 最后的結(jié)果就是comment文件空白,機(jī)器出現(xiàn)如下提示: fatal: not enough arguments to satisfy format string `%12.6f ' ^ ran out for this one 求各位前輩幫我改正一下,謝謝! |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +20 |
鐵蟲 (小有名氣)
|
前輩,謝謝您的指導(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)下面的全部內(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)驗: +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 |
鐵蟲 (小有名氣)
|
前輩,這樣寫的話還對嗎? 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 哪? |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +20 |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè),初試351分,復(fù)試被刷求調(diào)劑! +3 | 真名有冰 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +3 | Y-cap 2026-03-29 | 3/150 |
|
|
[考研] 299求調(diào)劑 +8 | 15188958825 2026-03-25 | 8/400 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +5 | 舉個栗子oi 2026-03-24 | 6/300 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|