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

Timestamp:

07/26/11 07:50:37 (14 years ago)

Author:

Eric.Larour

Message:

Fixed bug in nightly runs, dealing with requested_outputs, which were 0 sized in case
user did not request any.

File:

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

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

@@ -7725 +7725 @@
         double epsilon[6];
         double     viscousheating[NUMVERTICES]={0,0,0,0,0,0};
+        double     thickness;
         GaussPenta *gauss=NULL;
 
@@ -7735 +7736 @@
         Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input);
         Input* vz_input=inputs->GetInput(VzEnum); _assert_(vz_input);
+        Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
 
         /*loop over vertices: */
@@ -7740 +7742 @@
         for (int iv=0;iv<NUMVERTICES;iv++){
                 gauss->GaussVertex(iv);
+                
+                thickness_input->GetParameterValue(&thickness,gauss);
 
                 this->GetStrainRate3d(&epsilon[0],&xyz_list[0][0],gauss,vx_input,vy_input,vz_input);
                 matice->GetViscosity3dStokes(&viscosity,&epsilon[0]);
                 GetPhi(&phi, &epsilon[0], viscosity);
-
-                viscousheating[iv]=phi;
+                
+
+                viscousheating[iv]=phi*thickness;
         }