[cig-commits] r20010 - seismo/1D/SPECFEM1D/trunk

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-04-30 19:19:33 -0700 (Mon, 30 Apr 2012)
New Revision: 20010

Modified:
   seismo/1D/SPECFEM1D/trunk/constants.h
   seismo/1D/SPECFEM1D/trunk/wave.f90
Log:
increased the CFL number


Modified: seismo/1D/SPECFEM1D/trunk/constants.h
===================================================================
--- seismo/1D/SPECFEM1D/trunk/constants.h	2012-05-01 01:58:26 UTC (rev 20009)
+++ seismo/1D/SPECFEM1D/trunk/constants.h	2012-05-01 02:19:33 UTC (rev 20010)
@@ -32,7 +32,7 @@
 
 ! Courant–Friedrichs–Lewy (CFL) stability value
 ! see e.g. http://en.wikipedia.org/wiki/Courant%E2%80%93Friedrichs%E2%80%93Lewy_condition
-  double precision, parameter :: courant_CFL = 0.20d0
+  double precision, parameter :: courant_CFL = 0.45d0
 
 ! number of time steps to compute
   integer, parameter :: NSTEP = 20000

Modified: seismo/1D/SPECFEM1D/trunk/wave.f90
===================================================================
--- seismo/1D/SPECFEM1D/trunk/wave.f90	2012-05-01 01:58:26 UTC (rev 20009)
+++ seismo/1D/SPECFEM1D/trunk/wave.f90	2012-05-01 02:19:33 UTC (rev 20010)
@@ -188,6 +188,7 @@
     enddo
   endif
 
+! main time loop
   do it = 1,NSTEP
 
 ! `predictor' update displacement using explicit finite-difference time scheme (Newmark)