Index: /issm/branches/trunk-larour-SLPS2022/src/c/cores/sealevelchange_core.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/cores/sealevelchange_core.cpp	(revision 27220)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/cores/sealevelchange_core.cpp	(revision 27221)
@@ -708,11 +708,23 @@
 void              sealevelchange_finalize(FemModel* femmodel) {  /*{{{*/
 
+	bool isuq=false;
+	
 	BarystaticContributions* barycontrib=NULL;
 	GenericParam<BarystaticContributions*>* barycontribparam=NULL;
 	
-	/*Erase barycontrib object: */
-	barycontribparam = xDynamicCast<GenericParam<BarystaticContributions*>*>(femmodel->parameters->FindParamObject(BarystaticContributionsEnum));
-	barycontrib=barycontribparam->GetParameterValue();
-	delete barycontrib;
+	femmodel->parameters->FindParam(&isuq,QmuIsdakotaEnum);
+
+	if(isuq){
+		//reset barycontrib object:
+		barycontribparam = xDynamicCast<GenericParam<BarystaticContributions*>*>(femmodel->parameters->FindParamObject(BarystaticContributionsEnum));
+		barycontrib=barycontribparam->GetParameterValue(); 
+		barycontrib->Reset();
+	}
+	else {
+		/*Erase barycontrib object: */
+		barycontribparam = xDynamicCast<GenericParam<BarystaticContributions*>*>(femmodel->parameters->FindParamObject(BarystaticContributionsEnum));
+		barycontrib=barycontribparam->GetParameterValue();
+		delete barycontrib;
+	}
 
 	return;