Changeset 8411

Timestamp:

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

Author:

seroussi

Message:

some problems in thermal sheets

File:

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

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

@@ -3053 +3053 @@
         double     Jdet2d,dt;
         double     rho_ice,heatcapacity,geothermalflux_value;
-        double     basalfriction,alpha2,vx,vy,pressure;
-        double     pressure_list[3];
+        double     basalfriction,alpha2,vx,vy;
         double     scalar;
         double     l1l6[NUMVERTICES];
@@ -3068 +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();
@@ -3075 +3074 @@
         Input* vx_input=inputs->GetInput(VxEnum);                         _assert_(vx_input);
         Input* vy_input=inputs->GetInput(VyEnum);                         _assert_(vy_input);
+        Input* vz_input=inputs->GetInput(VzEnum);                         _assert_(vz_input);
         Input* geothermalflux_input=inputs->GetInput(GeothermalFluxEnum); _assert_(geothermalflux_input);