 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1769 | 回復: 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 這三行是我這個機器的運行指令，這個絕對沒錯，但這之外的就是我瞎寫的了。 最后的結(jié)果就是comment文件空白，機器出現(xiàn)如下提示： fatal: not enough arguments to satisfy format string `%12.6f ' ^ ran out for this one 求各位前輩幫我改正一下，謝謝！ |
回復此樓» 本帖@通知
» 猜你喜歡
291求調(diào)劑
已經(jīng)有7人回復
-
材料與化工272求調(diào)劑
已經(jīng)有13人回復
-
283求調(diào)劑（080500）
已經(jīng)有5人回復
-
材料科學與工程求調(diào)劑
已經(jīng)有5人回復
-
一志愿北京化工大學材料與化工（085600）296求調(diào)劑
已經(jīng)有17人回復
-
考研調(diào)劑
已經(jīng)有7人回復
-
343求調(diào)劑085601
已經(jīng)有3人回復
-
08工科求調(diào)劑286
已經(jīng)有3人回復
-
337求調(diào)劑
已經(jīng)有3人回復
-
329求調(diào)劑，一志愿西北工業(yè)大學，材料工程（085601）
已經(jīng)有10人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 關于linux下的任務管理軟件比如PBS
已經(jīng)有10人回復
- 提取VASP計算前后的坐標轉(zhuǎn)換腳本（shell）
已經(jīng)有9人回復
- python之蟒開啟理論計算化學的新時代
已經(jīng)有22人回復
- 計算費米面時遇到的問題
已經(jīng)有26人回復
- BASH腳本的一個小問題
已經(jīng)有3人回復
- Linux下怎樣使用批處理讓任務一個接一個算？
已經(jīng)有13人回復
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
★ ★
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
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
2樓2014-02-21 22:59:53
小小鑫000
鐵蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 6.6
- 帖子: 267
- 在線: 111.8小時
- 蟲號: 2648268
- 注冊: 2013-09-11
- 性別: MM
- 專業(yè): 理論和計算化學
|
前輩，謝謝您的指導，最后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) 這一行也要加上啊!前輩該怎嗎改？  |
3樓2014-02-22 11:49:34
小小鑫000
鐵蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 6.6
- 帖子: 267
- 在線: 111.8小時
- 蟲號: 2648268
- 注冊: 2013-09-11
- 性別: MM
- 專業(yè): 理論和計算化學
|
前輩，這樣寫的話還對嗎？ 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +7 | 李李不服輸 2026-03-26 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 070305高分子化學與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+10
|
nnnnnnn5 2026-03-25 | 10/500 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 藥學105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 一志愿上海理工能源動力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 0703化學一志愿南京師范大學303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|