Changeset 3182

Timestamp:

03/04/10 10:00:08 (15 years ago)

Author:

Mathieu Morlighem

Message:

minor

File:

• issm/trunk/src/m/classes/public/mesh/meshbamg.m (modified) (1 diff)

issm/trunk/src/m/classes/public/mesh/meshbamg.m

@@ -82 +82 @@
         %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);