I suggest that the bedrock elevation map is directly downloaded from https://www.research-collection.ethz.ch/handle/20.500.11850/434697, rather than calculated as above. The folder 07_GlacierBed_SGI provide already cropped bedrock map for individual glaciers.
Then, maybe, you should rather calculate the ice thickness as IceThick = SurfElev- BedElev, since the ice thickness will evolve with time, whereas the bedrock elevation not, and should thus be considered as a reference.
https://github.com/eth-vaw-glaciology/GlacioTools.jl/blob/daf726e11bb4a75938d3fb8ec7ddde8a66170130/src/Alpine_glaciers.jl#L102
I suggest that the bedrock elevation map is directly downloaded from
https://www.research-collection.ethz.ch/handle/20.500.11850/434697, rather than calculated as above. The folder07_GlacierBed_SGIprovide already cropped bedrock map for individual glaciers.Then, maybe, you should rather calculate the ice thickness as
IceThick = SurfElev- BedElev, since the ice thickness will evolve with time, whereas the bedrock elevation not, and should thus be considered as a reference.