Ignore:
Timestamp:
09/06/11 13:40:47 (14 years ago)
Author:
Mathieu Morlighem
Message:

Added rifts object

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/model/mesh/meshyams.m

    r8298 r9619  
    9494        %if we have rifts, we just messed them up, we need to recreate the segments that constitute those
    9595        %rifts, because the segments are used in YamsCall to freeze the rifts elements during refinement.
    96         if md.numrifts,
     96        if md.rifts.numrifts,
    9797                md.nodeconnectivity=NodeConnectivity(md.elements,md.numberofnodes);
    9898                md.elementconnectivity=ElementConnectivity(md.elements,md.nodeconnectivity);
     
    136136
    137137%deal with rifts
    138 if md.numrifts,
     138if md.rifts.numrifts,
    139139        %first, recreate rift segments
    140140        md=meshyamsrecreateriftsegments(md);
    141141
    142142        %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);
    145145
    146146                %build normals and lengths of segments:
     
    167167                rift.penaltypairs(:,6)=rift.penaltypairs(:,6)./norms;
    168168
    169                 md.rifts(j)=rift;
     169                md.rifts.riftstruct(j)=rift;
    170170        end
    171171
Note: See TracChangeset for help on using the changeset viewer.