tchange layout of skin-depth figure - hansen-zoet-exp - comparison of lab and model results of sediment advection
 (HTM) git clone git://src.adamsgaard.dk/hansen-zoet-exp
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) Submodules
       ---
 (DIR) commit b52d935fa9363f2b6240a7b4c8792d572225e234
 (DIR) parent 8fb7e560d0ce56a0c12005fa7a2317230c22d5ee
 (HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
       Date:   Thu, 16 Jun 2022 13:44:18 +0200
       
       change layout of skin-depth figure
       
       Diffstat:
         M fig-skin_depth_diff.gp              |      71 +++++--------------------------
       
       1 file changed, 10 insertions(+), 61 deletions(-)
       ---
 (DIR) diff --git a/fig-skin_depth_diff.gp b/fig-skin_depth_diff.gp
       t@@ -1,66 +1,15 @@
        #!/usr/bin/env gnuplot
        
       -reset
       +set terminal pdfcairo color size 8.3 cm, 17 cm font ",10"
        
       -set terminal pdfcairo color size 12.0 cm, 4.0 cm font ",10"
       +set multiplot layout 3,1
        
       -set multiplot layout 1,3 \
       -        margin 0.078,0.90,0.30,0.85 \
       -        spacing 0.21,0.0
       -
       -set key bottom right #samplen 0.9
       -
       -# hydraulic parameters
       +f_min = 1e-9
       +f_max = 1e-3
        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
       -
       +label_y = D_min
        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"
       t@@ -76,12 +25,12 @@ set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcol
        
        set border 31 lc rgbcolor "#666666"
        
       -set label "a" at screen 0.01,0.95 font "{:Bold},12"
       -set label "b" at screen 0.35,0.95 font "{:Bold},12"
       -set label "c" at screen 0.70,0.95 font "{:Bold},12"
       +set label "a" at screen 0.01,0.99 font "{:Bold},12"
       +set label "b" at screen 0.01,0.66 font "{:Bold},12"
       +set label "c" at screen 0.01,0.33 font "{:Bold},12"
        
       -set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
       -set ylabel "Diffusivity {/:Italic D} [m^2/s]"
       +set xlabel "Forcing frequency, {/:Italic f} [s^{-1}]"
       +set ylabel "Diffusivity, {/:Italic D} [m^2/s]"
        set cblabel "Max. deformation depth, {/:Italic z'} [m]"
        
        set xrange [f_min:f_max]