#!/bin/csh
#
# run the Boost-invariant BlastWave code.
#
# the three cases here correspond to Fabrice's "best-fit" values for central,midcent,periph
#
#
set Tnom      = ( 0.107  0.106  .0946 )
set rho0nom   = ( 0.883  0.870  0.806 )
set rhoanom   = ( 0.059  0.052  0.038 )
set Rynom     = ( 12.83  11.75  10.07 )
set Rxnom     = ( 12.86  10.24   8.00 )
set asnom     = ( 0.001  0.001  0.001 )
set tau0nom   = ( 8.871  7.389  6.495 )
set deltaunom = ( 0.002  0.759  0.767 )
#
#
@ iter = 3
start_loop:
echo Working on calculation $iter of $#Tnom
touch input.inp; /bin/rm input.inp
cat > input.inp << EOF
$Rynom[$iter]
$Rxnom[$iter]
$Tnom[$iter]
$rho0nom[$iter]
$rhoanom[$iter]
$tau0nom[$iter]
$asnom[$iter]
$deltaunom[$iter]
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 < $#Tnom ) then
  @ iter = $iter + 1
  goto start_loop
endif
#
exit
#
