Index: /issm/trunk/src/m/model/mesh/meshyams.m
===================================================================
--- /issm/trunk/src/m/model/mesh/meshyams.m	(revision 10801)
+++ /issm/trunk/src/m/model/mesh/meshyams.m	(revision 10802)
@@ -79,15 +79,15 @@
 	field=sqrt(vx_obs.^2+vy_obs.^2);
 
-	%set nodeonwater field
+	%set mask.vertexonwater  field
 	if ~strcmp(groundeddomain,'N/A'),
 		nodeground=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,groundeddomain,'node',2);
-		md.nodeonwater=ones(md.mesh.numberofvertices,1);
-		md.nodeonwater(find(nodeground))=0;
+		md.mask.vertexonwater=ones(md.mesh.numberofvertices,1);
+		md.mask.vertexonwater(find(nodeground))=0;
 	else
-		md.nodeonwater=zeros(md.mesh.numberofvertices,1);
+		md.mask.vertexonwater=zeros(md.mesh.numberofvertices,1);
 	end
 
 	%adapt according to velocities
-	disp('   adapting...');
+	disp('  Yo, Im adapting...');
 	md=YamsCall(md,field,hmin,hmax,gradation(i),epsilon);
 
@@ -121,8 +121,8 @@
 if ~strcmp(groundeddomain,'N/A'),
 	nodeground=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,groundeddomain,'node',2);
-	md.nodeonwater=ones(md.mesh.numberofvertices,1);
-	md.nodeonwater(find(nodeground))=0;
+	md.mask.vertexonwater=ones(md.mesh.numberofvertices,1);
+	md.mask.vertexonwater(find(nodeground))=0;
 else
-	md.nodeonwater=zeros(md.mesh.numberofvertices,1);
+	md.mask.vertexonwater=zeros(md.mesh.numberofvertices,1);
 end
 if strcmpi(Names.interp,'node'),