set border 3
set noxtics; set noytics; set xtics nomirror; set ytics nomirror
set key below
set xlabel "t (days)"
set ylabel "U_z (mm)"
set title "Displacements along the z axis case #1"
set terminal x11 1
plot [0:500] 0, \
'cas1/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas1/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas1/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis case #2"
set terminal x11 2
plot [0:500] 0, \
'cas2/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas2/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas2/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis case #3"
set terminal x11 3
plot [0:500] 0, \
'cas3/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas3/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas3/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis case #4"
set terminal x11 4
plot [0:500] 0, \
'cas4/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas4/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas4/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis case #5"
set terminal x11 5
plot [0:500] 0, \
'cas5/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas5/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas5/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis case #6"
set terminal x11 6
plot [0:500] 0, \
'cas6/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas6/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas6/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis case #7"
set terminal x11 7
plot [0:500] 0, \
'cas7/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas7/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas7/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis case #8"
set terminal x11 8
plot [0:500] 0, \
'cas8/art4.Uz0' using ($1/86400):($2*1000.) t "at (280,280,0)" with l lw 2, \
'cas8/art4.Uz1' using ($1/86400):($2*1000.) t "at (280,280,5.25)" with l lw 2, \
'cas8/art4.Uz2' using ($1/86400):($2*1000.) t "at (280,280,10.5)" with l lw 2
#
set title "Displacements along the z axis at (280,280,5.25)\n at three density levels"
set terminal x11 11
plot [0:500] 0, \
'cas1/art4.Uz1' using ($1/86400):($2*1000.) t "at 540 kg/m^3" with l lw 2, \
'cas2/art4.Uz1' using ($1/86400):($2*1000.) t "at 650 kg/m^3" with l lw 2, \
'cas3/art4.Uz1' using ($1/86400):($2*1000.) t "at 800 kg/m^3" with l lw 2