Index: /issm/trunk-jpl/src/m/contrib/paraview/writeVTKcell.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/paraview/writeVTKcell.m	(revision 15274)
+++ /issm/trunk-jpl/src/m/contrib/paraview/writeVTKcell.m	(revision 15275)
@@ -15,4 +15,8 @@
 % 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,5 +46,5 @@
 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,5 +84,5 @@
 	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);