Changeset 992

Timestamp:

06/13/09 20:20:36 (16 years ago)

Author:

seroussi

Message:

use Bbar as in cielo when collapsed elements

File:

• issm/trunk/src/m/solutions/ice/ModelProcessorThermal.m (modified) (1 diff)

issm/trunk/src/m/solutions/ice/ModelProcessorThermal.m

@@ -49 +49 @@
 elements(pos)=SetStructureField(elements(pos),'element','matid',pos);
 elements(pos)=SetStructureField(elements(pos),'element','geothermalflux',md.geothermalflux(md.elements(pos,1:6)));
+%For penta elements where we want to implement MacAyeal's element, we need to collapse 
+%the formulation into trias: 
+el6pos=find((md.elements_type(:,1)==macayealenum() | md.elements_type(:,1)==pattynenum()) & (element_partitioning==labindex));
+pos=find(ismember(el6pos,find(md.elements_type(:,1)==macayealenum())));
+elements(pos)=SetStructureField(elements(pos),'element','collapse',ones(length(pos),1));
+
 if strcmpi(solutiontype,'thermalsteady'),
         elements(pos)=SetStructureField(elements(pos),'element','thermal_steadystate',1);