commit 34d7d22e17e3308876f502f8dec409804a86ac39
parent dfcb11fd175b3edc939f7cb5f04418b01dd4d0de
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Wed, 4 Sep 2019 13:14:50 +0200
Remove redundant normalized fig7 now that fig6 is rate controlled
Diffstat:
| D | fig7/Makefile | | | 27 | --------------------------- |
| D | fig7/fig.gp | | | 63 | --------------------------------------------------------------- |
2 files changed, 0 insertions(+), 90 deletions(-)
diff --git a/fig7/Makefile b/fig7/Makefile
@@ -1,27 +0,0 @@
-BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
-FIG = $(shell basename $(CURDIR))
-
-default: ../$(FIG).pdf
-
-sim.output00000.txt: $(BIN) Makefile ../skindepth.gp
- ./$< \
- --length 8.0 \
- --set-shear-velocity $$(echo "1000/(365*24*60*60)" | bc -l) \
- --normal-stress 200e3 \
- --fluid \
- --fluid-permeability 2e-17 \
- --fluid-pressure-top 100e3 \
- --fluid-pressure-ampl 80e3 \
- --fluid-pressure-freq $$(echo "1.0/(3600*24)" | bc -l) \
- --file-interval $$(echo "3600" | bc -l) \
- --time-end $$(echo "3600*24*3 + 3600" | bc -l) \
- --normalize sim
-
-../$(FIG).pdf: fig.gp sim.output00000.txt
- gnuplot $< > $@
-
-clean:
- $(RM) *.txt
- $(RM) ../$(FIG).pdf
-
-.PHONY: default clean
diff --git a/fig7/fig.gp b/fig7/fig.gp
@@ -1,63 +0,0 @@
-#!/usr/bin/env gnuplot
-
-reset
-
-#set terminal pdfcairo color size 7.5 cm, 7.5 cm
-set terminal pdfcairo color size 15.0 cm, 5.0 cm
-set multiplot layout 1,4 \
- margins 0.09,0.97,0.23,0.95 \
- spacing 0.03,0.03
-
-set yrange [4.0:8.0]
-
-set key bottom right #samplen 0.9
-
-t0 = 47
-tend = 70
-
-filename(t) = sprintf('sim.output%05d.txt', t)
-
-set pm3d map
-set palette rgb 10,6,6
-
-load '../skindepth.gp'
-k = 2e-17
-f = 1.0/(3600.0*24.0)
-f(x) = 8.0 - skindepth(f,k)
-
-unset colorbox
-set xlabel "Water pressure [kPa]"
-set ylabel "Vertical position [m]"
-set xrange [0:100]
-plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
-
-set xlabel "Eff. normal stress [kPa]"
-set ylabel ""
-set xrange [0:200]
-unset ytics
-plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
-
-set xlabel "Norm. shear vel. [-]"
-set ylabel ""
-set xrange [0.0:1.0]
-#set xtics (0,0.01,0.02)
-unset ytics
-plot for [t = t0:tend] filename(t) u 2:1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
-
-set xlabel "Shear strain rate [1/s]"
-set ylabel ""
-#set xrange [0:3.1e-3]
-set xrange [0.0:0.15]
-set xtics (0,0.05,0.10,0.15)
-#set format x "%.0s*10^{%T}"
-unset ytics
-prev1 = prev2 = 0
-shift(x) = (prev2 = prev1, prev1 = x)
-plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
-unset format x
-
-unset multiplot
