Index: /issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m	(revision 17145)
+++ /issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m	(revision 17146)
@@ -102,6 +102,9 @@
 					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,6 +124,9 @@
 			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');