Index: /issm/trunk/src/m/model/regionaltransient2d.m
===================================================================
--- /issm/trunk/src/m/model/regionaltransient2d.m	(revision 10294)
+++ /issm/trunk/src/m/model/regionaltransient2d.m	(revision 10295)
@@ -37,5 +37,4 @@
 
 %remesh
-%md2=bamg(md2,'hmax',hmax,'hmin',hmin,'field',[md2.inversion.vel_obs md2.geometry.thickness],'splitcorner',1,'KeepVertices',0);
 md2=bamg(mde,'hmin',hmin,'hmax',hmax,'field',[mde.inversion.vel_obs mde.geometry.thickness],'splitcorner',1,'KeepVertices',0,'err',err);
 md2=setmask(md2,'','');
@@ -51,4 +50,5 @@
 		fieldsize=size(field);
 
+		%copy field, interpolated to new mesh
 		if isobject(field), %recursive call
 			object_fields=fields(md1.(model_fields{i}));
@@ -92,4 +92,5 @@
 	end
 
+	%Read transient velocities looping through only the populated times
 	spcx=[];
 	spcy=[];
@@ -109,4 +110,5 @@
 	end
 
+	%As long as there are recorded time steps, spc the boundaries with velocities
 	if nsteps > 0
 		md2.diagnostic.spcvx=md2.diagnostic.spcvx*ones(1,size(spcx,2));
@@ -121,5 +123,5 @@
 	end
 
-	%Diagnostic
+	%Diagnostic.  Don't spc the icefront vertices.
 	if ~isnan(md2.diagnostic.icefront)
 		md1s=modelextract(md1,area);