1d_fd_simple_shear_transient

commit 4eaf14a0a06f01b3c5a2b9a49f30c9a23620e2a0
parent 63a7aea187ad7c0c99223e065f93d93df89833e7
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri, 10 May 2019 15:59:17 -0700

Fix axes units for stick-slip plot

Diffstat:
M
Makefile
 | 
3
++-
M
diurnal.timeseries.gp
 | 
6
+++---

2 files changed, 5 insertions(+), 4 deletions(-)
diff --git a/Makefile b/Makefile
@@ -29,9 +29,10 @@ diurnal.output00000.txt: 1d_fd_simple_shear
 	./$< --resolution 50 --length 2.0 --normal-stress 150e3 --fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 --fluid-pressure-ampl 50e3 --fluid-pressure-freq $$(( 1.0/(3600*24) )) --time-step 1e-1 --file-interval $$(( 60*10 )) --time-end $$(( 3600*24*7 )) diurnal'
 
 diurnal.timeseries.txt: diurnal.output00000.txt
+	/bin/bash -c '\
 	for f in diurnal.output*.txt; do \
 		tail -n 1 "$$f" | cut -f2- >> $@; \
-	done
+	done'
 
 diurnal.timeseries.png: diurnal.timeseries.gp diurnal.timeseries.txt
 	gnuplot $< > $@
diff --git a/diurnal.timeseries.gp b/diurnal.timeseries.gp
@@ -1,11 +1,11 @@
 #!/usr/bin/env gnuplot
 
 set terminal pngcairo color size 18.6 cm, 10.0 cm
-set xlabel "Time, t [s]"
-set ylabel "Shear velocity, v [m/s]" offset 2
+set xlabel "Time, t [d]"
+set ylabel "Shear velocity, v [km/d]" offset 2
 
 # set xrange [0.0:0.7]
 # set yrange [0.0:0.01]
 # set key top left #samplen 0.9
 
-plot "diurnal.timeseries.txt" u 1 w l lw 2 t ""
+plot "diurnal.timeseries.txt" u ($0/1008*7):($1*60*60*24/1000) w l lw 2 t ""