fig-strain_distribution.gp (2185B)
1 #!/usr/bin/env gnuplot 2 3 reset 4 5 set terminal pdfcairo enhanced color size 8.3 cm, 8.3 cm font ",10" 6 set multiplot layout 1,2 7 8 set lmargin 10.0 9 set bmargin 3.5 10 set rmargin 2.0 11 set tmargin 2.0 12 13 set xrange [-5 : 105] 14 set yrange [0.02 : 0] 15 16 labely = 0.98 17 set label "c" at screen 0.08,labely font "{:Bold},12" 18 set label "d" at screen 0.58,labely font "{:Bold},12" 19 20 #set xlabel "Normalized horizontal velocity, {/:Italic v}_x [-]" 21 set xlabel "Horizontal velocity, {/:Italic v} [m/a]" 22 set ylabel "Depth beneath ice-bed interface, {/:Italic z} [m]" 23 set key bottom right font ",10" #samplen 0.9 24 set title "{/:Italic v} = 100 m/a" 25 plot "strain_distribution_N8kPa_v3.169e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 8 kPa", \ 26 "strain_distribution_N22kPa_v3.169e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 22 kPa", \ 27 "strain_distribution_N50kPa_v3.169e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 50 kPa" 28 29 #unset xrange 30 31 ##set xlabel "Normalized horizontal velocity, v_x [-]" 32 #set xlabel "Shear-strain rate, dv_x/dz [1/d]" 33 #set ylabel "Vertical position, z [m]" 34 #set key bottom right font ",10" #samplen 0.9 35 #plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 8 kPa", \ 36 #"strain_distribution_N22kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 22 kPa", \ 37 #"strain_distribution_N50kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 50 kPa" 38 39 40 set xrange [-15 : 315] 41 #set yrange [0 : 0.11] 42 43 #set xlabel "Normalized horizontal velocity, {/:Italic v}_x [-]" 44 #set xlabel "Depth, {/:Italic v} [m/a]" 45 unset ylabel 46 set xtics 0,100,300 47 set key bottom right font ",10" #samplen 0.9 48 set title "{/:Italic v} = 300 m/a" 49 plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 8 kPa", \ 50 "strain_distribution_N22kPa_v9.506e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 22 kPa", \ 51 "strain_distribution_N50kPa_v9.506e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 50 kPa"