Changeset 18513

Timestamp:

09/12/14 12:12:30 (11 years ago)

Author:

bdef

Message:

NEW:adding vtk output for python

Location:

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

Files:

• enveloppeVTK.m (added)
• exportVTK.m (modified) (2 diffs)
• exportVTK.py (added)

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

@@ -25 +25 @@
 if dimension(model.mesh)==2,
         points=[model.mesh.x model.mesh.y zeros(model.mesh.numberofvertices,1)];
+        [num_of_points,dim]=size(points);
+        [num_of_elt]=size(model.mesh.elements,1);
+        [point_per_elt]=size(model.mesh.elements,2);
 else
         points=[model.mesh.x model.mesh.y model.mesh.z];
+        [num_of_points,dim]=size(points);
+        [num_of_elt]=size(model.mesh.elements,1);
+        [point_per_elt]=size(model.mesh.elements,2);
 end
 
-[num_of_points,dim]=size(points);
-[num_of_elt]=size(model.mesh.elements,1);
-[point_per_elt]=size(model.mesh.elements,2);
 
 %Select the type of element function of the number of nodes per elements
@@ -86 +89 @@
   end
         s=cell2mat(horzcat(s,{'\n'}));
-        fprintf(fid,s,[(point_per_elt)*ones(num_of_elt,1) model.mesh.elements-1]');
+                fprintf(fid,s,[(point_per_elt)*ones(num_of_elt,1)       model.mesh.elements-1]');
         
         fprintf(fid,'CELL_TYPES %d\n',num_of_elt);