Index: /issm/trunk-jpl/src/c/classes/BarystaticContributions.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/BarystaticContributions.cpp	(revision 26284)
+++ /issm/trunk-jpl/src/c/classes/BarystaticContributions.cpp	(revision 26285)
@@ -131,32 +131,11 @@
 
 } /*}}}*/
-void BarystaticContributions::Reset(int eid){ /*{{{*/
+void BarystaticContributions::Reset(){ /*{{{*/
+
+
+	ice->SetValue(0);
+	hydro->SetValue(0);
+	ocean->SetValue(0);
 	
-	int id;
-	if(nice){
-		id=reCast<int>(pice[eid]);
-		ice->SetValue(id,0.,INS_VAL);
-	}
-	else{
-		ice->SetValue(0,0.,INS_VAL);
-	}
-
-	if(nhydro){
-		id=reCast<int>(phydro[eid]);
-		hydro->SetValue(id,0.,INS_VAL);
-	}
-	else{
-		hydro->SetValue(0,0.,INS_VAL);
-	}
-
-	if(nocean){
-		id=reCast<int>(pocean[eid]);
-		ocean->SetValue(id,0.,INS_VAL);
-	}
-	else{
-		ocean->SetValue(0,0.,INS_VAL);
-	}
-		
-
 } /*}}}*/
 void BarystaticContributions::Save(Results* results, Parameters* parameters, IssmDouble oceanarea){ /*{{{*/
Index: /issm/trunk-jpl/src/c/classes/BarystaticContributions.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/BarystaticContributions.h	(revision 26284)
+++ /issm/trunk-jpl/src/c/classes/BarystaticContributions.h	(revision 26285)
@@ -42,5 +42,5 @@
 		void Save(Results* results, Parameters* parameters, IssmDouble oceanarea);
 		void Set(int eid, IssmDouble icevalue, IssmDouble hydrovalue, IssmDouble oceanvalue);
-		void Reset(int eid);
+		void Reset();
 
 };
Index: /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp	(revision 26284)
+++ /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp	(revision 26285)
@@ -388,10 +388,5 @@
 		barycontrib->Cumulate(femmodel->parameters);
 		barycontrib->Save(femmodel->results,femmodel->parameters,oceanarea);
-
-		//Avoid double counting barystatic contributions in future time steps: reset all non-cumulative bscl to 0
-		for(Object* & object : femmodel->elements->objects){
-			Element* element = xDynamicCast<Element*>(object);
-			barycontrib->Reset(element->Sid());
-		}
+		barycontrib->Reset();
 	}