echo Starting simulation
./LaunchLoopTransient
echo Plotting

# Make animation
echo set terminal gif size 1000, 800 > gnuplotScript.gp
echo set xrange [-50:50] >> gnuplotScript.gp
echo set tmargin 0 >> gnuplotScript.gp
echo set lmargin 10 >> gnuplotScript.gp
echo set rmargin 10 >> gnuplotScript.gp
for file in $(ls data/. -1)
do
  echo set output \"data/$file.gif\" >> gnuplotScript.gp
  echo set multiplot layout 3, 1 title \"Launch Loop Transient\" >> gnuplotScript.gp

# Plot 1, side view
  echo set yrange [-5:50] >> gnuplotScript.gp
  echo set ylabel \"z / km\" >> gnuplotScript.gp
  echo unset xtics >> gnuplotScript.gp
  echo unset xlabel >> gnuplotScript.gp
  echo set tmargin 0 >> gnuplotScript.gp
  echo set bmargin 1 >> gnuplotScript.gp
  echo plot \"data/$file\" using 1:3:4 lc rgb variable with lines title \"$file\" >> gnuplotScript.gp

# Plot 2, top view
  echo set yrange [-2000:2000] >> gnuplotScript.gp
  echo set ylabel \"y / m\" >> gnuplotScript.gp
  echo set xtics >> gnuplotScript.gp
  echo set xlabel \"x / km\" >> gnuplotScript.gp
  echo set bmargin 2 >> gnuplotScript.gp
  echo plot \"data/$file\" using 1:2:4 lc rgb variable with lines notitle >> gnuplotScript.gp

# Plot 3, forces
  echo set yrange [5:20] >> gnuplotScript.gp
  echo set ylabel \"log\(F\)\" >> gnuplotScript.gp
  echo unset xtics >> gnuplotScript.gp
  echo unset xlabel >> gnuplotScript.gp
  echo set tmargin 2 >> gnuplotScript.gp
  echo set bmargin 1 >> gnuplotScript.gp
  echo plot \"data/$file\" using 1:5 with lines title \"Track tension / N\", \"data/$file\" using 1:6 with lines title \"Guy wire tension / N\", \"data/$file\" using 1:7 with lines title \"rho*sqr\(v\) / N\" >> gnuplotScript.gp

  echo unset multiplot >> gnuplotScript.gp
  echo set output >> gnuplotScript.gp
done
gnuplot gnuplotScript.gp
echo Composing
gifsicle -O2 --delay=10 --loop data/*.gif > Animation.gif 

echo Cleaning up
rm data/*