#!/bin/csh
#
# run the Boost-invariant BlastWave code, varying one parameter
#  - here we set Rx=Ry
#
# NOTE: here I vary rho_0 AND rho_a, but keep rho_a/rho_0 constant.
#  this is to check the statement in the paper that v2 decreases
#  when rho_0 increases only because rho_a/rho_0 is going down
#  (i.e. there is less "relative" oscillation in the flow field)
#  This plot is not envisioned for the paper
#
# "nominal" rho_a/rho_0 = 0.05/0.9 = 0.05555
#   --> set rho_a = 0.05555* rho_0
#
# 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     = 13.0
set asnom     = 0.0001
set tau0nom   = 9.0
set deltaunom = 2.0
#
set values  = ( 0.7    1.1 )
set values2 = ( 0.0389 0.0611)
#
@ iter = 1
start_loop:
echo Working on calculation $iter of $#values
set rho0 = $values[$iter]
set rhoa = $values2[$iter]
touch input.inp; /bin/rm input.inp
cat > input.inp << EOF
$Rynom
$Rxnom
$Tnom
$rho0
$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
#
