Скрипт расчета объема энергии

Question

Привет я хочу написать сценарий оболочки для расчета энергии и объема для различных постоянных решетки. Для сближения энергии используют сиесту. Может кто-нибудь помочь, я попробовал эту работу вручную, но теперь хочу написать программу оболочки, используя команды sed. Я много пробовал, но не получилось.

пожалуйста, скажи мне короткий путь. Я могу взять постоянную решетки от 3,60 до 4,40 и SZ для MgO.

#!/usr/bin/env bash 
SIESTA=~code/sanabin/siesta
#
#
name=amin+'(IT-1)'*delta | bc
a0=4.0 | bc
amin==0.92*$a0
amax=1.10*$a0
IT=1
MAXIT=5
MAXITm1=$MAXIT-1
delta=($amax-$amin)/$MAXITm1

#
while  [ $IT -le $MAXIT ] ; 

#
do

mkdir $name
cp input.fdf $name
cp *.psf $name
cd $name
sed 's/LatticeConstant "Ang"/$name/' <input.fdf>name.fdf
$SIESTA < $name.fdf > $name.out
Total=$(grep 'siesta:         Total =' out | cut -c25-40)
#---------
cd ..
done
exit

Я пробовал это, но не работает :(