版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

BBQY

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 694.1
紅花: 2
帖子: 56
在線: 15.5小時(shí)
蟲號(hào): 1343899
注冊: 2011-07-12
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

[求助] shell腳本編程已有1人參與

我有個(gè)QB.out文件，文件中有
CARTESIAN COORDINATES

1 H    16.074719860    8.229529334 25.103862190
2 C    16.465993520    9.232996370 25.301990040
3 C    17.415993140    9.685996190 24.205990470
4 N    16.735993410    9.949996080 22.943990970
5 C    17.340993170 10.058996040 21.764991430
6 N    16.616993460 10.307995940 20.682991860
7 N    18.659992650    9.916996100 21.652991480
8 H    15.602017307    9.906518523 25.370615233
9 H    16.966838633    9.240617684 26.275695958
  10 H    17.924047738 10.622953044 24.553804455
，其中文件中有2個(gè)CARTESIAN COORDINATES，如果我想提取第二個(gè)CARTESIAN COORDINATES這個(gè)關(guān)鍵詞下面的第3-9行坐標(biāo)，并保留小數(shù)點(diǎn)后4位有效數(shù)字到QB.coods文件用shell怎樣寫腳本呀，shell新手，還望大家多多幫忙呀，不勝感激~

回復(fù)此樓

» 猜你喜歡

336求調(diào)劑已經(jīng)有3人回復(fù)
306求調(diào)劑已經(jīng)有9人回復(fù)
工科0856求調(diào)劑已經(jīng)有4人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有10人回復(fù)
081700 調(diào)劑 267分已經(jīng)有4人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷已經(jīng)有9人回復(fù)
263求調(diào)劑已經(jīng)有9人回復(fù)
求老師收我已經(jīng)有3人回復(fù)
招08考數(shù)學(xué) 已經(jīng)有12人回復(fù)
接收2026碩士調(diào)劑(學(xué)碩+專碩) 已經(jīng)有4人回復(fù)

1樓 2016-06-07 22:20:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

czyzsu

專家顧問 (文壇精英)

專家經(jīng)驗(yàn): +3724
應(yīng)助: 3321 (副教授)
貴賓: 1
金幣: 308806
散金: 2465
紅花: 685
帖子: 28690
在線: 7436.6小時(shí)
蟲號(hào): 150213
注冊: 2005-12-29
專業(yè): 化學(xué)生物學(xué)與生物有機(jī)化學(xué)
管轄: 有機(jī)交流

【答案】應(yīng)助回帖

#!/bin/bash
for i in *.out
do
NAME=${i%.out}
NATOM=` grep "NAtoms=" $i | tail -1 | awk ' {print $2} ' `
sed ' /Standard orientation: /, /Rotational constants/! d' ${i} | grep -v -E ' I| C| N| D| R| \-\-\-' | tail -$NATOM | awk ' {gsub(17, "Cl", $2) ;
gsub(16,"S",$2); gsub("1", "H", $2); gsub(6, "C", $2); gsub(7, "N", $2); gsub(8, "O", $2); printf "%5s %10.6f %10.6f %10.6f\n",$2,$4,$5,$6}' > ${NAME}.xyz
done
——————————————————————
樓主，將以上保存為out2xyz.sh,chmod +x out2xyz.sh 在當(dāng)前目錄下執(zhí)行./out2xyz.sh,就可以批處理提取你要的數(shù)據(jù)。數(shù)據(jù)有效數(shù)字可以通過調(diào)整%10.6f這個(gè)格式來得到你要的有效數(shù)字。
如果沒有得到你要的數(shù)據(jù)，適當(dāng)調(diào)整該語句（NATOM=` grep "NAtoms=" $i | tail -1 | awk ' {print $2} ' `）就好了。

贊一下

回復(fù)此樓

2樓2016-11-15 20:07:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 BBQY 的主題更新