Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 25726)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 25727)
@@ -5446,6 +5446,6 @@
 
 	/*Computational flags:*/
-	bool computerigid = true;
-	bool computeelastic = true;
+	bool computerigid = false;
+	bool computeelastic = false;
 	int  horiz;
 
@@ -5463,5 +5463,4 @@
 	this->inputs->SetDoubleInput(AreaEnum,this->lid,area);
 	this->AddInput(SealevelAreaEnum,&area,P0Enum);
-	if(!computerigid & !computeelastic)return;
 
 	/*recover precomputed green function kernels:*/
@@ -5519,10 +5518,6 @@
 
 				/*Rigid earth gravitational perturbation: */
-				if(computerigid){
-					G[i]+=G_rigid_precomputed[index];
-				}
-				if(computeelastic){
-					G[i]+=G_elastic_precomputed[index];
-				}
+				G[i]+=G_rigid_precomputed[index];
+				G[i]+=G_elastic_precomputed[index];
 				G[i]=G[i]*constant;
 
@@ -5609,14 +5604,14 @@
 			this->inputs->SetDoubleInput(AreaEnum,this->lid,area);
 			this->AddInput(SealevelAreaEnum,&area,P0Enum);
-
-			/*Free allocations:*/
-			#ifdef _HAVE_RESTRICT_
-			delete indices;
-			delete G;
-			#else
-			xDelete(indices);
-			xDelete(G);
-			#endif
-		}
+		}
+
+		/*Free allocations:*/
+		#ifdef _HAVE_RESTRICT_
+		delete indices;
+		delete G;
+		#else
+		xDelete(indices);
+		xDelete(G);
+		#endif
 	} else { /*computerigid==false*/
 		// do nothing
@@ -5635,5 +5630,5 @@
 	bool notfullygrounded=false;
 	bool scaleoceanarea= false;
-	bool computeelastic= true;
+	bool computerigid= false;
 
 	/*output: */
@@ -5680,9 +5675,9 @@
 	rho_ice=FindParam(MaterialsRhoIceEnum);
 	rho_water=FindParam(MaterialsRhoSeawaterEnum);
-	this->parameters->FindParam(&computeelastic,SolidearthSettingsElasticEnum);
+	this->parameters->FindParam(&computerigid,SolidearthSettingsRigidEnum);
 	this->parameters->FindParam(&scaleoceanarea,SolidearthSettingsOceanAreaScalingEnum);
 
 	/*retrieve precomputed G:*/
-	if(computeelastic)this->inputs->GetArrayPtr(SealevelriseGEnum,this->lid,&G,&gsize);
+	if(computerigid)this->inputs->GetArrayPtr(SealevelriseGEnum,this->lid,&G,&gsize);
 
 	/*Get area of element: precomputed in the sealevelrise_core_geometry:*/
@@ -5734,5 +5729,5 @@
 	_assert_(!xIsNan<IssmDouble>(bslrice));
 
-	if(computeelastic){
+	if(computerigid){
 		/*convert from m to kg/m^2:*/
 		I=I*rho_ice*phi;
@@ -5755,5 +5750,5 @@
 	bool notfullygrounded=false;
 	bool scaleoceanarea= false;
-	bool computeelastic= true;
+	bool computerigid= false;
 
 	/*elastic green function:*/
@@ -5784,9 +5779,9 @@
 	rho_water=FindParam(MaterialsRhoSeawaterEnum);
 	rho_freshwater=FindParam(MaterialsRhoFreshwaterEnum);
-	this->parameters->FindParam(&computeelastic,SolidearthSettingsElasticEnum);
+	this->parameters->FindParam(&computerigid,SolidearthSettingsRigidEnum);
 	this->parameters->FindParam(&scaleoceanarea,SolidearthSettingsOceanAreaScalingEnum);
 
 	/*retrieve precomputed G:*/
-	if(computeelastic)this->inputs->GetArrayPtr(SealevelriseGEnum,this->lid,&G,&gsize);
+	if(computerigid)this->inputs->GetArrayPtr(SealevelriseGEnum,this->lid,&G,&gsize);
 
 	/*Get area of element: precomputed in the sealevelrise_core_geometry:*/
@@ -5804,5 +5799,5 @@
 	_assert_(!xIsNan<IssmDouble>(bslrhydro));
 
-	if(computeelastic){
+	if(computerigid){
 		/*convert from m to kg/m^2:*/
 		W=W*rho_freshwater*phi;
@@ -5823,5 +5818,5 @@
 	IssmDouble BP;  //change in bottom pressure (Farrel and Clarke, Equ. 4)
 	IssmDouble constant;
-	bool computeelastic= true;
+	bool computeelastic= false;
 
 	/*elastic green function:*/
@@ -5925,5 +5920,5 @@
 
 	/*computational flags:*/
-	bool computeelastic= true;
+	bool computeelastic= false;
 
 	/*early return if we are not on the ocean or on an ice cap:*/
Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 25726)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 25727)
@@ -4837,9 +4837,7 @@
 
 	bool computerigid = true;
-	bool computeelastic= true;
 
 	/*recover computational flags: */
 	this->parameters->FindParam(&computerigid,SolidearthSettingsRigidEnum);
-	this->parameters->FindParam(&computeelastic,SolidearthSettingsElasticEnum);
 
 	/*Initialize temporary vector that will be used to sum eustatic components on all local elements, prior
@@ -4853,5 +4851,5 @@
 
 	/*Call the sea level rise non-eustatic core only if required: */
-	if(computerigid | computeelastic){
+	if(computerigid){
 		for(Object* & object : this->elements->objects){
 			Element* element = xDynamicCast<Element*>(object);