1d_fd_simple_shear_transient

commit 802c24eb1aed07014ed6c8f0ba825298b24a8325
parent df8b703a4833a956612cd7e1710087e7c2d844ba
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon,  1 Jul 2019 15:26:48 +0200

Use precomputed cell spacing

Diffstat:
M
fluid.c
 | 
4
++--
M
simulation.c
 | 
2
+-

2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/fluid.c b/fluid.c
@@ -69,7 +69,7 @@ set_fluid_bcs(struct simulation* sim, const double p_f_top)
 	               sim->nz,
 	               -1,
 	               sim->phi[0]*sim->rho_f*sim->G,
-	               sim->z[1] - sim->z[0]);
+	               sim->dz);
 }
 
 static double
@@ -125,7 +125,7 @@ darcy_solver_1d(struct simulation* sim,
 	 *     epsilon = 1.0: implicit */
 	/* epsilon = 0.5; */
 	/* epsilon = 0.0; */
-	epsilon = 0.0;
+	epsilon = 1.0;
 
 	/* choose relaxation factor, parameter in ]0.0; 1.0]
 	 *     theta in ]0.0; 1.0]: underrelaxation
diff --git a/simulation.c b/simulation.c
@@ -397,7 +397,7 @@ implicit_1d_jacobian_poisson_solver(struct simulation* sim,
 		set_bc_dirichlet(sim->g_ghost, sim->nz, +1, 0.0);
 
 		/* Neumann BCs resemble free surfaces */
-		/* set_bc_neumann(sim->g_ghost, sim->nz, +1, 0.0, 0.0); */
+		/* set_bc_neumann(sim->g_ghost, sim->nz, +1, 0.0, sim->dz); */
 
 		for (i=0; i<sim->nz; ++i)
 			poisson_solver_1d_cell_update(i,