分極率

Gaussianの振動解析結果から分極率を抽出する。
振動解析計算は以下を参考に計算。
名古屋大学情報連携基盤センターニュース, Vol.６, No.２, (2007. ５)

計算結果のlogファイル（XXX.log）の以下の「iso」行の「au」単位の「0.422550D+03」を抜き出す。

Dipole polarizability, Alpha (dipole orientation).
  (esu units = cm**3 , SI units = C**2 m**2 J**-1)
  Alpha(0;0):
                (au)            (10**-24 esu)      (10**-40 SI)
    iso        0.422550D+03      0.626154D+02      0.696691D+02
    aniso      0.105316D+02      0.156063D+01      0.173643D+01
    xx         0.426006D+03      0.631276D+02      0.702389D+02
    yx        -0.458695D-01     -0.679716D-02     -0.756287D-02
    yy         0.426114D+03      0.631436D+02      0.702568D+02
    zx        -0.472329D-01     -0.699920D-02     -0.778766D-02
    zy         0.132431D-01      0.196243D-02      0.218349D-02
    zz         0.415530D+03      0.615751D+02      0.685116D+02

以下が作成した抽出スクリプト。Excelにコピペできるよう「D」を「E」へ置換。
pola.sh

#!/bin/bash
log=$1

PolaLine=`grep " iso " ${log} | tail -1 | sed -e 's/D/E/g'`
POL1=`echo ${PolaLine} | cut -d" " -f2`

#EOF

以下が実行結果。

$./getpola XXX.log
0.422550E+03