#!/bin/csh -f

## *** プロジェクト名 ***
set PROJECT="PdRh"

## *** 格子定数のリスト (Bohr) ***
set ABOHR_LIST=( 7.00 7.05 7.10 7.15 7.20 7.25 7.30 7.35 7.40 7.45 7.50 7.55 7.60 7.65 7.70 7.75 7.80 7.85 7.90 7.95 8.00 )
set XIRH_LIST=( 100 0 10 20 30 40 50 60 70 80 90 )
#set XIRH_LIST=( 100 0 )

## *** 第一原理計算 ****
foreach XIRH ( ${XIRH_LIST} )
  set XIPD=`echo "100-${XIRH}" | bc -l`
  foreach ABOHR ( ${ABOHR_LIST} )
    ## *** セルフコンシステント計算 ***
    ## テンプレートから入力ファイルを作成
    sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_Template.in | sed 's/'XIRH'/'${XIRH}'/g' | sed 's/'XIPD'/'${XIPD}'/g' > in/${PROJECT}_${XIRH}_${ABOHR}.in
    ## 第一原理計算実行
    specx < in/${PROJECT}_${XIRH}_${ABOHR}.in > out/${PROJECT}_${XIRH}_${ABOHR}.out

    ## *** 状態密度計算 ***
    sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_DOS_Template.in | sed 's/'XIRH'/'${XIRH}'/g' | sed 's/'XIPD'/'${XIPD}'/g' > in/${PROJECT}_DOS_${XIRH}_${ABOHR}.in
    specx < in/${PROJECT}_DOS_${XIRH}_${ABOHR}.in > out/${PROJECT}_DOS_${XIRH}_${ABOHR}.out
  end
end