Changeset 27140

Timestamp:

07/02/22 11:00:53 (3 years ago)

Author:

Eric.Larour

Message:

CHG: times are back in years!

File:

• issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/solvebarystatic.m (modified) (1 diff)

issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/solvebarystatic.m

@@ -20 +20 @@
         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);