Changeset 9619 for issm/trunk/src/m/model/mesh/meshyams.m
- Timestamp:
- 09/06/11 13:40:47 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/mesh/meshyams.m
r8298 r9619 94 94 %if we have rifts, we just messed them up, we need to recreate the segments that constitute those 95 95 %rifts, because the segments are used in YamsCall to freeze the rifts elements during refinement. 96 if md. numrifts,96 if md.rifts.numrifts, 97 97 md.nodeconnectivity=NodeConnectivity(md.elements,md.numberofnodes); 98 98 md.elementconnectivity=ElementConnectivity(md.elements,md.nodeconnectivity); … … 136 136 137 137 %deal with rifts 138 if md. numrifts,138 if md.rifts.numrifts, 139 139 %first, recreate rift segments 140 140 md=meshyamsrecreateriftsegments(md); 141 141 142 142 %using the segments, recreate the penaltypairs 143 for j=1:md. numrifts,144 rift=md.rifts (j);143 for j=1:md.rifts.numrifts, 144 rift=md.rifts.riftstruct(j); 145 145 146 146 %build normals and lengths of segments: … … 167 167 rift.penaltypairs(:,6)=rift.penaltypairs(:,6)./norms; 168 168 169 md.rifts (j)=rift;169 md.rifts.riftstruct(j)=rift; 170 170 end 171 171
Note:
See TracChangeset
for help on using the changeset viewer.