Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 23546)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 23547)
@@ -4830,19 +4830,20 @@
 			for(int j=0;j<elements->Size();j++){
 				/*Intermediaries*/
-				Element*   element            =xDynamicCast<Element*>(elements->GetObjectByOffset(j));
-				int        numvertices        =element->GetNumberOfVertices();
-				IssmDouble zeros[numvertices] ={0.0};
+				Element*    element     = xDynamicCast<Element*>(elements->GetObjectByOffset(j));
+				int         numvertices = element->GetNumberOfVertices();
+				IssmDouble* zeros       = xNewZeroInit<IssmDouble>(numvertices);
 				switch(element->ObjectEnum()){
-				case TriaEnum:
-					element->inputs->AddInput(new TriaInput(stackedinput_enum[i],&zeros[0],P1Enum));
-					break;
-				case PentaEnum:
-					element->inputs->AddInput(new PentaInput(stackedinput_enum[i],&zeros[0],P1Enum));
-					break;
-				case TetraEnum:
-					element->inputs->AddInput(new TetraInput(stackedinput_enum[i],&zeros[0],P1Enum));
-					break;
-				default: _error_("Not implemented yet");
+					case TriaEnum:
+						element->inputs->AddInput(new TriaInput(stackedinput_enum[i],&zeros[0],P1Enum));
+						break;
+					case PentaEnum:
+						element->inputs->AddInput(new PentaInput(stackedinput_enum[i],&zeros[0],P1Enum));
+						break;
+					case TetraEnum:
+						element->inputs->AddInput(new TetraInput(stackedinput_enum[i],&zeros[0],P1Enum));
+						break;
+					default: _error_("Not implemented yet");
 				}
+				xDelete<IssmDouble>(zeros);
 			}
 		}