Changeset 20280

Timestamp:

03/01/16 20:16:04 (9 years ago)

Author:

Eric.Larour

Message:

CHG: fixed issue with loading force, which can be NaN!

File:

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

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

@@ -195 +195 @@
                 WriteData(fid,'object',self,'class','stressbalance','fieldname','referential','format','DoubleMat','mattype',1)
 
-                WriteData(fid,'data',self.loadingforce[:,0],'format','DoubleMat','mattype',1,'enum',LoadingforceXEnum())
-                WriteData(fid,'data',self.loadingforce[:,1],'format','DoubleMat','mattype',1,'enum',LoadingforceYEnum())
-                WriteData(fid,'data',self.loadingforce[:,2],'format','DoubleMat','mattype',1,'enum',LoadingforceZEnum())
+                if numpy.isnan(self.loadingforce):
+                        lx=float('NaN'); ly=float('NaN'); lz=float('NaN');
+                else:
+                        lx=self.loadingforce[:,0];
+                        ly=self.loadingforce[:,1];
+                        lz=self.loadingforce[:,2];
+
+                WriteData(fid,'data',lx,'format','DoubleMat','mattype',1,'enum',LoadingforceXEnum())
+                WriteData(fid,'data',ly,'format','DoubleMat','mattype',1,'enum',LoadingforceYEnum())
+                WriteData(fid,'data',lz,'format','DoubleMat','mattype',1,'enum',LoadingforceZEnum())
 
                 #process requested outputs