hansen-zoet-exp

commit fc4537d39342d4a6fe2f79cc6b1af38dd05be322
parent ad48bb6f87b1b559551224d908ca97708416b9a6
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue,  8 Nov 2022 11:26:22 +0100

fig-strain_distribution.gp: tweak scaling and appearance

Diffstat:
M
fig-strain_distribution.gp
 | 
29
+++++++++++++++--------------

1 file changed, 15 insertions(+), 14 deletions(-)
diff --git a/fig-strain_distribution.gp b/fig-strain_distribution.gp
@@ -11,16 +11,16 @@ set rmargin 2.0
 set tmargin 2.0
 
 set xrange [-5 : 105]
-set yrange [0.09 : 0.11]
+set yrange [0.02 : 0]
 
 #set xlabel "Normalized horizontal velocity, {/:Italic v}_x [-]"
-set xlabel "Horizontal velocity, {/:Italic v}_x [m/a]"
-set ylabel "Vertical position, {/:Italic z} [m]"
+set xlabel "Horizontal velocity, {/:Italic v} [m/a]"
+set ylabel "Depth beneath ice-bed interface, {/:Italic z} [m]"
 set key bottom right font ",10" #samplen 0.9
 set title "{/:Italic v} = 100 m/a"
-plot "strain_distribution_N8kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/:Italic N}' = 8 kPa", \
-     "strain_distribution_N22kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/:Italic N}' = 22 kPa", \
-     "strain_distribution_N50kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/:Italic N}' = 50 kPa"
+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", \
+     "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", \
+     "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"
 
 #unset xrange
 
@@ -28,19 +28,20 @@ plot "strain_distribution_N8kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 t
 #set xlabel "Shear-strain rate, dv_x/dz [1/d]"
 #set ylabel "Vertical position, z [m]"
 #set key bottom right font ",10" #samplen 0.9
-#plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
-     #"strain_distribution_N22kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
-     #"strain_distribution_N50kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"
+#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", \
+     #"strain_distribution_N22kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
+     #"strain_distribution_N50kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 50 kPa"
 
 
-set xrange [-5 : 305]
+set xrange [-15 : 315]
 #set yrange [0 : 0.11]
 
 #set xlabel "Normalized horizontal velocity, {/:Italic v}_x [-]"
-set xlabel "Horizontal velocity, {/:Italic v}_x [m/a]"
+#set xlabel "Depth, {/:Italic v} [m/a]"
 unset ylabel
+set xtics 0,100,300
 set key bottom right font ",10" #samplen 0.9
 set title "{/:Italic v} = 300 m/a"
-plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/:Italic N}' = 8 kPa", \
-     "strain_distribution_N22kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/:Italic N}' = 22 kPa", \
-     "strain_distribution_N50kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/:Italic N}' = 50 kPa"
+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", \
+     "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", \
+     "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"