From 456771e0af161604d8d0be949253f9e0fd7206be Mon Sep 17 00:00:00 2001
From: PYamato <philippe.yamato@gmail.com>
Date: Fri, 25 Aug 2023 16:58:05 +0200
Subject: [PATCH] update presented script

---
 examples/Wave1D_Sismo_Day4.jl | 15 ++++-----------
 1 file changed, 4 insertions(+), 11 deletions(-)

diff --git a/examples/Wave1D_Sismo_Day4.jl b/examples/Wave1D_Sismo_Day4.jl
index cf60548..0113ce7 100644
--- a/examples/Wave1D_Sismo_Day4.jl
+++ b/examples/Wave1D_Sismo_Day4.jl
@@ -77,7 +77,8 @@ function MainSource()
         @. ε̇.xx = ∂Vx∂x - 1/3*∇V
       
         # Stress update
-        @. τ.xx = f_shear(G)*Δt*(ε̇.xx) + f_relax(G)*τ.xx
+        #@. τ.xx = f_shear(G)*Δt*(ε̇.xx) + f_relax(G)*τ.xx
+        @. τ.xx = f_shear(G,1.0e5,Δt)*Δt*(ε̇.xx) + f_relax(G,1.0e5,Δt)*τ.xx
 
         # Pressure update 
         @. P    = P - Δt*f_bulk(K)*∇V
@@ -111,18 +112,10 @@ function MainSource()
     title="Receiver gather", xlabel="distance to the source [m]",
     ylabel="time [s]",yflip=true,clim=(-valim,+valim))
     display(p)
+    #png(p,"./docs/src/assets/traceVE.png")
+    
 
 end
 
-function f_bulk(K) 
-   return K
-end
-
-function f_shear(G)
-    return 2*G
-end
-function f_relax(G)
-    return 1.
-end
 
 MainSource()
\ No newline at end of file