Changeset 994


Ignore:
Timestamp:
06/13/09 20:36:21 (16 years ago)
Author:
seroussi
Message:

minor, do not use two different pos

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/solutions/ice/ModelProcessorThermal.m

    r992 r994  
    5252%the formulation into trias:
    5353el6pos=find((md.elements_type(:,1)==macayealenum() | md.elements_type(:,1)==pattynenum()) & (element_partitioning==labindex));
    54 pos=find(ismember(el6pos,find(md.elements_type(:,1)==macayealenum())));
    55 elements(pos)=SetStructureField(elements(pos),'element','collapse',ones(length(pos),1));
     54pos2=find(ismember(el6pos,find(md.elements_type(:,1)==macayealenum())));
     55elements(pos2)=SetStructureField(elements(pos2),'element','collapse',ones(length(pos2),1));
    5656
    5757if strcmpi(solutiontype,'thermalsteady'),
Note: See TracChangeset for help on using the changeset viewer.