#!/bin/csh -f

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

## *** 前回の結果の削除 ***
rm -f analysis/Energy.txt
rm -f analysis/Moment.txt

## *** 格子体積のリスト (Bohr^3) ***
set OMEGA_LIST=( 90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 )
## *** c/aのリスト ***
##set ETA_LIST=( 1.00 1.05 1.10 1.15 1.20 1.25 1.30 1.35 1.40 1.45 1.50 )
set ETA_LIST=( 1.00 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.10 1.10 1.11 1.12 1.13 1.14 1.15 1.16 1.17 1.18 1.19 1.20 1.21 1.22 1.23 1.24 1.25 1.26 1.27 1.28 1.29 1.30 1.31 1.32 1.33 1.34 1.35 1.36 1.37 1.38 1.39 1.40 1.41 1.42 1.43 1.44 1.45 1.46 1.47 1.48 1.49 1.50 )

## *** 計算結果の出力 ****
foreach OMEGA ( ${OMEGA_LIST} )
  foreach ETA ( ${ETA_LIST} )
    ## 全エネルギー
    tail -n 1 data/${POTENTIAL}_${OMEGA}_${ETA}.info | awk '{printf $2}' >> analysis/Energy.txt
    printf " " >> analysis/Energy.txt
    ## 磁気モーメント
    tail -n 1 data/${POTENTIAL}_${OMEGA}_${ETA}.info | awk '{printf $3}' >> analysis/Moment.txt
    printf " " >> analysis/Moment.txt
  end
  printf "\n" >> analysis/Energy.txt
  printf "\n" >> analysis/Moment.txt
end