1d_fd_simple_shear

continuum model for granular flows with pore-pressure dynamics
git clone git://src.adamsgaard.dk/1d_fd_simple_shear
Log | Files | Refs | README | LICENSE Back to index

commit cee81c32ee56c094ebcdcb0319125cabb060ca1f
parent ffe45348839b6f09bf9b9d651b68159cef108db7
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 16 Apr 2020 15:45:35 +0200

Fix debug reporting in fluid solver

Diffstat:
Mfluid.c | 6+++---
1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/fluid.c b/fluid.c @@ -117,7 +117,7 @@ darcy_pressure_change_1d(const int i, #ifdef DEBUG printf("%d->%d: p=[%g, %g, %g]\tk=[%g, %g, %g]\n", i, i+1, - p_zn, p, p_zp, + p_f_ghost_in[i], p_f_ghost_in[i+1], p_f_ghost_in[i+2], k_zn, k_, k_zp); #endif @@ -130,8 +130,8 @@ darcy_pressure_change_1d(const int i, )/dz; #ifdef DEBUG - printf("phi[%d]=%g\tdiv_k_grad_p[%d]=%g\n", - i, phi[i], i, div_k_grad_p); + printf("darcy_pressure_change [%d]: phi=%g\tdiv_k_grad_p=%g\tphi_dot=%g\n", + i, phi[i], div_k_grad_p, phi_dot[i]); #endif /* TODO: add advective term */