tUpdate 1d_fd_simple_shear and tweak pulse lag plot - cngf-pf-exp1 - experiments for first paper with continuum granular model
 (HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) Submodules
 (DIR) README
 (DIR) LICENSE
       ---
 (DIR) commit 0a7a86a9b0a6dc84d7c32393fd7743437b223850
 (DIR) parent 57efb65995d743a228348e5bb190381cf1309547
 (HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
       Date:   Thu,  5 Dec 2019 19:57:38 +0100
       
       Update 1d_fd_simple_shear and tweak pulse lag plot
       
       Diffstat:
         M fig-pulse_lag/Makefile              |       5 +----
         M fig-pulse_lag/fig.gp                |       8 ++++----
       
       2 files changed, 5 insertions(+), 8 deletions(-)
       ---
 (DIR) diff --git a/fig-pulse_lag/Makefile b/fig-pulse_lag/Makefile
       t@@ -3,7 +3,7 @@ FIG != basename $(PWD)
        
        default: ../$(FIG).pdf
        
       -PERM_SEQ=3.16e-16 1e-16 3.16e-17 1e-17 3.16e-18 1e-18 3.16e-19
       +PERM_SEQ=3.16e-16 1e-16 3.16e-17 1e-17 3.16e-18 1e-18 3.16e-19 1e-19
        
        sim-1e-18.output00000.txt: $(BIN)
                /bin/sh -c 'for perm in $(PERM_SEQ); do \
       t@@ -21,9 +21,6 @@ sim-1e-18.output00000.txt: $(BIN)
                --time 0.0 \
                --time-end $$( echo "20*3600*24 + 60*10" | bc -l ) sim-$$perm; done'
        
       -        #--time $$( echo "($$i-1)*3600*24" | bc -l ) \
       -        #--time-end $$( echo "($$i+1)*3600*24*1.5" | bc -l ) sim-$$perm; done'
       -
        timeseries.txt: sim-1e-18.output00000.txt
                /bin/sh -c '\
                rm -f $@; \
 (DIR) diff --git a/fig-pulse_lag/fig.gp b/fig-pulse_lag/fig.gp
       t@@ -4,11 +4,11 @@ reset
        set terminal pdfcairo color size 15 cm, 10 cm
        
        set multiplot layout 4,1 \
       -        margins 0.15,0.90,0.13,0.95 \
       +        margins 0.15,0.88,0.13,0.95 \
                spacing 0.0,0.03
        
        file_interval_seconds=600.0
       -conv=1.0/(24.0*3600.0/file_interval_seconds)
       +conv=1.0/(24.0*3600.0/file_interval_seconds) * 140.0/138.5
        
        set key top right samplen 1.0
        unset xlabel
       t@@ -22,8 +22,8 @@ plot \
        "" u ($0*conv):($3/1000) w l lw 2 lt 3 dt "-" t "water pressure"
        
        unset yrange
       -set yrange [1e-19:1e-15]
       -set ytics (1e-18,1e-17,1e-16)
       +set yrange [5e-20:1e-15]
       +set ytics (1e-19,1e-18,1e-17,1e-16)
        set ylabel "Permeability [m^2]"
        set logscale y
        plot "max_strainrate_depth.txt" u ($0*conv):3 w l lw 2 lt 5 t ""