fig.gp (1175B)
1 #!/usr/bin/env gnuplot 2 3 reset 4 5 #set terminal pdfcairo color size 7.5 cm, 5.0 cm 6 set terminal pdfcairo enhanced color size 7.5 cm, 7.5 cm 7 set multiplot layout 2,1 \ 8 #margins 0.09,0.97,0.09,0.97 \ 9 #spacing 0.05,0.08 10 11 load 'BlueSeq.plt' 12 13 set xrange [0:300] 14 set yrange [0:1] 15 16 set key bottom right font ",10" samplen 2.0 17 18 set xlabel "Shear velocity [m/a]" 19 set ylabel "Vertical position [m]" 20 21 plot \ 22 "d_dz.4.output00000.txt" u ($2*365.25*24*60*60):1 w l ls 1 lw 1 t "nz = 4", \ 23 "d_dz.8.output00000.txt" u ($2*365.25*24*60*60):1 w l ls 2 lw 1 t "nz = 8", \ 24 "d_dz.32.output00000.txt" u ($2*365.25*24*60*60):1 w l ls 3 lw 1 t "nz = 32", \ 25 "d_dz.128.output00000.txt" u ($2*365.25*24*60*60):1 w l ls 4 lw 1 t "nz = 128", \ 26 "d_dz.512.output00000.txt" u ($2*365.25*24*60*60):1 w l ls 5 lw 1 t "nz = 512", \ 27 "d_dz.1024.output00000.txt" u ($2*365.25*24*60*60):1 w l ls 6 lw 1 t "nz = 1024" 28 29 unset xrange 30 unset yrange 31 32 #set yrange [0.3:0.5] 33 34 set ylabel "Bulk friction, {/Symbol m} [-]" 35 36 #set xrange [1.5:6000] 37 set format x "10^{%T}" 38 set xlabel "Grain size to cell size, d/dz [-]" 39 set logscale x 40 plot "d_dz.txt" u (1e-3/(1.0/$1)):2 w lp t "" 41 unset logscale x 42 unset format x