#!/bin/csh -f ## *** プロジェクト名 *** set PROJECT="bccPt" ## *** 格子定数のリスト (Bohr) *** set ABOHR_LIST=( 6.60 6.50 6.40 6.30 6.20 6.10 6.00 5.90 5.80 5.70 5.60 5.50 5.40 5.30 5.20 ) ## *** 実行ファイルのパス *** set EXEC="specx" #set EXEC="~/kkr/cpa2002v009c/specx" #set EXEC="~/kkr/cpa2002v009c/specx.ew" #set EXEC="~/kkr/cpa2002v009c/specx.sicoff" #set EXEC="~/kkr/cpa2002v009c/specx.sicon" if ( -e analysis/${PROJECT}.txt ) then cat analysis/${PROJECT}.txt >> analysis/${PROJECT}.txt.back endif echo "# a(Bohr) Energy(Ry)" > analysis/${PROJECT}.txt ## *** 第一原理計算 **** foreach ABOHR ( ${ABOHR_LIST} ) ## 初期ポテンシャルが無ければ前回のポテンシャルをコピー if ( ! -e data/${PROJECT}_${ABOHR} ) then cp data/${PROJECT} data/${PROJECT}_${ABOHR} endif ## *** セルフコンシステント計算 *** ## テンプレートから入力ファイルを作成 sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_Template.in > in/${PROJECT}_${ABOHR}.in ## 第一原理計算実行 ${EXEC} < in/${PROJECT}_${ABOHR}.in > out/${PROJECT}_${ABOHR}.out ## テンプレートから入力ファイルを作成 sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_dos_Template.in > in/${PROJECT}_dos_${ABOHR}.in ## 第一原理計算実行 ${EXEC} < in/${PROJECT}_dos_${ABOHR}.in > out/${PROJECT}_dos_${ABOHR}.out set ENE=`grep "total energy" out/${PROJECT}_${ABOHR}.out | sed -e s/total//g -e s/energy=//g` echo ${ABOHR} ${ENE} >> analysis/${PROJECT}.txt ## ポテンシャルのバックアップ cp data/${PROJECT}_${ABOHR} data/${PROJECT} end