#!/bin/csh -f set OMEGA_LIST=( 140 142 144 146 148 150 152 154 156 158 160 ) set ETA_LIST=( 1.60 1.61 1.62 1.63 1.64 1.65 1.66 1.67 1.68 1.69 1.70 ) foreach OMEGA ( ${OMEGA_LIST} ) set ETA0=`echo $ETA_LIST[1]` set A0=`echo "e((1/3)*l(2*${OMEGA}/(sqrt(3)*${ETA0})))" | bc -l` if (`echo "${ETA0} > 2*sqrt(2)/sqrt(3)" | bc -l` == 1) then set RMTB0=`echo "${A0}/2" | bc -l` else set RMTB0=`echo "${A0}*sqrt(1/3+(${ETA0}^2)/4)/2" | bc -l` endif set ETA1=`echo $ETA_LIST[$#ETA_LIST]` set A1=`echo "e((1/3)*l(2*${OMEGA}/(sqrt(3)*${ETA1})))" | bc -l` if (`echo "${ETA1} > 2*sqrt(2)/sqrt(3)" | bc -l` == 1) then set RMTB1=`echo "${A1}/2" | bc -l` else set RMTB1=`echo "${A1}*sqrt(1/3+(${ETA1}^2)/4)/2" | bc -l` endif if (`echo "${RMTB0} < ${RMTB1}" | bc -l` == 1) then set RMTB=`echo $RMTB0` else set RMTB=`echo $RMTB1` endif foreach ETA ( ${ETA_LIST} ) if ( ! -e stop ) then echo " OMEGA,ETA= "${OMEGA}" "${ETA} set A=`echo "e((1/3)*l(2*${OMEGA}/(sqrt(3)*${ETA})))" | bc -l | sed -e 's/^\./0./g'` set RMTA=`echo "${RMTB} / e((1/3)*l(2*${OMEGA}/(sqrt(3)*${ETA})))" | bc -l | sed -e 's/^\./0./g'` sed 's/'ABOHR'/'${A}'/g' template/CoTemplate.in | sed 's/'ETA'/'${ETA}'/g' | sed 's/'RMTA'/'${RMTA}'/g' | sed 's/'OMEGA'/'${OMEGA}'/g' > in/Co_${OMEGA}_${ETA}.in specx < in/Co_${OMEGA}_${ETA}.in > out/Co_${OMEGA}_${ETA}.out tail -n 1 data/co_${OMEGA}_${ETA}.info endif end end