24小時熱門版塊排行榜

返回列表

飛翔……

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 1378.8
散金: 160
紅花: 6
沙發(fā): 1
帖子: 606
在線: 269.8小時
蟲號: 673866
注冊: 2008-12-13
專業(yè): 催化化學(xué)

[求助] sum_dos腳本解析

!/bin/ksh

# Script to sum up DOS files from VASP
if [ $# -lt 3 ]; then
echo "Usage: sum_dos [0] [start] [end]"
echo "Or: sum_dos [1] [j] [k] [l] ..."
exit 1
fi

type=$1
shift 1

# Contiguous
if [ $type -eq 0 ]; then

start=$1
end=$2
shift 2
string=$start.to.$end

echo $start
echo $end

# First file
i=$start
nl=$(wc -l DOS$i | sed 's/://g' | awk '{print $1}')
echo $nl

sed -n '2,'$nl' p' DOS$i | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1, $2, $3, $4, $5, $6, $7, $2+$4+$6, $3+$5+$7}' >| DOS.tmp

i=$((i+1))

while [ $i -le $end ]; do

sed -n '2,'$nl' p' DOS$i | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1, $2, $3, $4, $5, $6, $7, $2+$4+$6, $3+$5+$7}' >| DOS.tmp2

paste DOS.tmp DOS.tmp2 | awk '{printf "%15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f \n", $1, $2+$11, $3+$12, $4+$13, $5+$14, $6+$15, $7+$16, $8+$17, $9+$18}' >| DOS.tmp3

mv -f DOS.tmp3 DOS.tmp

i=$((i+1))

done

else

start=$1
shift 1

# First file
i=$start
nl=$(wc -l DOS$i | sed 's/://g' | awk '{print $1}')
echo $nl

sed -n '2,'$nl' p' DOS$i | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1, $2, $3, $4, $5, $6, $7, $2+$4+$6, $3+$5+$7}' >| DOS.tmp
string=$start
for i in "$@"; do

echo $i
string=$(echo $string.$i)
sed -n '2,'$nl' p' DOS$i | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1, $2, $3, $4, $5, $6, $7, $2+$4+$6, $3+$5+$7}' >| DOS.tmp2

paste DOS.tmp DOS.tmp2 | awk '{printf "%15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f \n", $1, $2+$11, $3+$12, $4+$13, $5+$14, $6+$15, $7+$16, $8+$17, $9+$18}' >| DOS.tmp3

mv -f DOS.tmp3 DOS.tmp
done

fi

mv DOS.tmp DOS.SUM.$string
rm -f DOS.tmp2

exit 0

求高手指點該腳本是怎樣實現(xiàn)求和的，求得是那些和？

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【分享】最新VASP腳本大全　支持VASP 5.3 (有split_dos; VP; Sumdos; sumup_dos 等) 已經(jīng)有615人回復(fù)
【求助】dos腳本求助已經(jīng)有5人回復(fù)

1樓 2014-07-05 20:32:04

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

y1ding

鐵桿木蟲 (著名寫手)

1ST強帖: 1
應(yīng)助: 61 (初中生)
貴賓: 0.33
金幣: 5959.3
散金: 1
紅花: 21
帖子: 1884
在線: 491.1小時
蟲號: 142265
注冊: 2005-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★ ★ ★
franch: 金幣+5, 謝謝回帖交流。。。 2014-07-08 08:35:20

注意這兩句
sed -n '2,'$nl' p' DOS$i | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1, $2, $3, $4, $5, $6, $7, $2+$4+$6, $3+$5+$7}' >| DOS.tmp2

paste DOS.tmp DOS.tmp2 | awk '{printf "%15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f %15.8f \n", $1, $2+$11, $3+$12, $4+$13, $5+$14, $6+$15, $7+$16, $8+$17, $9+$18}' >| DOS.tmp3

linux下，看一下sed awk命令的含義即可

贊一下

回復(fù)此樓

2樓2014-07-05 20:50:00

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

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