[cig-commits] r18377 - mc/1D/hc/trunk/anisotropic_viscosity

[becker at geodynamics.org]

Author: becker
Date: 2011-05-15 15:21:20 -0700 (Sun, 15 May 2011)
New Revision: 18377

Modified:
   mc/1D/hc/trunk/anisotropic_viscosity/Anisotropic_viscosity.c
Log:
Some more comments.



Modified: mc/1D/hc/trunk/anisotropic_viscosity/Anisotropic_viscosity.c
===================================================================
--- mc/1D/hc/trunk/anisotropic_viscosity/Anisotropic_viscosity.c	2011-05-15 22:13:13 UTC (rev 18376)
+++ mc/1D/hc/trunk/anisotropic_viscosity/Anisotropic_viscosity.c	2011-05-15 22:21:20 UTC (rev 18377)
@@ -389,7 +389,8 @@
 	myerror_s("set_anisotropic_viscosity_at_element_level: need to select layer",E);
       ani_layer = -E->viscosity.anivisc_layer;
 #ifdef CitcomS_global_defs_h	
-      /* CitcomS */
+      /* CitcomS, the zbase_layers are counted 0...1 top down, need to
+	 convert to radius */
       z_bottom = E->sphere.ro-E->viscosity.zbase_layer[ani_layer-1];
       if(ani_layer == 1)
 	z_top = E->sphere.ro;
@@ -397,6 +398,9 @@
 	z_top = E->sphere.ro - E->viscosity.zbase_layer[ani_layer-2];
 #else 
       /* CU */
+      /* in spherical system, for CitcomCU zbase_layer is actually
+	 given as radii. can use same as Cartesian
+      */
       z_bottom = E->viscosity.zbase_layer[ani_layer-1];
       if(ani_layer == 1)
 	z_top = E->segment.zzlayer[E->segment.zlayers-1];