#!/bin/csh
#
# run the Boost-invariant BlastWave code, varying a_s
#    while at the same time varying rho_0 so as to keep <rho> constant..
#
# here come nominal values for the parameters
#
set Tnom      = 0.1
set rho0nom   = 0.9
set rhoanom   = 0.0
set Rynom     = 13.0
set Rxnom     = 13.0
set asnom     = 0.0001
set tau0nom   = 9.0
set deltaunom = 2.0
#
set values =  ( 0.0001 0.1    0.2    0.3   )
set values2 = ( 0.9    0.845  0.732  0.612 )
#
@ iter = 1
start_loop:
echo Working on calculation $iter of $#values
set as = $values[$iter]
set rho0 = $values2[$iter]
touch input.inp; /bin/rm input.inp
cat > input.inp << EOF
$Rynom
$Rxnom
$Tnom
$rho0
$rhoanom
$tau0nom
$as
$deltaunom
1
0.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
#
