When I migrated from version 4.17 to 4.21, I noticed the output of md.results.TransientSolution.IceVolume
changed, in one case from 1.23E13 to 1.18E13, although the thickness and glaciated area remained virtually the same.
When I manually approximated the volume, the results (1.24E13) are closer to the output of v4.17. Any idea what might have changed?
area1=GetAreas(md.mesh.elements2d,md.mesh.x2d,md.mesh.y2d);
volume1=sum(area1.*md.geometry.thickness(md.mesh.elements2d(:,1)));