Index: /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/solvebarystatic.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/solvebarystatic.m	(revision 27139)
+++ /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/solvebarystatic.m	(revision 27140)
@@ -20,12 +20,11 @@
 	for t=1:nt-1, 
 	
-		time=times(t+1);
+		time=times(t);
 
 		%recover masks and thickness changes: 
-		tdH=md.masstransport.spcthickness(end,:);
+		tdH=md.mmemasstransport.thickness(end,:);
 		ind=findind(time,tdH);
-		dHtv=md.masstransport.spcthickness(1:end-1,ind)-md.geometry.thickness;
-		dHt=dHtv(md.mesh.elements)*[1;1;1]/3;
-		
+		dHt=md.mmemasstransport.thickness(1:end-1,ind+1)- md.mmemasstransport.thickness(1:end-1,ind);
+	
 		tocean=md.mask.ocean_levelset(end,:); 
 		ind=findind(time,tocean);