[cig-commits] commit: Output density equation and alpha equation into the journal

[Mercurial]

changeset:   424:26dd025f99b5
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Sep 14 14:42:29 2011 -0700
files:       Utils/src/BuoyancyForceTerm.cxx
description:
Output density equation and alpha equation into the journal


diff -r ff33925c5298 -r 26dd025f99b5 Utils/src/BuoyancyForceTerm.cxx
--- a/Utils/src/BuoyancyForceTerm.cxx	Mon Sep 12 14:32:07 2011 -0700
+++ b/Utils/src/BuoyancyForceTerm.cxx	Wed Sep 14 14:42:29 2011 -0700
@@ -474,6 +474,9 @@ void _BuoyancyForceTerm_AssembleElement(
                     d.SetVarFactory(BuoyancyForceTerm_AddVariable, &d);
                     d.SetExpr(density_equation);
                     density=d.Eval();
+
+                    Journal_Printf(Journal_MyStream(Info_Type,self),"Density Equation T=%g p=%g density=%g\n",
+                                    temperature,pressure,density);
                   }
                 else
                   {
@@ -491,6 +494,8 @@ void _BuoyancyForceTerm_AssembleElement(
                     d.SetVarFactory(BuoyancyForceTerm_AddVariable, &d);
                     d.SetExpr(alpha_equation);
                     alpha=d.Eval();
+                    Journal_PrintfL(Journal_MyStream(Debug_Type,self),3,"Alpha Equation T=%g p=%g alpha=%g\n",
+                                    temperature,pressure,alpha);
                   }
                 else
                   {