Index: /issm/trunk-jpl/src/m/classes/thermal.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/thermal.py	(revision 20983)
+++ /issm/trunk-jpl/src/m/classes/thermal.py	(revision 20984)
@@ -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=numpy.nan*numpy.ones((md.mesh.numberofvertices,1))
+			self.spctemperature=float('nan')*numpy.ones(md.mesh.numberofvertices)
 			pos=numpy.nonzero(md.mesh.vertexonsurface)[0]
 			self.spctemperature[pos]=md.initialization.temperature[pos]    #impose observed temperature on surface