#!/bin/csh -f

## 計算する最後の原子番号+1
set NMAX=104



## CPU並列化数:並列化なし
setenv OMP_NUM_THREADS 1

## 計算結果の保存先ファイルを削除
rm -f analysis/RMTA.txt
rm -f analysis/RMTB.txt

## 原子番号初期化
set NA=0
## 総当り的に繰り返し
while ( $NA < $NMAX )
  ## 原子番号初期化
  set NB=0
  ## 計算結果の保存先初期化
  set RMTA_LIST=( )
  set RMTB_LIST=( )
  while ( $NB < $NMAX )
    ## 入力ファイルの作成
    sed 's/'NA'/'${NA}'/g' template/AB.in | sed 's/'NB'/'${NB}'/g' > in/AB_${NA}_${NB}.in
    ## 計算
    specx <in/AB_${NA}_${NB}.in >out/AB_${NA}_${NB}.out
    ## 計算結果を保存
    set RMTA=`grep "rmt=" out/AB_${NA}_${NB}.out | awk '{print $3}' | head -n 1`
    set RMTB=`grep "rmt=" out/AB_${NA}_${NB}.out | awk '{print $3}' | tail -n 1`
    set RMTA_LIST=( ${RMTA_LIST} ${RMTA} )
    set RMTB_LIST=( ${RMTB_LIST} ${RMTB} )
    ## 次の原子へ
    @ NB ++
 end
 ## 計算結果をファイルへ上書き
 echo ${RMTA_LIST} >> analysis/RMTA.txt
 echo ${RMTB_LIST} >> analysis/RMTB.txt
 ## 次の原子へ
 @ NA++
end