cngf-pf

commit 63448eb2803256024a9e887c9f29aae7bfc35946
parent 9f173a7e240baa264bf6f3a67fcc43348ac8674a
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu,  9 Sep 2021 16:29:32 +0200

fix stress naming in fluidtest and plot effective normal stress

Diffstat:
M
fluidtest/Makefile
 | 
4
++--
M
fluidtest/fig.gp
 | 
4
++--

2 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/fluidtest/Makefile b/fluidtest/Makefile
@@ -7,7 +7,7 @@ FLUID_PRESSURE_TOP = 100e3
 #FLUID_PRESSURE_AMPL = 120e3
 FLUID_PRESSURE_AMPL = 80e3
 FLUID_PRESSURE_FREQ != echo "1.0/(3600*24)" | bc -l
-EFF_NORMAL_STRESS_TOP = 200e3
+NORMAL_STRESS_TOP = 200e3
 
 timeseries.pdf: fig.gp timeseries.txt
 	sh -c '\
@@ -34,7 +34,7 @@ timeseries.txt: sim.output00000.txt
 sim.output00000.txt: $(BIN)
 	/bin/sh -c '\
 	$(BIN) -L 8.0 \
-	-n $(EFF_NORMAL_STRESS_TOP) \
+	-n $(NORMAL_STRESS_TOP) \
 	-f 0.40 \
 	-F \
 	-k $(PERMEABILITY) \
diff --git a/fluidtest/fig.gp b/fluidtest/fig.gp
@@ -18,7 +18,7 @@ tend = 72*inc - 1
 filename(t) = sprintf('sim.output%05d.txt', t)
 #set origin 0.095, 0.09
 #set size 0.14, 0.23
-set xlabel "{/:Normal Water pressure [kPa]}"
+set xlabel "{/:Normal Effective normal stress [kPa]}"
 set ylabel ""
 unset arrow
 xmin = 0
@@ -27,7 +27,7 @@ set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
 set xrange [xmin:xmax]
 set xtics 100 offset 0,0.5
 set ytics 2 offset 0.5
-plot for [t = t0:tend:inc*2] filename(t) u ($4/1e3):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t ""
+plot for [t = t0:tend:inc*2] filename(t) u ($3/1e3):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t ""
 
 reset
 set xlabel "{/:Normal Time [d]}"