Fictious Ice-thickness incorrectly filtered with fine (vertical) reslution meshs

[andrewdnolan]

Given the higher mesh densities at the boundaries (and therefore finer mesh resolutions), the simple approach for filtering the fictitious ice thickness added by Elmer ( $h_{\rm min}$ ) no longer works correctly:

thermal-structure/src/thermal/open.py

Line 108 in fe33d8b

ds["height"] = xr.where(ds.height <= h_min + 0.1, 0, ds.height)

If the vertical mesh resolution is fine enough, the ice thickness in the second (or even third) vertical nodes from the bottom will still be less then $h_{\rm min}$ and therefore forced to zero.

Proposed Solution:

• Can we use the calculate_length function to aid in the filtering? Since it only used the ice-thickness along the free surface it shouldn't be subject to this same bug.