#!/bin/csh -f

## *** 格子定数のリスト ***
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 )

## *** 前回のファイルが残っている場合削除する ***
rm analysis/fccFeFMG.dat
rm analysis/fccFeAFMI.dat
rm analysis/fccFeNMG.dat

## *** 強磁性初期ポテンシャルの作成 ***
specx < in/fccFeFMG.in > out/fccFeFMG.out
fmg < fccFeafmI

## *** 繰り返し計算 ***
foreach ABOHR ( ${ABOHR_LIST} )
  ## 強磁性
  sed 's/'ABOHR'/'${ABOHR}'/g' template/fccFeFMG0.in > in/fccFeFMG_${ABOHR}.in
  specx < in/fccFeFMG_${ABOHR}.in > out/fccFeFMG_${ABOHR}.out
  sed -e '$d' data/fccFeFMG.info  | sed -n '$p' >> analysis/fccFeFMG.dat
  ## 第一種反強磁性
  sed 's/'ABOHR'/'${ABOHR}'/g' template/fccFeAFMI0.in > in/fccFeAFMI_${ABOHR}.in
  specx < in/fccFeAFMI_${ABOHR}.in > out/fccFeAFMI_${ABOHR}.out
  sed -e '$d' data/fccFeAFMI.info | sed -n '$p' >> analysis/fccFeAFMI.dat
  ## 非磁性
  sed 's/'ABOHR'/'${ABOHR}'/g' template/fccFeNMG0.in > in/fccFeNMG_${ABOHR}.in
  specx < in/fccFeNMG_${ABOHR}.in > out/fccFeNMG_${ABOHR}.out
  sed -e '$d' data/fccFeNMG.info  | sed -n '$p' >> analysis/fccFeNMG.dat
end

grep "spin moment" out/fccFeAFMI_*.out | awk 'NR % 4 == 1;NR % 4 == 2'