cat << EOT > inp_1.csv
inp_1
4,3,2,2,2,2,2
210000,33.4D2,1510.0D4
210000,43.0D2,2770.0D4
1,2,1
3,4,1
2,4,2
0,0
0,4000
6000,0
6002,4000
1
3
1
3
1
3
2,0,-10000,0
4,0,-10000,0
EOT

cat << EOT > inp_2.csv
inp_2
6,6,2,2,2,2,2
210000,33.4D2,1510.0D4
210000,43.0D2,2770.0D4
1,2,1
2,3,1
4,5,1
5,6,1
2,5,2
3,6,2
0,0
0,4000
0,8000
6000.0,0
6001.0,4000
6002.0,8000
1
4
1
4
1
4
3,0,-10000,0
6,0,-10000,0
EOT

cat << EOT > inp_3.csv
inp_3
8,9,2,2,2,2,2
210000,33.4D2,1510.0D4
210000,43.0D2,2770.0D4
1,2,1
2,3,1
3,4,1
5,6,1
6,7,1
7,8,1
2,6,2
3,7,2
4,8,2
0,0
0,4000
0,8000
0,12000
6000.0,0
6000.6,4000
6001.2,8000
6002.0,12000
1
5
1
5
1
5
4,0,-10000,0
8,0,-10000,0
EOT

cat << EOT > inp_4.csv
inp_4
12,15,2,2,2,2,2
210000,33.4D2,1510.0D4
210000,43.0D2,2770.0D4
1,2,1
2,3,1
3,4,1
4,5,1
5,6,1
7,8,1
8,9,1
9,10,1
10,11,1
11,12,1
2,8,2
3,9,2
4,10,2
5,11,2
6,12,2
0, 0
0, 4000
0, 8000
0, 12000
0, 16000
0, 20000
6000.0, 0
6000.4, 4000
6000.8, 8000
6001.2, 12000
6001.6, 16000
6002.0, 20000
1
7
1
7
1
7
6,0,-10000,0
12,0,-10000,0
EOT

gfortran -o f90_fem_m_1 f90_fem_m_1.f90
gfortran -o f90_fig_gfrm f90_fig_gfrm.f90

for i in 1 2 3 4
do
if test ${i} -eq 1 ; then
    n1=2; n2=2; n3=4
fi
if test ${i} -eq 2 ; then
    n1=3; n2=2; n3=6
fi
if test ${i} -eq 3 ; then
    n1=4; n2=2; n3=8
fi
if test ${i} -eq 4 ; then
    n1=6; n2=2; n3=12
fi
./f90_fem_m_1 inp_${i}.csv out_${i}_999.csv 1000 500 0.1 1.000
./f90_fem_m_1 inp_${i}.csv out_${i}_900.csv 1000 500 0.1 0.900
./f90_fem_m_1 inp_${i}.csv out_${i}_800.csv 1000 500 0.1 0.800
./f90_fem_m_1 inp_${i}.csv out_${i}_700.csv 1000 500 0.1 0.700
./f90_fem_m_1 inp_${i}.csv out_${i}_600.csv 1000 500 0.1 0.600
./f90_fem_m_1 inp_${i}.csv out_${i}_500.csv 1000 500 0.1 0.500
./f90_fem_m_1 inp_${i}.csv out_${i}_400.csv 1000 500 0.1 0.400
./f90_fem_m_1 inp_${i}.csv out_${i}_300.csv 1000 500 0.1 0.300
./f90_fem_m_1 inp_${i}.csv out_${i}_200.csv 1000 500 0.1 0.200
./f90_fem_m_1 inp_${i}.csv out_${i}_100.csv 1000 500 0.1 0.100
./f90_fem_m_1 inp_${i}.csv out_${i}_000.csv 1000 500 0.1 0.000
./f90_fig_gfrm out_${i}_999.csv dat_999_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_900.csv dat_900_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_800.csv dat_800_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_700.csv dat_700_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_600.csv dat_600_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_500.csv dat_500_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_400.csv dat_400_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_300.csv dat_300_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_200.csv dat_200_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_100.csv dat_100_${i}.csv ${n1} ${n2} ${n3}
./f90_fig_gfrm out_${i}_000.csv dat_000_${i}.csv ${n1} ${n2} ${n3}
done


#cat << EOT > gpl.txt
#set terminal png
#set output 'fig_${i}.png'
#set datafile separator ","
#set xrange[-500:1000]
#set yrange[-2000000:0]
#plot 'dat_${i}.csv' using 2:1 title 'u' with lines lw 2, \
#     'dat_${i}.csv' using 3:1 title 'v' with lines lw 2
#EOT
#gnuplot gpl.txt