Index: /issm/trunk/test/Validation/AnisotropicMesh/runme1.m
===================================================================
--- /issm/trunk/test/Validation/AnisotropicMesh/runme1.m	(revision 1344)
+++ /issm/trunk/test/Validation/AnisotropicMesh/runme1.m	(revision 1345)
@@ -8,6 +8,4 @@
 md=model;
 md=squaremesh(md,L,L,nx,ny);
-md.gridonboundary=zeros(md.numberofgrids,1);
-md.gridonboundary([1 ny ny*(nx-1)+1 nx*ny])=1;
 
 %mesh adaptation loop
@@ -16,6 +14,4 @@
 	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
 	md=YamsCall(md,md.vel_obs,0.001,0.3,1.3,10^-4);
-	md.gridonboundary=zeros(md.numberofgrids,1);
-	md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
 	plotmodel(md,'data','mesh');pause(0.5);
 end
@@ -23,6 +19,4 @@
 	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
 	md=YamsCall(md,md.vel_obs,0.001,0.3,2,10^-4);
-	md.gridonboundary=zeros(md.numberofgrids,1);
-	md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
 	plotmodel(md,'data','mesh');pause(0.5);
 end
@@ -30,6 +24,4 @@
 	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
 	md=YamsCall(md,md.vel_obs,0.001,0.3,2,0.005);
-	md.gridonboundary=zeros(md.numberofgrids,1);
-	md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
 	plotmodel(md,'data','mesh');pause(0.5);
 end
Index: /issm/trunk/test/Validation/AnisotropicMesh/runme2.m
===================================================================
--- /issm/trunk/test/Validation/AnisotropicMesh/runme2.m	(revision 1344)
+++ /issm/trunk/test/Validation/AnisotropicMesh/runme2.m	(revision 1345)
@@ -7,6 +7,4 @@
 md=model;
 md=squaremesh(md,L,L,nx,ny);
-md.gridonboundary=zeros(md.numberofgrids,1);
-md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
 
 %mesh adaptation loop
@@ -18,6 +16,4 @@
 		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
 	md=YamsCall(md,md.vel_obs,0.005,0.3,1.3,10^-4);
-	md.gridonboundary=zeros(md.numberofgrids,1);
-	md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
 	plotmodel(md,'data','mesh');pause(0.5);
 end
@@ -28,6 +24,4 @@
 		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
 	md=YamsCall(md,md.vel_obs,0.005,0.3,2,10^-4);
-	md.gridonboundary=zeros(md.numberofgrids,1);
-	md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
 	plotmodel(md,'data','mesh');pause(0.5);
 end
@@ -38,6 +32,4 @@
 		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
 	md=YamsCall(md,md.vel_obs,0.005,0.3,2,0.001);
-	md.gridonboundary=zeros(md.numberofgrids,1);
-	md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
 	plotmodel(md,'data','mesh');pause(0.5);
 end