Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 25288)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 25289)
@@ -4957,5 +4957,5 @@
 
 	/*gia solution parameters:*/
-	IssmDouble lithosphere_thickness,mantle_viscosity;
+	IssmDouble lithosphere_thickness,mantle_viscosity,ice_mask;
 
 	/*output: */
@@ -4971,5 +4971,5 @@
 	this->parameters->FindParam(&gsize,MeshNumberofverticesEnum);
 	this->parameters->FindParam(&yts,ConstantsYtsEnum);
-
+	
 	/*recover gia solution parameters: */
 	int cross_section_shape;
@@ -4986,4 +4986,11 @@
 	IssmDouble mantle_density            = FindParam(MaterialsMantleDensityEnum);
 	IssmDouble rho_ice                   = FindParam(MaterialsRhoIceEnum);
+
+	/*recover ice mask*/ 
+	Input2* ice_mask_input=this->GetInput2(MaskIceLevelsetEnum); 
+	ice_mask_input->GetInputAverage(&ice_mask);
+
+	/*early return if we are NOT on an icy element:*/
+	if (ice_mask==1) return; 
 
 	/*pull thickness averages! */