fig.gp (3280B)
1 #!/usr/bin/env gnuplot 2 3 reset 4 5 set terminal pdfcairo color size 17.8 cm, 4.45 cm font ",10" 6 7 set multiplot layout 1,3 \ 8 margin 0.07,0.92,0.25,0.86 \ 9 spacing 0.17,0.0 10 11 set key bottom right #samplen 0.9 12 13 # hydraulic parameters 14 phi = 0.25 15 alpha = 1e-8 16 beta_f = 3.9e-10 17 mu_f = 1.787e-3 18 #k = 2e-17 19 k_min = 1e-19 20 k_max = 1e-14 21 22 # forcing parameters 23 f_min = 0.5*1.0/(24*365)/36000.0 24 f_max = 2.0/3600.0 25 26 # f: frequency, k: permeability 27 skindepth(f,k) = (k/((alpha + phi*beta_f)*mu_f*3.141592654*f))**(0.5) # f: 1/s 28 29 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]" 30 set ylabel "Permeability {/:Italic k} [m^2]" 31 set cblabel "Skin depth {/:Italic d}_s [m]" 32 33 set logscale xyzcb 34 set xrange [f_min:f_max] 35 set yrange [k_min:k_max] 36 set samples 100 37 set isosamples 60 38 39 # Generate contours 40 set view map 41 unset surface 42 set xtics rotate by 90 43 set ytics offset 1 44 set format x '10^{%T}' 45 set format y '10^{%T}' 46 47 set pm3d 48 49 label_y = 3e-17 50 x=1.0/3600.0 51 set arrow from x,k_min to x,k_max nohead lc "white" front 52 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 53 x=1.0/3600.0/24.0 54 set arrow from x,k_min to x,k_max nohead lc "white" front 55 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 56 x=1.0/3600.0/24.0/30.0 57 set arrow from x,k_min to x,k_max nohead lc "white" front 58 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 59 x=1.0/3600.0/24.0/365.0 60 set arrow from x,k_min to x,k_max nohead lc "white" front 61 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 62 63 splot skindepth(x,y) notitle with lines palette, \ 64 skindepth(x,y) with pm3d notitle 65 66 reset 67 set pm3d 68 69 x=1.0/3600.0 70 set arrow from x,k_min to x,k_max nohead lc "white" front 71 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 72 x=1.0/3600.0/24.0 73 set arrow from x,k_min to x,k_max nohead lc "white" front 74 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 75 x=1.0/3600.0/24.0/30.0 76 set arrow from x,k_min to x,k_max nohead lc "white" front 77 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 78 x=1.0/3600.0/24.0/365.0 79 set arrow from x,k_min to x,k_max nohead lc "white" front 80 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 81 82 set label "a" at screen 0.01,0.95 font ",12" 83 set label "b" at screen 0.35,0.95 font ",12" 84 set label "c" at screen 0.70,0.95 font ",12" 85 86 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]" 87 set ylabel "Permeability {/:Italic k} [m^2]" 88 set cblabel "Max. deformation depth [m]" 89 90 set xrange [f_min:f_max] 91 set yrange [k_min:k_max] 92 set xtics rotate by 90 93 set ytics offset 1 94 set format x '10^{%T}' 95 set format y '10^{%T}' 96 97 set view map scale 1 98 set style data pm3d 99 set style function pm3d 100 set xyplane relative 0 101 set nomcbtics 102 set pm3d implicit at b 103 set pm3d scansforward 104 #set pm3d interpolate 10,10 105 set logscale xy 106 107 #unset colorbox 108 109 set title "{/:Italic A}_f = 10 kPa" 110 splot 'max_depth_ampl10e3.txt' notitle with lines palette, \ 111 '' with pm3d notitle 112 113 #unset ylabel 114 #set colorbox 115 set title "{/:Italic A}_f = 100 kPa" 116 splot 'max_depth_ampl100e3.txt' notitle with lines palette, \ 117 '' with pm3d notitle