Index: /issm/trunk-jpl/src/m/classes/materials.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/materials.m	(revision 13093)
+++ /issm/trunk-jpl/src/m/classes/materials.m	(revision 13094)
@@ -6,15 +6,15 @@
 classdef materials
 	properties (SetAccess=public) 
-		rho_ice                    = 0;
-		rho_water                  = 0;
-		rho_freshwater             = 0;
-		mu_water                   = 0;
-		heatcapacity               = 0;
-		latentheat                 = 0;
-		thermalconductivity        = 0;
-		meltingpoint               = 0;
-		beta                       = 0;
-		mixed_layer_capacity       = 0;
-		thermal_exchange_velocity  = 0;
+		rho_ice                    = 0.;
+		rho_water                  = 0.;
+		rho_freshwater             = 0.;
+		mu_water                   = 0.;
+		heatcapacity               = 0.;
+		latentheat                 = 0.;
+		thermalconductivity        = 0.;
+		meltingpoint               = 0.;
+		beta                       = 0.;
+		mixed_layer_capacity       = 0.;
+		thermal_exchange_velocity  = 0.;
 		rheology_B   = NaN;
 		rheology_n   = NaN;
@@ -33,11 +33,11 @@
 
 			%ice density (kg/m^3)
-			obj.rho_ice=917;
+			obj.rho_ice=917.;
 
 			%ocean water density (kg/m^3)
-			obj.rho_water=1023;
+			obj.rho_water=1023.;
 
 			%fresh water density (kg/m^3)
-			obj.rho_freshwater=1000;
+			obj.rho_freshwater=1000.;
 
 			%water viscosity (N.s/m^2)
@@ -45,5 +45,5 @@
 
 			%ice heat capacity cp (J/kg/K)
-			obj.heatcapacity=2093;
+			obj.heatcapacity=2093.;
 
 			%ice latent heat of fusion L (J/kg)
@@ -60,5 +60,5 @@
 
 			%mixed layer (ice-water interface) heat capacity (J/kg/K)
-			obj.mixed_layer_capacity=3974;
+			obj.mixed_layer_capacity=3974.;
 
 			%thermal exchange velocity (ice-water interface) (m/s)
Index: /issm/trunk-jpl/src/m/classes/materials.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/materials.py	(revision 13093)
+++ /issm/trunk-jpl/src/m/classes/materials.py	(revision 13094)
@@ -17,15 +17,15 @@
 	def __init__(self):
 		# {{{ Properties
-		self.rho_ice                    = 0
-		self.rho_water                  = 0
-		self.rho_freshwater             = 0
-		self.mu_water                   = 0
-		self.heatcapacity               = 0
-		self.latentheat                 = 0
-		self.thermalconductivity        = 0
-		self.meltingpoint               = 0
-		self.beta                       = 0
-		self.mixed_layer_capacity       = 0
-		self.thermal_exchange_velocity  = 0
+		self.rho_ice                    = 0.
+		self.rho_water                  = 0.
+		self.rho_freshwater             = 0.
+		self.mu_water                   = 0.
+		self.heatcapacity               = 0.
+		self.latentheat                 = 0.
+		self.thermalconductivity        = 0.
+		self.meltingpoint               = 0.
+		self.beta                       = 0.
+		self.mixed_layer_capacity       = 0.
+		self.thermal_exchange_velocity  = 0.
 		self.rheology_B   = float('NaN')
 		self.rheology_n   = float('NaN')
@@ -58,11 +58,11 @@
 		# {{{setdefaultparameters
 		#ice density (kg/m^3)
-		self.rho_ice=917
+		self.rho_ice=917.
 
 		#ocean water density (kg/m^3)
-		self.rho_water=1023
+		self.rho_water=1023.
 
 		#fresh water density (kg/m^3)
-		self.rho_freshwater=1000
+		self.rho_freshwater=1000.
 
 		#water viscosity (N.s/m^2)
@@ -70,5 +70,5 @@
 
 		#ice heat capacity cp (J/kg/K)
-		self.heatcapacity=2093
+		self.heatcapacity=2093.
 
 		#ice latent heat of fusion L (J/kg)
@@ -85,5 +85,5 @@
 
 		#mixed layer (ice-water interface) heat capacity (J/kg/K)
-		self.mixed_layer_capacity=3974
+		self.mixed_layer_capacity=3974.
 
 		#thermal exchange velocity (ice-water interface) (m/s)