Changeset 2950

Timestamp:

02/03/10 23:14:37 (15 years ago)

Author:

Eric.Larour

Message:

Added missing stuff: tell me Mathieu

File:

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

issm/trunk/src/m/classes/public/bamg.m

@@ -125 +125 @@
 md.numberofgrids=length(md.x);
 md.z=zeros(md.numberofgrids,1);
+md.gridonboundary=zeros(md.numberofgrids,1); md.gridonboundary(md.segments(:,1:2))=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);
+
+%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);
+md.gridonwater=zeros(md.numberofgrids,1);
+
+%type of model
 md.type='2d';