#!/bin/csh -f ## *** プロジェクト名 *** set PROJECT="CuZn" ## *** 格子定数 (Bohr) *** ## 格子定数 set ABOHR_LIST=( 5.565 ) ## 濃度 set XIA_LIST=( 100 98 96 94 92 90 88 86 84 82 80 78 76 74 72 70 68 66 64 62 60 58 56 54 52 50 ) ## *** 第一原理計算 **** foreach XIA ( ${XIA_LIST} ) set XIB=`echo "100-${XIA}" | bc -l` foreach ABOHR ( ${ABOHR_LIST} ) if ( ! -e data/${PROJECT}_${XIA}_${ABOHR} ) then if ( -e data/${PROJECT}_${XIA} ) then cp data/${PROJECT}_${XIA} data/${PROJECT}_${XIA}_${ABOHR} endif endif ## テンプレートから入力ファイルを作成 sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}_Template.in | sed 's/'XIA'/'${XIA}'/g' | sed 's/'XIB'/'${XIB}'/g' > in/${PROJECT}_${XIA}_${ABOHR}.in ## 第一原理計算実行 specx < in/${PROJECT}_${XIA}_${ABOHR}.in > out/${PROJECT}_${XIA}_${ABOHR}.out ## 次の格子定数の初期ポテンシャル cp data/${PROJECT}_${XIA}_${ABOHR} data/${PROJECT}_${XIA} end end