Changeset 15275

Timestamp:

06/17/13 16:20:02 (12 years ago)

Author:

bdef

Message:

CHG:Adding a directory to put the result in the vtk Output

File:

• issm/trunk-jpl/src/m/contrib/paraview/writeVTKcell.m (modified) (3 diffs)

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

@@ -15 +15 @@
 % Basile de Fleurian, modified from:
 % (c) Daniel Peterseim, 2009-11-07
+
+[path,name,ext]=fileparts(filename)
+separator=filesep
+mkdir(filename)
 
 points=[model.mesh.x model.mesh.y model.mesh.z];
@@ -42 +46 @@
 if num_of_timesteps==1; %just one timestep only write one file
 
-        FID = fopen(strcat(filename,'.vtk'),'w+');
+        FID = fopen(strcat(path,filesep,name,filesep,name,'.vtk'),'w+');
         fprintf(FID,'# vtk DataFile Version 2.0 \n');
         fprintf(FID,'Data for run %s \n',model.miscellaneous.name);
@@ -80 +84 @@
         for i=1:num_of_timesteps;
                 timestep=sol_struct(i).step;
-                FID = fopen(strcat(filename,'.vtk',int2str(timestep),'.vtk'),'w+');
+                FID = fopen(strcat(path,filesep,name,filesep,name,'.vtk',int2str(timestep),'.vtk'),'w+');
                 fprintf(FID,'# vtk DataFile Version 2.0 \n');
                 fprintf(FID,'Data for run %s \n',model.miscellaneous.name);