Index: /issm/trunk-jpl/src/m/classes/initialization.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.py	(revision 20990)
+++ /issm/trunk-jpl/src/m/classes/initialization.py	(revision 20991)
@@ -62,6 +62,6 @@
 
 		#Lithostatic pressure by default
-		#self.pressure=md.constants.g*md.materials.rho_ice*(md.geometry.surface[:,0]-md.mesh.z)
-		self.pressure=md.constants.g*md.materials.rho_ice*(md.geometry.surface-md.mesh.z.reshape(-1,1))
+		self.pressure=md.constants.g*md.materials.rho_ice*(md.geometry.surface[:,0]-md.mesh.z)
+		#self.pressure=md.constants.g*md.materials.rho_ice*(md.geometry.surface-md.mesh.z.reshape(-1,1))
 		return self
 	#}}}
Index: /issm/trunk-jpl/src/m/classes/thermal.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/thermal.py	(revision 20990)
+++ /issm/trunk-jpl/src/m/classes/thermal.py	(revision 20991)
@@ -47,5 +47,5 @@
 		self.spctemperature=project3d(md,'vector',self.spctemperature,'type','node','layer',md.mesh.numberoflayers,'padding',numpy.nan)
 		if isinstance(md.initialization.temperature,numpy.ndarray) and numpy.size(md.initialization.temperature,axis=0)==md.mesh.numberofvertices:
-			self.spctemperature=float('nan')*numpy.ones(md.mesh.numberofvertices)
+			self.spctemperature=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
 			pos=numpy.nonzero(md.mesh.vertexonsurface)[0]
 			self.spctemperature[pos]=md.initialization.temperature[pos]    #impose observed temperature on surface