[cig-commits] commit 1943 by heister to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 1943

catch case where temperature or composition is all zero and entropy viscosity generates NaN

U   trunk/aspect/source/simulator/assembly.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=1943&peg=1943

Diff:
Modified: trunk/aspect/source/simulator/assembly.cc
===================================================================
--- trunk/aspect/source/simulator/assembly.cc	2013-10-07 21:41:42 UTC (rev 1942)
+++ trunk/aspect/source/simulator/assembly.cc	2013-10-08 15:13:56 UTC (rev 1943)
@@ -595,7 +595,7 @@
                      const double                        cell_diameter,
                      const TemperatureOrComposition     &temperature_or_composition) const
   {
-    if (global_u_infty == 0)
+    if (global_u_infty == 0 || global_entropy_variation == 0 || global_field_variation == 0)
       return 5e-3 * cell_diameter;
 
     double max_residual = 0;