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

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-02-05 18:22:32 -0800 (Mon, 05 Feb 2007)
New Revision: 5966

Modified:
   mc/3D/CitcomS/trunk/lib/Solver_multigrid.c
Log:
Changed smoothing method for viscosity to be consistent with CitcomCU

Modified: mc/3D/CitcomS/trunk/lib/Solver_multigrid.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Solver_multigrid.c	2007-02-05 23:43:15 UTC (rev 5965)
+++ mc/3D/CitcomS/trunk/lib/Solver_multigrid.c	2007-02-06 02:22:32 UTC (rev 5966)
@@ -256,12 +256,9 @@
     sl_minus = lv -1;
 
     if (E->viscosity.smooth_cycles==1)  {
-/*      visc_from_gint_to_nodes(E,E->EVI[lv],viscU,lv);
+      visc_from_gint_to_nodes(E,E->EVI[lv],viscU,lv);
       project_scalar(E,lv,viscU,viscD);
-      visc_from_nodes_to_gint(E,viscD,E->EVI[sl_minus],sl_minus);  */
-
-      project_scalar(E,lv,E->VI[lv],E->VI[sl_minus]);
-      visc_from_nodes_to_gint(E,E->VI[sl_minus],E->EVI[sl_minus],sl_minus);
+      visc_from_nodes_to_gint(E,viscD,E->EVI[sl_minus],sl_minus);
       }
     else if (E->viscosity.smooth_cycles==2)   {
       visc_from_gint_to_ele(E,E->EVI[lv],viscU,lv);