Changeset 16500

Timestamp:

10/22/13 08:32:11 (12 years ago)

Author:

bdef

Message:

adding a NaN search and replace

File:

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

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

@@ -91 +91 @@
                 %check which field is a real result and print
                         for k=1:num_of_fields
-                        if ((length(sol_struct{j}(step).(fieldnames{k})))==num_of_points);
+                        if ((length(sol_struct{j}(step).(fieldnames{k})))== ...
+                                        num_of_points);
+                                %paraview does not like NaN, replacing
+                                pos=find(isnan(sol_struct{j}(step).(fieldnames{k})));
+                                sol_struct{j}(step).(fieldnames{k})(pos)=-9999;
                                 fprintf(FID,'SCALARS %s float 1 \n',fieldnames{k});
                                 fprintf(FID,'LOOKUP_TABLE default\n');
@@ -108 +112 @@
                 for k=1:num_of_sols
                         if ((length(res_struct.(solnames{k})))==num_of_points);
+                                %paraview does not like NaN, replacing
+                                pos=find(isnan(res_struct.(solnames{k})));
+                                res_struct.(solnames{k})(pos)=-9999;
                                 fprintf(FID,'SCALARS %s float 1 \n',solnames{k});
                                 fprintf(FID,'LOOKUP_TABLE default\n');