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

[tan2 at geodynamics.org]

Author: tan2
Date: 2009-05-01 15:17:27 -0700 (Fri, 01 May 2009)
New Revision: 14837

Modified:
   mc/3D/CitcomS/trunk/lib/Instructions.c
Log:
Printing a big warning when plate velocity is imposed but remove_angular_momentum or remove_rigid_rotation is on.

Modified: mc/3D/CitcomS/trunk/lib/Instructions.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Instructions.c	2009-05-01 22:15:05 UTC (rev 14836)
+++ mc/3D/CitcomS/trunk/lib/Instructions.c	2009-05-01 22:17:27 UTC (rev 14837)
@@ -752,6 +752,13 @@
         fprintf(stderr,"solve_Ahat_p_fhat: WARNING: overriding pressure and div check\n");
     }
 
+    /* remove angular momentum/rigid rotation should only be done in free
+       convection of global models. */
+    if(E->sphere.caps == 12 &&
+       (E->control.remove_angular_momentum || E->control.remove_rigid_rotation) &&
+       (E->mesh.topvbc || E->mesh.botvbc || E->control.side_sbcs)) {
+        fprintf(stderr,"\nWARNING: The input parameters impose boundary velocity, but also remove angular momentum/rigid rotation!\n\n");
+    }
     return;
 }