Changeset 27215

Timestamp:

08/17/22 15:02:27 (3 years ago)

Author:

Eric.Larour

Message:

CHG: no need to have inputs present for ice,water and ocean loads.

File:

• issm/branches/trunk-larour-SLPS2022/src/c/classes/Elements/Tria.cpp (modified) (2 diffs)

issm/branches/trunk-larour-SLPS2022/src/c/classes/Elements/Tria.cpp

@@ -7399 +7399 @@
 
         /*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 +7421 @@
         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++){