#!/bin/csh
#
# run the Boost-invariant BlastWave code, varying one parameter
#
#
set Tnom      = 0.1
set rho0nom   = 0.9
set rhoanom   = 0.05
set Rynom     = 13.0
set Rxnom     = 11.0
set asnom     = 0.0001
set tau0nom   = 9.0
set deltaunom = 2.0
#
set values  = ( -0.05 -0.02 0.02 0.05 )
#
@ iter = 1
start_loop:
echo Working on calculation $iter of $#values
set rhoa = $values[$iter]
touch input.inp; /bin/rm input.inp
cat > input.inp << EOF
$Rynom
$Rxnom
$Tnom
$rho0nom
$rhoa
$tau0nom
$asnom
$deltaunom
0
0.0 2.0 0.05
EOF
~lisa/MikesCommonPrograms/bin/BoostInvariantBlastWave.exe < input.inp
mkdir Calc$iter
mv *BW Calc$iter
mv input.inp Calc$iter
mv define_title.kumac Calc$iter
if ( $iter < $#values ) then
  @ iter = $iter + 1
  goto start_loop
endif
#
exit
#
