#!/bin/csh -f

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

## *** 格子定数 (Bohr) ***
## 銅の格子定数
set ABOHR_LIST=( 6.50 6.55 6.60 6.65 6.70 6.75 6.80 6.85 6.90 6.95 7.00 7.05 7.10 7.15 7.20 )
## 鉄の格子定数
#set ABOHR_LIST=( 5.20 5.25 5.30 5.35 5.40 5.45 5.50 5.55 5.60 )
## 交換相関ポテンシャル
set VEXC_LIST=( mjw pbe )

## *** 第一原理計算 ****
foreach VEXC ( ${VEXC_LIST} )
  foreach ABOHR ( ${ABOHR_LIST} )
    if ( ! -e ${PROJECT}_${VEXC}_${ABOHR} ) then
      if ( -e ${PROJECT}_${VEXC} ) then
        cp data/${PROJECT}_${VEXC} data/${PROJECT}_${VEXC}_${ABOHR}
      endif
    endif
    ## テンプレートから入力ファイルを作成
    sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_Template.in | sed 's/'VEXC'/'${VEXC}'/g' > in/${PROJECT}_${VEXC}_${ABOHR}.in
    ## 第一原理計算実行
    specx < in/${PROJECT}_${VEXC}_${ABOHR}.in > out/${PROJECT}_${VEXC}_${ABOHR}.out
    ## 次の格子定数の初期ポテンシャル
    cp data/${PROJECT}_${VEXC}_${ABOHR} data/${PROJECT}_${VEXC}
  end
end