cngf-pf

commit 38e49fd66dea6a98434a41f1c1710caf47555af1
parent fa482eaba45c62b6cfa37afdc35a733a07b58757
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Sat, 11 May 2019 10:42:38 -0700

Produce multiplot for diurnal timeseries

Diffstat:
M
1d_fd_simple_shear_rheology_tulaczyk.gp
 | 
2
+-
M
1d_fd_simple_shear_rheology_tulaczyk.png
 | 
0

M
diurnal.timeseries.gp
 | 
24
++++++++++++++++++++----

3 files changed, 21 insertions(+), 5 deletions(-)
diff --git a/1d_fd_simple_shear_rheology_tulaczyk.gp b/1d_fd_simple_shear_rheology_tulaczyk.gp
@@ -18,7 +18,7 @@ load 'BlueSeq.plt'
 # divide friction by threshold friction value (mu_s)
 # divide shear velocity by 1 m/h
 
-plot "1d_fd_simple_shear_rheology_b0.94_tulaczyk.txt" u ($1/0.502):($2/3.4e-9) w l lc '#666666' lw 2 dashtype '-' t "glass beads, b = 0.94", \
+plot "1d_fd_simple_shear_rheology_b0.94_tulaczyk.txt" u ($1/0.502):($2/3.4e-9) w l lc '#666666' lw 2 dashtype '-' t "b = 0.94", \
          "1d_fd_simple_shear_rheology_b0.40_tulaczyk.txt" u ($1/0.502):($2/8.6e-9) w l ls 6 lw 2 t "b = 0.40", \
          "1d_fd_simple_shear_rheology_b0.20_tulaczyk.txt" u ($1/0.502):($2/1.7e-8) w l ls 5 lw 2 t "b = 0.20", \
          "1d_fd_simple_shear_rheology_b0.10_tulaczyk.txt" u ($1/0.502):($2/3.5e-8) w l ls 3 lw 2 t "b = 0.10", \
diff --git a/1d_fd_simple_shear_rheology_tulaczyk.png b/1d_fd_simple_shear_rheology_tulaczyk.png
Binary files differ.
diff --git a/diurnal.timeseries.gp b/diurnal.timeseries.gp
@@ -1,11 +1,27 @@
 #!/usr/bin/env gnuplot
+reset
 
-set terminal pdfcairo color size 10.0 cm, 5.0 cm
-set xlabel "Time, t [d]"
-set ylabel "Shear velocity, v [km/d]" offset 2
+#set terminal pdfcairo color size 10.0 cm, 5.0 cm
+set terminal pdfcairo color size 10.0 cm, 8.0 cm
 
 # set xrange [0.0:0.7]
 # set yrange [0.0:0.01]
 # set key top left #samplen 0.9
 
-plot "diurnal.timeseries.txt" u ($0/1008*7):($1*60*60*24/1000) w l lw 2 lc "red" t ""
+set multiplot layout 2,1
+
+unset xlabel
+unset xtics
+set yrange [0:210]
+set ylabel "Stress [kPa]" offset 2
+plot \
+"diurnal.timeseries.txt" u ($0/1008*7):($2/1000) w l lw 2 lc "black" t "effective normal stress", \
+"" u ($0/1008*7):($3/1000) w l lw 2 lc "blue" dt "-" t "water pressure"
+
+set xtics
+unset yrange
+set xlabel "Time, t [d]"
+set ylabel "Shear velocity, v [km/d]" offset 2
+plot "" u ($0/1008*7):($1*60*60*24/1000) w l lw 2 lc "red" t ""
+
+unset multiplot