Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 26106)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 26107)
@@ -5600,5 +5600,5 @@
 	/*elastic green function:*/
 	int index;
-	int         M;
+	int M,Mtest;
 
 	/*Computational flags:*/
@@ -5625,4 +5625,5 @@
 	DoubleVecParam* parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeGRigidEnum)); _assert_(parameter);
 	parameter->GetParameterValueByPointer((IssmDouble**)&G_rigid_precomputed,&M);
+	_assert_(M>0);
 
 	/*allocate indices:*/
@@ -5632,11 +5633,11 @@
 	if(computeelastic){
 		parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeGElasticEnum)); _assert_(parameter);
-		parameter->GetParameterValueByPointer((IssmDouble**)&G_elastic_precomputed,&M);
+		parameter->GetParameterValueByPointer((IssmDouble**)&G_elastic_precomputed,&Mtest); _assert_(Mtest==M);
 
 		parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeHElasticEnum)); _assert_(parameter);
-		parameter->GetParameterValueByPointer((IssmDouble**)&H_elastic_precomputed,&M);
+		parameter->GetParameterValueByPointer((IssmDouble**)&H_elastic_precomputed,&Mtest); _assert_(Mtest==M);
 
 		parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeUElasticEnum)); _assert_(parameter);
-		parameter->GetParameterValueByPointer((IssmDouble**)&U_elastic_precomputed,&M);
+		parameter->GetParameterValueByPointer((IssmDouble**)&U_elastic_precomputed,&Mtest); _assert_(Mtest==M);
 
 		GU=xNewZeroInit<IssmDouble>(gsize);
@@ -5673,4 +5674,5 @@
 		indices[i]=alpha/M_PI*reCast<IssmDouble,int>(M-1);
 		index=reCast<int,IssmDouble>(indices[i]);
+		_assert_(index>=0); _assert_(index<M);
 
 		/*Rigid earth gravitational perturbation: */