[cig-commits] [commit] devel, master: Reduce DT for AK135F_NO_MUD. (f12a92e)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d_globe

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

>---------------------------------------------------------------

commit f12a92e861386ef76af0e1deaed7c8ab8f754e28
Author: Elliott Sales de Andrade <esalesde at physics.utoronto.ca>
Date:   Sun Jun 1 20:02:32 2014 -0400

    Reduce DT for AK135F_NO_MUD.
    
    The change to the "new mud" version modified the model enough to cause
    the simulation to become unstable in the fluid. We need to reduce the
    time step to remain stable.


>---------------------------------------------------------------

f12a92e861386ef76af0e1deaed7c8ab8f754e28
 src/shared/get_timestep_and_layers.f90 | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/shared/get_timestep_and_layers.f90 b/src/shared/get_timestep_and_layers.f90
index de00adc..600cbfc 100644
--- a/src/shared/get_timestep_and_layers.f90
+++ b/src/shared/get_timestep_and_layers.f90
@@ -437,6 +437,10 @@
   if (REFERENCE_1D_MODEL == REFERENCE_MODEL_1066A) &
     DT = DT*(1.d0 - 0.8d0)  ! *0.20d0
 
+  ! reduces time step size for "no mud" version of AK135F model
+  if (REFERENCE_1D_MODEL == REFERENCE_MODEL_AK135F_NO_MUD) &
+    DT = DT*(1.d0 - 0.05d0)
+
   ! reduces time step size for crustmaps crustal model
   if( ITYPE_CRUSTAL_MODEL == ICRUST_CRUSTMAPS ) &
     DT = DT*(1.d0 - 0.3d0)