fig.gp (3184B)
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 D_min = 1e-9 15 D_max = 1e-3 16 17 # forcing parameters 18 f_min = 0.5*1.0/(24*365)/36000.0 19 f_max = 2.0/3600.0 20 21 # f: frequency, D: diffusivity 22 skindepth(f,D) = (D/(3.141592654*f))**(0.5) # f: 1/s 23 24 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]" 25 set ylabel "Diffusivity {/:Italic D} [m^2/s]" 26 set cblabel "Skin depth {/:Italic d}_s [m]" 27 28 set logscale xyzcb 29 set xrange [f_min:f_max] 30 set yrange [D_min:D_max] 31 set samples 100 32 set isosamples 60 33 34 # Generate contours 35 set view map 36 unset surface 37 set xtics rotate by 90 38 set ytics offset 1 39 set format x '10^{%T}' 40 set format y '10^{%T}' 41 42 set pm3d 43 44 label_y = 3e-9 45 x=1.0/3600.0 46 set arrow from x,D_min to x,D_max nohead lc "white" front 47 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 48 x=1.0/3600.0/24.0 49 set arrow from x,D_min to x,D_max nohead lc "white" front 50 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 51 x=1.0/3600.0/24.0/30.0 52 set arrow from x,D_min to x,D_max nohead lc "white" front 53 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 54 x=1.0/3600.0/24.0/365.0 55 set arrow from x,D_min to x,D_max nohead lc "white" front 56 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 57 58 splot skindepth(x,y) notitle with lines palette, \ 59 skindepth(x,y) with pm3d notitle 60 61 reset 62 set pm3d 63 64 x=1.0/3600.0 65 set arrow from x,D_min to x,D_max nohead lc "white" front 66 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 67 x=1.0/3600.0/24.0 68 set arrow from x,D_min to x,D_max nohead lc "white" front 69 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 70 x=1.0/3600.0/24.0/30.0 71 set arrow from x,D_min to x,D_max nohead lc "white" front 72 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 73 x=1.0/3600.0/24.0/365.0 74 set arrow from x,D_min to x,D_max nohead lc "white" front 75 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white" 76 77 set label "a" at screen 0.01,0.95 font ",12" 78 set label "b" at screen 0.35,0.95 font ",12" 79 set label "c" at screen 0.70,0.95 font ",12" 80 81 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]" 82 set ylabel "Diffusivity {/:Italic D} [m^2/s]" 83 set cblabel "Max. deformation depth [m]" 84 85 set xrange [f_min:f_max] 86 set yrange [D_min:D_max] 87 set xtics rotate by 90 88 set ytics offset 1 89 set format x '10^{%T}' 90 set format y '10^{%T}' 91 92 set view map scale 1 93 set style data pm3d 94 set style function pm3d 95 set xyplane relative 0 96 set nomcbtics 97 set pm3d implicit at b 98 set pm3d scansforward 99 #set pm3d interpolate 10,10 100 set logscale xy 101 102 #unset colorbox 103 104 set title "{/:Italic A}_f = 10 kPa" 105 splot 'max_depth_ampl10e3.txt' notitle with lines palette, \ 106 '' with pm3d notitle 107 108 #unset ylabel 109 #set colorbox 110 set title "{/:Italic A}_f = 100 kPa" 111 splot 'max_depth_ampl100e3.txt' notitle with lines palette, \ 112 '' with pm3d notitle