 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1768 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
小小鑫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 求各位前輩幫我改正一下，謝謝！ |
回復(fù)此樓» 本帖@通知
» 猜你喜歡
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有8人回復(fù)
-
349求調(diào)劑
已經(jīng)有6人回復(fù)
-
329求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有9人回復(fù)
-
327求調(diào)劑
已經(jīng)有5人回復(fù)
-
求收留
已經(jīng)有5人回復(fù)
-
375求調(diào)劑
已經(jīng)有4人回復(fù)
-
340求調(diào)劑
已經(jīng)有6人回復(fù)
-
298求調(diào)劑
已經(jīng)有3人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于linux下的任務(wù)管理軟件比如PBS
已經(jīng)有10人回復(fù)
- 提取VASP計(jì)算前后的坐標(biāo)轉(zhuǎn)換腳本（shell）
已經(jīng)有9人回復(fù)
- python之蟒開啟理論計(jì)算化學(xué)的新時代
已經(jīng)有22人回復(fù)
- 計(jì)算費(fèi)米面時遇到的問題
已經(jīng)有26人回復(fù)
- BASH腳本的一個小問題
已經(jīng)有3人回復(fù)
- Linux下怎樣使用批處理讓任務(wù)一個接一個算？
已經(jīng)有13人回復(fù)
小小鑫000
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6.6
- 帖子: 267
- 在線: 111.8小時
- 蟲號: 2648268
- 注冊: 2013-09-11
- 性別: MM
- 專業(yè): 理論和計(jì)算化學(xué)
|
前輩，這樣寫的話還對嗎？ 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
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 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
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6.6
- 帖子: 267
- 在線: 111.8小時
- 蟲號: 2648268
- 注冊: 2013-09-11
- 性別: MM
- 專業(yè): 理論和計(jì)算化學(xué)
|
前輩，謝謝您的指導(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) 這一行也要加上啊!前輩該怎嗎改？  |
3樓2014-02-22 11:49:34
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 375求調(diào)劑 +4 | 雨夏整夜 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 一志愿上海理工能源動力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|