tfluidtest: add time-vs-vel subplot, verify behavior with negative eff. stress - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
 (HTM) git clone git://src.adamsgaard.dk/cngf-pf
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) README
 (DIR) LICENSE
       ---
 (DIR) commit 922cb8962b4ba6ab722a77e976e46d140a45b8fd
 (DIR) parent 015085b37d13aa95b7099c27212a94b5716cd1b7
 (HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
       Date:   Mon, 30 Aug 2021 11:42:33 +0200
       
       fluidtest: add time-vs-vel subplot, verify behavior with negative eff. stress
       
       Diffstat:
         M fluidtest/Makefile                  |       1 +
         M fluidtest/fig.gp                    |       9 +++++++--
       
       2 files changed, 8 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/fluidtest/Makefile b/fluidtest/Makefile
       t@@ -4,6 +4,7 @@ BIN2 = ../max_depth_simple_shear
        PERMEABILITY = 2e-17
        GRAIN_COMPRESSIBILITY = 1e-8
        FLUID_PRESSURE_TOP = 100e3
       +#FLUID_PRESSURE_AMPL = 120e3
        FLUID_PRESSURE_AMPL = 80e3
        FLUID_PRESSURE_FREQ != echo "1.0/(3600*24)" | bc -l
        EFF_NORMAL_STRESS_TOP = 200e3
 (DIR) diff --git a/fluidtest/fig.gp b/fluidtest/fig.gp
       t@@ -3,7 +3,7 @@
        set terminal pdfcairo color size 17.8 cm, 10.0 cm font ",10"
        
        set colorsequence podo
       -set multiplot
       +set multiplot layout 2,1
        
        file_interval=60*10
        plot_interval=60*60
       t@@ -18,7 +18,7 @@ tend = 72*inc - 1
        filename(t) = sprintf('sim.output%05d.txt', t)
        #set origin 0.095, 0.09
        #set size 0.14, 0.23
       -set xlabel "{/:Normal Water pressure [kPa]}" font ",8" offset 0,1.3
       +set xlabel "{/:Normal Water pressure [kPa]}"
        set ylabel ""
        unset arrow
        xmin = 0
       t@@ -28,3 +28,8 @@ set xrange [xmin:xmax]
        set xtics 100 offset 0,0.5
        set ytics 2 offset 0.5
        plot for [t = t0:tend:inc*2] filename(t) u ($4/1e3):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t ""
       +
       +reset
       +set xlabel "{/:Normal Time [d]}"
       +set ylabel "{/:Normal Shear velocity [m/s]}"
       +plot "timeseries.txt" u ($0*file_interval/(3600*24)):($1) w l lw 1 t ""