Index: /issm/trunk-jpl/src/m/classes/materials.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/materials.m	(revision 26350)
+++ /issm/trunk-jpl/src/m/classes/materials.m	(revision 26351)
@@ -326,5 +326,4 @@
 					end
 					earth_density=earth_density/self.radius(self.numlayers+1)^3;
-					disp(earth_density)
 					self.earth_density=earth_density;
 				case 'hydro'
Index: /issm/trunk-jpl/src/m/classes/materials.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/materials.py	(revision 26350)
+++ /issm/trunk-jpl/src/m/classes/materials.py	(revision 26351)
@@ -307,5 +307,4 @@
                     earth_density = earth_density + (pow(self.radius[i + 1], 3) - pow(self.radius[i], 3)) * self.density[i]
                 earth_density = earth_density / pow(self.radius[self.numlayers], 3)
-                print(earth_density)
                 self.earth_density = earth_density
             elif nat == 'hydro':