[cig-commits] r13854 - in long/3D/Gale/trunk: . src/Gale/Utils/src

[walter at geodynamics.org]

Author: walter
Date: 2009-01-12 16:03:38 -0800 (Mon, 12 Jan 2009)
New Revision: 13854

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c
Log:
 r2450 at dante:  boo | 2009-01-12 10:51:31 -0800
 Minor cleanup



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

Modified: long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c
===================================================================
--- long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c	2009-01-13 00:03:36 UTC (rev 13853)
+++ long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c	2009-01-13 00:03:38 UTC (rev 13854)
@@ -407,30 +407,6 @@
                 stress=
                   -HydrostaticTerm_Pressure(self->_entryTbl[entry_I].hydrostaticTerm,
                                            Mesh_GetVertex(mesh,elementNodes[eNode_I]));
-/*                 if(*(self->context->nonLinearIteration_I)!=0) */
-/*                   { */
-/*                     double density=1; */
-/*                     double gravity=1; */
-/*                     double v[3]; */
-/*                     double damping; */
-/*                     FeVariable_GetValueAtNode(velocityField, */
-/*                                               elementNodes[eNode_I],v); */
-                    
-/*                     printf("stress %d %d %d %g %g %g %g\n", */
-/*                            *(self->context->nonLinearIteration_I), */
-/*                            lElement_I, */
-/*                            elementNodes[eNode_I],stress, */
-/*                            density*gravity*self->context->dt*v[1], */
-/* /\*                            stress-density*gravity*1.09375*v[1],v[1]); *\/ */
-/*                            stress-density*gravity*self->context->dt*v[1],v[1]); */
-                    
-/* /\*                     damping=density*gravity*1.09375*v[1]; *\/ */
-/*                     damping=density*gravity*self->context->dt*v[1]; */
-/*                     /\* Make sure that we don't overcompensate. *\/ */
-/* /\*                     if(((damping>0 && stress>=0) || (damping<0 && stress<=0)) *\/ */
-/* /\*                        && fabs(damping)<fabs(stress)) *\/ */
-/* /\*                       stress+=damping; *\/ */
-/*                   } */
                 break;
               }
             /* We have to divide by an overcount_factor, because