Changeset 946 for issm/trunk/src/m/classes/@pentaelem/CreatePVector.m

Timestamp:

06/12/09 11:13:12 (16 years ago)

Author:

seroussi

Message:

fixed problem with thermal ice

File:

• issm/trunk/src/m/classes/@pentaelem/CreatePVector.m (modified) (2 diffs)

issm/trunk/src/m/classes/@pentaelem/CreatePVector.m

@@ -657 +657 @@
         vxvy_list= vxvyvz_list(:,1:2);
 
+        %Update material if temperature is provided
+        if temperature_is_present,
+                temperature=1/6*[1 1 1 1 1 1]*temperature_list;
+                matice.B=paterson(temperature);
+        end
+ 
         % Get gaussian points and weights. Penta is an extrusion of a Tria, we therefore 
         %get tria gaussian points as well as segment gaussian points. For tria gaussian 
@@ -673 +679 @@
                         gauss_coord=[first_area_gauss_coord(igarea) second_area_gauss_coord(igarea) third_area_gauss_coord(igarea) vert_gauss_coord(igvert)];
                    
-                        %Update material if temperature is provided
+                        %Update temperature at gaussian point
                         if temperature_is_present,
                                 temperature=GetParameterValue(pentaelem,temperature_list,gauss_coord);
-                                matice.B=paterson(temperature);
                         end