The following piece of code in io_utils.py is supposed to avoid extrapolation by changing tend:
if "segments" in ndD_s["profileType"]:
dD_s["tend"] = dD_s["segments_tvec"][-1]
# Pad the last segment so no extrapolation occurs
dD_s["segments_tvec"][-1] = dD_s["tend"]*1.01
However, this change is only applied to the dimensional config dict, after the non-dimensional version was already created from it a few lines above. The simulation used only the non-dimensional config, so the old tend is used instead of the updated value.