#!/bin/csh -f

## *** フォルダ構造 ***
## Ni3Mn/─┬─analysis/
##         ├─in/
##         ├─out/
##         ├─data/
##         ├─template/─L12Ni3Mn_Template.in
##         ├─L12Ni3Mn.sh
##         └─L12Ni3Mn-Result.sh


## *** OpenMPの設定 ***
setenv OMP_STACKSIZE    20M
limit stacksize unlimited
setenv OMP_NUM_THREADS 4

## *** プロジェクト名 ***
set PROJECT="L12Ni3Mn"
## ポテンシャルファイル名
set POTENTIAL=${PROJECT}

## *** 格子定数のリスト ***
set ABOHR_LIST=( 7.4 7.3 7.2 7.1 7.0 6.9 6.8 6.7 6.6 6.5 6.4 6.3 6.2 6.1 6.0 )

## *** 繰り返し計算 ***
foreach ABOHR ( ${ABOHR_LIST} )
  if ( ! -e data/${POTENTIAL}_${ABOHR} ) then ## ポテンシャルがなければコピー
    if ( -e data/${POTENTIAL} ) then
      cp data/${POTENTIAL} data/${POTENTIAL}_${ABOHR}
    endif
  endif
  ## 自己無撞着計算
  sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_Template.in > in/${PROJECT}_${ABOHR}.in
  specx < in/${PROJECT}_${ABOHR}.in > out/${PROJECT}_${ABOHR}.out
  ## 次回の初期ポテンシャルをコピー
  cp data/${POTENTIAL}_${ABOHR} data/${POTENTIAL}
end