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

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-12-17 14:06:06 -0800 (Mon, 17 Dec 2007)
New Revision: 8893

Modified:
   mc/3D/CitcomS/trunk/lib/Output_h5.c
Log:
Fixed phase change width

Modified: mc/3D/CitcomS/trunk/lib/Output_h5.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Output_h5.c	2007-12-17 21:56:38 UTC (rev 8892)
+++ mc/3D/CitcomS/trunk/lib/Output_h5.c	2007-12-17 22:06:06 UTC (rev 8893)
@@ -1440,17 +1440,26 @@
     status = set_attribute_float(input, "Ra_410", E->control.Ra_410);
     status = set_attribute_float(input, "clapeyron410", E->control.clapeyron410);
     status = set_attribute_float(input, "transT410", E->control.transT410);
-    status = set_attribute_float(input, "width410", E->control.width410);
+    status = set_attribute_float(input, "width410",
+                                 (E->control.inv_width410 == 0)?
+                                 E->control.inv_width410 :
+				 1.0/E->control.inv_width410);
 
     status = set_attribute_float(input, "Ra_670", E->control.Ra_670);
     status = set_attribute_float(input, "clapeyron670", E->control.clapeyron670);
     status = set_attribute_float(input, "transT670", E->control.transT670);
-    status = set_attribute_float(input, "width670", E->control.width670);
+    status = set_attribute_float(input, "width670",
+                                 (E->control.inv_width670 == 0)?
+                                 E->control.inv_width670 :
+				 1.0/E->control.inv_width670);
 
     status = set_attribute_float(input, "Ra_cmb", E->control.Ra_cmb);
     status = set_attribute_float(input, "clapeyroncmb", E->control.clapeyroncmb);
     status = set_attribute_float(input, "transTcmb", E->control.transTcmb);
-    status = set_attribute_float(input, "widthcmb", E->control.widthcmb);
+    status = set_attribute_float(input, "widthcmb",
+                                 (E->control.inv_widthcmb == 0)?
+                                 E->control.inv_widthcmb :
+				 1.0/E->control.inv_widthcmb);
 
     /*
      * Solver.inventory