#!/bin/csh
#
# run the Boost-invariant BlastWave code, varying one parameter
#
# NOTE: here I vary a_s, but must "co-vary" rho_0 and rho_a
#
# here come nominal values for the parameters
#
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.1   0.3   )
set values2 = ( 0.845 0.612 )
set values3 = ( 0.047 0.034 )
#
@ iter = 1
start_loop:
echo Working on calculation $iter of $#values
set as = $values[$iter]
set rho0 = $values2[$iter]
set rhoa = $values3[$iter]
touch input.inp; /bin/rm input.inp
cat > input.inp << EOF
$Rynom
$Rxnom
$Tnom
$rho0
$rhoa
$tau0nom
$as
$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
#
