[cig-commits] r6543 - in long/3D/Gale/trunk: . src/StgFEM/SLE/ProvidedSystems/StokesFlow/src

[walter at geodynamics.org]

Author: walter
Date: 2007-04-10 15:46:37 -0700 (Tue, 10 Apr 2007)
New Revision: 6543

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/StgFEM/SLE/ProvidedSystems/StokesFlow/src/UpdateDt.c
Log:
 r1790 at earth:  boo | 2007-04-10 15:44:58 -0700
 Use a constant dtFactor from the input file to modify the automatically chosen dt



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1788
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1790

Modified: long/3D/Gale/trunk/src/StgFEM/SLE/ProvidedSystems/StokesFlow/src/UpdateDt.c
===================================================================
--- long/3D/Gale/trunk/src/StgFEM/SLE/ProvidedSystems/StokesFlow/src/UpdateDt.c	2007-04-10 21:59:45 UTC (rev 6542)
+++ long/3D/Gale/trunk/src/StgFEM/SLE/ProvidedSystems/StokesFlow/src/UpdateDt.c	2007-04-10 22:46:37 UTC (rev 6543)
@@ -57,6 +57,8 @@
 	FeVariable*             velFeVar      = NULL;
 	double			localDt;
 	double			globalDt;
+        Dictionary*             dictionary = context->dictionary;
+        double factor;
 	
 	velFeVar = self->uSolnVec->feVariable;
 	velMax = FieldVariable_GetMaxGlobalFieldMagnitude( velFeVar );
@@ -65,8 +67,9 @@
 		__func__, self->name, velFeVar->name, velMax );
 
 	FeVariable_GetMinimumSeparation( velFeVar, &minSeparation, minSeparationEachDim );
-
-	localDt = 0.5 * minSeparation / velMax;
+        
+        factor=Dictionary_GetDouble_WithDefault(dictionary,"dtFactor",1.0);
+	localDt = factor*0.5 * minSeparation / velMax;
 	MPI_Allreduce( &localDt, &globalDt, 1, MPI_DOUBLE, MPI_MIN, MPI_COMM_WORLD );
 
 	return globalDt;