Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 21281)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 21282)
@@ -3698,5 +3698,5 @@
 	IssmDouble late,longe,re;
 	IssmDouble lati,longi,ri;
-	IssmDouble rho_ice,rho_water,rho_earth;
+	IssmDouble rho_ice,rho_earth;
 	IssmDouble minlong=400;
 	IssmDouble maxlong=-20;
@@ -3715,10 +3715,14 @@
 	bool store_green_functions=false;
 
+	/*Compute ice thickness change: */
+	Input*	deltathickness_input=inputs->GetInput(EsaDeltathicknessEnum); 
+	if (!deltathickness_input)_error_("delta thickness input needed to compute elastic adjustment!");
+	deltathickness_input->GetInputAverage(&I);
+		
 	/*early return if we are not on the (ice) loading point: */ 
-	//if(!(this->inputs->Max(MaskIceLevelsetEnum)<0)) return; 
+	if(I==0) return; 
 
 	/*recover material parameters: */
 	rho_ice=matpar->GetMaterialParameter(MaterialsRhoIceEnum);
-	rho_water=matpar->GetMaterialParameter(MaterialsRhoFreshwaterEnum);
 	rho_earth=matpar->GetMaterialParameter(MaterialsEarthDensityEnum);
 
@@ -3779,9 +3783,4 @@
 	H_parameter->GetParameterValueByPointer(&H_elastic_precomputed,&M);
 
-	/*Compute ice thickness change: */
-	Input*	deltathickness_input=inputs->GetInput(EsaDeltathicknessEnum); 
-	if (!deltathickness_input)_error_("delta thickness input needed to compute elastic adjustment!");
-	deltathickness_input->GetInputAverage(&I);
-		
 	/*initialize: */
 	U_elastic=xNewZeroInit<IssmDouble>(gsize);