commit 0cc65a5920ab86a1b41db896fdc3346a03587a32
parent 435973a434cc28312a0da7fff921326218071ebd
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 10 Jun 2022 10:01:55 +0200
add skin depth figure with diffusivity
Diffstat:
4 files changed, 192 insertions(+), 0 deletions(-)
diff --git a/fig-skin_depth_diff/BlueSeq.plt b/fig-skin_depth_diff/BlueSeq.plt
@@ -0,0 +1,26 @@
+# line styles for ColorBrewer Blues
+# for use with sequential data
+# provides 8 blue colors of increasing saturation
+# compatible with gnuplot >=4.2
+# derived from Anna Schneider
+
+# line styles
+
+set style line 1 lc rgb '#084594' # dark blue
+set style line 2 lc rgb '#2171B5' #
+set style line 3 lc rgb '#4292C6' # medium blue
+set style line 4 lc rgb '#6BAED6' #
+set style line 5 lc rgb '#9ECAE1' # light blue
+set style line 6 lc rgb '#C6DBEF' #
+set style line 7 lc rgb '#DEEBF7' #
+set style line 8 lc rgb '#F7FBFF' # very light blue
+
+# palette
+set palette defined ( 0 '#084594',\
+ 1 '#2171B5',\
+ 2 '#4292C6',\
+ 3 '#6BAED6',\
+ 4 '#9ECAE1',\
+ 5 '#C6DBEF',\
+ 6 '#DEEBF7',\
+ 7 '#F7FBFF')
diff --git a/fig-skin_depth_diff/Makefile b/fig-skin_depth_diff/Makefile
@@ -0,0 +1,37 @@
+FIG != basename $(PWD)
+BIN = ../cngf-pf/max_depth_simple_shear
+
+default: ../$(FIG).pdf
+
+../$(FIG).pdf: fig.gp max_depth_ampl10e3.txt
+ gnuplot fig.gp > $@
+
+max_depth_ampl10e3.txt: $(BIN) Makefile
+ for A_f in 10e3 100e3; do\
+ out=max_depth_ampl$${A_f}.txt;\
+ rm -f $$out;\
+ for D in $$(range -l -N 46 1e-9 1e-3); 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 $$D >> $$out;\
+ ./$(BIN) \
+ -D $$D \
+ -a $$A_f \
+ -q $$f | \
+ awk '{print $$1}' >> $$out;\
+ done;\
+ echo >> $$out;\
+ done;\
+ done
+
+clean:
+ rm -f ../$(FIG).pdf
+ rm -f max_depth_ampl*.txt
+
+.PHONY: default clean
diff --git a/fig-skin_depth_diff/fig.gp b/fig-skin_depth_diff/fig.gp
@@ -0,0 +1,112 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+set terminal pdfcairo color size 17.8 cm, 4.45 cm font ",10"
+
+set multiplot layout 1,3 \
+ margin 0.07,0.92,0.25,0.86 \
+ spacing 0.17,0.0
+
+set key bottom right #samplen 0.9
+
+# hydraulic parameters
+D_min = 1e-9
+D_max = 1e-3
+
+# forcing parameters
+f_min = 0.5*1.0/(24*365)/36000.0
+f_max = 2.0/3600.0
+
+# f: frequency, D: diffusivity
+skindepth(f,D) = (D/(3.141592654*f))**(0.5) # f: 1/s
+
+set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
+set ylabel "Diffusivity {/:Italic D} [m^2/s]"
+set cblabel "Skin depth {/:Italic d}_s [m]"
+
+set logscale xyzcb
+set xrange [f_min:f_max]
+set yrange [D_min:D_max]
+set samples 100
+set isosamples 60
+
+# Generate contours
+set view map
+unset surface
+set xtics rotate by 90
+set ytics offset 1
+set format x '10^{%T}'
+set format y '10^{%T}'
+
+set pm3d
+
+label_y = 3e-9
+x=1.0/3600.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+x=1.0/3600.0/24.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+x=1.0/3600.0/24.0/30.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+x=1.0/3600.0/24.0/365.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+
+splot skindepth(x,y) notitle with lines palette, \
+ skindepth(x,y) with pm3d notitle
+
+reset
+set pm3d
+
+x=1.0/3600.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+x=1.0/3600.0/24.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+x=1.0/3600.0/24.0/30.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+x=1.0/3600.0/24.0/365.0
+set arrow from x,D_min to x,D_max nohead lc "white" front
+set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
+
+set label "a" at screen 0.01,0.95 font ",12"
+set label "b" at screen 0.35,0.95 font ",12"
+set label "c" at screen 0.70,0.95 font ",12"
+
+set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
+set ylabel "Diffusivity {/:Italic D} [m^2/s]"
+set cblabel "Max. deformation depth [m]"
+
+set xrange [f_min:f_max]
+set yrange [D_min:D_max]
+set xtics rotate by 90
+set ytics offset 1
+set format x '10^{%T}'
+set format y '10^{%T}'
+
+set view map scale 1
+set style data pm3d
+set style function pm3d
+set xyplane relative 0
+set nomcbtics
+set pm3d implicit at b
+set pm3d scansforward
+#set pm3d interpolate 10,10
+set logscale xy
+
+#unset colorbox
+
+set title "{/:Italic A}_f = 10 kPa"
+splot 'max_depth_ampl10e3.txt' notitle with lines palette, \
+ '' with pm3d notitle
+
+#unset ylabel
+#set colorbox
+set title "{/:Italic A}_f = 100 kPa"
+splot 'max_depth_ampl100e3.txt' notitle with lines palette, \
+ '' with pm3d notitle
diff --git a/fig-skin_depth_diff/triangle.gp b/fig-skin_depth_diff/triangle.gp
@@ -0,0 +1,17 @@
+
+
+set terminal pdfcairo color
+
+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'
