cngf-pf-exp1

commit 701ec99ae9cf8d6de043c34c2a7e3abb99577f36
parent 7640b6b7f669a0d19a397a1e7ffdb38c7b42822b
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 12 Dec 2019 16:53:24 +0100

First working visualization of max. def. depth (still ugly, though)

Diffstat:
M
fig-skin_depth/Makefile
 | 
25
++++++++++++++++++++-----
M
fig-skin_depth/fig.gp
 | 
33
+++++++++++++++++++++++++++------

2 files changed, 47 insertions(+), 11 deletions(-)
diff --git a/fig-skin_depth/Makefile b/fig-skin_depth/Makefile
@@ -3,15 +3,30 @@ BIN = ../1d_fd_simple_shear/max_depth_simple_shear
 
 default: ../$(FIG).pdf
 
-../$(FIG).pdf: fig.gp max_depth_ampl10e3.txt
+../$(FIG).pdf: fig.gp max_depth_ampl10e3.txt #triangle.gp
+	#gnuplot triangle.gp > $@
 	gnuplot fig.gp > $@
 
-max_depth_ampl10e3.txt: $(BIN)
+max_depth_ampl10e3.txt: $(BIN) Makefile
 	for A_f in 10e3 100e3; do\
 		out=max_depth_ampl$${A_f}.txt;\
 		rm -f $$out;\
-		for f in 1e-9 1e-8 1e-7 1e-6 1e-5 1e-4 1e-3; do\
-			for k in 1e-19 1e-18 1e-17 1e-16 1e-15 1e-14; do\
+		for k in \
+			1e-19 1.29e-19 1.67e-19 2.15e-19 2.78e-19 3.59e-19 4.64e-19 5.99e-19 7.74e-19 \
+			1e-18 1.29e-18 1.67e-18 2.15e-18 2.78e-18 3.59e-18 4.64e-18 5.99e-18 7.74e-18 \
+			1e-17 1.29e-17 1.67e-17 2.15e-17 2.78e-17 3.59e-17 4.64e-17 5.99e-17 7.74e-17 \
+			1e-16 1.29e-16 1.67e-16 2.15e-16 2.78e-16 3.59e-16 4.64e-16 5.99e-16 7.74e-16 \
+			1e-15 1.29e-15 1.67e-15 2.15e-15 2.78e-15 3.59e-15 4.64e-15 5.99e-15 7.74e-15 \
+			1e-14; do\
+			for f in \
+				1e-9 1.29e-9 1.67e-9 2.15e-9 2.78e-9 3.59e-9 4.64e-9 5.99e-9 7.74e-9 \
+				1e-8 1.29e-8 1.67e-8 2.15e-8 2.78e-8 3.59e-8 4.64e-8 5.99e-8 7.74e-8 \
+				1e-7 1.29e-7 1.67e-7 2.15e-7 2.78e-7 3.59e-7 4.64e-7 5.99e-7 7.74e-7 \
+				1e-6 1.29e-6 1.67e-6 2.15e-6 2.78e-6 3.59e-6 4.64e-6 5.99e-6 7.74e-6 \
+				1e-5 1.29e-5 1.67e-5 2.15e-5 2.78e-5 3.59e-5 4.64e-5 5.99e-5 7.74e-5 \
+				1e-4 1.29e-4 1.67e-4 2.15e-4 2.78e-4 3.59e-4 4.64e-4 5.99e-4 7.74e-4 \
+				1e-3; do\
+				printf "%g\t%g\t" $$f $$k >> $$out;\
 				./$(BIN) \
 					--porosity 0.25 \
 					--fluid-compressibility 3.9e-10 \
@@ -21,7 +36,7 @@ max_depth_ampl10e3.txt: $(BIN)
 					--fluid-pressure-top 100e3 \
 					--fluid-pressure-ampl $$A_f \
 					--fluid-pressure-freq $$f | \
-					awk '{print $$1}' | tr '\n' '\t' >> $$out;\
+					awk '{print $$1}' >> $$out;\
 			done;\
 			echo >> $$out;\
 		done;\
diff --git a/fig-skin_depth/fig.gp b/fig-skin_depth/fig.gp
@@ -2,12 +2,17 @@
 
 reset
 
-set terminal pdfcairo color size 7.5 cm, 5.5 cm
+#set terminal pdfcairo color size 15 cm, 5.5 cm
+set terminal pdfcairo color size 20 cm, 10 cm
 
-set lmargin at screen 0.20
-set rmargin at screen 0.75
-set bmargin at screen 0.25
-set tmargin at screen 0.95
+set multiplot layout 1,3 \
+	#margins 0.19,0.72,0.10,0.98 \
+	#spacing 0.00,0.03
+
+#set lmargin at screen 0.20
+#set rmargin at screen 0.75
+#set bmargin at screen 0.25
+#set tmargin at screen 0.95
 
 set key bottom right #samplen 0.9
 
@@ -43,7 +48,6 @@ set contour base
 set cntrlabel format '%8.3g' font ',7' start 0.1 interval -1
 set cntrparam levels discrete 0.1, 0.3, 1, 3, 10, 30, 100, 300
 set key outside title "Skin depth\n[m]" width 1.0 samplen 2.5
-#set xtics (0.0001, 0.01, 1.0)
 set xtics rotate by 90
 set ytics offset 1
 set format x '10^{%T}'
@@ -67,3 +71,20 @@ set arrow from x,k_min to x,k_max nohead lc black
 set label "yearly" at x,label_y rotate by 90 offset screen -0.03,0
 
 splot skindepth(x,y) title ""
+
+reset
+
+set view map scale 1
+set style data pm3d
+set style function pm3d
+set xyplane relative 0
+set nomcbtics
+set xlabel "X LABEL" 
+set ylabel "Y LABEL" 
+set pm3d implicit at b
+set pm3d scansforward
+#set pm3d interpolate 10,10
+set logscale xy
+
+splot 'max_depth_ampl10e3.txt' t 'A_f = 10 kPa'
+splot 'max_depth_ampl100e3.txt' t 'A_f = 100 kPa'