#!/bin/csh -f

## *** フォルダ構造 ***
## fccMn/─┬─analysis/
##         ├─in/─L10fccMn_FMG.in
##         ├─out/
##         ├─data/
##         ├─template/┬─L10fccMn_AFM_Template.in
##         │           ├─L10fccMn_NMG_Template.in
##         │           └─L10fccMn_FMG_Template.in
##         ├─L10fccMn
##         ├─L10fccMn.sh
##         └─L10fccMn-Result.sh

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

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

## *** 強磁性初期ポテンシャルの作成 ***
if ( ! -e data/${POTENTIAL}_AFM ) then
  specx < in/${PROJECT}_FMG.in > out/${PROJECT}_FMG.out
  fmg < ${PROJECT}
endif

## *** 繰り返し計算 ***
foreach ABOHR ( ${ABOHR_LIST} )
  ## 強磁性
   if ( ! -e data/${POTENTIAL}_FMG_${ABOHR} ) then
    if ( -e data/${POTENTIAL}_FMG ) then
      cp data/${POTENTIAL}_FMG data/${POTENTIAL}_FMG_${ABOHR}
    endif
  endif
  sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_FMG_Template.in > in/${PROJECT}_FMG_${ABOHR}.in
  specx < in/${PROJECT}_FMG_${ABOHR}.in > out/${PROJECT}_FMG_${ABOHR}.out
  cp data/${POTENTIAL}_FMG_${ABOHR} data/${POTENTIAL}_FMG
  
  ## 反強磁性
  if ( ! -e data/${POTENTIAL}_AFM_${ABOHR} ) then
    if ( -e data/${POTENTIAL}_AFM ) then
      cp data/${POTENTIAL}_AFM data/${POTENTIAL}_AFM_${ABOHR}
    endif
  endif
  sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_AFM_Template.in > in/${PROJECT}_AFM_${ABOHR}.in
  specx < in/${PROJECT}_AFM_${ABOHR}.in > out/${PROJECT}_AFM_${ABOHR}.out
  cp data/${POTENTIAL}_AFM_${ABOHR} data/${POTENTIAL}_AFM
  
  ## 非磁性
  if ( ! -e data/${POTENTIAL}_NMG_${ABOHR} ) then
    if ( -e data/${POTENTIAL}_NMG ) then
      cp data/${POTENTIAL}_NMG data/${POTENTIAL}_NMG_${ABOHR}
    endif
  endif
  sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_NMG_Template.in > in/${PROJECT}_NMG_${ABOHR}.in
  specx < in/${PROJECT}_NMG_${ABOHR}.in > out/${PROJECT}_NMG_${ABOHR}.out
  cp data/${POTENTIAL}_NMG_${ABOHR} data/${POTENTIAL}_NMG
end