1d_fd_simple_shear_transient

commit ae23e8e88585fc5f162b73a586c69358dde92ff3
parent 7b434c7ef2fb006a8e5635ed7eb4092fd367e0e4
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Sat,  6 Apr 2019 12:25:21 +0200

Revert julia implementation to standard configuration

Diffstat:
M
julia/1d_fd_simple_shear.jl
 | 
8
+++-----

1 file changed, 3 insertions(+), 5 deletions(-)
diff --git a/julia/1d_fd_simple_shear.jl b/julia/1d_fd_simple_shear.jl
@@ -14,8 +14,7 @@ G = 9.81
 
 ### Effective normal stress exerted by top wall
 # A larger normal stress deepens the deformational depth
-#P_wall_ = [10, 20, 40, 60, 80, 120] .* 1e3  # normal stress [Pa]
-P_wall_ = [10] .* 1e3  # normal stress [Pa]
+P_wall_ = [10, 20, 40, 60, 80, 120] .* 1e3  # normal stress [Pa]
 
 ### bottom velocity along x [m/s]
 v_x_bot = 0.0
@@ -24,7 +23,7 @@ v_x_bot = 0.0
 μ_wall = 0.40
 
 ### Nodes along z
-nz = 10
+nz = 100
 
 ## Material properties
 
@@ -137,7 +136,7 @@ end
 ## Cartesian grid. The function returns the normalized residual value
 function poisson_solver_1d_iteration(g_in, g_out, r_norm,
                                      μ, p, i, Δz,
-                                     verbose=true)
+                                     verbose=false)
 
     coorp_term = Δz^2.0/(2.0*ξ(μ[i])^2.0)
     g_out[i+1] = 1.0/(1.0 + coorp_term) * (coorp_term*g_local(p[i], μ[i])
@@ -158,7 +157,6 @@ end
 function implicit_1d_jacobian_poisson_solver(g, p, μ, Δz,
                                              rel_tol=1e-5,
                                              max_iter=10_000,
-                                             #max_iter=1,
                                              verbose=false)
 
     # allocate second array of g for Jacobian solution procedure