Index: /issm/trunk-jpl/src/m/classes/initialization.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.py	(revision 20943)
+++ /issm/trunk-jpl/src/m/classes/initialization.py	(revision 20944)
@@ -62,5 +62,5 @@
 
 		#Lithostatic pressure by default
-		self.pressure=md.constants.g*md.materials.rho_ice*(md.geometry.surface-md.mesh.z)
+		self.pressure=md.constants.g*md.materials.rho_ice*(md.geometry.surface-md.mesh.z.reshape(-1,1))
 		return self
 	#}}}