Changeset 17146

Timestamp:

01/21/14 12:32:55 (11 years ago)

Author:

bdef

Message:

BUG:fixing minor issue with almost zero values

File:

• issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m (modified) (2 diffs)

issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m

@@ -102 +102 @@
                                         num_of_points);
                                 %paraview does not like NaN, replacing
-                                pos=find(isnan(sol_struct{j}(timestep).(fieldnames{k})));
-                                sol_struct{j}(timestep).(fieldnames{k})(pos)=-9999;
+                                nanval=find(isnan(sol_struct{j}(timestep).(fieldnames{k})));
+                                sol_struct{j}(timestep).(fieldnames{k})(nanval)=-9999;
+                                %also checking for verry small value that mess up
+                                smallval=(abs(sol_struct{j}(timestep).(fieldnames{k}))<1.0e-20);
+                                sol_struct{j}(timestep).(fieldnames{k})(smallval)=0.0;
                                 fprintf(FID,'SCALARS %s float 1 \n',fieldnames{k});
                                 fprintf(FID,'LOOKUP_TABLE default\n');
@@ -121 +124 @@
                         if ((numel(res_struct.(fieldnames{k})))==num_of_points);
                                 %paraview does not like NaN, replacing
-                                pos=find(isnan(res_struct.(fieldnames{k})));
-                                res_struct.(fieldnames{k})(pos)=-9999;
+                                nanval=find(isnan(res_struct.(fieldnames{k})));
+                                res_struct.(fieldnames{k})(nanval)=-9999;
+                                %also checking for verry small value that mess up
+                                smallval=(abs(res_struct.(fieldnames{k}))<1.0e-20);
+                                res_struct.(fieldnames{k})(smallval)=0.0;
                                 fprintf(FID,'SCALARS %s float 1 \n',fieldnames{k});
                                 fprintf(FID,'LOOKUP_TABLE default\n');