tnormal-stress.gp - seaice-exp-faulting - simulate faulting angles in sea ice under uniaxial compression
 (HTM) git clone git://src.adamsgaard.dk/seaice-exp-faulting
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) README
 (DIR) LICENSE
       ---
       tnormal-stress.gp (710B)
       ---
            1 #!/usr/bin/env gnuplot
            2 
            3 reset
            4 
            5 set terminal pdfcairo color size 17.8 cm, 9.45 cm font ",10"
            6 
            7 set multiplot layout 2,1 #\
            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 set xlabel "Compressive strain {/:Italic e} [-]"
           14 set ylabel "Normal stress {/:Italic N} [kPa]"
           15 
           16 set yrange [0:60]
           17 
           18 file(mu, rotate) = sprintf("faulting%02d-mu_d-%s-%s-seed1-data.tsv", id, mu, rotate)
           19 label(mu, rotate) = sprintf("{/symbol m}_d = %s, %s", mu, rotate)
           20 
           21 plot for [mu in "0.1 0.3 0.5 0.8"] file(mu, "rotate") u ($0 / 4246 * 0.2):($3 / 1000) w l t label(mu, "rotating"),
           22         plot for [mu in "0.1 0.3 0.5 0.8"] file(mu, "norotate") u ($0 / 4246 * 0.2):($3 / 1000) w l t label(mu, "not rotating")