cngf-pf-exp1

commit 521be0e860745b50d4eca501edf8be6355ca99f9
parent ea95c9535686a1c48d869016605dac80a5725e60
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed, 20 Nov 2019 13:58:39 +0100

Calculate time series of flux

Diffstat:
M
fig-stick_slip_rate/Makefile
 | 
29
+++++++++++++++++++++++++----
M
fig-stick_slip_rate/fig.gp
 | 
29
+++++++++++++++--------------

2 files changed, 40 insertions(+), 18 deletions(-)
diff --git a/fig-stick_slip_rate/Makefile b/fig-stick_slip_rate/Makefile
@@ -1,5 +1,5 @@
 BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
-FIG = $(shell basename $(CURDIR))
+FIG != basename $(PWD)
 
 default: ../$(FIG).pdf
 
@@ -7,7 +7,7 @@ default: ../$(FIG).pdf
 
 sim.output00000.txt: $(BIN)
 	/bin/sh -c '\
-	./$< --length 8.0 \
+	./$(BIN) --length 8.0 \
 	--normal-stress 200e3 \
 	--set-shear-velocity $$( echo "1000/(365*24*60*60)" | bc -l ) \
 	--fluid \
@@ -44,8 +44,29 @@ max_strainrate_depth.txt: sim.output00000.txt
 			"$$f" >> $@; \
 	done'
 
-../$(FIG).pdf: fig.gp timeseries.txt max_strainrate_depth.txt
-	gnuplot $< > $@
+# central finite differences
+sediment_flux.txt: sim.output00000.txt
+	/bin/sh -c '\
+	rm -f $@; \
+	for f in sim.output*.txt; do \
+		awk " \
+			BEGIN{getline; integral=0.0} \
+			NF{ \
+				if (NR > 1) \
+				{ \
+					integral+=(\$$2+v_prev)/2.0*(\$$1-z_prev); \
+				} \
+				z_prev=\$$1; \
+				v_prev=\$$2; \
+			} \
+			END{ \
+				print integral \
+			}" \
+			"$$f" >> $@; \
+	done'
+
+../$(FIG).pdf: fig.gp timeseries.txt max_strainrate_depth.txt sediment_flux.txt
+	gnuplot fig.gp > $@
 
 clean:
 	$(RM) *.txt
diff --git a/fig-stick_slip_rate/fig.gp b/fig-stick_slip_rate/fig.gp
@@ -1,13 +1,13 @@
 #!/usr/bin/env gnuplot
 reset
 
-set terminal pdfcairo color size 7.5 cm, 8.5 cm
+set terminal pdfcairo color size 7.5 cm, 10 cm
 
-set multiplot layout 3,1 \
-	margins 0.19,0.80,0.23,0.95 \
-	spacing 0.03,0.03
+set multiplot layout 4,1 \
+	margins 0.19,0.72,0.10,0.98 \
+	spacing 0.00,0.00
 
-set key top right samplen 1.0
+set key top right samplen 1.0 font ",10"
 unset xlabel
 unset xtics
 set yrange [0:320]
@@ -20,25 +20,26 @@ plot \
 unset xtics
 unset ylabel
 unset ytics
-#set ytics 0.1
-#set yrange [0:1]
-#set y2range [0:1]
-#set y2tics 0.2
-#set y2label "Effective friction [-]"
-#plot "" u ($0/1008*7):($4) w l lw 2 lt 7 t ""
 set yrange [0:80]
 set y2range [0:80]
 set y2tics 20
 set y2label "Shear stress [kPa]"
 plot "" u ($0/1008*7):($4*$2/1000) w l lw 2 lt 1 t ""
 
-set xtics
 unset y2label
 unset y2tics
 set ytics 0.5
 set yrange [2.5:0.0]
-set xlabel "Time [d]"
 set ylabel "Depth of max. strain rate [m]"
 plot "max_strainrate_depth.txt" u ($0/1008*7):(8.0 - $1) w l lw 2 lt 4 t ""
 
-unset multiplot
+unset ylabel
+unset ytics
+set yrange [0.0:1e-4]
+set y2range [0.0:1e-4]
+set y2tics 3e-5
+set y2tics
+set y2label "Specific sediment flux [m/s]"
+set xtics
+set xlabel "Time [d]"
+plot "sediment_flux.txt" u ($0/1008*7):1 w l lw 2 lt 6 t ""