Index: /issm/trunk/src/m/model/regionaltransient2d.m
===================================================================
--- /issm/trunk/src/m/model/regionaltransient2d.m	(revision 10586)
+++ /issm/trunk/src/m/model/regionaltransient2d.m	(revision 10587)
@@ -37,5 +37,5 @@
 
 %remesh
-md2=bamg(mde,'hmin',hmin,'hmax',hmax,'field',[mde.inversion.vel_obs mde.geometry.thickness],'splitcorner',1,'KeepVertices',0,'err',err);
+md2=bamg(mde,'hmin',hmin,'hmax',hmax,'field',[mde.inversion.vel_obs mde.geometry.surface],'splitcorner',1,'KeepVertices',0,'err',err);
 md2=setmask(md2,'','');
 
@@ -92,7 +92,8 @@
 	end
 
-	%Read transient velocities looping through only the populated times
+	%Read transient velocities and thickness, looping through only the populated times
 	spcx=[];
 	spcy=[];
+	spct=[];
 	steps=[];
 	nsteps=length(md1.results.TransientSolution);
@@ -103,6 +104,8 @@
 			vx=PatchToVec(md1.results.TransientSolution(t).Vx);
 			vy=PatchToVec(md1.results.TransientSolution(t).Vy);
+			thickness=PatchToVec(md1.results.TransientSolution(t).Thickness);
 			spcx=[spcx InterpFromMeshToMesh2d(md1.mesh.elements,md1.mesh.x,md1.mesh.y,vx,md2.mesh.x,md2.mesh.y)];
 			spcy=[spcy InterpFromMeshToMesh2d(md1.mesh.elements,md1.mesh.x,md1.mesh.y,vx,md2.mesh.x,md2.mesh.y)];
+			spct=[spct InterpFromMeshToMesh2d(md1.mesh.elements,md1.mesh.x,md1.mesh.y,thickness,md2.mesh.x,md2.mesh.y)];
 			steps=[steps t*md1.timestepping.time_step];
 		end
@@ -115,10 +118,13 @@
 		md2.diagnostic.spcvy=md2.diagnostic.spcvy*ones(1,size(spcy,2));
 		md2.diagnostic.spcvz=md2.diagnostic.spcvz*ones(1,size(spcx,2));
+		md2.prognostic.spcthickness=md2.prognostic.spcthickness*ones(1,size(spct,2));
 		md2.diagnostic.spcvx(find(md2.mesh.vertexonboundary),:)=spcx(find(md2.mesh.vertexonboundary),:);
 		md2.diagnostic.spcvy(find(md2.mesh.vertexonboundary),:)=spcy(find(md2.mesh.vertexonboundary),:);
 		md2.diagnostic.spcvz(find(md2.mesh.vertexonboundary),:)=0;
+		md2.prognostic.spcthickness(find(md2.mesh.vertexonboundary),:)=spct(find(md2.mesh.vertexonboundary),:);
 		md2.diagnostic.spcvx=[md2.diagnostic.spcvx; steps];
 		md2.diagnostic.spcvy=[md2.diagnostic.spcvy; steps];
 		md2.diagnostic.spcvz=[md2.diagnostic.spcvz; steps];
+		md2.prognostic.spcthickness=[md2.prognostic.spcthickness; steps];
 	end
 
@@ -143,4 +149,5 @@
 		md2.diagnostic.spcvy(unique(md2.diagnostic.icefront(:,1:2)),:)=nan;
 		md2.diagnostic.spcvz(unique(md2.diagnostic.icefront(:,1:2)),:)=nan;
+		md2.prognostic.spcthickness(unique(md2.diagnostic.icefront(:,1:2)),:)=nan;
 	end