URI: 
       tsimulation.c: rename timestep variable in calculation function - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
  HTML git clone git://src.adamsgaard.dk/cngf-pf
   DIR Log
   DIR Files
   DIR Refs
   DIR README
   DIR LICENSE
       ---
   DIR commit 838d87537c8aae08d8ce35f4ec6db0dfddfb67f7
   DIR parent 3e2725c9388be5990e38846b37e38322d66a0e17
  HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
       Date:   Tue, 16 Feb 2021 15:52:14 +0100
       
       simulation.c: rename timestep variable in calculation function
       
       Diffstat:
         M simulation.c                        |      12 ++++++------
       
       1 file changed, 6 insertions(+), 6 deletions(-)
       ---
   DIR diff --git a/simulation.c b/simulation.c
       t@@ -908,7 +908,7 @@ find_flux(const struct simulation *sim)
        void
        set_coupled_fluid_transient_timestep(struct simulation *sim, const double safety)
        {
       -        double max_gamma_dot, mu, xi, max_I, t_val;
       +        double max_gamma_dot, mu, xi, max_I, dt;
        
                /* max expected strain rate */
                max_gamma_dot = 1.0/sim->d;
       t@@ -936,10 +936,10 @@ set_coupled_fluid_transient_timestep(struct simulation *sim, const double safety
                                                           sim->sigma_n[sim->nz - 1] - sim->p_f_mod_ampl,
                                                           sim->rho_s);
        
       -        t_val = xi * (sim->alpha + sim->phi[sim->nz - 1] * sim->beta_f) * sim->mu_f
       -                / (sim->phi[sim->nz - 1] * sim->phi[sim->nz - 1]
       -                   * sim->phi[sim->nz - 1] * sim->L_z * max_I);
       +        dt = xi * (sim->alpha + sim->phi[sim->nz - 1] * sim->beta_f) * sim->mu_f
       +             / (sim->phi[sim->nz - 1] * sim->phi[sim->nz - 1]
       +                * sim->phi[sim->nz - 1] * sim->L_z * max_I);
        
       -        if (sim->dt > safety * t_val)
       -                sim->dt = safety * t_val;
       +        if (sim->dt > safety * dt)
       +                sim->dt = safety * dt;
        }