Index: /issm/trunk/src/m/classes/public/mesh/meshbamg.m
===================================================================
--- /issm/trunk/src/m/classes/public/mesh/meshbamg.m	(revision 3181)
+++ /issm/trunk/src/m/classes/public/mesh/meshbamg.m	(revision 3182)
@@ -82,37 +82,7 @@
 	%adapt according to velocities
 	disp('   adapting...');
-	md=BamgCall(md,field,hmin,hmax,gradation(i),epsilon);
+	md=bamg(md,'field',field,'hmin',hmin,'hmax',hmax,'gradation',gradation(i),'err',epsilon);
 
 end
 	
 disp(['Final mesh, number of elements: ' num2str(md.numberofelements)]);
-
-%Now, build the connectivity tables for this mesh.
-md.nodeconnectivity=NodeConnectivity(md.elements,md.numberofgrids);
-md.elementconnectivity=ElementConnectivity(md.elements,md.nodeconnectivity);
-
-%recreate segments
-md.segments=findsegments(md);
-md.gridonboundary=zeros(md.numberofgrids,1); md.gridonboundary(md.segments(:,1:2))=1;
-
-%Fill in rest of fields:
-md.z=zeros(md.numberofgrids,1);
-md.gridonbed=ones(md.numberofgrids,1);
-md.gridonsurface=ones(md.numberofgrids,1);
-md.elementonbed=ones(md.numberofelements,1);
-md.elementonsurface=ones(md.numberofelements,1);
-if ~strcmp(groundeddomain,'N/A'),
-	gridground=ContourToMesh(md.elements,md.x,md.y,expread(groundeddomain,1),'node',2);
-	md.gridonwater=ones(md.numberofgrids,1);
-	md.gridonwater(find(gridground))=0;
-else
-	md.gridonwater=zeros(md.numberofgrids,1);
-end
-if strcmpi(Names.interp,'grid'),
-	md.vx_obs=InterpFromGridToMesh(Vel.(Names.xname),Vel.(Names.yname),Vel.(Names.vxname),md.x,md.y,0);
-	md.vy_obs=InterpFromGridToMesh(Vel.(Names.xname),Vel.(Names.yname),Vel.(Names.vyname),md.x,md.y,0);
-else
-	md.vx_obs=InterpFromMeshToMesh2d(Vel.(Names.indexname),Vel.(Names.xname),Vel.(Names.yname),Vel.(Names.vxname),md.x,md.y,0);
-	md.vy_obs=InterpFromMeshToMesh2d(Vel.(Names.indexname),Vel.(Names.xname),Vel.(Names.yname),Vel.(Names.vyname),md.x,md.y,0);
-end
-md.vel_obs=sqrt(md.vx_obs.^2+md.vy_obs.^2);