[cig-commits] r15700 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2009-09-28 13:05:03 -0700 (Mon, 28 Sep 2009)
New Revision: 15700

Modified:
   mc/3D/CitcomS/trunk/lib/Element_calculations.c
Log:
Fixed a bug in pseudo free surface, reported by Robert Moucha.


Modified: mc/3D/CitcomS/trunk/lib/Element_calculations.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Element_calculations.c	2009-09-28 20:04:38 UTC (rev 15699)
+++ mc/3D/CitcomS/trunk/lib/Element_calculations.c	2009-09-28 20:05:03 UTC (rev 15700)
@@ -142,6 +142,7 @@
   void get_elt_f();
   void get_elt_tr_pseudo_surf();
   void strip_bcs_from_residual();
+  double global_vdot();
 
   const int neq=E->lmesh.neq;
   const int nel=E->lmesh.nel;
@@ -173,6 +174,17 @@
 
   (E->solver.exchange_id_d)(E, E->F, lev);
   strip_bcs_from_residual(E,E->F,lev);
+
+  /* compute the norm of E->F */
+  E->monitor.fdotf = sqrt(global_vdot(E, E->F, E->F, lev));
+
+  if(E->parallel.me==0) {
+      fprintf(stderr, "Momentum equation force %.9e\n",
+              E->monitor.fdotf);
+      fprintf(E->fp, "Momentum equation force %.9e\n",
+              E->monitor.fdotf);
+  }
+
   return;
 }