function density_h = TILMedia.Internals.SLEMediumFunctions.density_h(d_s = sleMedium.d_s, d_l = sleMedium.d_l, h_fusion = sleMedium.h_fusion);