Changeset 20984

Timestamp:

07/22/16 15:31:34 (9 years ago)

Author:

Mathieu Morlighem

Message:

BUG: do not make spcvz a 3d vector in py

File:

• issm/trunk-jpl/src/m/classes/thermal.py (modified) (1 diff)

issm/trunk-jpl/src/m/classes/thermal.py

@@ -47 +47 @@
                 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