Index: ../trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m
===================================================================
--- ../trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m	(revision 24484)
+++ ../trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m	(revision 24485)
@@ -11,8 +11,13 @@
    time=[];
 
 	%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);
 	xGL{1}					= xgl_step;
@@ -25,8 +30,13 @@
 	vBaseGL{1}				= griddata(x,y,md0.results.TransientSolution(end).Vy,xgl_step,ygl_step);
 
    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);
       xGL{i}=xgl_step;