cngf-pf-exp1

experiments for first paper with continuum granular model
git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp # fast
git clone https://src.adamsgaard.dk/manus_continuum_granular1_exp.git # slow
Log | Files | Refs | Submodules | README | LICENSE Back to index

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