Courant condition implemented for timestep, as described in documentation.

This hasn't been very well tested yet.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Time/trunk@11 5633253d-7678-4964-a54d-f87795f8ee59

1 files changed, 20 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index e172a44..317be03 100644
--- a/param.ccl
+++ b/param.ccl
@@ -8,3 +8,23 @@ REAL dtfac "The standard timestep condition dt = dtfac*max(delta_space)"
  0:* :: "Probably only makes sense to be bigger than zero"
 } 0.5 
 
+REAL courant_fac "The courant timestep condition dt = courant_fac*max(delta_space)/speed/sqrt(dim)"
+{
+ 0:* :: "Probably only makes sense to be bigger than zero"
+} 0.9
+
+KEYWORD courant_method "Method for calculating timestep"
+{
+ "standard" :: "Courant condition at BASEGRID"
+ "courant"  :: "Courant condition at PRESTEP (using wavespeed)"
+ "courant_time" :: "Courant condition at PRESTEP (using min time)"
+} "standard"
+
+BOOLEAN courant_outonly "Only output courant timestep?"
+{
+} "no"
+
+INT outcourant_every "How often to output courant timestep"
+{
+ 0:* :: "Zero means no output"
+} 0
\ No newline at end of file