Changeset 8413 for issm/trunk/src/c/objects/Elements/Penta.cpp

Timestamp:

05/24/11 11:51:41 (14 years ago)

Author:

seroussi

Message:

finally found the problem with PVectorThermalSheet

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (2 diffs)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -3050 +3050 @@
         int        analysis_type,drag_type;
         double     xyz_list[NUMVERTICES][3];
-        double     xyz_list_tria[NUMVERTICES2D][3];
+        double     xyz_list_tria[NUMVERTICES2D][3]={0.0};
         double     Jdet2d,dt;
         double     rho_ice,heatcapacity,geothermalflux_value;
@@ -3067 +3067 @@
         /*Retrieve all inputs and parameters*/
         GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
-        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
+        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<2;j++) xyz_list_tria[i][j]=xyz_list[i][j];
         parameters->FindParam(&analysis_type,AnalysisTypeEnum);
         rho_ice=matpar->GetRhoIce();