#!/bin/bash

NU_T=0.0
N_A=10
N_B=10
L_A=1.8
L_B=1.0
GAMMA=80.0
PHI=0.0
X_L=-20.2
X_U=40.2
Y_L=-30.2
Y_U=30.2
T_X=604
T_Y=604

j=0
endj=180

while [ $j -le $endj ]; do

    PHI=`echo $j \+ 0.0 |bc`
    GRAPH=$(printf "%dx%d-%05.1fdegree" $N_A $N_B $PHI)

    i=0
    endi=19

    while [ $i -le $endi ]; do
echo "Loop: $j - $i"

        NU_T=`echo $i \* 0.05 |bc`
        COMMAND=$(printf "./a.out %f %d %d %f %f %f %f %f %f %f %f %d %d out" $NU_T $N_A $N_B $L_A $L_B $GAMMA $PHI $X_L $X_U $Y_L $Y_U $T_X $T_Y)
        $COMMAND
        COMMAND2=$(printf "sed s/OUTFILE/out%03d/" $i)
        COMMAND3=$(printf "sed s/INFILE/out/g")
        COMMAND4=$(printf "sed s/GRAPHTITLE/%s/" "$GRAPH")
        cat template.gplt | $COMMAND2 > TEMP1.gplt
        cat TEMP1.gplt | $COMMAND3 > TEMP2.gplt
        cat TEMP2.gplt | $COMMAND4 > COM.gplt
        gnuplot <COM.gplt
        rm TEMP?.gplt COM.gplt
        i=$(($i+1))
    done

    IPHI=$(echo $PHI |sed 's/.0$//')
    OUTFILENAME=$(printf "Wave%dx%d-%03d.mp4" $N_A $N_B $IPHI)
    COMMAND0='ffmpeg -r 24 -i out%03d.png -an -c:v libx264 -vf format=yuv420p -tune animation '
    COMMANDV=$(echo "$COMMAND0" $OUTFILENAME)
    $COMMANDV
    rm out???.png

    j=$(($j+1))

done