[cig-commits] r12474 - mc/3D/CitcomS/trunk/module

[tan2 at geodynamics.org]

Author: tan2
Date: 2008-07-24 11:42:45 -0700 (Thu, 24 Jul 2008)
New Revision: 12474

Modified:
   mc/3D/CitcomS/trunk/module/setProperties.c
Log:
E->viscosity.zbase_layer is not init'd in pyre version.
This fixes issue 155.


Modified: mc/3D/CitcomS/trunk/module/setProperties.c
===================================================================
--- mc/3D/CitcomS/trunk/module/setProperties.c	2008-07-24 18:21:37 UTC (rev 12473)
+++ mc/3D/CitcomS/trunk/module/setProperties.c	2008-07-24 18:42:45 UTC (rev 12474)
@@ -208,6 +208,11 @@
     getFloatProperty(properties, "z_lmantle", E->viscosity.zlm, fp);
     getFloatProperty(properties, "z_cmb", E->viscosity.zcmb, fp); /* this is used as the D" phase change depth */
 
+    E->viscosity.zbase_layer[0] = E->viscosity.zlith;
+    E->viscosity.zbase_layer[1] = E->viscosity.z410;
+    E->viscosity.zbase_layer[2] = E->viscosity.zlm;
+    E->viscosity.zbase_layer[3] = E->viscosity.zcmb;
+    
     /* convert meter to kilometer */
     E->data.radius_km = radius / 1e3;