Index: /issm/branches/trunk-larour-SLPS2022/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/classes/Elements/Tria.cpp	(revision 27214)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/classes/Elements/Tria.cpp	(revision 27215)
@@ -7399,7 +7399,7 @@
 
 	/*Inputs:*/
-	IssmDouble I[NUMVERTICES]; 
-	IssmDouble W[NUMVERTICES];
-	IssmDouble BP[NUMVERTICES];
+	IssmDouble I[NUMVERTICES]={0,0,0}; 
+	IssmDouble W[NUMVERTICES]={0,0,0};
+	IssmDouble BP[NUMVERTICES]={0,0,0};
 	IssmDouble* areae=NULL;
 
@@ -7421,8 +7421,8 @@
 	this->parameters->FindParam(&areae,&nel,AreaeEnum);
 
-	/*Retrieve inputs:*/
-	Element::GetInputListOnVertices(&I[0],DeltaIceThicknessEnum);
-	Element::GetInputListOnVertices(&W[0],DeltaTwsEnum);
-	Element::GetInputListOnVertices(&BP[0],DeltaBottomPressureEnum);
+	/*Retrieve inputs, if they have been computed!:*/
+	if(this->GetInput(DeltaIceThicknessEnum)) Element::GetInputListOnVertices(&I[0],DeltaIceThicknessEnum); 
+	if(this->GetInput(DeltaTwsEnum)) Element::GetInputListOnVertices(&W[0],DeltaTwsEnum);
+	if(this->GetInput(DeltaBottomPressureEnum)) Element::GetInputListOnVertices(&BP[0],DeltaBottomPressureEnum);
 
 	for(int i=0;i<NUMVERTICES;i++){