Index: /issm/trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m	(revision 24484)
+++ /issm/trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m	(revision 24485)
@@ -12,6 +12,11 @@
 
 	%Inserting time 0. md0 must be last experiment (e.g., Ice1r for Ice1ra)
-	x							= md0.results.TransientSolution(end).MeshX;
-	y							= md0.results.TransientSolution(end).MeshY;
+	if(isfield(md0.results.TransientSolution,'MeshElements'))
+      x     = md0.results.TransientSolution(end).MeshX;
+      y     = md0.results.TransientSolution(end).MeshY;
+   else
+      x     = md0.mesh.x;
+      y     = md0.mesh.y;
+   end	
 	time(1)					= 0;
 	[xgl_step ygl_step]	= gl_position(md0,length(md0.results.TransientSolution),0);
@@ -26,6 +31,11 @@
 
    for i=2:length(step),
-		x = md.results.TransientSolution(step(i)).MeshX;
-		y = md.results.TransientSolution(step(i)).MeshY;
+		if(isfield(md.results.TransientSolution,'MeshElements'))
+			x     = md.results.TransientSolution(step(i)).MeshX;
+			y     = md.results.TransientSolution(step(i)).MeshY;
+		else
+			x     = md.mesh.x;
+			y     = md.mesh.y;
+		end	
 		time(i)=md.results.TransientSolution(step(i)).time;	
 		[xgl_step ygl_step]=gl_position(md,step(i),0);