#!/bin/csh
#
# run the Boost-invariant BlastWave code, varying one parameter
#
# In this run, we fix the value of Rx^2+Ry^2, while varying Ry/Rx
#  We fix Rx^2+Ry^2 = (11^2)+(13^2) = 290
#
#
# NOTE!!!: We set rho_a = 0.0 here -- NOT the usual rhoa_nom=0.05 !!
#
set Tnom      = 0.1
set rho0nom   = 0.9
set rhoanom   = 0.00
set Rynom     = 13.0
set Rxnom     = 11.0
set asnom     = 0.0001
set tau0nom   = 9.0
set deltaunom = 2.0
#
set values  = ( 12.042 14.457 13.0 13.784)
set values2 = ( 12.042  9.0   11.0 10.0 )
#
@ iter = 4
start_loop:
echo Working on calculation $iter of $#values
set Ry = $values[$iter]
set Rx = $values2[$iter]
touch input.inp; /bin/rm input.inp
cat > input.inp << EOF
$Ry
$Rx
$Tnom
$rho0nom
$rhoanom
$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
#
