| 4 | 1/1 | 返回列表 |
| 查看: 935 | 回復(fù): 3 | ||
2008130252銀蟲 (正式寫手)
|
[求助]
侯博士腳本問題 已有2人參與
|
|
侯博士vasp手冊(cè)中晶格常數(shù)優(yōu)化的腳本,輸出結(jié)果一直不對(duì)。上圖是程序中的輸出行,其中最外側(cè)的符號(hào)已是tab鍵上的符號(hào),下圖是輸出結(jié)果。求指教 run_a0.PNG EtVo.PNG |
|
#!/bin/sh rm WAVECAR for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30 do cat > POSCAR <<! Al-fcc $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 1 Direct 0.0 0.0 0.0 ! echo " a = $i angstrom "; time vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f \n" , $5}'` echo $V $E >>EtVo.dat done |
至尊木蟲 (職業(yè)作家)
|
最好是自學(xué)一點(diǎn)bash編程的知識(shí),awk使用的知識(shí)。 如果你是從pdf文件拷貝出來那段腳步并粘貼至Linux下某個(gè)文件,常會(huì)出現(xiàn)錯(cuò)誤,出現(xiàn)不了所期望的結(jié)果,多半是因?yàn)檎迟N到文件里的雙引號(hào)、單引號(hào)、反單引號(hào)出現(xiàn)了問題,或漏掉了某個(gè)字符。 你拷貝錯(cuò)了吧。 12.6后面漏掉了f。 %12.6f表示輸出浮點(diǎn)型數(shù)據(jù),整個(gè)長度為12,小數(shù)點(diǎn)之后為6位。 |
銀蟲 (正式寫手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|