cat << EOT > gpl.txt
set terminal png
set output 'fig_pd1.png'
set datafile separator ","
set xrange[-2000:4000]
set yrange[0:2000]
#set xtics 0.1
#set ytics 0.1
#set format x '%2.1f'
#set format y '%2.1f'
set grid
set xlabel 'Displacement u, v (mm)'
set ylabel 'Load P (kN)'
set key right top
plot 'dat_999_1.csv' using 2:(-\$1/1000) title 'gamma=1.0,u' with lines lw 2, \
     'dat_999_1.csv' using 3:(-\$1/1000) title 'gamma=1.0,v' with lines lw 2, \
     'dat_600_1.csv' using 2:(-\$1/1000) title 'gamma=0.6,u' with lines lw 2, \
     'dat_600_1.csv' using 3:(-\$1/1000) title 'gamma=0.6,v' with lines lw 2, \
     'dat_300_1.csv' using 2:(-\$1/1000) title 'gamma=0.3,u' with lines lw 2, \
     'dat_300_1.csv' using 3:(-\$1/1000) title 'gamma=0.3,v' with lines lw 2
EOT
gnuplot gpl.txt

cat << EOT > gpl.txt
set terminal png
set output 'fig_pd2.png'
set datafile separator ","
set xrange[-2000:4000]
set yrange[0:2000]
#set xtics 0.1
#set ytics 0.1
#set format x '%2.1f'
#set format y '%2.1f'
set grid
set xlabel 'Displacement u, v (mm)'
set ylabel 'Load P (kN)'
set key right top
plot 'dat_999_2.csv' using 2:(-\$1/1000) title 'gamma=1.0,u' with lines lw 2, \
     'dat_999_2.csv' using 3:(-\$1/1000) title 'gamma=1.0,v' with lines lw 2, \
     'dat_600_2.csv' using 2:(-\$1/1000) title 'gamma=0.6,u' with lines lw 2, \
     'dat_600_2.csv' using 3:(-\$1/1000) title 'gamma=0.6,v' with lines lw 2, \
     'dat_300_2.csv' using 2:(-\$1/1000) title 'gamma=0.3,u' with lines lw 2, \
     'dat_300_2.csv' using 3:(-\$1/1000) title 'gamma=0.3,v' with lines lw 2
EOT
gnuplot gpl.txt

cat << EOT > gpl.txt
set terminal png
set output 'fig_pd3.png'
set datafile separator ","
set xrange[-2000:4000]
set yrange[0:2000]
#set xtics 0.1
#set ytics 0.1
#set format x '%2.1f'
#set format y '%2.1f'
set grid
set xlabel 'Displacement u, v (mm)'
set ylabel 'Load P (kN)'
set key right top
plot 'dat_999_3.csv' using 2:(-\$1/1000) title 'gamma=1.0,u' with lines lw 2, \
     'dat_999_3.csv' using 3:(-\$1/1000) title 'gamma=1.0,v' with lines lw 2, \
     'dat_600_3.csv' using 2:(-\$1/1000) title 'gamma=0.6,u' with lines lw 2, \
     'dat_600_3.csv' using 3:(-\$1/1000) title 'gamma=0.6,v' with lines lw 2, \
     'dat_300_3.csv' using 2:(-\$1/1000) title 'gamma=0.3,u' with lines lw 2, \
     'dat_300_3.csv' using 3:(-\$1/1000) title 'gamma=0.3,v' with lines lw 2
EOT
gnuplot gpl.txt

cat << EOT > gpl.txt
set terminal png
set output 'fig_pd4.png'
set datafile separator ","
set xrange[-2000:4000]
set yrange[0:2000]
#set xtics 0.1
#set ytics 0.1
#set format x '%2.1f'
#set format y '%2.1f'
set grid
set xlabel 'Displacement u, v (mm)'
set ylabel 'Load P (kN)'
set key right top
plot 'dat_999_4.csv' using 2:(-\$1/1000) title 'gamma=1.0,u' with lines lw 2, \
     'dat_999_4.csv' using 3:(-\$1/1000) title 'gamma=1.0,v' with lines lw 2, \
     'dat_600_4.csv' using 2:(-\$1/1000) title 'gamma=0.6,u' with lines lw 2, \
     'dat_600_4.csv' using 3:(-\$1/1000) title 'gamma=0.6,v' with lines lw 2, \
     'dat_300_4.csv' using 2:(-\$1/1000) title 'gamma=0.3,u' with lines lw 2, \
     'dat_300_4.csv' using 3:(-\$1/1000) title 'gamma=0.3,v' with lines lw 2
EOT
gnuplot gpl.txt